SEARCH · Search NASA

Results for “Compositional simulation”

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.

SourceResearch connectionFrom yearThrough year

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

At least 847 records · Page 47

All sources NASA NTRS DOE OSTI

DOE OSTI2015

Effective thermal conductivity of metal and non-metal particulate composites with interfacial thermal resistance at high volume fraction of nano to macro-sized spheres

In this study, we propose a theoretical model to compute the effective thermal conductivity of metal and dielectric spherical particle reinforced composites with interfacial thermal resistance. We consider a wide range of filler volume fraction with sizes ranging from nano- to macro-scale. The model, based on the differential effective medium theory, accounts for particle interactions through two sets of volume fraction corrections. The first correction accounts for a finite volume of composite and the second correction introduces a self-crowding factor that allows us to develop an accurate model for particle interaction even for high volume fraction of fillers. The model is examined to other published models, experiments, and numerical simulations for different types of composites. We observe an excellent agreement between the model and published datasets over a wide range of particle volume fractions and material properties of the composite constituents.

36 MATERIALS SCIENCE↗

NASA NTRS1987

Mechanical and electrical properties of graphite fiber-epoxy composites made from pristine and bromine intercalated fibers

The mechanical and electrical properties of pristine and bromine intercalated graphite fiber-epoxy composites were compared. The two types of composite were similar in terms of tensile modulus, tensile strength, and Poisson's Ratio. However, the interlaminar shear strength of the brominated composite was 18 percent greater than its pristine counterpart. Only slight differences were observed in flexural properties. A five-fold decrease was observed in the electrical resistivity of the brominated composite parallel to the axis of the fibers, resulting in a unidirectional resistivity of about 90 microOmega/cm. Transverse resistivity was unaffected. Both types of composite were subjected to a simulated lightning strike of 10 KJ (at a peak current of 150 kA), and the composite with the intercalated graphite exhibited less damage.

Jaworske, Donald A.↗

NASA NTRS2004

Testing of Composite Fan Vanes With Erosion-Resistant Coating Accelerated

The high-cycle fatigue of composite stator vanes provided an accelerated life-state prior to insertion in a test stand engine. The accelerated testing was performed in the Structural Dynamics Laboratory at the NASA Glenn Research Center under the guidance of Structural Mechanics and Dynamics Branch personnel. Previous research on fixturing and test procedures developed at Glenn determined that engine vibratory conditions could be simulated for polymer matrix composite vanes by using the excitation of a combined slip table and electrodynamic shaker in Glenn's Structural Dynamics Laboratory. Bench-top testing gave researchers the confidence to test the coated vanes in a full-scale engine test.

Bowman, Cheryl L.↗

NASA NTRS2006

Micro-Energy Rates for Damage Tolerance and Durability of Composite Structures

In this paper, the adhesive bond strength of lap-jointed graphite/aluminum composites is examined by computational simulation. Computed micro-stress level energy release rates are used to identify the damage mechanisms associated with the corresponding acoustic emission (AE) signals. Computed damage regions are similarly correlated with ultrasonically scanned damage regions. Results show that computational simulation can be used with suitable NDE methods for credible in-service monitoring of composites.

Chamis, Christos C.↗

NASA NTRS2013

Skipped Stage Modeling and Testing of the CPAS Main Parachutes

The Capsule Parachute Assembly System (CPAS) has undergone the transition from modeling a skipped stage event using a simulation that treats a cluster of parachutes as a single composite canopy to the capability of simulating each parachute individually. This capability along with data obtained from skipped stage flight tests has been crucial in modeling the behavior of a skipping canopy as well as the crowding effect on non-skipping ("lagging") neighbors. For the finite mass inflation of CPAS Main parachutes, the cluster is assumed to inflate nominally through the nominal fill time, at which point the skipping parachute continues inflating. This sub-phase modeling method was used to reconstruct three flight tests involving skipped stages. Best fit inflation parameters were determined for both the skipping and lagging canopies.

Varela, Jose G.↗

NASA NTRS2018

An Engineering Approach to Analyze Damage Initiation Modes in Tapered Composite Structures

In this work an engineering approach is demonstrated for analyzing damage initiation modes in tapered composite structures. The analysis methodology includes simulation of the non-linear static response of tapered composite structures under static tension loads to predict the location of interfacial delamination initiation. Furthermore, the developed methodology provides a strength-based criterion to assess whether damage initiation will occur in the inter-laminar delamination or intra-laminar matrix cracking mode. Based on the results of the analysis, a tapered composite structure is fabricated and tested under displacement-controlled quasi-static tension loading. The damage initiation location captured experimentally is compared with the analysis towards achieving preliminary qualitative validation. The linear stiffness of the tapered composite structure is predicted within 15% of the experimental average thereby achieving preliminary quantitative validation.

Rao, Prabhakar M.↗

NASA NTRSDate not supplied

An Engineering Approach to Analyze Damage Initiation Modes in Tapered Composite Structures

Prabhakar M Rao↗

DOE OSTI2025

Modeling Equilibrium Solid–Liquid Interfaces under Effective Constant Chemical Potential Using Machine Learning Interatomic Potentials

The chemical potential (μ) of species in solution is essential for understanding various chemical processes at interfaces. Molecular dynamics (MD) simulations, constrained by fixed compositions, cannot maintain constant chemical potential with reference to a targeted concentration or chemical potential under nonequilibrium or dynamic conditions, as solute species can migrate to the interface and deplete (or enrich) the bulk due to solute-interface interactions. In this study, we introduce a simple and computationally efficient approach named iterative quasi-constant chemical potential molecular dynamics (iqCμMD) simulation, which helps simulate targeted molar concentrations of species in solution. iqCμMD overcomes the limitations of conventional MD by adjusting the number of species in the solution to reach a target bulk concentration (chemical potential), which allows simulation of the interface under the bulk conditions comparable to experiment. We demonstrate our approach using machine learning interatomic potential (MLIP)-based MD simulations of the Na 2 SO 4,aq –graphene interface, and to show the transferability of our approach, we also perform classical force field-based MD simulations of NaCl aq –air and NaCl aq –graphite interfaces, which produce comparable results to previous CμMD simulations. Our results also show that the iqCμMD approach efficiently achieves the desired bulk ion concentration within two iterations, and by utilizing MLIPs, we can achieve converged results using relatively small-scale simulations compared to previous CμMD simulations. By combining iqCμMD with MLIP-driven simulations, solid–liquid interfaces can be modeled under an effective constant chemical potential with DFT-level accuracy. Here, we show that iqCμMD offers a robust and simple computational framework for constant chemical potential simulations, as its only requirement is to be able to converge interfacial simulations with a measurable bulk region.

Chemical structure↗

DOE OSTI2020

Unlocking multiphysics design guidelines on Si/C composite nanostructures for high-energy-density and robust lithium-ion battery anode

In general, current material fabrication guidance for novel designs of Si/C composite particle materials focuses on electrochemical behavior and redox reactions at the nano/micro level. However, such guidance cannot provide detailed information for predicting mechanical deformations of the composite particles, especially when the mechanical field coupled with electrochemical and thermal fields. Here, we establish an electro-chemo-mechanical model and implement it to quantitatively analyze the multiphysics behavior of five representative Si/C composite nanostructures. Modeling and computation discover that yolk-shell and dual-shell structures are more robust in terms of particle fractures. When considering electrochemical performance, the yolk-shell structure is the best among the compared five Si/C composites. Finally, we map design guidance to further illustrate quantitative structure-property relationships. This study provides novel insights on Si/C composite nanostructure anode material design and additional powerful design tools for next-generation high-energy-density lithium-ion batteries.

25 ENERGY STORAGE↗

NASA NTRS2010

Modification of a Macromechanical Finite-Element Based Model for Impact Analysis of Triaxially-Braided Composites

A macro level finite element-based model has been developed to simulate the mechanical and impact response of triaxially-braided polymer matrix composites. In the analytical model, the triaxial braid architecture is simulated by using four parallel shell elements, each of which is modeled as a laminated composite. For the current analytical approach, each shell element is considered to be a smeared homogeneous material. The commercial transient dynamic finite element code LS-DYNA is used to conduct the simulations, and a continuum damage mechanics model internal to LS-DYNA is used as the material constitutive model. The constitutive model requires stiffness and strength properties of an equivalent unidirectional composite. Simplified micromechanics methods are used to determine the equivalent stiffness properties, and results from coupon level tests on the braided composite are utilized to back out the required strength properties. Simulations of quasi-static coupon tests of several representative braided composites are conducted to demonstrate the correlation of the model. Impact simulations of a represented braided composites are conducted to demonstrate the capability of the model to predict the penetration velocity and damage patterns obtained experimentally.

Goldberg, Robert K.↗

NASA NTRS2015

Development and Characterization of a Rate-Dependent Three-Dimensional Macroscopic Plasticity Model Suitable for Use in Composite Impact Problems

Several key capabilities have been identified by the aerospace community as lacking in the material/models for composite materials currently available within commercial transient dynamic finite element codes such as LS-DYNA. Some of the specific desired features that have been identified include the incorporation of both plasticity and damage within the material model, the capability of using the material model to analyze the response of both three-dimensional solid elements and two dimensional shell elements, and the ability to simulate the response of composites composed with a variety of composite architectures, including laminates, weaves and braids. In addition, a need has been expressed to have a material model that utilizes tabulated experimentally based input to define the evolution of plasticity and damage as opposed to utilizing discrete input parameters (such as modulus and strength) and analytical functions based on curve fitting. To begin to address these needs, an orthotropic macroscopic plasticity based model suitable for implementation within LS-DYNA has been developed. Specifically, the Tsai-Wu composite failure model has been generalized and extended to a strain-hardening based orthotropic plasticity model with a non-associative flow rule. The coefficients in the yield function are determined based on tabulated stress-strain curves in the various normal and shear directions, along with selected off-axis curves. Incorporating rate dependence into the yield function is achieved by using a series of tabluated input curves, each at a different constant strain rate. The non-associative flow-rule is used to compute the evolution of the effective plastic strain. Systematic procedures have been developed to determine the values of the various coefficients in the yield function and the flow rule based on the tabulated input data. An algorithm based on the radial return method has been developed to facilitate the numerical implementation of the material model. The presented paper will present in detail the development of the orthotropic plasticity model and the procedures used to obtain the required material parameters. Methods in which a combination of actual testing and selective numerical testing can be combined to yield the appropriate input data for the model will be described. A specific laminated polymer matrix composite will be examined to demonstrate the application of the model.

Impact↗

NASA NTRS1993

On predicting contamination levels of HALOE optics aboard UARS using direct simulation Monte Carlo

A three-dimensional version of the direct simulation Monte Carlo method is adapted to assess the contamination environment surrounding a highly detailed model of the Upper Atmosphere Research Satellite. Emphasis is placed on simulating a realistic, worst-case set of flowfield and surface conditions and geometric orientations in order to estimate an upper limit for the cumulative level of volatile organic molecular deposits at the aperture of the Halogen Occultation Experiment. Problems resolving species outgassing and vent flux rates that varied over many orders of magnitude were handled using species weighting factors. Results relating to contaminant cloud structure, cloud composition, and statistics of simulated molecules impinging on the target surface are presented, along with data related to code performance. Using procedures developed in standard contamination analyses, the cumulative level of volatile organic deposits on HALOE's aperture over the instrument's 35-month nominal data collection period is estimated to be about 2700A.

Woronowicz, Michael S.↗

NASA NTRS2015

Using Neutron Spectroscopy to Constrain the Composition and Provenance of Phobos and Deimos

The origin of the Martian moons Phobos and Deimos is obscure and enigmatic. Hypotheses include the capture of asteroids originally from the outer main belt or beyond, residual material left over from Mars' formation, and accreted ejecta from a large impact on Mars, among others. Measurements of reflectance spectra indicate a similarity to dark, red D-type asteroids, but could indicate a highly space-weathered veneer. Here we suggest a way of constraining the near-surface composition of the two moons, for comparison to known meteoritic compositions. Neutron spectroscopy, particularly the thermal and epithermal neutron flux, distinguishes clearly between various classes of meteorites and varying hydrogen (water) abundances. Perhaps most surprising of all, a rendezvous with Phobos or Deimos is not necessary to achieve this. A low-cost mission based on the LADEE spacecraft design in an eccentric orbit around Mars can encounter Phobos every 2 weeks. As few as five flyby encounters at speeds of 2.3 kilometers per second and closest-approach distance of 3 kilometers provide sufficient data to distinguish between ordinary chondrite, water-bearing carbonaceous chondrite, ureilite, Mars surface, and aubrite compositions. A one-Earth year mission design includes many more flybys at lower speeds and closer approach distances, as well as similar multiple flybys at Deimos in the second mission phase, as described in the Phobos And Deimos Mars Environment (PADME) mission concept. This presentation will describe the expected thermal and epithermal neutron fluxes based on MCNP6 (Monte Carlo N (i.e. Neutron)-Particle transport code (version 6) simulations of different meteorite compositions and their uncertainties.

thtermal and epithermal neutron fluxes↗

DOE OSTI2026

Global energy balance simulations of shattered pellet injection

The dynamics of shattered pellet injection (SPI) shutdowns are simulated using a time-dependent global energy balance model, based on a modification of the KPRAD framework. The new SPI particle source in the model calculates the ablation of individual pellet fragments that enter the plasma as a temporally resolved plume, thus capturing the effects of earlier fragments on the ablation of those that follow, which has a significant impact on the overall assimilation. Despite the reduced physics and the global averaging of all quantities, results from a large number of DIII-D, KSTAR, and JET experiments are well reproduced, including the plasma cooling timescales, particle assimilation, and current quench (CQ) rates. Cooling timescales and CQ rates are in good agreement for pellets containing as little as ∼1% neon, while particle assimilations are most accurate for neon fractions above ∼15% by number of atoms. Below this, the assimilation tends to be overestimated due to the lack of radial transport in the particle balance, which becomes important in the low-Z limit. Predictive simulations of mixed-composition dual-SPI shutdowns in ITER are compared against those with the 3D non-linear magnetohydrodynamic code JOREK and are found to reproduce overall trends observed in the higher-fidelity modeling across a range of injection scenarios. The general success of the model points to the critical role of energy balance in determining SPI particle assimilation and the subsequent disruption dynamics and highlights the value of these simulations for experimental interpretation and for optimizing the deployment of computationally expensive, higher-fidelity models.

energy balance model↗

DOE OSTI2018

Liquidus temperature in the spinel primary phase field: A comparison between optical and crystal fraction methods

Liquidus temperature (TL) was measured for simulated high-level waste borosilicate glasses covering a Hanford composition region, using an optical method and a crystal-fraction extrapolation method with X-ray diffraction data from isothermal heat treatments. The 38 glasses encompassed a one-component-at-a-time variation of a 16-component matrix. The TL values ranged from 1006°C to 1603°C. First-order polynomial models were fitted to data to obtain coefficients in terms of the effect of 1 mass% component addition on the TL: TL-increasing oxides Cr2O3 (264°C), “Others” (minor components, 163°C), oxides of noble metals (137°C), NiO (91°C), Al2O3 and Fe2O3 (~19–21°C); TL-decreasing oxides K2O (-26°C), Na2O (-41°C), and Li2O (-68°C); oxides of little effect MnO, P2O5, ZrO2, F, Bi2O3, SiO2, B2O3, and CaO (9 to -12°C). Also presented are temperatures (T1%) at which glasses contain 1 vol% of spinel as these values are considered relevant to the Hanford Tank Waste Treatment and Immobilization Plant. The measured and estimated values are compared and contrasted and the effect of TL and T1% on glass formulation is discussed.

crystal fraction↗

DOE OSTI2022

Identification of Distribution of Elastic Constants of SiC/SiC Braided Tubes

There has been considerable interest in chemical vapor infiltration (CVI) manufactured silicon carbide fiber and silicon carbide matrix (SiCf/SiCm) composite tubes due to their superior mechanical properties. There are many parameters in the manufacturing of SiCf/SiCm composite tubes such as elastic constants of constituent SiC fiber and SiC matrix, braiding angle, porosity, and others whose values vary over a range. A sensitivity analysis will be carried out to determine the dominant parameters among the abovementioned ones. We intend to determine the variability in elastic constants of SiCf/SiCm composite tubes considering the variability in dominant parameters. A chosen number of combinations of dominant variables will be sampled through Latin hypercube sampling (LHS) method. The elastic constants of composite tube namely Young's moduli in circumferential and longitudinal direction, Poisson's ratio and shear modulus will be calculated at sampled points through finite element (FE) analysis. A polynomial response surface (PRS) will be built for each of the elastic constants to be used a surrogate to FE analysis. once the predictive accuracy of PRS is verified at the validation set of data points, the PRS will be invoked 105 times in Monte Carlo simulations (MCS). The variability can be calculated based on MCS-generated data which will be presented at the conference.

42 ENGINEERING↗

DOE OSTI2019

Preliminary Report on the Study of Beam-Induced Background Effects at a Muon Collider

Physics at a multi-TeV muon collider needs a change of perspective for the detector design due to the large amount of background induced by muon beam decays. Preliminary studies, based on simulated data, on the composition and the characteristics of the particles originated from the muon decays and reaching the detectors are presented here. The reconstruction performance of the physics processes $H\to b\bar b$ and $Z\to b\bar b$ has been investigated for the time being without the effect of the machine induced background. A preliminary study of the environment hazard due to the radiation induced by neutrino interactions with the matter is presented using the FLUKA simulation program.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS2015

Incorporation of Plasticity and Damage Into an Orthotropic Three-Dimensional Model with Tabulated Input Suitable for Use in Composite Impact Problems

The need for accurate material models to simulate the deformation, damage and failure of polymer matrix composites under impact conditions is becoming critical as these materials are gaining increased usage in the aerospace and automotive industries. While there are several composite material models currently available within commercial transient dynamic finite element codes, several features have been identified as being lacking in the currently available material models that could substantially enhance the predictive capability of the impact simulations. A specific desired feature pertains to the incorporation of both plasticity and damage within the material model. Another desired feature relates to using experimentally based tabulated stress-strain input to define the evolution of plasticity and damage as opposed to specifying discrete input properties (such as modulus and strength) and employing analytical functions to track the response of the material. To begin to address these needs, a combined plasticity and damage model suitable for use with both solid and shell elements is being developed for implementation within the commercial code LS-DYNA. The plasticity model is based on extending the Tsai-Wu composite failure model into a strain-hardening based orthotropic plasticity model with a non-associative flow rule. The evolution of the yield surface is determined based on tabulated stress-strain curves in the various normal and shear directions and is tracked using the effective plastic strain. The effective plastic strain is computed by using the non-associative flow rule in combination with appropriate numerical methods. To compute the evolution of damage, a strain equivalent semi-coupled formulation is used, in which a load in one direction results in a stiffness reduction in multiple coordinate directions. A specific laminated composite is examined to demonstrate the process of characterizing and analyzing the response of a composite using the developed model.

Impact↗