Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “composite models”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 505 records · Page 28

Processing and mechanical characterization of short carbon fiber-reinforced epoxy composites for material extrusion additive manufacturing

Fiber-reinforced polymer composites have been extensively utilized in recent years as feedstock materials for material extrusion additive manufacturing (AM) processes to improve strength, stiffness, and functionality of printed parts over unfilled printed polymers. However, the widespread adoption of AM of fiber-reinforced polymer composites requires a deeper understanding of the process-structure-property relationships in printed components, and such relationships are not well understood yet. Fiber length is critically important to the mechanical performance of short fiber composites, but very few studies to-date have focused on how the fiber length distribution (FLD) evolves during processing of composite feedstocks and how this evolution affects printing behavior and mechanical properties in 3D-printed composites. Here, FLD is measured for carbon fiber reinforced epoxy composites over a wide range of ink compositions and shear mixing times, and the distributions are fit with a Weibull-type distribution function. The effects of FLD on the tradeoff between ink processability, ink rheology, printing behavior and mechanical properties are investigated. Furthermore, the effects of printing parameters (nozzle size and print speed) on mechanical anisotropy and fiber orientation distribution (FOD) in printed composites are explored. Mechanical properties of printed composites are characterized via 3 pt-flexural testing, and microstructure is investigated using optical and scanning electron microscopy (SEM), and x-ray computed tomography. Finally, the fitted Weibull parameters are fed into a composite model that incorporates FLD and FOD, and model predictions are found to be in excellent agreement with experimental observations.

3D printing↗

Oxidation of Boron Nitride in Composites

Boron nitride (BN) is a prime candidate for fiber coatings in silicon carbide (SiC) fiber-reinforced SiC matrix composites. The properties of BN allow the fiber to impart beneficial composite properties to the matrix, even at elevated temperatures. The problem with BN is that it is readily attacked by oxygen. Although BN is an internal component of the composite, a matrix crack or pore can create a path for hot oxygen to attack the BN. This type of attack is not well understood. A variety of phenomena have been observed. These include borosilicate glass formation, volatilization of the BN, and under some conditions, preservation of the BN. In this study at the NASA Lewis Research Center, a series of BN materials and BN-containing model composites were methodically examined to understand the various issues dealing with the oxidation of BN in composites. Initial studies were done with a series of monolithic BN materials prepared by hot pressing and chemical vapor deposition (CVD). From these studies, we found that BN showed a strong orientation effect in oxidation and was extremely sensitive to the presence of water vapor in the environment. In addition, CVD material deposited at a high temperature showed much better oxidation behavior than CVD material deposited at a lower temperature.

Jacobson, Nathan S.↗

Preliminary IHLW Formulation Algorithm Description

This report documents the initial algorithm that could be used by the Waste Treatment and Immobilization Plant (WTP) in batching high-level waste (HLW) and glass-forming chemicals (GFCs) in the HLW melter feed preparation vessel (MFPV) (HFP-VSL-00001 and -00005). Not all Hanford tank waste can be accommodated by the models developed for this report and significant expansion of the model boundaries could be achievable to reduce the WTP mission life and total canister production count. The immobilized HLW (IHLW) must meet a series of constraints to be acceptable for disposal in the Monitored Geologic Repository, which are contained in the Specification 1 of the Contract (DOE 2000), the Waste Acceptance Product Specifications (WAPS, DOE 1996), and the Waste Acceptance System Requirements Document (WASRD, DOE 2007). The IHLW Waste Form Compliance Plan (WCP, 24590-HLW-PL-RT-07-0001, Rev 3) specifies that the formulation algorithm will be developed and used to comply with the constraints associated with glass composition and properties. This report is not an engineering calculation, does not provide design input, and is not an engineering study. Algorithm inputs include the chemical analyses of the blended HLW in the HLW blend vessel (HBV) (HLP-VSL-00028, the volume and composition of the MFPV heel, the volume and composition of the MFPV after waste addition, the volume and composition of MFPV batch after GFC addition, the compositions of individual GFCs, and the mass of glass in each canister. In addition to these inputs, uncertainties in the HLW composition and processing parameters are included in the algorithm. Using the above inputs, the algorithm calculates the following outputs: 1) the volume of HLW to be transferred from the HBV to the MFPV, 2) the mass of each GFC for addition to the MFPV, 3) the composition of the glass that will be produced along with uncertainties, and 4) the predicted properties, with associated uncertainties, of the resulting IHLW. The algorithm uses the property-composition models to calculate properties with associated uncertainties and compares them with various constraints to ensure that a processable feed is formulated and a compliant IHLW is produced. The GFC additions are determined using an optimization approach to provide high confidence that the HLW glass will meet all product quality requirements and key processing constraints. For most HLW batches there are many possible glass compositions that meet all constraints. In these cases, the glass composition is optimized for a series of target component concentrations and target property values. The algorithm also incorporates process measurement and product quality uncertainties, based on the work of Piepel et al. (2005). Estimates of the various process and measurement uncertainties that affect glass compositions and predicted glass properties have been previously reported (Piepel et al. 2005, 2006) and the impacts of these estimated uncertainties on the IHLW composition envelope that meets product quality and processing-related properties with sufficient confidence were evaluated. The details of work performed to date to develop this initial GFC addition and batching algorithm are summarized in Sections 4 and 5. An example data set is used to illustrate the calculations of the algorithm summarized in Section 6. Finally, in Section 7, there is a statement of the required work to achieve a final operational IHLW formulation control algorithm. This report is not an engineering calculation, does not provide design input, and is not an engineering study.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Interlaminar R-Curve Effects in Thermoplastic Composites: Experiments and Model Correlations

The interlaminar fracture behavior of thermoplastic composites under Mode-I and Mode-II static loading is investigated in this research. The material system under consideration is a carbon fiber reinforced low melt semi-crystalline resin TC1225 LMPAEK T700G (T700/LMPAEK)†. Characterization experiments were conducted to measure the load-displacement response, peak load, and propagation of delamination from a pre-implanted Kapton® insert. The Mode-I response was characterized using a double cantilever beam (DCB) test according to the ASTM D5528-22 standard. It was observed that the propagation exhibited significant resistance to delamination growth, i.e., an R-Curve effect. The test data indicated that the steady-state fracture toughness, G SS I , was approximately twice the mean initiation fracture toughness, G Ic . The Mode-II response of an end notch flexure (ENF) specimen tested according to the ASTM D7905-19 standard exhibited an unexpected stable behavior upon attaining peak load. Prior to attaining peak load, the Mode-II response exhibited some non-linearity, which may be attributed to nonlinear bulk material response, e.g., hyper-elasticity, viscoelasticity, plasticity, etc., and potentially friction at the load/support rollers and R-Curve effects. Finite element models based on the Cohesive Zone Method (CZM) and the Virtual Crack Closure Technique (VCCT) were used to predict the DCB and ENF specimen responses. DCB analysis results obtained using G SS I predicted propagation consistent with experimental data. However, results obtained with the as characterized G Ic were inaccurate in predicting delamination growth onset and peak load. ENF analysis results obtained with the as-characterized Mode-II fracture toughness, G IIc , were inaccurate in predicting the peak load and propagation. Hence, R-Curve effects were included in both DCB and ENF simulations, which enabled improved predictions. The experimentally determined Mode-I R-Curve was used for the DCB simulations, whereas an R-Curve was derived for Mode-II assuming that the change in compliance of the specimen was associated with crack growth only.

Prabhakar M Rao↗

Interlaminar R-Curve Effects in Thermoplastic Composites: Experiments and Model Correlations

The interlaminar fracture behavior of thermoplastic composites under Mode-I and Mode-II static loading is investigated in this research. The material system under consideration is a carbon fiber reinforced low melt semi-crystalline resin TC1225 LMPAEK T700G (T700/LMPAEK)†. Characterization experiments were conducted to measure the load-displacement response, peak load, and propagation of delamination from a pre-implanted Kapton® insert. The Mode-I response was characterized using a double cantilever beam (DCB) test according to the ASTM D5528-22 standard. It was observed that the propagation exhibited significant resistance to delamination growth, i.e., an R-Curve effect. The test data indicated that the steady-state fracture toughness, G SS I , was approximately twice the mean initiation fracture toughness, G Ic . The Mode-II response of an end notch flexure (ENF) specimen tested according to the ASTM D7905-19 standard exhibited an unexpected stable behavior upon attaining peak load. Prior to attaining peak load, the Mode-II response exhibited some non-linearity, which may be attributed to nonlinear bulk material response, e.g., hyper-elasticity, viscoelasticity, plasticity, etc., and potentially friction at the load/support rollers and R-Curve effects. Finite element models based on the Cohesive Zone Method (CZM) and the Virtual Crack Closure Technique (VCCT) were used to predict the DCB and ENF specimen responses. DCB analysis results obtained using G SS I predicted propagation consistent with experimental data. However, results obtained with the as characterized G Ic were inaccurate in predicting delamination growth onset and peak load. ENF analysis results obtained with the as-characterized Mode-II fracture toughness, G IIc , were inaccurate in predicting the peak load and propagation. Hence, R-Curve effects were included in both DCB and ENF simulations, which enabled improved predictions. The experimentally determined Mode-I R-Curve was used for the DCB simulations, whereas an R-Curve was derived for Mode-II assuming that the change in compliance of the specimen was associated with crack growth only.

Prabhakar Rao↗

Micromechanical models for textile structural composites

The objective is to develop micromechanical models for predicting the stiffness and strength properties of textile composite materials. Two models are presented to predict the homogeneous elastic constants and coefficients of thermal expansion of a textile composite. The first model is based on rigorous finite element analysis of the textile composite unit-cell. Periodic boundary conditions are enforced between opposite faces of the unit-cell to simulate deformations accurately. The second model implements the selective averaging method (SAM), which is based on a judicious combination of stiffness and compliance averaging. For thin textile composites, both models can predict the plate stiffness coefficients and plate thermal coefficients. The finite element procedure is extended to compute the thermal residual microstresses, and to estimate the initial failure envelope for textile composites.

Marrey, Ramesh V.↗

Modeling creep behavior of fiber composites

A micromechanical model for the creep behavior of fiber composites is developed based on a typical cell consisting of a fiber and the surrounding matrix. The fiber is assumed to be linearly elastic and the matrix nonlinearly viscous. The creep strain rate in the matrix is assumed to be a function of stress. The nominal stress-strain relations are derived in the form of differential equations which are solved numerically for off-axis specimens under uniaxial loading. A potential function and the associated effective stress and effective creep strain rates are introduced to simplify the orthotropic relations.

Chen, J. L.↗

Nonmethane hydrocarbon chemistry in the remote marine boundary layer

A photochemical model of the remote marine boundary layer (MBL) is presented, with focus placed on the role of reactive nonmethane hydrocarbons (NMHC). A wide range of NMHC air-sea fluxes with various relative distributions of NMHC regions are considered. In particular, the flux magnitude at which NMHC emissions become significant, and then dominant, players in MBL chemistry is identified. Emphasis is placed on diurnal variability, diurnal ozone variations and sensitivity to NMHC emission fluxes, to CO, O3, H2O, and UV light, and to kinetics and isometric composition. Model runs indicate that, in the range consistent with current observations, the NMHCs may either dominate MBL chemistry, or simply be contributors at the 10-percent level. These model runs also show that existing observations of NMHCs in ocean water find them to scarce for fluxes from bulk-flux air-sea gas exchange models to be consistent with the fluxes needed in the proposed model to maintain the lowest observed MBL NMHC.

Donahue, Neil M.↗

Improved interfacial li-ion transport in composite polymer electrolytes via surface modification of LLZO

Composite polymer electrolytes that incorporate ceramic fillers in a polymer matrix offer mechanical strength and flexibility as solid electrolytes for lithium metal batteries. However, fast Li + transport between polymer and Li + -conductive filler phases is not a simple achievement due to high barriers for Li + exchange across the interphase. This study demonstrates how modification of Li 7 La 3 Zr 2 O 12 (LLZO) nanofiller surfaces with silane chemistries influences Li + transport at local and global electrolyte scales. Anhydrous reactions covalently link amine-functionalized silanes [(3-aminopropyl)triethoxysilane (APTES)] to LLZO nanoparticles, which protects LLZO in air. APTES functionalization lowers the poly (ethylene oxide) (PEO)-LLZO interphase resistance to half that of unmodified LLZO and increases effective Li + transference number, while insulating Al 2 O 3 completely blocks ion exchange and lowers transference number and conductivity in PEO-lithium bis(trifluoromethanesulfonyl)imide (LiTFSI)-LLZO composites. Modeling an inner resistive interphase between LLZO and PEO surrounded by an outer conductive interphase explains non-linear conductivity trends. Solid-state 7 Li & 6 Li nuclear magnetic resonance shows Li + only exchanges between PEO-LiTFSI and some LLZO interphase, with no appreciable Li + transport through bulk LLZO. Surface functionalization is a promising path toward lowering the polymer-ceramic interphase resistance. This work demonstrates that local changes in Li + transport affect macroscopic performance, highlighting the intricate relationships between all interfaces in inherently heterogeneous composite polymer electrolytes.

25 ENERGY STORAGE↗

A Shell/3D Modeling Technique for Delaminations in Composite Laminates

A shell/3D modeling technique was developed for which a local solid finite element model is used only in the immediate vicinity of the delamination front. The goal was to combine the accuracy of the full three-dimensional solution with the computational efficiency of a plate or shell finite element model. Multi-point constraints provide a kinematically compatible interface between the local 3D model and the global structural model which has been meshed with plate or shell finite elements. For simple double cantilever beam (DCB), end notched flexure (ENF), and single leg bending (SLB) specimens, mixed mode energy release rate distributions were computed across the width from nonlinear finite element analyses using the virtual crack closure technique. The analyses served to test the accuracy of the shell/3D technique for the pure mode I case (DCB), mode II case (ENF) and a mixed mode I/II case (SLB). Specimens with a unidirectional layup where the delamination is located between two 0 plies, as well as a multidirectional layup where the delamination is located between two non-zero degree plies, were simulated. For a local 3D model extending to a minimum of about three specimen thicknesses in front of and behind the delamination front, the results were in good agreement with mixed mode strain energy release rates obtained from computations where the entire specimen had been modeled with solid elements. For large built-up composite structures modeled with plate elements, the shell/3D modeling technique offers a great potential, since only a relatively small section in the vicinity of the delamination front needs to be modeled with solid elements.

Krueger, Ronald↗

Composite fracture using the normal stress ratio theory

Experiments were performed on center-notched specimens of unidirectional AS4/3501-6 graphite-epoxy subjected to varied far-field loading conditions. Specimens tested include on-axis and off-axis tensile coupons and Iosipescu shear specimens, allowing the study of normal, shear, and combined far-field loading conditions. From the tests, the direction of initial crack extension and the far-field stresses at crack initiation and fracture are obtained. The normal stress ratio theory is subsequently applied within a macroscopic-level homogeneous anisotropic stress solution to predict the direction of crack extension, critical stresses, and other crack growth behavior for the experiments. Comparison between experimental results and theoretical predictions indicates that application of the normal stress ratio theory within a macroscopic-level stress analysis has significant potential as an approach to modeling composite crack growth.

Beuth, J. L., Jr.↗

High Pressure Melting Curve of Fe‐Si: Implication for the Thermal Properties in Mercury's Core

The motion of liquid iron (Fe) alloy materials in the outer core drives the dynamo, which generates Mercury's magnetic field. The assessment of core models requires laboratory measurements of the melting temperature of Fe alloys at high pressure. Here, we experimentally determined the melting curve of Fe9wt%Si and Fe17wt%Si up to 17 GPa using in situ and ex situ measurements of intermetallic fast diffusion that serves as the melting criterion in a large-volume press. Our determined melting slopes are comparable with previous studies up to about 17 GPa. However, when extrapolated, our melting slopes significantly deviate from previous studies at higher pressures. For Mercury's core with a model composition of Fe9wt%Si, the melting temperature-depth profile determined in our study is lower by ∼150–250 K when compared with theoretical calculations. Using the new melting curve of Fe9wt%Si and the electrical resistivity values from a previous study of Fe8.5wt%Si, we estimate that the electronic thermal conductivity of liquid Fe9wt%Si is 30 Wm −1 K −1 at the Mercury's CMB pressure of 5 GPa and 37 Wm −1 K −1 at an assumed ICB of 21 GPa, corresponding to heat flux values of 23 mWm −2 and 32 mWm −2 , respectively. These values provide new constraints on the core models.

58 GEOSCIENCES↗

Global Geodetic Signatures of the Antarctic Ice Sheet

Four scenarios of present day Antarctic ice sheet mass change are developed from comprehensive reviews of the available glaciological and oceanographic evidence. The gridded scenarios predict widely varying contributions to secular sea level change xi ranging from -1.1 to 0.45 mm/yr, and predict polar motion m and time-varying low-degree gravitational coefficients J(sub l) that differ significantly from earlier estimates. A reasonably linear relationship between the rate of sea level change from Antarctica xi(sub A) and the predicted Antarctic J(sub l) is found for the four scenarios. This linearity permits a series of forward models to be constructed that incorporate the effects of ice mass changes in Antarctica, Greenland, and distributed smaller glaciers, as well as postglacial rebound (assuming the ICE-3G deglaciation history), with the goal of obtaining optimum reconciliation between observed constraints on J(sub l) and sea level rise xi. Numerous viable combinations of lower mantle viscosity and hydrologic sources are found that safely "observed" in the range of 1 to 2-2.5 mm/yr and observed J(sub l) for degrees 2, 3, and 4. In contrast, rates of global sea level rise above 2.5 mm/yr are inconsistent with available J(sub l) observations. The successful composite models feature a pair of lower mantle viscosity solutions arising from the sensitivity of J(sub l) to glacial rebound. The paired values are well separated at xi = 1 mm/yr, but move closer together as xi is increased, and, in fact, merge around xi = 2 - 2.5 mm/yr, revealing an intimate relation between xi and preferred lower mantle viscosity. This general pattern is quite robust and persists for different J(sub l) solutions, for variations in source assumptions, and for different styles of lower mantle viscosity stratification. Tighter J(sub l) constraints for l greater than 2 may allow some viscosity stratification schemes and source assumptions to be excluded in the future. For a given total observed xi the sea level rise from Antarctica xi(sub A) is tightly constrained and ranges from 0 to + 1 mr/yr (corresponding to an ablating ice sheet) as estimates of xi are raised from 1 to 2.5 mm/yr. However, when the degree 3 zonal harmonic constraint is removed, the solutions show little sensitivity to Antarctic mass balance, emphasizing the need for a well determined odd-degree secular zonal harmonic for determining polar ice mass balance.

James, Thomas S.↗

Development of a Rubber-Based Product Using a Mixture Experiment: A Challenging Case Study

Many products used in daily life are made by blending two or more components. The properties of such products typically depend on the relative proportions of the components. Experimental design, modeling, and data analysis methods for mixture experiments provide for efficiently determining the component proportions that will yield a product with desired properties. This article presents a case study of the work performed to develop a new rubber formulation for an o-ring (a circular gasket) with requirements specified on 10 product properties. Each step of the study is discussed, including: 1) identifying the objective of the study and requirements for properties of the o-ring, 2) selecting the components to vary and specifying the component constraints, 3) constructing a mixture experiment design, 4) measuring the responses and assessing the data, 5) developing property-composition models, 6) selecting the new product formulation, and 7) confirming the selected formulation in manufacturing. The case study includes some challenging and new aspects, which are discussed in the article.

Full-cubic mixture model↗

A composition and thickness model for lunar impact crater and basin deposits

A model was developed for determining the composition and thicknesses of continuous deposits of large lunar impact craters and basins. Results of a photogeologic study and topographic analysis of continuous deposits of the lunar crater Delisle, which support the model, show that numerous secondary craters and crater chains with concentric dunes on their uprange rims occur well within the mapped deposits. At any given radius from Delisle, the secondary craters are of equal freshness and not vastly different in size, but with increasing radius they become better defined. The upper surface of the continuous deposits reaches the level of the pre-existing mare surface well within the mapped boundary of the deposits.

Morrison, R. H.↗

Toward the Development of a Fundamentally Based Chemical Model for Cyclopentanone: High-Pressure-Limit Rate Constants for H Atom Abstraction and Fuel Radical Decomposition

Theoretical aspects of the development of a chemical kinetic model for the pyrolysis and combustion of a cyclic ketone, cyclopentanone, are considered. We present calculated thermodynamic and kinetic data for the first time for the principal species including 2- and 3-oxo-cyclopentyl radicals, which are in reasonable agreement with the literature. Furthermore, these radicals can be formed via H atom abstraction reactions by H and Ö atoms and OH, HO 2 , and CH 3 radicals, the rate constants of which have been calculated. Abstraction from the β-hydrogen atom is the dominant process when OH is involved, but the reverse holds true for HO 2 radicals. We also determined the subsequent β-scission of the radicals formed, and it is shown that recent tunable VUV photoionization mass spectrometry experiments can be interpreted in this light. The bulk of the calculations used the composite model chemistry G4, which was benchmarked in the simplest case with a coupled cluster treatment, CCSD(T), in the complete basis set limit.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Multiscale Modeling of Ceramic Matrix Composites

Results of multiscale modeling simulations of the nonlinear response of SiC/SiC ceramic matrix composites are reported, wherein the microstructure of the ceramic matrix is captured. This micro scale architecture, which contains free Si material as well as the SiC ceramic, is responsible for residual stresses that play an important role in the subsequent thermo-mechanical behavior of the SiC/SiC composite. Using the novel Multiscale Generalized Method of Cells recursive micromechanics theory, the microstructure of the matrix, as well as the microstructure of the composite (fiber and matrix) can be captured.

ceramic matrix composites↗

Waste Glass Property Database and Data Qualification Plan

The U.S. Department of Energy vitrification facilities currently employ glass property-composition models to ensure processability of waste streams and acceptability of the final glass products. The efficiency (e.g., process flexibility, reduced down time, and reduced cost) of waste vitrification is directly tied to the size of the waste processing envelope. To increase the size of the processing envelope, a larger database was developed. The database will improve the prediction accuracy of glass properties, leading to higher waste loadings, higher waste throughput, broader process flexibility, and reduced mission life. Some of the data in this expanded database has insufficient quality assurance (QA) rigor for use in nuclear facility operation and glass qualification. A plan was developed to qualify high-value data under the appropriate QA rigor for use at the Hanford Waste Treatment and Immobilization Plant’s High-Level Waste (HLW) Facility. This qualification plan first compares the data coverage over the composition region of interest to Hanford HLW by property. Those studies that significantly improve the data coverage in the Hanford HLW glass composition region were prioritized for qualification activities. Qualification methods including QA equivalence, peer review, data corroboration, and confirmatory testing have been assigned to each high-priority study.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗