Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Material Point Method”

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 523 records · Page 29

Physical and Biological Properties of Cobalt- and Copper-Doped Calcium Phosphates as Bone Substitute Materials

Arthritis, osteoporosis, and other bone diseases and defects are common medical issues worldwide. By utilizing ceramic biomaterials as bone substitutes to treat these diseases, bone regeneration can be promoted and toxicity from the substitute material can be limited. Calcium phosphate (CaP) ceramics have become popular as bone substitutes due to their biocompatibility and similarity in composition to natural bone. The purpose of this research was to investigate the physical properties and biological responses of CaP-based bone substitutes, doped with metal ions. Some metal ions are present in natural bone in small amounts, and recent studies show they affect the biological responses of CaPs significantly. Dopants such as magnesium (Mg), strontium (Sr), silicon (Si) and iron (Fe) alter the phase composition, mechanical properties, osteogenesis and angiogenesis properties of CaPs, depending on their concentration and method of addition. In the present study, the effects of cobalt (Co) and copper (Cu) on the physical and biological properties of two main CaP materials, brushite cement (BrC) and tricalcium phosphate (TCP), have been investigated. Different concentrations of dopants in forms of oxide and/or chloride were selected and their effect on phase composition, sintering behavior, density, setting time, compressive strength and in vitro interaction with osteosarcoma and osteoblast cells were studied. Different techniques of sample preparation, challenges during cement preparation and compact sample pressing, sintering, and methods of dopant addition are discussed. The presence of Cu in tricalcium phosphate was found to increase thermal stability of the material and decrease the impact of sintering temperature on porosity. In addition, Cu caused a reduction in expression of inflammatory gene markers by human osteoblast cells and an increased expression at early time points due of osteoinductive markers.The addition of Cu and Co to brushite cement caused an increase in thermal stability of the material, and the addition of Cu caused significant increase in setting time. Small dopant amounts of Co caused decrease in setting time, but higher amounts caused an increase. In addition, Cu dopant in small amounts resulted in an increase of compressive strength. Both Cu and Co led to a decrease in expression of inflammatory gene markers by osteoblast and osteosarcoma cells respectively, and Cu also caused an increase in osteoinductive marker expression. Thus, these results indicate that Cu and Co have a positive impact on the physical, mechanical and biological properties of calcium phosphate biomaterials.

Cummings, Haley V.↗

Constrained model calibration of grain structure dependent spall dynamics in shock-loaded tantalum

Here, we perform a gas gun experiment by shock loading tantalum samples of varying grain structures to assess the suitability of a numerical model for simulating spall behavior. The observed differences in spall strength, as well spallation and re-compression history, are not captured in uncalibrated hydrodynamic simulations. An optimization is performed on the Johnson spall model to determine the best parameters that fit the observed trends. Linear stability analysis is employed to motivate bounds on those parameters. Herein, optimized simulations agree well with the experimental results, reproducing pullback depth and recompression timescales across the different samples tested. Further, the observed pullback time of the single crystal sample was found to imply, via the stability analysis, a percolation threshold in good agreement with the theoretical value for a body centered cubic lattice. Therefore, the combined linear stability and percolation analysis shows promise and may be applied to other materials with diverse microstructures. Collectively, the findings demonstrate that the model is suitable for reproducing spall-induced free surface behavior across various microstructures, but also points to caution in using model coefficients for uncalibrated microstructures.

58 GEOSCIENCES↗

Apparatus and method for heating a material in a transparent ampoule

An improved process for heating a material within a fused silica ampoule by radiation through the wall of the ampoule, while simultaneously passing a cooling gas around the ampoule is described. The radiation passes through a screen of fused silica so as to remove those components capable of directly heating the silica, therby increasing the temperature of the material within the ampoule above the strain point of the ampoule, while maintaining the exterior of the ampoule cool enough to prevent rupturing the amp.

Holland, L. R.↗

Emulation of radiation transport in 3D stochastic media using 1D planar Monte Carlo stochastic media radiation transport algorithms

A subset of stochastic media radiation transport problems involves those in which radiation is incident on a thin slab of stochastic material. Particle tracking in 3D for such problems is expensive, and 1D planar models lack accuracy because they only allow the material to change in one dimension. Therefore, we propose dimensional emulation, which through a slight modification allows existing 1D planar geometry stochastic media radiation transport models to reproduce results from the equivalent 3D models by allowing the material to change in all three dimensions, reproducing the fidelity of the 3D model for the low computational cost of the 1D planar model. In this work, we apply dimensional emulation to three Monte Carlo stochastic media radiation transport models: Chord Length Sampling (CLS), the Local Realization Preserving method (LRP), and a variant of Conditional Point Sampling (CoPS). For a common Markovian benchmark set, the 3D emulation variants of these algorithms are numerically verified to reproduce the results of the 3D variants within statistics while running 1.3 to 2 times faster in the implementation within Sandia National Laboratories open-source research code PlaybookMC. The 3D emulation variants are also shown to yield a 72%–92% reduction in error for the thin slab problems in comparison to the 1D benchmark. As a result, the 3D emulation variant of CLS and CoPS-1 are shown to reproduce 3D CLS results that were used to approximate results for a 3D spherical inclusion geometry benchmark set.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Measurement of crack length in width tapered beam experiments

The width tapered beam method for measuring fracture toughness has been contemplated for use in measuring the adhesion of photovoltaic materials, is promoted by it being a viable method for use directly on modules rather than to engineered test specimens. However, precise determination of crack length is often difficult as failure often involves tendril formation, cavitation voids, and other plastic deformation at the crack front in addition to obfuscation by opaque materials. Because fracture energy varies as the square of crack length, imprecise crack tip length results in large measurement uncertainties. In this work, we develop a method to estimate the crack tip length using the relationship between the applied force and the displacement of the load frame. The crack tip is determined for every data point and can be defined from a calibration curve related to the amount of energy put into the beam, or from a theoretically idealized estimation. Once a calibration curve is set up and the equations programmed, the fracture energy as a function of position on the newly formed crack interface can be accurately estimated. The test additionally provides much easier analysis and insight into the joint failure through more precise knowledge of fracture toughness at each de-bonded area.

42 ENGINEERING↗

Integrated Process Model Utilization and Development for Inertial Fusion Energy

This CRADA between LLNS and Longview Fusion Energy Systems utilized a modernized LLNL Integrated Process Model (IPM) to assess technoeconomic viability and power plant configurations for Inertial Fusion Energy (IFE). During the collaboration period, a streamlined IPM model was produced by LLNL, which consolidated a combination of 3 models under one file. Additionally, obsolete parameters were removed and parameters available for trade space analysis were organized in a user-friendly fashion on the front interface. Gain scaling curves were discussed and analyzed in the framework of the IPM. Additional updates were made to account for cost scaling to today's dollars using an average inflation rate. Heat transfer material costing was updated using a ground up approach with public vendor data. Further levelized cost of electricity (LCOE) methods and models were discussed and reviewed for applicability to the IPM. Prior published literature on optical scaling was also discussed. Moreover, as part of this CRADA, the Participant developed a fusion technology development roadmap and point designs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Integrated Process Model Utilization and Development for Inertial Fusion Energy

This CRADA between LLNS and Longview Fusion Energy Systems utilized a modernized LLNL Integrated Process Model (IPM) to assess technoeconomic viability and power plant configurations for Inertial Fusion Energy (IFE). During the collaboration period, a streamlined IPM model was produced by LLNL, which consolidated a combination of 3 models under one file. Additionally, obsolete parameters were removed and parameters available for trade space analysis were organized in a user-friendly fashion on the front interface. Gain scaling curves were discussed and analyzed in the framework of the IPM. Additional updates were made to account for cost scaling to today's dollars using an average inflation rate. Heat transfer material costing was updated using a ground up approach with public vendor data. Further levelized cost of electricity (LCOE) methods and models were discussed and reviewed for applicability to the IPM. Prior published literature on optical scaling was also discussed. Moreover, as part of this CRADA, the Participant developed a fusion technology development roadmap and point designs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Factors influencing the ultrasonic stress wave factor evaluation of composite material structures

To demonstrate that the finite-element model can be used to investigate some of the factors influencing the ultrasonic stress wave evaluation of materials, a hypothetical case was studied in which classical vibration theory was used. Vibration analysis and experiments for the undamaged case were conducted on an isotropic aluminum plate and a unidirectional graphite-epoxy plate, using a point source to excite the plates. The finite-element solution correlated within eight percent with the exact method. The frequencies predicted by the finite-element model were observed in the experiments in both plates, although in the composite plate, additional frequencies were observed which could not be accounted for. Damaged isotropic plates were also considered. The effects of increasing damage severity with constant damage area, and increasing damage area with constant severity on the resonant frequencies were analyzed.

Rebello, C. J.↗

Strain energy release rate analysis of delamination in a tapered laminate subjected to tension load

A tapered composite laminate subjected to tension load was analyzed using the finite-element method. The glass/epoxy laminate has a (+ or - 45)sub 3 group of plies dropped in three distinct steps, each 20 ply-thicknesses apart, thus forming a taper angle of 5.71 degrees. Steep gradients of interlaminar normal and shear stress on a potential delamination interface suggest the existence of stress singularities at the points of material and geometric discontinuities created by the internal plydrops. The delamination was assumed to initiate at the thin end of the taper on a -45/+45 interface and the delamination growth was simulated in both directions, i.e., along the taper and into the thin region. The strain-energy-release rate for a delamination growing into the thin laminate consisted predominantly of mode I (opening) component. For a delamination growing along the tapered region, the strain-energy-release rate was initially all mode I, but the proportion of mode I decreased with increase in delamination size until eventually total G was all mode II. The total G for both delamination tips increased with increase in delamination size, indicating that a delamination initiating at the end of the taper will grow unstably along the taper and into the thin laminate simultaneously.

Salpekar, S. A.↗

Real-time measurement of materials properties at high temperatures by laser produced plasmas

Determination of elemental composition and thermophysical properties of materials at high temperatures, as visualized in the context of containerless materials processing in a microgravity environment, presents a variety of unusual requirements owing to the thermal hazards and interferences from electromagnetic control fields. In addition, such information is intended for process control applications and thus the measurements must be real time in nature. A new technique is described which was developed for real time, in-situ determination of the elemental composition of molten metallic alloys such as specialty steel. The technique is based on time-resolved spectroscopy of a laser produced plasma (LPP) plume resulting from the interaction of a giant laser pulse with a material target. The sensitivity and precision were demonstrated to be comparable to, or better than, the conventional methods of analysis which are applicable only to post-mortem specimens sampled from a molten metal pool. The LPP technique can be applied widely to other materials composition analysis applications. The LPP technique is extremely information rich and therefore provides opportunities for extracting other physical properties in addition to the materials composition. The case in point is that it is possible to determine thermophysical properties of the target materials at high temperatures by monitoring generation and transport of acoustic pulses as well as a number of other fluid-dynamic processes triggered by the LPP event. By manipulation of the scaling properties of the laser-matter interaction, many different kinds of flow events, ranging from shock waves to surface waves to flow induced instabilities, can be generated in a controllable manner. Time-resolved detection of these events can lead to such thermophysical quantities as volume and shear viscosities, thermal conductivity, specific heat, mass density, and others.

Kim, Yong W.↗

An Overview of Prognosis Health Management Research at GRC for Gas Turbine Engine Structures With Special Emphasis on Deformation and Damage Modeling

Herein a general, multimechanism, physics-based viscoelastoplastic model is presented in the context of an integrated diagnosis and prognosis methodology which is proposed for structural health monitoring, with particular applicability to gas turbine engine structures. In this methodology, diagnostics and prognostics will be linked through state awareness variable(s). Key technologies which comprise the proposed integrated approach include 1) diagnostic/detection methodology, 2) prognosis/lifing methodology, 3) diagnostic/prognosis linkage, 4) experimental validation and 5) material data information management system. A specific prognosis lifing methodology, experimental characterization and validation and data information management are the focal point of current activities being pursued within this integrated approach. The prognostic lifing methodology is based on an advanced multi-mechanism viscoelastoplastic model which accounts for both stiffness and/or strength reduction damage variables. Methods to characterize both the reversible and irreversible portions of the model are discussed. Once the multiscale model is validated the intent is to link it to appropriate diagnostic methods to provide a full-featured structural health monitoring system.

Arnold, Steven M.↗

Hercules Single-Stage Reusable Vehicle (HSRV) Operating Base

Conceptual design for the layout of lunar-planetary surface support systems remains an important area needing further master planning. This paper explores a structured approach to organize the layout of a Mars-based site equipped for routinely flying a human-scale reusable taxi system. The proposed Hercules Transportation System requires a surface support capability to sustain its routine, affordable, and dependable operation. The approach organizes a conceptual Hercules operating base through functional station sets. The station set approach will allow follow-on work to trade design approaches and consider technologies for more efficient flow of material, energy, and information at future Mars bases and settlements. The station set requirements at a Mars site point to specific capabilities needed. By drawing from specific Hercules design characteristics, the technology requirements for surface-based systems will come into greater focus. This paper begins a comprehensive process for documenting functional needs, architectural design methods, and analysis techniques necessary for follow-on concept studies.

Surface Support Architectur↗

Synthesis of Layered Gold Tellurides AuSbTe and Au 2 Te 3 and Their Semiconducting and Metallic Behavior

Previous studies on natural samples of pampaloite (AuSbTe) revealed the crystal structure of a potentially cleavable and/or exfoliable material, while studies on natural and synthetic montbrayite (Sb-containing Au 2 Te 3 ) claimed various chemical compositions for this low-symmetry compound. Few investigations of synthetic samples have been reported for both materials, leaving much of their chemical, thermal, and electronic characteristics unknown. Here, in this study, we investigate the stability, electronic properties, and synthesis of the gold antimony tellurides AuSbTe and Au 1.9 Sb 0.46 Te 2.64 (montbrayite). Differential thermal analysis and in situ powder X-ray diffraction revealed that AuSbTe is incongruently melting, while Au 1.9 Sb 0.46 Te 2.64 is congruently melting. Calculations of the band structures and four-point resistivity measurements showed that AuSbTe is a semiconductor and Au 1.9 Sb 0.46 Te 2.64 a metal. Various synthesis attempts confirmed the limited stable chemical composition of Au 1.9 Sb 0.46 Te 2.64 , identified successful methods to synthesize both compounds, and highlighted the challenges associated with single-crystal synthesis of AuSbTe.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

X-ray Mapping of Terrestrial and Extraterrestrial Materials Using the Electron Microprobe

Lunar samples returned from the Apollo program motivated development of the Bence-Albee algorithm for the rapid and accurate analysis of lunar materials, and established interlaboratory comparability through its common use. In the analysis of mineral and rock fragments it became necessary to combine micro- and macroscopic analysis by coupling electron-probe microanalysis (EPMA) with automated stage point counting. A coarse grid that included several thousand points was used, and initially wavelength-dispersive (WDS) and later energydispersive (EDS) data were acquired at discrete stage points using approx. 5 sec count times. A approx 50 micrometer beam diameter was used for WDS and up to 500 micrometer beam diameter for EDS analysis. Average analyses of discretely sampled phases were coupled with the point count data to calculate the bulk composition using matrix algebra. Use of a defocused beam resulted in a contribution from multiple phases to each analytical point, and the analytical data were deconvolved relative to end-member phase chemistry on the fly. Impressive agreement was obtained between WDS and EDS measurements as well as comparison with bulk chemistry obtained by other methods. In the 30 years since these methods were developed, significant improvements in EPMA automation and computer processing have taken place. Digital beam control allows routine collection of x-ray maps by EDS, and stage mapping for WDS is conducted continuously at slew speed and incrementally by sampling at discrete points. Digital pulse processing in EDS systems has significantly increased the throughput for EDS mapping, and the ongoing development of Si-drift detector systems promises mapping capabilities rivaling WDS systems. Spectrum imaging allows a data cube of EDS spectra to be acquired and sophisticated processing of the original data is possible using matrix algebra techniques. The study of lunar and meteoritic materials includes the need to conveniently: (1) Characterize the sample at microscopic and macroscopic scales with relatively high sensitivity, (2) Determine the modal abundance of minerals, and (3) Identify and relocate discrete features of interest in terms of size and chemistry. The coupled substitution of cations in minerals can result in significant variation in mineral chemistry, but at similar average Z, leading to poor backscattered-electron (BSE) contrast discrimination of mineralogy. It is necessary to discriminate phase chemistry at both the trace element level and the major element level. To date, the WDS of microprobe systems is preferred for mapping due to high throughput and the ability to obtain the necessary intensity to discriminate phases at both trace and major element concentrations. It is desirable to produce fully quantitative compositional maps of geological materials, which requires the acquisition of k-ratio maps that are background and dead-time corrected, and which have been corrected by phi(delta z> or an equivalent algorithm at each pixel. To date, turnkey systems do not allow the acquisition of k-ratio maps and the rigorous correction in this manner. X-ray maps of a chondrule from the Ourique meteorite, and a comb-layered xenolith from the San Francisco volcanic field, have been analyzed and processed to extract phase information. The Ourique meteorite presents a challenge due to relatively low BSE contrast, and has been studied using spectrum imaging. X-ray maps for Si, Mg, and FeK(alpha) were used to produce RGB images. The xenolith sample contains sector-zoned augite, olivine, plagioclase, and basaltic glass. X-ray maps were processed using Lispix and ImageJ software to produce mineral phase maps. The x-ray maps for Mg, Ca, and Ti were used with traceback to generate binary images that were converted to RGB images. These approaches are successful in discriminating phases, but it is desirable to achieve the methods that were used on lunar samples 30 years ago on current microprobe systems. Curnt research includes x-ray mapping analysis of the Dalgety Downs chondrite by micro x-ray fluorescence and spectrum imaging, in collaboration with Kenny Witherspoon of IXRF Systems and Dale Newbury of NIST.

Carpenter, P.↗

High-throughput combinatorial approach expedites the synthesis of a lead-free relaxor ferroelectric system

Developing novel lead-free ferroelectric materials is crucial for next-generation microelectronic technologies that are energy efficient and environment friendly. However, materials discovery and property optimization are typically time-consuming due to the limited throughput of traditional synthesis methods. In this work, we use a high-throughput combinatorial synthesis approach to fabricate lead-free ferroelectric superlattices and solid solutions of (Ba 0.7 Ca 0.3 )TiO 3 (BCT) and Ba(Zr 0.2 Ti 0.8 )O 3 (BZT) phases with continuous variation of composition and layer thickness. High-resolution x-ray diffraction (XRD) and analytical scanning transmission electron microscopy (STEM) demonstrate high film quality and well-controlled compositional gradients. Ferroelectric and dielectric property measurements identify the “optimal property point” achieved at the composition of 48BZT–52BCT. Displacement vector maps reveal that ferroelectric domain sizes are tunable by varying {BCT–BZT} N superlattice geometry. This high-throughput synthesis approach can be applied to many other material systems to expedite new materials discovery and properties optimization, allowing for the exploration of a large area of phase space within a single growth.

36 MATERIALS SCIENCE↗

Impact of the adhesion-force lever-arm “a” on the Rock ‘n’ Roll resuspension model and how to compute it from contact mechanics

The paper proposes a methodology to fill a gap in the context of mechanistic particle resuspension models related to the adhesion and lift force lever arms ( a ): to date, only one experimental observation exists from a centrifuge experiment with a specific particle/surface pair which found r P / a = 100 and assumed that the contact points were symmetric with respect to the particle center. However, the Reeks and Hall’s (2001) model is highly sensitive to this parameter, which therefore cannot be used as-is for different particle diameters or materials. This paper proposes a methodology to compute the adhesion and lift force lever arms ( a A and a L , respectively) based on Atomic Force Microcopy (AFM) surface roughness measurements and contact mechanics arguments. In practice, a spherical, perfectly smooth particle is brought in contact (computationally) with the measured surface roughness and is then rotated by the point of contact to determine the two or more asperities over which the particle would sit on the surface. A method to compute the asperity radius of curvature is devised to compute the asperity deformation and the correct partitioning of the total deformation between the particle and the asperity is presented. For the glass surface under consideration, the asperity deformation was negligible, i.e., it was below the van der Waals radius, which is the limit of applicability of contact mechanics. The Reeks and Hall’s (2001) model is then extended to account for the joint probability of a A and a L . Significant differences in particle resuspended after a day in common atmospheric conditions are found between the approach that uses the lever-arm distributions or their mean values. Also, important differences are found when comparing the lever-arm distribution and r P / a = 100. Hence, one should be cautious of using r P / a = 100 for surface/particle pairs differing from those in the original study. Finally, limits of applicability of the theories used are presented together with considerations regarding how to extend the method to non-spherical particles and different materials.

42 ENGINEERING↗

Molecular dynamics simulations of radiation cascade evolution near cellular dislocation structures in additively manufactured stainless steels

The dislocation dense cellular structures observed in additively manufactured (AM) stainless steels can allow these materials to have enhanced strength and ductility compared to conventionally manufactured materials. Effective design of radiation resistant materials often necessitates creating microstructures containing of a high density of dislocations that act as sinks for radiation induced defects. This work uses molecular dynamics simulations to study the impact of pre-existing defects on radiation damage in stainless steel 316L fabricated by the laser powder bed fusion process. The evolution of synthetically generated dislocations in a crystal system in response to multiple radiation collision cascades suggest that dislocation dense regions reduce the probability of surviving point defects forming clusters. This finding seems to result from the dislocation cores’ ability to absorb point defects, thereby limiting the formation of clusters during the recombination phase. However, after successive cascades, the ability of the dislocation entanglement to prevent defect formation is diminished and defect cluster formation trends to an equilibrium. Accordingly, it is hypothesized that in an experimental setting the dislocation cells will initially act as neutral sinks for point defects, potentially delaying the onset of radiation damage effects; but, these structures will be degraded due to radiation-enhanced diffusion and lose their effectiveness at doses above 1.5–2 displacements per atom. Nevertheless, the ability to produce complex and feature-specific microstructures suggests promise for the design of radiation tolerant materials using AM methods.

36 MATERIALS SCIENCE↗

Score-based denoising for atomic structure identification

We propose an effective method for removing thermal vibrations that complicate the task of analyzing complex dynamics in atomistic simulation of condensed matter. Our method iteratively subtracts thermal noises or perturbations in atomic positions using a denoising score function trained on synthetically noised but otherwise perfect crystal lattices. The resulting denoised structures clearly reveal underlying crystal order while retaining disorder associated with crystal defects. Purely geometric, agnostic to interatomic potentials, and trained without inputs from explicit simulations, our denoiser can be applied to simulation data generated from vastly different interatomic interactions. The denoiser is shown to improve existing classification methods, such as common neighbor analysis and polyhedral template matching, reaching perfect classification accuracy on a recent benchmark dataset of thermally perturbed structures up to the melting point. Demonstrated here in a wide variety of atomistic simulation contexts, the denoiser is general, robust, and readily extendable to delineate order from disorder in structurally and chemically complex materials.

36 MATERIALS SCIENCE↗