SEARCH · Engineering Papers
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.
Determine how these plastic current collectors work and can they obviate the battery design burden for achieving passive propagation resistance?
Selective laser melting is receiving increasing interest as an additive manufacturing technique. Residual stresses induced by the large temperature gradients and inhomogeneous cooling process can favour the generation of cracks. In this work, a crystal plasticity finite element model is developed to simulate the formation of residual stresses and to understand the correlation between plastic deformation, grain orientation and residual stresses in the additive manufacturing process. The temperature profile and grain structure from thermal-fluid flow and grain growth simulations are implemented into the crystal plasticity model. An element elimination and reactivation method is proposed to model the melting and solidification and to reinitialize state variables, such as the plastic deformation, in the reactivated elements. The accuracy of this method is judged against previous method based on the stiffness degradation of liquid regions by comparing the plastic deformation as a function of time induced by thermal stresses. The method is used to investigate residual stresses parallel and perpendicular to the laser scan direction, and the correlation with the maximum Schmid factor of the grains along those directions. The magnitude of the residual stress can be predicted as a function of the depth, grain orientation and position with respect to the molten pool. The simulation results are directly comparable to X-ray diffraction experiments and stress–strain curves.
Classically, crystal plasticity modeling has used a range of constitutive equations, in which the incorporation of additional physics-based relationships typically results in additional model parameters. These additional parameters need to be reliably calibrated, which often necessitates the use of a range of experimental data acquired at multiple length scales. In this work, a crystal plasticity based finite element (CPFE) model for a dual-phase Titanium alloy, Ti–6Al–4V, is developed. The α and β phases of the microstructure are explicitly modeled. The model is calibrated using a systematic optimization routine and experimental data that consist of macroscopic stress-strain curves coupled with lattice strains on different crystallographic planes for the two phases. These experimental data were obtained from in situ high energy X-ray diffraction experiments for multiple material pedigrees, with varying crystallographic orientation distribution and β volume fractions. Depending on the thermomechanical-processing route and the heat treatment used to manufacture the alloy, Ti–6Al–4V can exist in a wide number of microstructural forms, which often results in the α and β phases either having well aligned slip systems (following the Burgers orientation relationship (BOR)) or possessing no alignment of the slip systems across the interphase boundary (not following the BOR). In this study, the fully-calibrated CPFE model is used to gain a comprehensive understanding of the deformation behavior of Ti–6Al–4V, specifically, the effect of microstructures that follow the BOR (or not) on time-dependent cyclic loading (including the effects of dwell hold times).
Here, ionic liquids (ILs) nanostructuring at electrified interfaces is of both fundamental and practical interest as these materials are increasingly gaining prominence in energy storage and conversion processes. However, much remains unresolved about IL potential-controlled (re)organization under highly polarized interfaces, mostly due to the difficulty of selectively probing both the distal and proximal surface layers of adsorbed ions. In this work, the structural dynamics of the innermost layer (<10 nm from the surface) were independently interrogated from that of the ionic layers in the sub-surface region (>100 nm from the surface), using an infrared (IR) spectroscopy approach. By tuning the metal fill factor of gold films deposited on conductive metal oxide-modified IR internal reflection elements, the charge-driven (re)structuring of the inner and distal layers of 1-butyl-1-methylpyrrolidinium trifluoromethanesulfonate is unveiled. Within a relatively wide potential region (~±1 V) bounding the potential of zero charges, the ionic liquid is shown to undergo a reversible (i.e., soft) reorganization whereby the innermost layer of anions (cations) is exchanged by a layer of cations (anions). Kinetically unhindered changes in the number density of constituent cations and anions largely follow electrostatic expectations in the subsurface region, whereas the innermost layer exhibits a pronounced hysteresis and very slow relaxation. Under larger negative potential bias, IL restructuring is characterized by a highly irreversible (i.e., hard) and intense interfacial densification of the BMPy + cations, consistent with the formation of nanoscale segregated liquids. The outcomes of this work reveal a plastic IL nanostructuring under a strong electric field.
The paper reviews recent work on fundamentals of elastic-plastic finite-element analysis and its applications to the mechanics of crack opening and growth in ductile solids. The presentation begins with a precise formulation of incremental equilibrium equations and their finite-element forms in a manner valid for deformations of arbitrary magnitude. Special features of computational procedures are outlined for accuracy in view of the near-incompressibility of elastic-plastic response. Applications to crack mechanics include the analysis of large plastic deformations at a progressively opening crack tip, the determination of J integral values and of limitations to J characterizations of the intensity of the crack tip field, and the determination of crack tip fields in stable crack growth.
A formulation for a dimensionless coefficient c b is derived that represents a scale of microstructural softening for alloys which follow Kocks-Mecking (K-M) work hardening behavior. The variation of the true plastic strain ε p between the proportional limit σ y and the strength σ u at the instability is determined using the Considère criterion. Parameterization of the model is limited to variables expressly measured within tensile experiments. Further development is now made for the softening factor c bi through the individual and sequential stages 3 and 4 of plastic deformation during K-M work hardening Q behavior. Application is shown for tensile test results of Ti-6Al-4V made by different additively manufactured (AM) processes. It is found that the variation in plastic strain as a function of c bi produces a continuous curve representative of the alloy system. The results of data analysis indicate that Θ o3 and c b3 increase, while Θ o4 and c b4 decrease, as the total plastic strain ε p increases. Furthermore, formulations derived for K-M stages 3 and 4 enable the evaluation of other material parameters such as the activation volume ν* for the onset of plastic deformation. This activation volume is found to be near constant, at a ν*-value of 0.353±0.036 nm 3 as computed using a strain-rate sensitivity of strength exponent m of 0.014, irrespective of the AM method used to produce the Ti-6Al-4V alloy.
We report at present, the available Ramberg–Osgood (R–O) parameters for different metals (e.g., in ASME code and other literatures) are static (generally based on a tensile curve). These static R–O parameters cannot accurately model the cyclic plasticity behavior. This work presents the cyclic R–O material hardening parameters for 316 stainless steel similar metal welds. The parameters were estimated under various conditions (in-air at room temperature, 300 °C in-air, and in-air at primary water conditions for a pressurized water reactor (PWR)). It is anticipated that the reported results would be useful for computational mechanics-based shakedown analysis and fatigue life estimation of PWR components.
Synergizing wire arc-directed energy deposition (WA-DED) additive manufacturing (AM) with particle-strengthened creep strength-enhanced ferritic (CSEF) steels enables fabrication and repair of critical power-plant components. Investigations focused on fusion-welded particle-strengthened CSEF steels, such as Grade 91 steel, have linked microstructurally heterogeneous regions—forming due to heat affected zones (HAZ)—with premature failure during elevated temperature service. Fusion-based AM, including WA-DED, likewise generates microstructurally and plastically heterogeneous regions due to spatiotemporally varying thermokinetics during deposition. However, works investigating such microstructural heterogeneities, their implications for mechanical behavior, and strategies to mitigate their formation remain scarce. This work identifies microstructurally and plastically heterogeneous regions within the WA-DED-processed Grade 91 steel. Spatial microhardness variations in the as-fabricated specimen correlate with the variation in the attributes of grain, martensitic microstructure, and precipitates across the fusion zone and HAZ. Digital image correlation-enabled tensile tests performed at 500 °C revealed pronounced deformation localization and a wave-like strain distribution, with wavelength close to the melt pool depth, indicating susceptibility of the as-fabricated components to premature creep failure. Such heterogeneity in microstructural and mechanical behavior was attributed to recurring solid-state phase transformations. Subsequently, an interlayer temperature control strategy was implemented, wherein maintaining interlayer temperature above the martensitic start temperature mitigated the heterogeneous microstructural and plastic response in the as-fabricated condition. Findings open pathways to achieving deformation-localization- and creep-resistant microstructures in WA-DED fabricated particle-strengthened CSEF steel components, reducing reliance on post-welding heat treatments—conventionally required to enhance creep resistance—and enabling on-demand, short lead-time fabrication of next-generation power-plant components.
This work describes three diffuse-interface methods for the simulation of immiscible, compressible multiphase fluid flows and elastic-plastic deformation in solids. The first method is the localized-artificial-diffusivity approach of Cook, Subramaniam et al., and Adler and Lele, in which artificial diffusion terms are added to the individual phase mass fraction transport equations and are coupled with the other conservation equations. The second method is the gradient-form approach that is based on the quasi-conservative method of Shukla et al., in which the diffusion and sharpening terms (together called regularization terms) are added to the individual phase volume fraction transport equations and are coupled with the other conservation equations. The third approach is the divergence-form approach that is based on the fully conservative method of Jain et al., in which the regularization terms are added to the individual phase volume fraction transport equations and are coupled with the other conservation equations. In the present study, all three diffuse-interface methods are used in conjunction with a four-equation, multicomponent mixture model, in which pressure and temperature equilibria are assumed among the various phases. The primary objective of this work is to compare these three methods in terms of their ability to: maintain constant interface thickness throughout the simulation; conserve mass, momentum, and energy; and maintain accurate interface shape for long-time integration. The second objective of this work is to consistently extend these methods to model interfaces between solid materials with strength. To assess and compare the methods, they are used to simulate a wide variety of problems, including (1) advection of an air bubble in water, (2) shock interaction with a helium bubble in air, (3) shock interaction and the collapse of an air bubble in water, and (4) Richtmyer–Meshkov instability of a copper–aluminum interface. The current work focuses on comparing these methods in the limit of relatively coarse grid resolution, which illustrates the true performance of these methods. In conclusion, this is because it is rarely practical to use hundreds of grid points to resolve a single bubble or drop in large-scale simulations of engineering interest.
This work quantifies the plastic deformation in terms of dislocation density for a series of cold sprayed Al-Cu coatings with copper contents from 2 mass% to 5 mass%. The success of the deposition and consolidation of the feedstock powders during the cold spray process rely on the occurrence of significant plastic deformation. Inert gas atomized Al-Cu alloy powders were sprayed onto substrates made of an Al-Cu-Mg-Mn alloy (AA2024) to produce dense coatings using a low pressure cold spray system with helium as the carrier gas. X-ray diffraction patterns were obtained from the Al-Cu feedstock powder material and the cold sprayed coatings using a monochromatic X-ray source, and the dislocation density was determined via an X-ray whole profile analysis. Increasing the Cu alloy content (from 2 mass% to 5 mass%) systematically increased the dislocation density in the Al-Cu coatings from (4.3 ± 0.5) × 10 14 m -2 to (7.5 ± 0.8) × 10 14 m -2 . The dislocation densities in the feedstock powders ranging from (0.4 ± 0.1) × 10 14 m -2 to (1.8 ± 0.2) × 10 14 m -2 were all lower than the dislocation densities in the corresponding coatings. The increasing deformation level in the Al-Cu coatings with Cu additions were confirmed by the classic and modified Williamson-Hall analyses of X-ray diffraction data, and peak breadth measurements from neutron diffraction data. Finally, a high density of dislocations was also observed in these coatings via electron backscatter diffraction and transmission electron microscopy.
Solid phase processing (SPP) is a promising alloy fabrication technique to produce fine and homogeneous grain structures for high-performance alloys. However, there is very limited modeling capability to understand and predict the grain refinement during SPP. In this work, the crystal plasticity theory was used to study elastic-plastic deformation in polycrystalline aluminums under large shear deformation. Two approaches, kernel averaged misorientation (KAM) and grain reference orientation deviation (GROD), were used to assess the grain misorientations. The geometrically necessary dislocation (GND) density was computed with the plastic strain rate. The deformation simulations were carried out under two loading conditions to investigate the effect of loading paths on the evolutions of grain misorientation and GND density. The results show that the regions with high misorientation and GND density first appear near grain boundaries. These regions then extend toward interior grains. The loading path affects dislocation system activation and dislocation recovery, hence dislocation evolution and misorientation. In conclusion, both two- and three-dimensional simulations showed that the spatial and temporal evolutions of GROD, KAM, and GND density in are closely correlated, which indicates they all can be used as criteria of grain refinement or recrystallization.
The advantages offered by additive manufacturing over traditional processes has driven a great deal of industrial and academic interest in recent years. However, the process is relatively new and requires additional investigation to become sufficiently mature for wide scale industrial adoption. Electron beam melting powder bed fusion is one technology that has shown promise for fabricating high temperature resistant materials such as nickel based superalloys. The resulting microstructures typically exhibit a strong fiber texture in the build direction giving rise to anisotropic time-dependent deformation behavior. In order to accelerate the qualification of these materials for industrial adoption accurate numerical models are needed for simulating their behavior. In this work a crystal plasticity model including non-Schmid effects is presented for capturing creep anisotropy observed in additively manufactured IN738LC. The model is calibrated via a probabilistic framework where model parameters are treated as random variables. An iterative sequential design strategy is utilized to efficiently identify the probability density of the unknown model parameters. As a case study the model is utilized to investigate the behavior of randomly oriented equiaxed grain clusters sometimes observed embedded in the additively manufactured columnar structure. A synthetic realization is simulated and uncertainty is propagated through to the full-field response. Results indicate that these features are the source of significant creep relaxation and strain accumulation which partially explains observed grain boundary decohesion at these locations.
This work utilized crystal plasticity simulations and experimental thermo-mechanical processing to examine how Zn and Ca content affects texture evolution in Mg–Zn–Ca alloys during simulated hot rolling. Four different compositions were studied: unalloyed Mg and three Mg alloys of ZX0p50 (Mg- 0.5 wt% Zn- 0.1 wt% Ca), ZX30 (Mg-3 wt% Zn- 0.1 wt% Ca), and ZX31 (Mg-3 wt% Zn- 0.3 wt% Ca). Multi-pass Gleeble experiments simulated hot rolling using plane strain conditions. Unalloyed Mg developed a strong basal texture, while the addition of Zn and Ca significantly weakened the texture in ZX30 and ZX31; this reduction was not seen in ZX0p50. Notably, ZX31 exhibited a split basal texture aligned with the rolling direction. Crystal plasticity simulations revealed that the weaker basal texture in ZX30 and ZX31 resulted from reduced activation of basal and twinning modes compared to unalloyed Mg and Zx0p50. Increased Zn and Ca content raised the ratios of basal to pyramidal critical resolved sheer stress (CRSS) and twin to pyramidal CRSS, enhancing pyramidal slip activity and weakening the basal texture. Extension twinning played a critical role in texture development during the first deformation pass. In unalloyed Mg and ZX0p50, twinning led to a strong basal texture, while ZX30 and ZX31 had weak basal textures that increased only slightly with further deformation. The split basal poles simulated in ZX31 were consistent with the experimental findings, highlighting the interplay between alloy composition and texture evolution.
The implementation of soil health-promoting practices, such as cover cropping and compost application, has important implications for nutrient cycling and management in agroecosystems. At the same time, plant belowground carbon (C) allocation patterns can influence nutrient cycling and availability in soil through changes to the microbial community, but the effects may depend on the crop genotype and management practices in place. We evaluated belowground C allocation patterns using 13 C labeling and root architecture in two genotypes of winter wheat (Triticum aestivum) with different levels of exudation and belowground allocation strategies in soils with contrasting compost amendment legacy (108.7 Mg ha -1 every 2 years over 10 years vs. no compost). We also measured microbial community structure and function in the rhizosphere and quantified uptake of residue-derived N from 15 N-labelled cover crop residues. We found an interactive effect between soil management and genotype, where in the no-compost soil, the high-exudation genotype (Snowmass) increased exudation by over 4-fold, while the low-exudate genotype (Byrd) increased only 2-fold. While we did not observe genotype differences in rhizosphere enzyme activity or dissolved N pools, residue N uptake was 1.8 times greater for Snowmass in the compost-amended soil. There were more rhizosphere microbial taxa associated with the high-exudate genotype (Snowmass); nine bacterial and seven fungal families were indicative of Snowmass, versus one bacterial and four fungal families for Byrd. Our results suggest that the high-exudation strategy can influence the rhizosphere microbial community, and lead to greater short-term residue N uptake in high SOM soil. By directly linking root architecture, exudation, microbial communities, and N mineralization and uptake dynamics, this work demonstrates that plasticity in root C allocation is genotype-specific and influences microbial communities and nutrient cycling depending on the soil health context.
Inspired by the parallelism and efficiency of the brain, several candidates for artificial synapse devices have been developed for neuromorphic computing, yet a nonlinear and asymmetric synaptic response curve precludes their use for backpropagation, the foundation of modern supervised learning. Spintronic devices—which benefit from high endurance, low power consumption, low latency, and CMOS compatibility—are a promising technology for memory, and domain-wall magnetic tunnel junction (DW-MTJ) devices have been shown to implement synaptic functions such as long-term potentiation and spike-timing dependent plasticity. In this work, we propose a notched DW-MTJ synapse as a candidate for supervised learning. Using micromagnetic simulations at room temperature, we show that notched synapses ensure the non-volatility of the synaptic weight and allow for highly linear, symmetric, and reproducible weight updates using either spin transfer torque (STT) or spin–orbit torque (SOT) mechanisms of DW propagation. We use lookup tables constructed from micromagnetics simulations to model the training of neural networks built with DW-MTJ synapses on both the MNIST and Fashion-MNIST image classification tasks. Accounting for thermal noise and realistic process variations, the DW-MTJ devices achieve classification accuracy close to ideal floating-point updates using both STT and SOT devices at room temperature and at 400 K. Our work establishes the basis for a magnetic artificial synapse that can eventually lead to hardware neural networks with fully spintronic matrix operations implementing machine learning.
Proposed apparatus cuts and polishes specimens while preventing contamination of outside environment or of subsequent specimens processed in it. Designed for use in zero gravity but also includes features useful in cutting and polishing of toxic or otherwise hazardous materials on Earth. Includes remote manipulator for handling specimens, cutting and polishing wire, inlets for gas and liquid, and outlets for waste liquid and gas. Replaceable plastic liner surrounds working space.
Abstract In this work we tailor the response of trusses loaded beyond their yield limit. The truss structures are modeled using finite strain theory and rate-independent elasto-plasticity. We design trusses with a tailored mechanical response that is between “elastic” and “elastic-ideal-plastic” subject to the volume constraint. The design updates are generated by the gradient-based Method of Moving Asymptotes (MMA) solver and the sensitivities of the response functions are computed using a path-dependent adjoint sensitivity analysis. The computations are performed in Matlab.
In this work, we present a model parameter calibration procedure for a physics-based crystal plasticity model. The calibration process utilizes a powerful statistics-based Bayesian calibration method. Calibration of the crystal plasticity parameters makes use of experimentally-measured data, i.e. compressive stress–strain response, from $\langle$100$\rangle$ and $\langle$123$\rangle$ single crystal copper dynamically loaded via Kolsky bar tests. The calibration of damage parameters is achieved using experimentally-measured free-surface velocity history data from plate impact test on $\langle$100$\rangle$ and $\langle$110$\rangle$ single crystal copper, which generates shock compression followed by dynamic tensile failure. A validation assessment is then carried out by comparing the calibrated model predictions and experimental measurements of the dynamic tensile damage generated in an impacted bicrystal copper plate. Lastly, a model-informed rationale for the experimentally-observed dependence of the spatial distribution of ductile damage (porosity) on crystallography is provided.