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 379 records · Page 21

Specialty functions singularity mechanics problems

The focus is in the development of more accurate and efficient advanced methods for solution of singular problems encountered in mechanics. At present, finite element methods in conjunction with special functions, boolean sum and blending interpolations are being considered. In dealing with systems which contain a singularity, special finite elements are being formulated to be used in singular regions. Further, special transition elements are being formulated to couple the special element to the mesh that models the rest of the system, and to be used in conjunction with 1-D, 2-D and 3-D elements within the same mesh. Computational simulation with a least squares fit is being utilized to construct special elements, if there is an unknown singularity in the system. A novel approach is taken in formulation of the elements in that: (1) the material properties are modified to include time, temperature, coordinate and stress dependant behavior within the element; (2) material properties vary at nodal points of the elements; (3) a hidden-symbolic computation scheme is developed and utilized in formulating the elements; and (4) special functions and boolean sum are utilized in order to interpolate the field variables and their derivatives along the boundary of the elements. It may be noted that the proposed methods are also applicable to fluids and coupled problems.

Sarigul, Nesrin↗

International Interlaboratory Compilation of Trace Element Concentrations in the CUP-2 Uranium Ore Concentrate Standard

The accurate and precise analysis of impurities in uranium ore concentrates is an important aspect of nuclear forensics. Defensibility of analytical data requires the use of analytical reference materials as part of quality control. This report presents a compilation of the results of trace element concentration measurements of the CUP-2 Uranium Ore Concentration Standard from 11 different laboratories. Various dissolution methods, instrumental platforms, and analytical methods were employed. The resulting data set contains concentration of 66 impurities with up to 139 individual data points for each impurity. Consensus values have been assigned to each impurity following a statistical analysis of the data set.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Modeling the Effective Elasticity of Anisotropic Porous Materials

The development and optimization of composite materials designed for thermal protection of NASA’s spacecraft require understanding their physical response to high-enthalpy environments. To predict their macro-scale properties and behavior, high-fidelity 3D simulations are performed at the microscale on realistic representations of these composites. The digital microstructures are generated either synthetically or through X-ray micro-computed tomography reconstructions. One of the main challenges in the prediction of the structural response of heatshield materials is the computation of the effective elasticity of the fibrous composite, as well as the understanding of the deformation and stresses generated at the microscale. These are driven by the fiber layout within the microstructure and the distribution of the infused matrix. In this effort, the micro-mechanical linear elastic behavior of fibrous ablators is modeled using a numerical method based on the Multi-Point Stress Approximation (MPSA) finite volume scheme, a generalization of the more commonly used Multi-Point Flux Approximation (MPFA) that was presented at the 10th Ablation Workshop. To predict the behavior of fibrous and woven architectures, algorithms that compute the local fiber orientation are used. The implementation of the MPSA was verified using analytical solutions, engineering test cases, and compared against legacy Finite Element Analysis (FEA) software. The stress analysis models were then applied to real geometries used by NASA in thermal protection systems such as fibrous preforms and woven materials and the results were compared to experimental data.

Elasticity↗

Determination of the Fracture Parameters in a Stiffened Composite Panel

A modified J-integral, namely the equivalent domain integral, is derived for a three-dimensional anisotropic cracked solid to evaluate the stress intensity factor along the crack front using the finite element method. Based on the equivalent domain integral method with auxiliary fields, an interaction integral is also derived to extract the second fracture parameter, the T-stress, from the finite element results. The auxiliary fields are the two-dimensional plane strain solutions of monoclinic materials with the plane of symmetry at x(sub 3) = 0 under point loads applied at the crack tip. These solutions are expressed in a compact form based on the Stroh formalism. Both integrals can be implemented into a single numerical procedure to determine the distributions of stress intensity factor and T-stress components, T11, T13, and thus T33, along a three-dimensional crack front. The effects of plate thickness and crack length on the variation of the stress intensity factor and T-stresses through the thickness are investigated in detail for through-thickness center-cracked plates (isotropic and orthotropic) and orthotropic stiffened panels under pure mode-I loading conditions. For all the cases studied, T11 remains negative. For plates with the same dimensions, a larger size of crack yields larger magnitude of the normalized stress intensity factor and normalized T-stresses. The results in orthotropic stiffened panels exhibit an opposite trend in general. As expected, for the thicker panels, the fracture parameters evaluated through the thickness, except the region near the free surfaces, approach two-dimensional plane strain solutions. In summary, the numerical methods presented in this research demonstrate their high computational effectiveness and good numerical accuracy in extracting these fracture parameters from the finite element results in three-dimensional cracked solids.

Lin, Chung-Yi↗

RCSB Protein Data Bank: supporting research and education worldwide through explorations of experimentally determined and computationally predicted atomic level 3D biostructures

The Protein Data Bank (PDB) was established as the first open-access digital data resource in biology and medicine in 1971 with seven X-ray crystal structures of proteins. Today, the PDB houses >210 000 experimentally determined, atomic level, 3D structures of proteins and nucleic acids as well as their complexes with one another and small molecules ( e.g. approved drugs, enzyme cofactors). These data provide insights into fundamental biology, biomedicine, bioenergy and biotechnology. They proved particularly important for understanding the SARS-CoV-2 global pandemic. The US-funded Research Collaboratory for Structural Bioinformatics Protein Data Bank (RCSB PDB) and other members of the Worldwide Protein Data Bank (wwPDB) partnership jointly manage the PDB archive and support >60 000 `data depositors' (structural biologists) around the world. wwPDB ensures the quality and integrity of the data in the ever-expanding PDB archive and supports global open access without limitations on data usage. The RCSB PDB research-focused web portal at https://www.rcsb.org/ (RCSB.org) supports millions of users worldwide, representing a broad range of expertise and interests. In addition to retrieving 3D structure data, PDB `data consumers' access comparative data and external annotations, such as information about disease-causing point mutations and genetic variations. RCSB.org also provides access to >1 000 000 computed structure models (CSMs) generated using artificial intelligence/machine-learning methods. To avoid doubt, the provenance and reliability of experimentally determined PDB structures and CSMs are identified. Related training materials are available to support users in their RCSB.org explorations.

59 BASIC BIOLOGICAL SCIENCES↗

Kinetic isolation tether experiment

Progress was made on the analysis of tether damping and on experimentation of the control system on the laboratory simulator. The damping analysis considers the dynamics of a long tether connecting two spacecraft in Earth orbit, one of the spacecraft having dominant mass. In particular, it considers the material damping of the tether. The results show that, with properly chosen tether material and braiding structure, longitudinal vibration of the tethered system is well damped. A particularly effective method of implementing attitude control for tethered satellites is to use the tether tension force to generate control torques by moving the tether attach point relative to the satellite center of mass. A scaled, one dimensional laboratory simulation of the KITE mission was built and preliminary experiments of the proposed attitude control system were performed. The simulator was built to verify theoretical predictions of attitude controllability, and to investigate the technological requirements in order to implement this concept. A detailed description of the laboratory apparatus is provided, and in addition, the results of the preliminary experiments are presented and discussed.

Powell, J. David↗

Measuring Thermal Conductivities Of Rough Specimens

Thermal conductivity of irregularly surfaced specimen measured more accurately and reliably by use of heat-transfer adhesive establishing contact between specimen and reference standards. Method useful when specimen made up of fragments of material. Adhesive ensures good thermal contact to pieces not matching others in height. In rough-surfaced or unequally sized specimens, heat-flow pattern becomes distorted, and temperatures at measurement points do not accurately represent those throughout reference standards and specimen. Filling voids with conductive adhesive restores nearly uniform distribution heat flow.

Striepens, Aloys H.↗

Experimental Investigation of the Distribution of Shock Effects in Regolith Impact Ejecta Using an Ejecta Recovery Chamber

Because the mass-flux of solar system meteoroids is concentrated in the approx. 200 microns size range, small-scale impacts play a key role in driving the space weathering of regoliths on airless bodies. Quantifying this role requires improved data linking the mass, density and velocity of the incoming impactors to the nature of the shock effects produced, with particular emphasis on effects, such as production of impact melt and vapor, that drive the optical changes seen in space weathered regoliths. Of particular importance with regard to space weathering is understanding not only the composition of the shock melt created in small-scale impacts, but also how it is partitioned volumetrically between the local impact site and more widely distributed ejecta. To improve the ability of hypervelocity impact experiments to obtain this type of information, we have developed an enclosed sample target chamber with multiple-geometry interior capture cells for in-situ retention of ejecta from granular targets. A key design objective was to select and test capture cell materials that could meet three requirements: 1) Capture ejecta fragments traveling at various trajectories and velocities away from the impact point, while inducing minimal additional damage relative to the primary shock effects; 2) facilitate follow-up characterization of the ejecta either on or in the cell material by analytical SEM, or ex-situ by microprobe, TEM and other methods; and 3) enable the trajectories of the captured and characterized ejecta to be reconstructed relative to the target.

Christoffersen, R.↗

Integration of MAC/GMC into CalculiX, an Open Source Finite Element Code

An analysis framework is presented that makes available multiscale analysis of composite structures using the open-source FEA solver package CalculiX CrunchiX (CCX). At the center of this framework is the coupling of the Finite Element Analysis - Micromechanics Analysis Code (FEAMAC) library from NASA's Micromechanics Analysis Code with Generalized Method of Cells (MAC/GMC) to allow micromechanics analysis. The results show that the proposed coupling can be used with appropriate care for multiscale FEA simulations of composite materials. The largest error reported in this validation was in a four-point bend test specimen with an error of less than 1% difference in the maximum deflection of the beam.

Yapor Genao, Francisco A.↗

Validation of Floating Node Method Using Three-Point Bend Doubler Under Quasi-Static Loading

The NASA Advanced Composite Project (ACP), an industry/government/university partnership, has embarked upon the task of developing technology that can aid in reducing the time line for structural certification of aircraft composite parts using a combination of technologies, one of which is high fidelity damage progression computational methods. Phase II of this project included a task for validating an approach based on the Floating Node Method combined with Directional Cohesive Elements (FNM-DCZE). This paper discusses predicted damage onset and growth in a three-point bend doubler specimen compared to experimental results. Sensitivity of the simulations to mesh refinement as well as key material properties and thermal effects are studied and reported. Overall, qualitative results suggest the main aspects of the damage progression have been captured, with the simulated damage morphology and sequence of events resembling closely what was observed experimentally. Quantitatively, the first load-peak is predicted. However, the re-loading observed in the experiments, after the first load peak, is not captured numerically, suggesting further investigation may be worth pursuing.

Selvarathinam, Alex S.↗

Optoelectronic property comparison for isostructural Cu 2 BaGeSe 4 and Cu 2 BaSnS 4 solar absorbers

To target mitigation of anti-site defect formation in Cu 2 ZnSnS 4–x Se x , a new class of chalcogenides, for which Ba or Sr (group 2) replace Zn (group 12), has recently been introduced for prospective solar absorber application. Cu 2 BaGeSe 4 (CBGSe) and Cu 2 BaSnS 4 (CBTS) are two such compounds, which share a common trigonal crystal structure (P3 1 space group) and similar quasi-direct band gap (~2 eV). While CBTS-based films have already been studied, there are no reports yet on films and solar cells based on related CBGSe. To identify key differences and similarities in the electronic properties between these two materials, electronic characteristics (e.g., carrier concentration, mobility, electron affinity, defect levels, recombination, and charge carrier kinetics) of vacuum-deposited CBGSe and CBTS films are compared using a variety of characterization methods. Hall effect measurements reveal that CBGSe films have relatively higher hole carrier concentration and lower mobility (3 × 10 15 cm –3 , 0.6 cm 2 V –1 s –1 ) compared to CBTS (5 × 10 12 cm –3 , 3.5 cm 2 V –1 s –1 ). Photoelectron spectroscopy yields low electron affinity values for both CBGSe (3.7 eV) and CBTS (3.3 eV), pointing to the necessity of pursuing low electron affinity buffer materials for both types of absorbers. At low temperatures, CBGSe films show free-exciton photoluminescence, as well as pronounced deep-level emission at ~1.4 eV, while CBTS films exhibit a strong bound-exciton signal with noticeably less intense deep-level emission than for CBGSe. Charge carrier kinetics, transport, and recombination properties of both types of films are also analyzed using optical-pump terahertz-probe spectroscopy and time-resolved microwave conductivity. The first CBGSe prototype solar cells (using chemical bath deposited CdS as a buffer layer) show a maximum of 1.5% efficiency with ~0.62 V open-circuit voltage. Furthermore, the measured properties point to possible limiting factors for CBGSe and related films for PV and optoelectronics and provide insights on possible approaches for improvement within this multinary chalcogenide family.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The effects of space radiation on thin films of YBa2Cu3O(7-x)

This investigation had two objectives: (1) to determine the effects of space radiation on superconductor parameters that are most important in space applications; and (2) to determine whether this effect can be simulated with Co-60 gamma rays, the standard test method for space materials. Thin films of yttrium barium copper oxide (YBCO) were formed by coevaporation of Y, BaF2, and Cu and post-annealing in wet oxygen at 850 C for 3.5 h. The substrate used was (100) silicon with an evaporated zirconia buffer layer. The samples were characterized by four point probe electrical measurements as a function of temperature. The parameters measured were the zero resistance transition temperature T(sub c) and the room temperature resistance. The samples were then exposed to Co-60 gamma-rays in air and in pure nitrogen, and to 780 keV electrons, in air. The parameters were then remeasured. The results are summarized. The results indicate little or no degradation in the parameters measured for samples exposed up to 10 Mrads of gamma-rays in nitrogen. However, complete degradation is preliminarily attributed to the high level of ozone generated in the chamber by the gamma-ray interaction with air. It can be concluded that: (1) the electron component of space radiation does not degrade the critical temperature of the YBCO films described, at least for energies around 800 keV and doses similar to those received by surface materials on spacecraft in typical remote sensing missions; and (2) for qualifying this and other superconducting materials against the space-radiation threat the standard test method used in the aerospace industry, namely, exposure to Co-60 gamma-rays in air, may require some further investigation. As a minimum, the sample must be either in vacuum or in positive nitrogen pressure.

Herschitz, R.↗

Atomistic modeling of meso -timescale processes with $\mathrm{SEAKMC}$: A perspective and recent developments

On-the-fly kinetic Monte Carlo (kMC) methods have recently garnered significant attentions after successful applications to various atomic-scale problems using a timescale outside the reach of classical molecular dynamics. These methods play a critical role in modeling atomistic meso-timescale processes, and it is therefore essential to further improve their capabilities. In this report we review one of the on-the-fly kMC methods, Self-Evolving Atomistic kinetic Monte Carlo (SEAKMC) and propose two schemes that considerably enhance the efficiency of saddle point searches (SPSs) during the simulations. The performance of these schemes is tested using the diffusion of point defects in bcc Fe. In addition, we discuss approaches to significantly mitigate limitations of these schemes, which further improves their efficiencies. Importantly, these schemes improve the SPS efficiency not only for SEAKMC but also for other on-the-fly kMC methods, broadening the applications of on-the-fly kMC simulations to complex meso-timescale problems.

36 MATERIALS SCIENCE↗

Evolution of microstructures in radiation fields using a coupled binary-collision Monte Carlo phase field approach

The simulation of radiation effects in materials broadly falls into two categories. At short time and length scales lies the modeling of primary radiation damage, such as point defect creation, energy deposition, and ballistic mixing. This is followed by the modeling at longer time scales of thermally activated microstructure evolution and defect reactions, such as recombination, clustering, and coarsening. The binary collision Monte Carlo method is an established, numerically efficient method for the computation of primary radiation damage. Conversely, the phase field method is a state of the art method for microstructure evolution on longer time and length scales. Here we present a concurrent coupling of these two methods, overcoming the difference between the discrete object Monte Carlo paradigm for primary radiation damage and the continuum field variable approach for microstructure evolution. The coupling is bidirectional, in which the microstructure evolution in the MOOSE finite element frame- work provides the spatial scattering data set for the charged particle trans- port and receives point defect, mass transport, and heat source terms from the simulated collision cascades that contribute to the field variable evolution. The concurrent coupling scheme is implemented in the code Magpie and demonstrated by investigating patterning for a model irradiated immiscible binary alloy. The results from the coupled binary collision Monte Carlo/phase field simulations reproduce the results of analytical models for phase separation, phase mixing, and patterning, supporting the approach and indicating its utility for modeling real materials systems.

36 MATERIALS SCIENCE↗

Geometrically Nonlinear Static Analysis of 3D Trusses Using the Arc-Length Method

Rigorous analysis of geometrically nonlinear structures demands creating mathematical models that accurately include loading and support conditions and, more importantly, model the stiffness and response of the structure. Nonlinear geometric structures often contain critical points with snap-through behavior during the response to large loads. Studying the post buckling behavior during a portion of a structure's unstable load history may be necessary. Primary structures made from ductile materials will stretch enough prior to failure for loads to redistribute producing sudden and often catastrophic collapses that are difficult to predict. The responses and redistribution of the internal loads during collapses and possible sharp snap-back of structures have frequently caused numerical difficulties in analysis procedures. The presence of critical stability points and unstable equilibrium paths are major difficulties that numerical solutions must pass to fully capture the nonlinear response. Some hurdles still exist in finding nonlinear responses of structures under large geometric changes. Predicting snap-through and snap-back of certain structures has been difficult and time consuming. Also difficult is finding how much load a structure may still carry safely. Highly geometrically nonlinear responses of structures exhibiting complex snap-back behavior are presented and analyzed with a finite element approach. The arc-length method will be reviewed and shown to predict the proper response and follow the nonlinear equilibrium path through limit points.

Hrinda, Glenn A.↗

Multireference diffusion Monte Carlo reaches 2D materials

Abstract Quantum confinement in 2D materials strongly enhances electronic correlation effects. Therefore, predicting the properties of these unique materials, with both a high level of accuracy and computational efficiency, without relying on adjustable parameters or functionals, remains an outstanding theoretical challenge. The majority of theoretical studies are based on the approximations of density functional theory (DFT). The reliability of DFT predictions are heavily dependent on the choice of an approximated exchange-correlation functional. Here, we estimate the magnitude of impact of correlation on the total energy for the quintessential 2D material, graphene, by performing and comparing state-of-the-art selected CI and quantum Monte Carlo extrapolated calculations for a single unit cell at the$$\Gamma$$point. We demonstrate that Self-Healing Diffusion Monte Carlo (SHDMC) obtains a very compact, but high-quality wavefunction for this system that lacks the strong basis set dependence displayed by state of the art quantum chemistry methods. The SHDMC wavefunction is of higher quality compared to that obtained from sCI, in the same orbital basis, while being$$\sim$$ 1000 times smaller in terms of determinant count compared to sCI. We also demonstrate that extrapolating SHDMC results to the infinite determinant limit compares extremely well with complete basis set extrapolated sCI. Our work paves the way for future validation of SHDMC applied to challenging 2D materials.

Science & Technology - Other Topics↗

The effects of space radiation on thin films of YBa2Cu3O(sub 7-x)

This investigation had two objectives: (1) to determine the effects of space radiation on superconductor parameters that are most important in space applications; and (2) to determine whether this effect can be simulated with Co-60 gamma rays, the standard test method for space materials. Thin films of yttrium barium copper oxide (YBCO) were formed by coevaporation of Y, BaF2, and Cu and post-annealing in wet oxygen at 850 C for 3.5 h. The substrate used was (100) silicon with an evaporated zirconia buffer layer. The samples were characterized by four point probe electrical measurements as a function of temperature. The parameters measured were the zero resistance transition temperature (T sub c) and the room temperature resistance. The samples were then exposed to Co-60 gamma-rays in air and in pure nitrogen, and to 780 keV electrons, in air. The parameters were then remeasured. The results are summarized. The results indicate little or no degradation in the parameters measured for samples exposed up to 10 Mrads of gamma-rays in nitrogen. However, complete degradation of samples exposed to 10-Mrad in air was observed. This degradation is preliminarily attributed to the high level of ozone generated in the chamber by the gamma-ray interaction with air. It can be concluded that: (1) the electron component of space radiation does not degrade the critical temperature of the YBCO films described, at least for energies around 800 keV and doses similar to those received by surface materials on spacecraft in typical remote sensing missions; and (2) for qualifying this and other superconducting materials against the space-radiation threat the standard test method in the aerospace industry, namely, exposure to Co-60 gamma-rays in air, may require some further investigation. As a minimum, the sample must be either in vacuum or in positive nitrogen pressure.

Herschitz, R.↗

Parallel projection—An improved return mapping algorithm for finite element modeling of shape memory alloys

Here, we present a novel finite element analysis of inelastic structures containing Shape Memory Alloys (SMAs). Phenomenological constitutive models for SMAs lead to material nonlinearities, that require substantial computational effort to resolve. Finite element analysis methods, which rely on Gauss quadrature integration schemes, must solve two sets of coupled differential equations: one at the global level and the other at the local, i.e. Gauss point level. In contrast to the conventional return mapping algorithm, which solves these two sets of coupled differential equations separately using a nested Newton procedure, we propose a scheme to solve the local and global differential equations simultaneously. In the process we also derive closed-form expressions used to update the internal/constitutive state variables, and unify the popular closest-point and cutting plane methods with our formulas. Numerical testing indicates that our method allows for larger thermomechanical loading steps and provides increased computational efficiency, over the standard return mapping algorithm.

42 ENGINEERING↗