MOLECULAR MODELING OF EPON 862-DETDA / CARBON COMPOSITES

The thermoset epoxy resin EPON 862, coupled with the DETDA hardening agent, are utilized as the polymer matrix component in many graphite (carbon fiber) composites. Because it is difficult to experimentally characterize the interfacial region, computational molecular modeling is a necessary tool for understanding the influence of the interfacial molecular structure on bulk-level material properties. The purpose of this research is to investigate the many possible variables that may influence the interfacial structure and the effect they will have on the mechanical behavior of the bulk level composite. Molecular models are established for EPON 862-DETDA polymer in the presence of a graphite surface. Material characteristics such as polymer mass-density, residual stresses, and molecular potential energy are investigated near the polymer/fiber interface. Because the exact degree of crosslinking in these thermoset systems is not known, many different crosslink densities (degrees of curing) are investigated. It is determined that a region exists near the carbon fiber surface in which the polymer mass density is different than that of the bulk mass density. These surface effects extend ~10 Å into the polymer from the center of the outermost graphite layer. Early simulations predict polymer residual stress levels to be higher near the graphite surface. It is also seen that the molecular potential energy in polymer atoms decreases with increasing crosslink density. New models are then established in order to investigate the interface between EPON 862-DETDA polymer and graphene nanoplatelets (GNPs) of various atomic thicknesses. Mechanical properties are extracted from the models using Molecular Dynamics techniques. These properties are then implemented into micromechanics software that utilizes the generalized method of cells to create representations of macro-scale composites. Micromechanics models are created representing GNP doped epoxy with varying number of graphene layers and interfacial polymer crosslink densities. The initial micromechanics results for the GNP doped epoxy are then taken to represent the matrix component and are re-run through the micromechanics software with the addition of a carbon fiber to simulate a GNP doped epoxy/carbon fiber composite. Micromechanics results agree well with experimental data, and indicate GNPs of 1 to 2 atomic layers to be highly favorable. The effect of oxygen bonded to the surface of the GNPs is lastly investigated. Molecular Models are created for systems with varying graphene atomic thickness, along with different amounts of oxygen species attached to them. Models are created for graphene containing hydroxyl groups only, epoxide groups only, and a combination of epoxide and hydroxyl groups. Results show models of oxidized graphene to decrease in both tensile and shear modulus. Attaching only epoxide groups gives the best results for mechanical properties, though pristine graphene is still favored.

Characterization of thermal and mechanical properties of polypropylene-based composites for fuel cell bipolar plates and development of educational tools in hydrogen and fuel cell technologies

In this project we developed conductive thermoplastic resins by adding varying amounts of three different carbon fillers: carbon black (CB), synthetic graphite (SG) and multi-walled carbon nanotubes (CNT) to a polypropylene matrix for application as fuel cell bipolar plates. This component of fuel cells provides mechanical support to the stack, circulates the gases that participate in the electrochemical reaction within the fuel cell and allows for removal of the excess heat from the system. The materials fabricated in this work were tested to determine their mechanical and thermal properties. These materials were produced by adding varying amounts of single carbon fillers to a polypropylene matrix (2.5 to 15 wt.% Ketjenblack EC-600 JD carbon black, 10 to 80 wt.% Asbury Carbon's Thermocarb TC-300 synthetic graphite, and 2.5 to 15 wt.% of Hyperion Catalysis International's FIBRILTM multi-walled carbon nanotubes) In addition, composite materials containing combinations of these three fillers were produced. The thermal conductivity results showed an increase in both through-plane and in-plane thermal conductivities, with the largest increase observed for synthetic graphite. The Department of Energy (DOE) had previously set a thermal conductivity goal of 20 W/m·K, which was surpassed by formulations containing 75 wt.% and 80 wt.% SG, yielding in-plane thermal conductivity values of 24.4 W/m·K and 33.6 W/m·K, respectively. In addition, composites containing 2.5 wt.% CB, 65 wt.% SG, and 6 wt.% CNT in PP had an in–plane thermal conductivity of 37 W/m·K. Flexural and tensile tests were conducted. All composite formulations exceeded the flexural strength target of 25 MPa set by DOE. The tensile and flexural modulus of the composites increased with higher concentration of carbon fillers. Carbon black and synthetic graphite caused a decrease in the tensile and flexural strengths of the composites. However, carbon nanotubes increased the composite tensile and flexural strengths. Mathematical models were applied to estimate through-plane and in-plane thermal conductivities of single and multiple filler formulations, and tensile modulus of single-filler formulations. For thermal conductivity, Nielsen's model yielded accurate thermal conductivity values when compared to experimental results obtained through the Flash method. For prediction of tensile modulus Nielsen's model yielded the smallest error between the predicted and experimental values. The second part of this project consisted of the development of a curriculum in Fuel Cell and Hydrogen Technologies to address different educational barriers identified by the Department of Energy. By the creation of new courses and enterprise programs in the areas of fuel cells and the use of hydrogen as an energy carrier, we introduced engineering students to the new technologies, policies and challenges present with this alternative energy. Feedback provided by students participating in these courses and enterprise programs indicate positive acceptance of the different educational tools. Results obtained from a survey applied to students after participating in these courses showed an increase in the knowledge and awareness of energy fundamentals, which indicates the modules developed in this project are effective in introducing students to alternative energy sources.

SYNTHESIS OF GRAPHENE AND ITS APPLICATIONS FOR DYE-SENSITIZED SOLAR CELLS

Graphene, which is a two-dimensional carbon material, exhibits unique properties that promise its potential applications in photovoltaic devices. Dye-sensitized solar cell (DSSC) is a representative of the third generation photovoltaic devices. Therefore, it is important to synthesize graphene with special structures, which possess excellent properties for dye-sensitized solar cells. This dissertation research was focused on (1) the effect of oxygen content on the structure of graphite oxide, (2) the stability of graphene oxide solution, (3) the application of graphene precipitate from graphene oxide solution as counter electrode for DSSCs, (4) the development of a novel synthesis method for the three-dimensional graphene with honeycomb-like structure, and (5) the exploration of honeycomb structured graphene (HSG) as counter electrodes for DSSCs. Graphite oxide is a crucial precursor to synthesize graphene sheets via chemical exfoliation method. The relationship between the oxygen content and the structures of graphite oxides was still not explored. In this research, the oxygen content of graphite oxide is tuned by changing the oxidation time and the effect of oxygen content on the structure of graphite oxide was evaluated. It has been found that the saturated ratio of oxygen to carbon is 0.47. The types of functional groups in graphite oxides, which are epoxy, hydroxyl, and carboxylgroups, are independent of oxygen content. However, the interplanar space and BET surface area of graphite oxide linearly increases with increasing O/C ratio. Graphene oxide (GO) can easily dissolve in water to form a stable homogeneous solution, which can be used to fabricate graphene films and graphene based composites. This work is the first research to evaluate the stability of graphene oxide solution. It has been found that the introduction of strong electrolytes (HCl, LiOH, LiCl) into GO solution can cause GO precipitation. This indicates that the electrostatic repulsion plays a critical role in stabilizing aqueous GO solution. Furthermore, the HCl-induced GO precipitation is a feasible approach to deposit GO sheets on a substrate as a Pt-free counter electrode for a dye-sensitized solar cell (DSSC), which exhibited 1.65% of power conversion efficiency. To explore broad and practical applications, large-scale synthesis with controllable integration of individual graphene sheets is essential. A novel strategy for the synthesis of graphene sheets with three-dimensional (3D) Honeycomb-like structure has been invented in this project based on a simple and novel chemical reaction (Li2O and CO to graphene and Li2CO3). The simultaneous formation of Li2CO3 with graphene not only can isolate graphene sheets from each other to prevent graphite formation during the process, but also determine the locally curved shape of graphene sheets. After removing Li2CO3, 3D graphene sheets with a honeycomb-like structure were obtained. This would be the first approach to synthesize 3D graphene sheets with a controllable shape. Furthermore, it has been demonstrated that the 3D Honeycomb-Structured Graphene (HSG) possesses excellent electrical conductivity and high catalytic activity. As a result, DSSCs with HSG counter electrodes exhibit energy conversion efficiency as high as 7.8%, which is comparable to that of an expensive noble Pt electrode.

CHARACTERIZATION OF THERMAL AND MECHANICAL PROPERTIES OF POLYPROPYLENE–BASED COMPOSITES FOR FUEL CELL BIPOLAR PLATES AND DEVELOPMENT OF EDUCATIONAL TOOLS IN HYDROGEN AND FUEL CELL TECHNOLOGIES

In this project we developed conductive thermoplastic resins by adding varying amounts of three different carbon fillers: carbon black (CB), synthetic graphite (SG) and multi–walled carbon nanotubes (CNT) to a polypropylene matrix for application as fuel cell bipolar plates. This component of fuel cells provides mechanical support to the stack, circulates the gases that participate in the electrochemical reaction within the fuel cell and allows for removal of the excess heat from the system. The materials fabricated in this work were tested to determine their mechanical and thermal properties. These materials were produced by adding varying amounts of single carbon fillers to a polypropylene matrix (2.5 to 15 wt.% Ketjenblack EC-600 JD carbon black, 10 to 80 wt.% Asbury Carbons’ Thermocarb TC-300 synthetic graphite, and 2.5 to 15 wt.% of Hyperion Catalysis International’s FIBRILTM multi-walled carbon nanotubes) In addition, composite materials containing combinations of these three fillers were produced. The thermal conductivity results showed an increase in both through–plane and in–plane thermal conductivities, with the largest increase observed for synthetic graphite. The Department of Energy (DOE) had previously set a thermal conductivity goal of 20 W/m·K, which was surpassed by formulations containing 75 wt.% and 80 wt.% SG, yielding in–plane thermal conductivity values of 24.4 W/m·K and 33.6 W/m·K, respectively. In addition, composites containing 2.5 wt.% CB, 65 wt.% SG, and 6 wt.% CNT in PP had an in–plane thermal conductivity of 37 W/m·K. Flexural and tensile tests were conducted. All composite formulations exceeded the flexural strength target of 25 MPa set by DOE. The tensile and flexural modulus of the composites increased with higher concentration of carbon fillers. Carbon black and synthetic graphite caused a decrease in the tensile and flexural strengths of the composites. However, carbon nanotubes increased the composite tensile and flexural strengths. Mathematical models were applied to estimate through–plane and in–plane thermal conductivities of single and multiple filler formulations, and tensile modulus of single–filler formulations. For thermal conductivity, Nielsen’s model yielded accurate thermal conductivity values when compared to experimental results obtained through the Flash method. For prediction of tensile modulus Nielsen’s model yielded the smallest error between the predicted and experimental values. The second part of this project consisted of the development of a curriculum in Fuel Cell and Hydrogen Technologies to address different educational barriers identified by the Department of Energy. By the creation of new courses and enterprise programs in the areas of fuel cells and the use of hydrogen as an energy carrier, we introduced engineering students to the new technologies, policies and challenges present with this alternative energy. Feedback provided by students participating in these courses and enterprise programs indicate positive acceptance of the different educational tools. Results obtained from a survey applied to students after participating in these courses showed an increase in the knowledge and awareness of energy fundamentals, which indicates the modules developed in this project are effective in introducing students to alternative energy sources.

In-plane thermal conductivity modeling of carbon filled liquid crystal polymer based resins

Adding conductive carbon fillers to insulating thermoplastic resins increases composite electrical and thermal conductivity. Often, as much of a single type of carbon filler is added to achieve the desired conductivity, while still allowing the material to be molded into a bipolar plate for a fuel cell. In this study, varying amounts of three different carbons (carbon black, synthetic graphite particles, and carbon fiber) were added to Vectra A950RX Liquid Crystal Polymer. The in-plane thermal conductivity of the resulting single filler composites were tested. The results showed that adding synthetic graphite particles caused the largest increase in the in-plane thermal conductivity of the composite. The composites were modeled using ellipsoidal inclusion problems to predict the effective in-plane thermal conductivities at varying volume fractions with only physical property data of constituents. The synthetic graphite and carbon black were modeled using the average field approximation with ellipsoidal inclusions and the model showed good agreement with the experimental data. The carbon fiber polymer composite was modeled using an assemblage of coated ellipsoids and the model showed good agreement with the experimental data.

Molecular modeling of EPON 862-DETDA polymer

EPON 862 is an epoxy resin which is cured with the hardening agent DETDA to form a crosslinked epoxy polymer and is used as a component in modern aircraft structures. These crosslinked polymers are often exposed to prolonged periods of temperatures below glass transition range which cause physical aging to occur. Because physical aging can compromise the performance of epoxies and their composites and because experimental techniques cannot provide all of the necessary physical insight that is needed to fully understand physical aging, efficient computational approaches to predict the effects of physical aging on thermo-mechanical properties are needed. In this study, Molecular Dynamics and Molecular Minimization simulations are being used to establish well-equilibrated, validated molecular models of the EPON 862-DETDA epoxy system with a range of crosslink densities using a united-atom force field. These simulations are subsequently used to predict the glass transition temperature, thermal expansion coefficients, and elastic properties of each of the crosslinked systems for validation of the modeling techniques. The results indicate that glass transition temperature and elastic properties increase with increasing levels of crosslink density and the thermal expansion coefficient decreases with crosslink density, both above and below the glass transition temperature. The results also indicate that there may be an upper limit to crosslink density that can be realistically achieved in epoxy systems. After evaluation of the thermo-mechanical properties, a method is developed to efficiently establish molecular models of epoxy resins that represent the corresponding real molecular structure at specific aging times. Although this approach does not model the physical aging process, it is useful in establishing a molecular model that resembles the physically-aged state for further use in predicting thermo-mechanical properties as a function of aging time. An equation has been predicted based on the results which directly correlate aging time to aged volume of the molecular model. This equation can be helpful for modelers who want to study properties of epoxy resins at different levels of aging but have little information about volume shrinkage occurring during physical aging.

Graphite exfoliation by supercritical carbon dioxide extraction

Supercritical carbon dioxide is used to exfoliate graphite, producing a small, several-layer graphitic flake. The supercritical conditions of 2000, 2500, and 3000 psi and temperatures of 40°, 50°, and 60°C, have been used to study the effect of critical density on the sizes and zeta potentials of the treated flakes. Photon Correlation Spectroscopy (PCS), Brunauer-Emmett-Teller (BET) surface area measurement, field emission scanning electron microscopy (FE-SEM), and atomic force microscopy (AFM) are used to observe the features of the flakes. N-methyl-2-pyrrolidinone (NMP), dimethylformamide (DMF), and isopropanol are used as co-solvents to enhance the supercritical carbon dioxide treatment. As a result, the PCS results show that the flakes obtained from high critical density treatment (low temperature and high pressure) are more stable due to more negative charges of zeta potential, but have smaller sizes than those from low critical density (high temperature and low pressure). However, when an additional 1-hour sonication is applied, the size of the flakes from low critical density treatment becomes smaller than those from high critical density treatment. This is probably due to more CO2 molecules stacked between the layers of the graphitic flakes. The zeta potentials of the sonicated samples were slightly more negative than nonsonicated samples. NMP and DMF co-solvents maintain stability and prevented reaggregation of the flakes better than isopropanol. The flakes tend to be larger and more stable as the treatment time increases since larger flat area of graphite is exfoliated. In these experiments, the temperature has more impact on the flakes than pressure. The BET surface area resultsshow that CO2 penetrates the graphite layers more than N2. Moreover, the negative surface area of the treated graphite indicates that the CO2 molecules may be adsorbed between the graphite layers during supercritical treatment. The FE-SEM and AFM images show that the flakes have various shapes and sizes. The effects of surfactants can be observed on the FE-SEM images of the samples in one percent by weight solution of SDBS in water since the sodium dodecylbenzene sulfonate (SDBS) residue covers all of the remaining flakes. The AFM images show that the vertical thickness of the graphitic flakes can ranges from several nanometers (less than ten layers thick), to more than a hundred nanometers. In conclusion, supercritical carbon dioxide treatment is a promising step compared to mechanical and chemical exfoliation techniques in the large scale production of thin graphitic flake, breaking down the graphite flakes into flakes only a fewer graphene layers thick.

Comparative study of anchorage strengths of epoxy coated hooked bars

Hooked reinforcing bars (rebar) are used frequently to carry the tension forces developed in beams and transferred to columns. Research into epoxy coated hooked bars has only been minimally performed and no research has been carried out incorporating the coating process found in ASTM A934. This research program compares hooked rebar that are uncoated, coated by ASTM A775, and coated by ASTM A934. In total, forty-two full size beam-column specimens were created, instrumented and tested to failure. The program was carried out in three phases. The first phase was used to refine the test setup and procedures. Phase two explored the spacing of column ties within the joint region. Phase three explored the three coating types found above. Each specimen included two hooked rebar which were loaded and measured independently for relative rebar slip. The load and displacement of the hooked rebar were analyzed, focusing on behavior at the levels of 30 ksi, 42 ksi and 60 ksi of rebar stress. Statistical and general comparisons were made using the coating types, tie spacing, and rebar stress level. Many of the parameters composing the rebar and concrete were also tested to characterize the components and specimens. All rebar tested met ASTM standards for tensile strength, but the newer ASTM A934 method seemed to produce slightly lower yield strengths. The A934 method also produced coating thicknesses that were very inconsistent and were higher than ASTM maximum limits in many locations. Continuity of coating surfaces was found to be less than 100% for both A775 and A934 rebar, but for different reasons. The many comparisons made did not always produce clear conclusions. The data suggests that the ACI Code (318-05) parameter of 1.2 for including epoxy coating on hooked rebar may need to be raised, possibly to 2.5, but more testing needs to be performed before such a large value change is set forth. This is particularly important as variables were identified which may have a larger influence on rebar capacity than the Development Length, of which the current 1.2 factor modifies. Many suggestions for future work are included throughout the thesis to help guide other researchers in carrying out successful and productive programs which will further the highly understudied topic of hooked rebar.

Photodegradation and photostabilization of weathered wood flour filled polyethylene composites

Wood plastic composites (WPCs) have gained popularity as building materials because of their usefulness in replacing solid wood in a variety of applications. These composites are promoted as being low-maintenance, high-durability products. However, it has been shown that WPCs exposed to weathering may experience a color change and/or loss in mechanical properties. An important requirement for building materials used in outdoor applications is the retention of their aesthetic qualities and mechanical properties during service life. Therefore, it is critical to understand the photodegradation mechanisms of WPCs exposed to UV radiation and to develop approaches to stabilize these composites (both unstabilized and stabilized) as well as the effect of weathering on the color fade and the retention of mechanical properties were characterized. Since different methods of manufacturing WPCs lead to different surface characteristics, which can influence weathering, the effect of manufacturing method on the photodegradation of WPCs was investigated first. Wood flour (WF) filled high-density polyethylene (HDPE) composite samples were either injection molded, extruded, or extruded and then planed. Fourier transform infrared (FTIR) spectroscopy was used to monitor the surface chemistry of the manufactured composites. The spectra showed that the surface of planed samples had more wood component than extruded and injection molded samples, respectively. After weathering, the samples were analyzed for color fade, and loss of flexural properties. The final lightness of the composites was not dependent upon the manufacturing method. However the mechanical property loss was dependent upon manufacturing method. The samples with more wood component at the surface (planed samples) experienced a larger percentage of total loss in flexural properties after weathering due to a greater effect of moisture on the samples. The change in surface chemistry of HDPE and WF/HDPE composites after weathering was studied using spectroscopic techniques. X-ray photoelectron spectroscopy (XPS) was used to characterize the occurrence of surface oxidation whereas FTIR spectroscopy was used to monitor the development of degradation products, such as carbonyl groups and vinyl groups, and to determine changes in HDPE crystallinity. Surface oxidation occurred immediately after exposure for both the neat HDPE and WF/HDPE composites. After weathering, the surface of the WF/HDPE composites was oxidized to a greater extent than the neat HDPE after weathering. This suggests that photodegradation is exacerbated by the addition of the carbonyl functional groups of the wood fibers within the HDPE atrix during composite manufacturing. While neat HDPE may undergo cross-linking in the initial stages of accelerated weathering, the WF may physically hinder the ability of the HDPE to cross-link resulting in the potential for HDPE chain scission to dominate in the initial weathering stages of the WF/HDPE composites. To determine which photostabilizers are most effective for WF/HDPE composites, factorial experimental designes were used to determine the effects of adding two hindered amine light stabilizers, an ultraviolet absorber, and a pigment on the color made and mechanical properties of both unweathered and UV weathered samples. Both the pigment and ultraviolet absorber were more effective photostabilizers for WF/HDPE composites than hinder amine light stabilizers. The ineffectiveness of hindered amine light stabilizers in protecting WPCs against UV radiation was attribuated to the acid/base reactions occurring between the WF and hindered amine light stabilizer. The efficiency of an ultraviolet absorber and/or pigment was also examined by incorporating different concentration of an ultraviolet absorber and/or pigment into WF/HDPE composites. Color change and flexural properties were determined after accelerated UV weathering. The lightness of the composite after weathering was influenced by the concentration of both the ultraviolet absorber by masking the bleaching wood component as well as blocking UV light. Flexural MOE loss was influenced by an increase in ultraviolet absorber concentration, but increasing pigment concentration from 1 to 2% had little influence on MOE loss. However, increasing both ultraviolet absorber and pigment concentration resulted in improved strength properties over the unstabilized composites after 3000 h of weather. Finally, the change in surface chemistry due to weathering of WF/HDPE composites that were either unstabilized or stabilized with an ultraviolet absorber and/or pigment was analyzed using FTIR spectroscopy. The samples were tested for loss in modulus of elasticity, carbonyl and vinyl group formation at the surface, and change in HDPE crystallinity. It was concluded that structural changes in the samples; carbonyl group formation, terminal vinyl group formation, and crystallinity changes cannot reliably be used to predict changes in modulus of elasticity using a simple linear relationship. The effect of cross-linking, chain scission, and crystallinity changes due to ultraviolet exposure as well as the interfacial degradation due to moisture exposure are inter-related factors when weathering HDPE and WF/HDPE composites.

Nanomechanics of cellulose crystals and cellulose-based polymer composites

Cellulose-polymer composites have potential applications in aerospace and transportation areas where lightweight materials with high mechanical properties are needed. In addition, these economical and biodegradable composites have been shown to be useful as polymer electrolytes, packaging structures, optoelectronic devices, and medical implants such as wound dressing and bone scaffolds. In spite of the above mentioned advantages and potential applications, due to the difficulties associated with synthesis and processing techniques, application of cellulose crystals (micro and nano sized) for preparation of new composite systems is limited. Cellulose is hydrophilic and polar as opposed to most of common thermoplastics, which are non-polar. This results in complications in addition of cellulose crystals to polymer matrices, and as a result in achieving sufficient dispersion levels, which directly affects the mechanical properties of the composites. As in other composite materials, the properties of cellulose-polymer composites depend on the volume fraction and the properties of individual phases (the reinforcement and the polymer matrix), the dispersion quality of the reinforcement through the matrix and the interaction between CNCs themselves and CNC and the matrix (interphase). In order to develop economical cellulose-polymer composites with superior qualities, the properties of individual cellulose crystals, as well as the effect of dispersion of reinforcements and the interphase on the properties of the final composites should be understood. In this research, the mechanical properties of CNC polymer composites were characterized at the macro and nano scales. A direct correlation was made between: Dispersion quality and macro-mechanical properties Nanomechanical properties at the surface and tensile properties CNC diameter and interphase thickness Lastly, individual CNCs from different sources were characterized and for the first time size-scale effect on their nanomechanical properties were reported. Then the effect of CNC surface modification on the mechanical properties was studied and correlated to the crystalline structure of these materials.

COMPUTATIONAL STUDY OF MICROSTRUCTURE-PROPERTYMECHANISM RELATIONS IN FERROIC COMPOSITES

Ferroic materials, as notable members of smart materials, have been widely used in applications that perform sensing, actuation and control. The macroscopic property change of ferroic materials may become remarkably large during ferroic phase transition, leading to the fact that the macroscopic properties can be tuned by carefully applying a suitable external field (electric, magnetic, stress). To obtain an enhancement in physical and/or mechanical properties, different kinds of ferroic composites have been fabricated. The properties of a ferroic composite are determined not only by the properties and relative amounts of the constituent phases, but also by the microstructure of individual phase such as the phase connectivity, phase size, shape and spatial arrangement. This dissertation mainly focuses on the computational study of microstructure – property – mechanism relations in two representative ferroic composites, i.e., two-phase particulate magnetoelectric (ME) composite and polymer matrix ferroelectric composite. The former is a great example of ferroic composite exhibiting a new property and functionality that neither of the constituent phases possesses individually. The latter well represents the kind of ferroic composites having property combinations that are better than the existing materials. Phase field modeling was employed as the computing tool, and the required models for ferroic composites were developed based on existing models for monolithic materials. Extensive computational simulations were performed to investigate the microstructure-property relations and the underlying mechanism in ferroic composites. In particulate, it is found that for ME composite 0-3 connectivity (isolated magnetostrictive phase) is necessary to exhibit ME effect, and small but finite electrical conductivity of isolated magnetic phase can beneficially enhance ME effect. It is revealed that longitudinal and transverse ME coefficients of isotropic 0-3 particulate composites can be effectively tailored by controlling magnetic domain structures without resort to anisotropic two-phase microstructures. Simulations also show that the macroscopic properties of the ferroelectricpolymer composites critically depend on the ferroelectric phase connectivity while are not sensitive to the sizes and internal grain structures of the ceramic particles. Texturing is found critical to exploit the paraelectric«ferroelectric phase transition and nonlinear polarization behavior in paraelectric polycrystal and its polymer matrix composite. Additionally, a Diffuse Interface Field model was developed to simulate packing and motion in liquid phase which is promising for studying the fabrication of particulatepolymer composites.