Search NASA⌕ Search

SEARCH · Search NASA

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 235 records · Page 13

New Tool Creates a Big Stir

A new self-adjusting, retractable pin tool for friction stir welding is now used in the manufacturing of components for NASA Space Shuttles. Friction stir welding is a process that makes straight-line welds without bringing the parent material to a liquid state. This is accomplished through high-speed rotation, which generates frictional heat between the welding tool and the piece being welded. This heat causes the material to soften to the point of plasticity without allowing it to melt. The plasticized material is then transferred from the front edge of the welding tool to the trail edge, where it joins the pieces being welded. However, a major flaw of this method is its reliance on a single-piece pin tool. The weld is left unfinished and a hole remains where the pin was inserted. The hole must be covered with a rivet in order to preserve the integrity of the weld. The NASA-developed pin tool, however, eliminates the need for this finishing step, as its retraction allows continuous rewelding at lesser depths, until the hole is completely closed. With this NASA technology, welding of higher strength alloys, as well as non-planer and variable thickness structures can be achieved.

Source record↗

The intact capture of hypervelocity dust particles using underdense foams

The impact of a hypervelocity projectile (greater than 3 km/s) is a process that subjects both the impactor and the impacted material to a large transient pressure distribution. The resultant stresses cause a large degree of fragmentation, melting, vaporization, and ionization (for normal densities). The pressure regime magnitude, however, is directly related to the density relationship between the projectile and target materials. As a consequence, a high-density impactor on a low-density target will experience the lowest level of damage. Historically, there have been three different approaches toward achieving the lowest possible target density. The first employs a projectile impinging on a foil or film of moderate density, but whose thickness is much less than the particle diameter. This results in the particle experiencing a pressure transient with both a short duration and a greatly reduced destructive effect. A succession of these films, spaced to allow nondestructive energy dissipation between impacts, will reduce the impactor's kinetic energy without allowing its internal energy to rise to the point where destruction of the projectile mass will occur. An added advantage to this method is that it yields the possibility of regions within the captured particle where a minimum of thermal modification has taken place. Polymer foams have been employed as the primary method of capturing particles with minimum degradation. The manufacture of extremely low bulk density materials is usually achieved by the introduction of voids into the material base. It must be noted, however, that a foam structure only has a true bulk density of the mixture at sizes much larger than the cell size, since for impact processes this is of paramount importance. The scale at which the bulk density must still be close to that of the mixture is approximately equal to the impactor. When this density criterion is met, shock pressures during impact are minimized, which in turn maximizes the probability of survival for the impacting particle. The primary objectives of the experiment are to (1) Examine the morphology of primary and secondary hypervelocity impact craters. Primary attention will be paid to craters caused by ejecta during hypervelocity impacts of different substrates. (2) Determine the size distribution of ejecta by means of witness plates and collect fragments of ejecta from craters by means of momentum-sensitive mcropore foam. (3) Assess the directionality of the flux by means of penetration-hole alignment of thin films placed above the cells. (4) Capture intact the particles that perforated the thin film and entered the cell. Capture media consisted of both previously flight-tested micropore foams and aerogel. The foams had different latent heats of fusion and, accordingly, will capture particles over a range of momenta. Aerogel was incorporated into the cells to determine the minimum diameter than can be captured intact.

Maag, Carl R.↗

A finite difference informed random walker (FDiRW) solver for strongly inhomogeneous diffusion problems

In nature, many complex multi-physics coupling problems exhibit strong diffusivity inhomogeneity. For instance, in the context of radionuclide absorption by porous wasteform materials within a flowing waste stream, the difference of species’ diffusivity in solid and liquid phases spans by 3~8 orders of magnitude. To solve the diffusion equations with strongly inhomogeneous diffusivity, traditional discretization-based methods, such as the Finite Difference Method (FDM), require infinitesimally small time steps (<10 -10 ) as high spatial resolutions are employed in most microstructure evolution processes, leading to prohibitively high computational costs. Here, this work developed an integrated numerical approach (FDiRW: Finite Difference informed Random Walk) to tackle this challenge. The idea is that utilizing the Random Walk concept, the fast diffusion is modeled as a superposition of point source’s solution for a concentration distribution while FDM is used to obtain the point source’s solution at each node. A mesh-coarsening algorithm is developed to generate an exclusive coarse mesh for FDiRW approach to maximize its efficiency. The effectiveness of the coarse mesh-based FDiRW approach is validated by benchmarking Finite Difference solutions. Numerical results demonstrated that FDiRW achieves a remarkable 1000x computational efficiency improvement over FDM while preserving desired accuracy for a medium-sized model of 192 × 192 × 192 grids. Finally, as models scale up, a floating-point operations (PLOPs) analysis of the FDiRW algorithm reveals that its computational complexity grows quadratically in terms of the number of nodes employed in computation.

36 MATERIALS SCIENCE↗

Tetravalent chromium doped laser materials and NIR tunable lasers

A method is described to improve and produce purer Cr.sup.4+-doped laser materials and lasers with reduced co-incorporation of chromium in any other valence states, such as Cr.sup.3+, Cr.sup.2+, Cr.sup.5+, and Cr.sup.6+. The method includes: 1) certain crystals of olivine structure with large cation (Ca) in octahedral sites such as Cr.sup.4+:Ca.sub.2GeO.sub.4, Cr.sup.4+:Ca.sub.2SiO.sub.4, Cr.sup.4+:Ca.sub.2Ge.sub.xSi.sub.1-xO.sub.4 (where 0<x<1), and/or 2) high-temperature solution growth techniques that enable the growth of the crystals below the temperature of polymorphic transitions by using low melting point solvent based on oxide, fluoride and/or chloride compounds. Purer Cr.sup.4+-doped laser materials are characterized by a relatively high concentration of Cr.sup.4+-lasing ion in crystalline host that makes these materials suitable for compact high power (thin disk/wedge) NIR laser applications.

Alfano, Robert R.↗

Impact of Thermoplastic Composites: Testing and Modeling

Thermoplastic Composites (TPCs) are being increasingly considered for aerospace applications given their faster manufacturing cycles and lower cost. The NASA High-Rate Composite Aircraft Manufacturing (HiCAM) project aims to evaluate and mature thermoplastic manufacturing technologies to achieve a 4X – 6X increase in the production rate of Next Generation Single Aisle (NGSA) commercial aircraft. Within the purview of the NASA HiCAM project, this work investigates the response of thermoplastic composites under High Energy Dynamic Impact (HEDI) conditions. HEDI tests were conducted on panels fabricated with a carbon fiber reinforced low melt semi-crystalline resin TC1225 LMPAEK T700G (T700/LMPAEK) material system. A blunt metal projectile was chosen as the impactor which impacted the T700/LMPAEK panels over a range of impact velocities imparted by a gas gun test setup. As expected, lower velocity impacts caused the projectile to rebound, whereas higher velocity impacts resulted in the projectile penetrating the test panels. For the test cases in which the projectile rebounded, the predominant damage modes recorded with Ultrasonic (UT) scans included interlaminar delaminations which exhibited a “rotating fan” like structure when viewed in the through-thickness direction. Additionally, the dynamic deflection of the center point of the back face of the test panels and the velocity of the projectile was also measured during the test. Finite element models were developed to model a test case in which the projectile rebounded. These models were developed using LS-DYNA® wherein the primary objective was to evaluate the ability of the material model MAT299 to capture the dynamic response and damage modes in the T700/LMPAEK test panels. MAT299 is a Deformation Gradient Decomposition (DGD) based Continuum Damage Mechanics (CDM) material model. The interlaminar delaminations were modeled using the cohesive contact formulation available within LS-DYNA, wherein the mixed-mode fracture is captured via the Benzegaggh-Kenane (B-K) mode-dependent fracture energy interpolation law. The above-mentioned modeling methods have mostly been applied to thermoset composites and adapting them to TPCs involves addressing challenges associated with appropriately representing material behavior. Therefore, the proposed paper would discuss the systematic approach taken to adapt thermoset composites modeling practices for applications to TPCs, while appropriately addressing material behavior. The proposed paper would show that the predicted back face center point deflection correlated reasonably well with experimental data and the peak deflection was predicted to be within 5% of the experimental measurements. The projectile rebound velocity, while predicted to be higher than the experimental measurement, was within reasonable bounds. Additionally, it will be shown that the predicted delamination shapes were in reasonable agreement with experimental data.

Composite Materials↗

Electronic structure and exchange interactions in altermagnetic MnGeP 2 in the quasiparticle self-consistent G W approach

The quasiparticle self-consistent 𝐺⁡𝑊 method is used to study the electronic band structure, optical dielectric function, and exchange interactions in chalcopyrite, 𝐼⁢ ‾ 4 ⁢2⁢𝑑, structure MnGeP 2 . The material is found to be an antiferromagnetic semiconductor with a lowest direct gap of 2.44 eV at the Γ point and a lower indirect gap of 1.87 eV from Γ to 𝑀. The material is an altermagnet because the two magnetic atoms of opposite spin are related by a twofold rotation operation perpendicular to the main fourfold rotation inversion axis. The spin splittings along a low symmetry line like 𝑃⁢𝑁 is sizable while at k points on the diagonal mirror planes or on the twofold symmetry axes the spin splitting is zero. The exchange interactions are calculated using a linear response approach. The antiferromagnetic exchange interaction between nearest neighbors in the primitive unit cell is dominating and found to be slightly decreasing upon carrier doping but not sufficiently to change the interaction to become ferromagnetic. The bare (noninteracting) and interacting transverse spin susceptibilities, which provide interatomic site exchange interactions after averaging over the muffin-tin spheres, are calculated from the 𝐺⁡𝑊 band structure and wave functions. From these exchange interactions, the spin wave spectra are obtained along the high symmetry lines and the Néel temperature is calculated using the mean-field and Tyablikov estimations. The dielectric function and the optical absorption spectra are calculated including excitonic effects using the Bethe Salpeter equation. The exchange interactions around Mn Ge defect sites is also studied. While we find it can generate ferromagnetic interactions with neighboring spins, we did not find direct evidence of producing an overall ferromagnetic phase. First, if Mn antisites are introduced by exchanging Mn with a nearby Ge, the interactions stay largely antiferromagnetic. Second, when we add additional Mn, in other words in Mn-rich stoichiometry, the Mn Ge antisites produce a strong ferromagnetic interaction primarily with the Mn in the same basal plane but weaker ferromagnetic interaction with adjacent plane Mn. The interactions between regular lattice Mn stay antiferromagnetic as before and thus favor keeping the antiferromagnetic order along the [001] direction. Adding Mn antisites, however, does lead to a metallic band structure.

36 MATERIALS SCIENCE↗

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↗

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↗

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.↗

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.↗

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.↗