Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite Elements”

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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Mechanical Performance of the First Two Prototype 4.5 m Long Nb$_3$Sn Low-$β$ Quadrupole Magnets for the Hi-Lumi LHC Upgrade

The U.S. High-Luminosity LHC Accelerator Upgrade Project (HL-LHC AUP) team is collaborating with CERN in the design and fabrication of the first 4.5 m long MQXFA magnets, a 150 mm aperture high-field Nb$_3$Sn quadrupole magnet that uses the aluminum shell-based bladder-and-key technology. The first two prototype magnets, MQXFAP1 and MQXFAP2, were assem-bled and tested while the first pre-series structure (MQXFA03) was in fabrication. This paper summarizes the mechanical perfor-mance of these prototype structures based on the comparison of measured strain gauge data with finite element model analyses from all load steps to powering. The MQXFAP1 magnet almost reached ultimate current before a short to ground was detected and the test was stopped. The MQXFAP2 magnet experienced a low training performance due to a fractured aluminum shell. MQXFAP1b was rebuilt with a new replacement coil, but an old coil limited the magnet from achieving the ultimate current. The mitigations and analyses of these prototype magnets are discussed in the context of the transition to pre-series production.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Mitigation of spall fracture by evolving porosity

Ductile materials subject to impact loading conditions can undergo spall fracture when an incoming compressive stress wave reflects off interfaces and free surfaces as a tensile stress wave. Experimental observations suggest that in ductile materials, spall fracture is driven by the evolution of porosity. However, the presence of initial porosity in ductile materials also introduces plastic compressibility, which can attenuate the incoming compressive stress wave and, as a result, reduce the amplitude of the reflected tensile stress wave. This, in turn, can mitigate or delay spall fracture. In this work, we report on finite deformation finite element calculations that analyze the response of porous ductile materials subjected to impact loading conditions. Two sets of calculations are carried out, in the first set the material contains initial porosity values ranging from 0% to 5% while in the second set the material also undergoes stress-controlled porosity nucleation. Both sets of calculations are carried out for a wide range of imposed impact velocities. Here, our results show that porosity in ductile materials can, under certain circumstances, mitigate spall fracture by attenuating stress waves. Results correlating the effects of impact velocity, initial porosity, and porosity nucleation on spall fracture are presented and the underlying mechanisms are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2021

GPU algorithms for Efficient Exascale Discretizations

In this paper we describe the research and development activities in the Center for Efficient Exascale Discretization within the US Exascale Computing Project, targeting state-of-the-art high-order finite-element algorithms for high-order applications on GPU-accelerated platforms. Furthermore, we discuss the GPU developments in several components of the CEED software stack, including the libCEED, MAGMA, MFEM, libParanumal, and Nek projects. We report performance and capability improvements in several CEED-enabled applications on both NVIDIA and AMD GPU systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Lattice Design and Advanced Modeling to Guide the Design of High-Performance Lightweight Structural Materials

Lightweight structural materials are required to increase the mobility of fission batteries. The materials must feature a robust combination of mechanical properties to demonstrate structural resilience. The primary objective of this project is to produce lightweight structural materials whose strength-to-weight ratios exceed those of the current widely used structural materials such as 316L stainless steels (316L SS). To achieve this, advanced modeling and simulation tools were employed to design lattice structures with different lattice parameters and different lattice types. A process was successfully developed for transforming lattice-structures models into Multiphysics Object Oriented Simulation Environment (MOOSE) inputs. Finite element modeling (FEM) was used to simulate the uniaxial tensile testing of the lattice-structured parts to investigate the stress distribution at a given displacement. The preliminary results showed that the lattice-structured sample displayed a lower Young’s modulus in comparison with the solid material and that the unit cell size of the lattice had a minimal effect. The novelty here is to apply up-front modeling to determine the best structure for the application before actually producing the sample. The approach of using modeling as a guiding tool for preliminary material design can significantly save time and cost for material development.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Simulation-driven optimization of high-order meshes in ALE hydrodynamics

Here we propose tools for high-order mesh optimization and demonstrate their benefits in the context of multi-material Arbitrary Lagrangian-Eulerian (ALE) compressible shock hydrodynamic applications. The mesh optimization process is driven by information provided by the simulation which uses the optimized mesh, such as shock positions, material regions, known error estimates, etc. These simulation features are usually represented discretely, for instance, as finite element functions on the Lagrangian mesh. The discrete nature of the input is critical for the practical applicability of the algorithms we propose and distinguishes this work from approaches that strictly require analytical information. Our methods are based on node movement through a high-order extension of the Target-Matrix Optimization Paradigm (TMOP). The proposed formulation is fully algebraic and relies only on local Jacobian matrices, so it is applicable to all types of mesh elements, in 2D and 3D, and any order of the mesh. We discuss the notions of constructing adaptive target matrices and obtaining their derivatives, reconstructing discrete data in intermediate meshes, node limiting that enables improvement of global mesh quality while preserving space-dependent local mesh features, and appropriate normalization of the objective function. The adaptivity methods are combined with automatic ALE triggers that can provide robustness of the mesh evolution and avoid excessive remap procedures. The benefits of the new high-order TMOP technology are illustrated on several simulations performed in the high-order ALE application BLAST.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Intergranular ductile failure of materials with plastically heterogeneous grains

In several precipitation hardened alloys that are susceptible to intergranular ductile failure, precipitation does not always occur uniformly throughout the microstructure, and regions close to grain boundaries may remain precipitate-free. These precipitate-free zones (PFZs) in the material microstructures result in plastically heterogeneous grains, since PFZs are expected to have lower yield strength but higher strain-hardenability compared to the precipitate containing grain interior. Experimentally, the presence of PFZs in precipitation-hardened alloys has been associated with both increase and decrease in materials’ ductility, with or without significant change in the strength. Thus, to understand and rationalize the experimental observations, we have carried out extensive microstructure-based finite element calculations of intergranular ductile failure in materials under tensile loading conditions. In the calculations, both the grain boundaries and PFZs are discretely modeled, and a wide range of the values of yield strength, strain-hardenability and width of PFZs in the material microstructures are analyzed. Our results show that the effects of PFZs on the overall mechanical response of the material strongly depend on the values of yield strength and strain-hardenability of PFZs. Here, there exists an optimum combination of the values of these two parameters that can result in the overall ductility and tensile strength of the material microstructures with PFZs being greater than the microstructures without PFZs.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Full-field quantitative visualization of shock-driven pore collapse and failure modes in PMMA

The dynamic collapse of pores under shock loading is thought to be directly related to hot spot generation and material failure, which is critical to the performance of porous energetic and structural materials. However, the shock compression response of porous materials at the local, individual pore scale is not well understood. This study examines, quantitatively, the collapse phenomenon of a single spherical void in PMMA at shock stresses ranging from 0.4 to 1.0 GPa. Using a newly developed internal digital image correlation technique in conjunction with plate impact experiments, full-field quantitative deformation measurements are conducted in the material surrounding the collapsing pore for the first time. The experimental results reveal two failure mode transitions as shock stress is increased: (i) the first in situ evidence of shear localization via adiabatic shear banding and (ii) dynamic fracture initiation at the pore surface. Numerical simulations using thermo-viscoplastic dynamic finite element analysis provide insights into the formation of adiabatic shear bands (ASBs) and stresses at which failure mode transitions occur. Further numerical and theoretical modeling indicates the dynamic fracture to occur along the weakened material inside an adiabatic shear band. Finally, analysis of the evolution of pore asymmetry and models for ASB spacing elucidate the mechanisms for the shear band initiation sites, and elastostatic theory explains the experimentally observed ASB and fracture paths based on the directions of maximum shear.

42 ENGINEERING↗
DOE OSTI2024

Fabrication, Modeling, and Testing of a Prototype Thermal Energy Storage Containment

Increasing penetration of variable renewable energy resources requires the deployment of energy storage at a range of durations. Long-duration energy storage (LDES) technologies will fulfill the need to firm variable renewable energy resource output year round; lithium-ion batteries are uneconomical at these durations. Thermal energy storage (TES) is one promising technology for LDES applications because of its siting flexibility and ease of scaling. Particle-based TES systems use low-cost solid particles that have higher temperature limits than the molten salts used in traditional concentrated solar power systems. A key component in particle-based TES systems is the containment silo for the high-temperature (>1100 degrees C) particles. This study combined experimental testing and computational modeling methods to design and characterize the performance of a particle containment silo for LDES applications. A laboratory-scale silo prototype was built and validated the congruent transient finite element analysis (FEA) model. The performance of a commercial-scale silo was then characterized using the validated model. The commercial-scale model predicted a storage efficiency above 95% after 5 days of storage with a design storage temperature of 1200 degrees C. Insulation material and concrete temperature limits were considered as well. The validation of the methodology means the FEA model can simulate a range of scenarios for future applications. This work supports the development of a promising LDES technology with implications for grid-scale electrical energy storage, but also for thermal energy storage for industrial process heating applications.

clean energy↗
DOE OSTI2021

Toughening of interface networks through the introduction of weak links

The resistance of polycrystalline materials to intergranular fracture depends on the toughness of the interfaces between neighboring grains. Using finite element simulations on model, two-dimensional microstructures, we demonstrate that the steady-state crack growth resistance of a polycrystal that contains exclusively high toughness interfaces is suboptimal: its crack growth resistance may be increased by introducing a minority fraction of weak interfaces at random locations within the interface network. Furthermore, the weaker are these interfaces, the greater is their toughening effect. These weak interfaces toughen the material through crack deflection and branching, nucleation of secondary cracks, and inelastic dissipation within intact interfaces. Furthermore, our findings suggest that fracture resistance-oriented material development should not aim to eliminate all weak interfaces, but rather judiciously incorporate a low fraction of them.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Quantifying leakage fields at ionic grain boundaries using off-axis electron holography

The electrical properties of interfaces in semiconductors and ionic conductors are immensely important in a wide range of applications. Electron holography is ideally suited for the direct measurement of the electrostatic potential of such interfaces. A key challenge with this approach is the contribution of the leakage field from the sample to the observed electron phase shift. This leakage field cannot be a priori independently determined and can cause an overestimation of the phase shift. In this work, we use finite element simulations to compute the three-dimensional electrostatic potential in the vicinity of an interface associated with a given interfacial charge density distribution. We then evaluate the predicted phase shift and demonstrate that the leakage field strongly affects the recovery of the projected interface potential. From the difference between the true potential and uncorrected, recovered potential, we propose a method to correct for this effect. We then demonstrate the application of this methodology to the analysis of experimental off-axis electron holography data acquired from the grain boundaries in lightly doped ceria.

36 MATERIALS SCIENCE↗
DOE OSTI2020

In situ measurement and modelling of the growth and length scale of twins in α -uranium

Twinning is an important deformation mechanism in crystalline materials. Twins nucleate, then grow, forming a discrete pattern of deformation bands with a certain thickness. The mechanism for twin growth is not completely understood. We present a simple phase field model which captures the physics of twin growth and reproduces the details observed during in situ electron backscatter diffraction experiments on α-uranium. The interaction between residual dislocations at twin interfaces and mobile dislocations in untwinned regions increases the stress needed for twin growth. This is described by a nonlocal term in the proposed constitutive equations: the nucleation stress of a twin increases proportionally to the twin phase field in a neighborhood. Competition between slip and twinning favors the nucleation of a new twin rather than the growth of a preexisting twin. The phase field model is coupled with a crystal plasticity finite element solver that includes the plastic deformation due to twinning. Simulations are able to reproduce the number of twins and their thicknesses as a function of the strain observed during in situ electron backscatter diffraction experiments. The stress-strain curve is also reproduced. This comparison allows the values of the nucleation stress and the interaction strength between residual and mobile dislocations to be found. It also gives an estimation of the density of residual dislocations that is consistent with the observed twin thickness. Finally, this model can be applied to understand microstructural effects in materials using twinning as a strengthening mechanism.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Controlling superconductivity of CeIrIn 5 microstructures by substrate selection

Superconductor/metal interfaces are usually fabricated in heterostructures that join these dissimilar materials. A conceptually different approach has recently exploited the strain sensitivity of heavy-fermion superconductors, selectively transforming regions of the crystal into the metallic state by strain gradients. The strain is generated by differential thermal contraction between the sample and the substrate. Here, we present an improved finite-element model that reliably predicts the superconducting transition temperature in CeIrIn 5 even in complex structures. Different substrates are employed to tailor the strain field into the desired shapes. Using this approach, both highly complex and strained as well as strain-free microstructures are fabricated to validate the model. This enables a high degree of control over the microscopic strain fields and forms the basis for more advanced structuring of superconductors as in Josephson junctions yet also finds natural use cases in any material class in which a modulation of the physical properties on a chip is desirable.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Shear Band Formation in Thin-Film Multilayer Columns Under Compressive Loading: A Mechanistic Study

Micro-pillar compression is a popular experimental technique used for characterizing the mechanical behavior of nano- and micro-laminates. The compressive stress–strain response of the column-shaped thin-film composite can be measured, and the deformation and damage features can be revealed by post-test cross-section microscopy. The development of plastic instability in the form of localized strain concentration (shear bands), leading to eventual failure, is frequently observed. In the present study, a computational approach is used to illustrate the commonality of shear band formation from a continuum standpoint. Systematic finite element analyses are conducted, showing that the strain field tends to become localized once plastic yielding commences. Distinct shear offsets of the layered structure can be revealed from the numerical model, which is similar to those observed in experiments. The actual appearance of shear bands depends on the materials’ constitutive behavior and precise geometries. Post-yield strain hardening reduces the propensity of shear band formation, while strain softening enhances it. Imperfections such as the undulated layer geometry, as well as the frictional characteristics between the specimen and test apparatus, can also influence the shear band morphology and overall stress–strain response.

finite element modeling↗
DOE OSTI2020

Synchrotron CT imaging of lattice structures with engineered defects

Understanding mechanical failure, crack propagation, and compressive behavior at the micrometer scale is essential for tailoring material properties for structural performance in cellular materials. Typically, modeling of traditional polymer foam materials is clouded by the lack of control in material morphology and its inherent stochastic structure. Additive manufacturing with sub-micrometer resolution provides a direct path for experimenters to specifically tailor structures needed by modelers to explicitly probe mechanical performance. Using laboratory-based 3D X-ray computed tomographic imaging (CT), the examination of deformation and damage provides a critical path to understand how these soft materials behave. Additionally, synchrotron CT yields realistic information at higher strain rates to directly validate the robustness of our finite element modeling. For this study, nanolithographic printing was employed to generate a series of engineered lattices with increasing levels of defects through the random removal of ligaments. These structures were mechanically tested and imaged with laboratory-based microCT. Additionally, synchrotron experiments were conducted in which the structures were imaged in 3D at 14 Hz during compression at a 0.4 s –1 strain rate. These 3D images show the changes in the structure as the ligaments bend, buckle and fracture in real time. This technique provides a robust framework for developing our methodologies and future exploration of engineered structures.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Predicting extreme anisotropy and shape variations in impact testing of tantalum single crystals

We report recent Taylor cylinder impact tests carried out for Ta single crystals showed strong variations in dimensional changes for different crystallographic directions aligned with the cylindrical axis. In order to capture the effect of crystallography on the deformation characteristics and final shapes of the impacted cylinders, a single crystal material subroutine is adapted and embedded in the solid mechanics/dynamics Finite Element solver Abaqus to simulate the aforementioned single crystal Ta Taylor impact experiments. Details of the coupled model implementation, and insights on the role played by single crystal anisotropic flow on the deformation behavior across a broad range of strain rates and temperatures for different single crystal orientations are presented and discussed. We demonstrate the predictive capability of the adopted crystal plasticity model to capture the significant role played by crystal orientation-induced anisotropy, as well as strain hardening and adiabatic heating, on the dynamic deformation response of crystalline materials. This re-emphasizes the need of microstructure-aware models to improve the accuracy of simulations for high-consequence engineering design.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Nanoscopic Imaging of Self-Propelled Ultrasmall Catalytic Nanomotors

Ultrasmall nanomotors (<100 nm) are highly desirable nanomachines for their size-specific advantages over their larger counterparts in applications spanning nanomedicine, directed assembly, active sensing, and environmental remediation. While there are extensive studies on motors larger than 100 nm, the design and understanding of ultrasmall nanomotors have been scant due to the lack of high-resolution imaging of their propelled motions with orientation and shape details resolved. Here, we report the imaging of the propelled motions of catalytically powered ultrasmall nanomotors─hundreds of them─at the nanometer resolution using liquid-phase transmission electron microscopy. These nanomotors are Pt nanoparticles of asymmetric shapes (“tadpoles” and “boomerangs”), which are colloidally synthesized and observed to be fueled by the catalyzed decomposition of NaBH4 in solution. Statistical analysis of the orientation and position trajectories of fueled and unfueled motors, coupled with finite element simulation, reveals that the shape asymmetry alone is sufficient to induce local chemical concentration gradient and self-diffusiophoresis to act against random Brownian motion. Our work elucidates the colloidal design and fundamental forces involved in the motions of ultrasmall nanomotors, which hold promise as active nanomachines to perform tasks in confined environments such as drug delivery and chemical sensing.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Measurement of helicon waves with phase contrast imaging on DIII-D – A theoretical feasibility study

A DIII-D high-beta H-mode discharge, with I p = 850 kA, B t = 2.1 T n e = 4 10 19 m –3 , has been designed to validate full wave modeling of helicon waves by optimizing the expected response of the Phase Contrast Imaging diagnostic. Helicon waves have been predicted to have high current drive efficiency off-axis without facing the accessibility issues of lower-hybrid waves. To test these predictions experimentally, DIII-D has recently commissioned a high-power helicon antenna. To confidently predict the behavior of helicon waves in future devices, measurements of their fundamental properties and validation against models will be essential. Phase contrast imaging (PCI) is an absolutely calibrated internal reference interferometer able to measure density fluctuations with radial wavenumbers k R between 1.5 cm –1 and 20 cm –1 . For helicon waves 2 cm –1 < k R < 10 cm –1 is expected, allowing PCI to measure their envelope and wavenumber spectrum. This makes PCI a powerful tool for the validation of state-of-the-art models, like the AORSA full wave code. AORSA is used to compute the density perturbations measured by the PCI with 2D calculations corresponding to 11 different toroidal mode numbers combined to resolve the trajectory of the helicon wave in 3D. This is necessary because the waves travel 120 degrees toroidally from the antenna to the PCI. Here, a cold plasma finite element model (CPFEM) [4] is used to predict propagation through the scrape-off layer. The result of the CPFEM model is connected to AORSA by creating an artificial Gaussian antenna on the last closed flux surface. PCI shows best results for waves with small vertical wavenumbers k z . Modeling the helicon waves for several past DIII-D experiments shows that k z is minimized if the intersection of the helicon and the PCI laser beams occurs in the midplane. For such an optimized scenario the predicted signal level is two orders of magnitude larger than the background density fluctuations arising from broadband turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Hydrogen charging and desorption from microstructural viewpoint: A method for deconvoluting hydrogen desorption spectra and unveiling the hydrogen-microstructure interaction

Understanding the interaction of hydrogen with microstructural features in metallic materials is crucial for designing hydrogen-resistant alloys. Although thermal desorption spectroscopy (TDS) is widely used for investigating the hydrogen binding behavior of various microstructural features, its application to face-centered cubic (fcc) metals and alloys that exhibit low hydrogen diffusivity is limited due to the lumped TDS desorption signals. This paper shows that, by coupling a Sofronis–McMeeking type hydrogen transport model with a microstructure-informed finite-element model, TDS data can be deconvoluted to reveal the underlying adsorption–diffusion–desorption processes, hydrogen diffusivity, and trap-binding energies. In conclusion, the austenitic steel SS316L in solution-annealed condition is used as a demonstration material, and we focused on investigating the interaction of deuterium (hydrogen isotope) with grain boundaries, which is difficult to investigate from experiments alone but critical for design of alloys for hydrogen infrastructure.

Finite element simulation, Polycrystalline Microst↗