Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “elastic materials”

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 1,117 records · Page 62

Diffusion of Sn donors in β-Ga 2 O 3

Diffusion of the n-type dopant Sn in β-Ga 2 O 3 is studied using secondary-ion mass spectrometry combined with hybrid functional calculations. The diffusion of Sn from a Sn-doped bulk substrate with surface orientation (001) into an epitaxial layer is observed after heat treatments in the temperature range of 1050–1250 °C. Calculated formation energies of Sn-related and intrinsic defects show that the migration of Sn is mediated by Ga vacancies ( V Ga ) through the formation and dissociation of intermittent mobile V Ga Sn Ga complexes. The evolution of the Sn concentration vs depth profiles after heat treatments can be well described by a reaction–diffusion model. Using model parameters guided by the hybrid functional calculations, we extract a V Ga Sn Ga complex migration barrier of 3.0 ± 0.4 eV with a diffusion coefficient of 2 × 10 –1 cm 2 /s. The extracted migration barrier is consistent with our theoretical predictions using the nudged elastic band method, which shows migration barriers of 3.42, 3.15, and 3.37 eV for the [100], [010], and [001] directions, respectively.

36 MATERIALS SCIENCE↗

PyDislocDyn: A Python code for calculating dislocation drag and other crystal properties

PyDislocDyn is a suite of python programs designed to perform various calculations for dislocation dynamics in the continuum limit. In particular, one of its main purposes is to calculate dislocation drag from phonon wind. Additional features include the averaging of elastic constants for polycrystals, the calculation of the dislocation field including its limiting velocities, and the calculation of dislocation self-energy and line tension.

36 MATERIALS SCIENCE↗

A FFT-based mesoscale continuum dislocation mechanics with defect energy: Applications to composites and polycrystals

A crystal plasticity elastoviscoplastic FFT (fast Fourier transform) formulation with a mesoscale continuum field dislocation mechanics model is presented, which incorporates a defect energy density that depends on GND densities and an associated material length scale. This allows to thermodynamically derive internal length scale dependent intra-crystalline backstress and Peach–Koehler force acting on GND densities. The model considers GND density evolution through a filtered numerical spectral approach, which is coupled with stress equilibrium through the elastoviscoplastic FFT algorithm. The discrete Fourier transform (DFT) method together with finite difference (FD) schemes is applied to solve both the backstress tensor and the Fourier–Green operator. Numerical results are first reported for two-phase laminate composites with plastic single crystal channels and elastic precipitates for shear loadings. Channel size effects are simulated and analyzed on the overall and local hardening behaviors during monotonous loadings. In addition, the evolutions of GND densities and the role of their associated backstress on size effects are examined during reversible shear loading. In a second part, the role of the defect energy internal length scale on polycrystal’s hardening during tension–compression is discussed. The results are compared to those obtained using FFT-based continuum field dislocation mechanics without defect energy.

36 MATERIALS SCIENCE↗

First-Principles Studies on Sc 2 RuZ (Z = Si, Ge, Sn) Inverse Heusler Alloys: Structural, Electronic, and Transport Properties

The continuous demand for efficient, nontoxic, and thermally stable materials for room-temperature energy conversion motivates the exploration of novel thermoelectric systems beyond the traditional magnetic Heusler alloys. While full and half-Heusler compounds, especially Co-, Ni-, and Mn-based systems, have demonstrated promising thermoelectric properties, their typically high operating temperatures and magnetic complexities limit their applicability in ambient thermal management. In this context, we investigate whether Sc-based inverse Heusler alloys can offer a viable nonmagnetic alternative with competitive thermoelectric performance. In this work, we perform a systematic first-principles study of the inverse Heusler compounds Sc 2 RuZ (Z = Si, Ge, Sn), focusing on their structural, electronic, mechanical, and thermodynamic-thermoelectric properties. Density Functional Theory (DFT) was employed to compute optimized lattice structures and band dispersion, while dynamical stability was assessed via phonon calculations. Thermoelectric transport coefficients, including Seebeck coefficient, electrical conductivity, and thermal conductivity, were estimated using the semiclassical Boltzmann transport theory within the constant relaxation time approximation. Our results show that all Sc 2 RuZ compounds are thermodynamically stable semiconductors with indirect band gaps of 0.12–0.16 eV and exhibit high elastic moduli, especially Sc 2 RuSn, which demonstrates superior stiffness and incompressibility. Importantly, all compounds display promising room-temperature thermoelectric characteristics, including high Seebeck coefficients and power factors. These findings reveal that Sc 2 RuZ alloys represent a rare class of stable, nonmagnetic inverse Heusler semiconductors with intrinsic thermoelectric potential at room temperature, unlike many existing Heusler systems optimized for spintronics or high-temperature operation. This work expands the known design space for Heusler-based thermoelectrics and offers a theoretical basis for experimental realization of efficient, low-temperature, nonmagnetic thermoelectric materials.

alloys↗

Evolution of anisotropy with saturation and its implications for the elastoplastic responses of clay rocks

Many clay rocks have distinct bedding planes. Experimental studies have shown that their mechanical properties evolve with the degree of saturation (DOS), often with higher stiffness and strength after drying. For transversely isotropic rocks, the effects of saturation can differ between the bed-normal (BN) and bed-parallel (BP) directions, which gives rise to saturation-dependent stiffness and strength anisotropy. Accurate prediction of the mechanical behavior of clay rocks under partially saturated conditions requires numerical models that can capture the evolving elastic and plastic anisotropy with DOS. In this study, we present an anisotropy framework for coupled solid deformation-fluid flow in unsaturated elastoplastic media. Here we incorporate saturation-dependent strength anisotropy into an anisotropic modified Cam-Clay (MCC) model and consider the evolving anisotropy in both the elastic and plastic responses. The model was calibrated using experimental data from triaxial tests to demonstrate its capability in capturing strength anisotropy at various levels of saturation. Through numerical simulations, we demonstrate the role of evolving stiffness and strength anisotropy in the mechanical behavior of clay rocks. Plane strain simulations of triaxial compression tests were also conducted to demonstrate the impacts of material anisotropy and DOS on the mechanical and fluid flow responses.

36 MATERIALS SCIENCE↗

Microstructure-Sensitive modeling of surface roughness and notch effects on extreme value fatigue response

Fatigue life knockdown due to surface roughness induces local stress concentrations. Physical fatigue experiments are crucial to qualify materials for use but are limited in their capability to examine the mechanisms underlying the fatigue crack formation driving force due to these stress concentrations. Microstructure-sensitive computational models can aid in this regard and are employed in this work to examine the effects of surface roughness in Al 7075-T6. Digital microstructure models are subjected to crystal plasticity finite element method (CPFEM) simulations to examine the combined effects of microstructure and grain scale asperities on the driving force for fatigue crack formation. Following elastic–plastic shakedown, mesoscale volume-averaged fatigue indicator parameters (FIPs) are computed within fatigue damage process zones of grains. An increase in the intensity of realistic surface roughness profiles corresponds accordingly to larger FIPs but it is difficult to precisely interpret the mechanisms and effects of this intensification. We thus parametrically investigate surface asperities that couple with microstructure and lack of constraint on slip at the free surface to describe microstructure-sensitive surface roughness knockdown effects on fatigue resistance. The depth of a single valley (i.e., notch) due to surface roughness is found to be more detrimental to fatigue resistance than the radius (i.e., notch width, not to be confused with notch acuity). Finally, we also find that different crystallographic textures, grain morphologies, and number of sampled grains have a limited effect on the depth of influence of these notches.

36 MATERIALS SCIENCE↗

Novel ceramic capacitors with ultrahigh energy density and efficiency (Final Technical Report)

Antiferroelectric ceramics are a special class of material that have shown great potential as the dielectric in electrical capacitors due to their high energy- and power-density. During each charge-discharge cycle, the ceramic undergoes transformation to a ferroelectric phase and resumes its antiferroelectric phase. The hysteresis associated with the transitions leads to a mediocre energy efficiency and service lifetime of antiferroelectric capacitors and, hence, their almost absence in commercial products. Under the support of this research project, we first formulated a universal lattice-compatibility theory that included electrostatic polarization energy along with elastic energy and thermal energy to understand the origin of the hysteresis in antiferroelectric oxides. Guided by this compatibility theory, we conducted high-throughput density functional theory (DFT) calculations to assess chemical modifiers and their effect on crystal structures of 400+ PbZrO 3 -based compositions. Down-selected compositions were experimentally validated for their suppressed hysteresis and higher energy efficiency. The verified low-hysteresis compositions were then expanded to an antiferroelectric ceramic library with nearly 500 new compositions (more than 1,500 samples) using high-throughput experiments involving ceramic synthesis and property screening. The large quantity of data generated (theory and experimental) in these tasks were processed by machine-learning techniques and identified trends were fed to the next iteration. In the end, we successfully discovered four compositions with near-zero hysteresis, yielding a world-record energy efficiency of 98.2% at an energy density of 3.0 J/cm 3 . Furthermore, our antiferroelectric ceramic capacitor reaches 79.5 million charge-discharge cycles lifetime, a factor of 80 enhancement over previous antiferroelectric ceramics with large hysteresis. These research accomplishments have not only met the milestones set in the SOPO, but also led to two patent filings, three journal publications (one of them was in Advanced Materials, impact factor 29.4), and nine oral presentations at various venues. Through the course of the project, three postdocs, four Ph.D. students, and one M.S. student were trained. In short, our project established a new methodology in searching next-generation functional ceramics on the fundamental side and discovered several high-efficiency antiferroelectric compositions for capacitors on the applied side. Once fabricated into the multilayer form for commercial applications, these ceramic capacitors can potentially enable the high temperature high power density DC-link capacitors that are critical for the next generation inverters in electric vehicles. The project also significantly contributed to the nation’s workforce development in the STEM fields.

36 MATERIALS SCIENCE↗

The fate of liquid crystal topological defects on chemically patterned surfaces during phase transitions

Controlling topological defects in liquid crystals (LCs) is an essential element in the development of areas such as directed self-assembly and micropatterning materials. However, during the phase transition on confined patterned surfaces, how the morphologies in one liquid crystalline phase change from deformations or defects into another phase is much less known. Here, we examine the fate of defects in a LC confined on a patterned surface during smectic-A–nematic and nematic–isotropic phase transitions, using experiments and simulation analyses. Upon heating from smectic-A to nematic, a Toric focal conical domain (TFCD) melts into a +1 converging boojum defect, which then transitioned into a concentric configuration as temperature increases, attributed to a steeper decrease of the bend and twist modulus compared to splay modulus. During cooling, TFCDs are developed from two distinct pathways depending on the cooling rates. Our continuum simulation recapitulates these transformations and provides elastic constant-based explanations for the two pathways. Although the phase transition pathways of defects are independent of the geometry of the confined patterns, the arrangement of FCDs is highly dependent on the size and shape of the patterns. Finally, taken together, this simple approach offers promising opportunities for tuning the micro- or nano-patterning of topological defects in liquid crystals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Hydrogen Localization and Cluster Formation in α-Zr from First-Principles Investigations

We report alpha-zirconium (α-Zr) is a hydrogen atom host material used in nuclear energy and semiconducting applications due to Zr's excellent properties such as corrosion resistance, low thermal neutrons cross-section, and thermodynamically favorable interstitial sequestering of H atoms. Herein we use density functional theory simulations to systematically analyze the initial stages of H insertion in the Zr matrix, and determine optimal hydrogen localization, distribution, and clustering tendencies from a thermodynamics and kinetics perspective. Our results suggest that H occupation of tetrahedral sites is strongly favored, consistent with previous experimental and theoretical studies. Further, a cooperative effect in H localization is shown where H atoms favor clustering. The H pair is found to be the most stable by occupying tetrahedral sites with pairs at second nearest-neighbor positions (TTn2). Among the four H atom clusters, the most favorable configuration has two TTn2 pairs forming "box " or "zig-zag" shaped-cluster with the pairs separated by 3.03 or 3.25 Å, respectively. Nudged elastic band simulations suggest that the formation of TTn2 is facile and is limited by the small diffusion barrier of the monomer of less than 0.1 eV. Further, the formation routes for the 4H atom clusters are also presented. The current study can serve as a foundation to systematically model H atom clustering in metals and alloys.

36 MATERIALS SCIENCE↗

Prediction of carbon nanostructure mechanical properties and the role of defects using machine learning

Graphene-based nanostructures hold immense potential as strong and lightweight materials, however, their mechanical properties such as modulus and strength are difficult to fully exploit due to challenges in atomic-scale engineering. This study presents a database of over 2,000 pristine and defective nanoscale CNT bundles and other graphitic assemblies, inspired by microscopy, with associated stress–strain curves from reactive molecular dynamics (MD) simulations using the reactive INTERFACE force field (IFF-R). These 3D structures, containing up to 80,000 atoms, enable detailed analyses of structure-stiffness-failure relationships. By leveraging the database and physics- and chemistry-informed machine learning (ML), accurate predictions of elastic moduli and tensile strength are demonstrated at speeds 1,000 to 10,000 times faster than efficient MD simulations. Hierarchical Graph Neural Networks with Spatial Information (HS-GNNs) are introduced, which integrate chemistry knowledge. HS-GNNs as well as extreme gradient boosted trees (XGBoost) achieve forecasts of mechanical properties of arbitrary carbon nanostructures with only 3 to 6% mean relative error. The reliability equals experimental accuracy and is up to 20 times higher than other ML methods. Predictions maintain 8 to 18% accuracy for large CNT bundles, CNT junctions, and carbon fiber cross-sections outside the training distribution. The physics- and chemistry-informed HS-GNN works remarkably well for data outside the training range while XGBoost works well with limited training data inside the training range. The carbon nanostructure database is designed for integration with multimodal experimental and simulation data, scalable beyond 100 nm size, and extendable to chemically similar compounds and broader property ranges. The ML approaches have potential for applications in structural materials, nanoelectronics, and carbon-based catalysts.

Winetrout, Jordan J.↗

A phase-field model of quasi-brittle fracture for pressurized cracks: Application to UO 2 high-burnup microstructure fragmentation

Here, in this paper, we present a phase-field model of quasi-brittle fracture with pressurized cracks, with dedicated applications for polycrystalline materials. The model is formulated as a minimization problem within the variational framework. The external work done by pressure on the crack surfaces is included in the objective function. Several careful modeling choices lead to a regularization-length-independent critical strength. The model is constructed to give a softening response with an underlying linear traction-separation law. The pressure-dependent softening response and the regularization of the prescribed pressure are demonstrated with a (quasi) one-dimensional numerical analysis. In a two-dimensional numerical analysis under plane strain assumptions, the critical stress corresponding to crack propagation (as predicted by our quasi-brittle fracture model) is compared with linear elastic fracture mechanics (LEFM) analytical solutions. Our model is further utilized to simulate fission-gas-induced fragmentation of the UO 2 high-burnup structure (HBS). Simulation results show that pressurized bubbles can cause crack nucleation and propagation, and that the bubble size and the surrounding external pressure affect the critical pressure corresponding to crack nucleation. Simulations of a partial HBS at different recrystallization stages show that different grain structures (due to recrystallization) also influence crack paths and fragmentation morphology.

36 MATERIALS SCIENCE↗

Endowing low fatigue for elastocaloric effect by refined hierarchical microcomposite in additive manufactured NiTiCuCo alloy

Abstract NiTiCu-based shape memory alloys have been considered as ideal materials for solid-state refrigeration due to their superb cycling stability for elastocaloric effect. However, the embrittlement and deterioration caused by secondary phase and coarse grains restrict their applications, and it is still challenging since the geometric components are required. Here, bulk NiTiCuCo parts with excellent forming quality were fabricated by laser powder bed fusion (LPBF) technique. The as-fabricated alloy exhibits refined three-phases hierarchical microcomposite formed based on the rapid cooling mode of LPBF, composed of intricate dendritic Ti 2 Ni–NiTi composite and nano Ti 2 Cu embedded inside the NiTi-matrix. This configuration endows far superior elastocaloric stability compared to the as-cast counterpart. The low fatigue stems from the strong elastic coupling between the interphases with reversible martensite transformation, revealed byin-situsynchrotron high-energy x-ray diffraction. The fabrication of NiTiCuCo alloy via LPBF fills the bill of complex geometric structures for elastocaloric NiTiCu alloys. The understanding of interphase micro-coupling could provide the guide for designing LPBF fabricated shape memory-based composites, enabling their applications for special demands on other functionalities.

Engineering↗

Interfacial toughness: dependence on surface roughness and test temperature

Interfacial toughness quantifies resistance to crack growth along an interface and in this investigation the toughness of an aluminum/epoxy interface was measured as a function of surface roughness and test temperature. The large strain response of the relatively ductile epoxy adhesive used in this study was also characterized. This epoxy adhesive exhibits intrinsic strain-softening after initial compressive yield and then deforms plastically at a roughly constant flow stress until it rapidly hardens at large compressive strains. Here, we find that interface toughness scales as the product of the temperature dependent epoxy yield strength and a length scale that characterizes surface roughness. The proposed scaling is based upon dimensional considerations of a model problem that assumes that the characteristic length scale of both the roughness and the crack-tip yield zone is small relative to the region dominated by the linear elastic asymptotic crack-tip stress field. Furthermore, the model assumes that interfacial failure occurs only after the epoxy begins to harden at large strains. The proposed relationship is validated by our interfacial toughness measurements.

36 MATERIALS SCIENCE↗

Micromechanical Fields Associated with Irregular Deformation Twins in Magnesium

Understanding and controlling the development of deformation twins is paramount for engineering strong and stable hexagonal close-packed (HCP) Mg alloys. Actual twins are often irregular in boundary morphology and twin crystallography, deviating from the classical picture commonly used in theory and simulation. In this work, the elastic strains and stresses around irregular twins are examined both experimentally and computationally to gain insight into how twins develop and the microstructural features that influence their development. A nanoprecession electron diffraction (N-PED) technique is used to measure the elastic strains within and around a $\left\{ {10\overline{1}2} \right\}$ tensile twin in AZ31B Mg alloy with nm scale resolution. A full-field elasto-viscoplastic fast Fourier transform (EVP-FFT) crystal plasticity model of the same sub-grain and irregular twin structure is employed to understand and interpret the measured elastic strain fields. The calculations predict spatially resolved elastic strain fields in good agreement with the measurement, as well as all the stress components and the dislocation density fields generated by the twin, which are not easily obtainable from the experiment. The model calculations find that neighboring twins, several twin thicknesses apart, have little influence on the twin-tip micromechanical fields. Furthermore, this work reveals that irregularity in the twin-tip shape has a negligible effect on the development of the elastic strains around and inside the twin. Importantly, the major contributor to these micromechanical fields is the alignment of the twinning shear direction with the twin boundary.

36 MATERIALS SCIENCE↗

Boron nitride neutron detector with the ability for detecting both thermal and fast neutrons

The detection of fast neutrons is regarded technically challenging because the interaction probability of fast neutron with matter is extremely low. Based on our recent development of hexagonal boron nitride (BN) semiconductor thermal neutron detectors with a record high efficiency of 59%, we report here the feasibility studies of BN detectors for detecting fast neutrons. A BN detector with a detection area of 2.1 cm 2 was fabricated from a 90 $μ$m thick BN epilayer. In the presence of a bare Cf-252 source emitting fast neutrons ranging from 1 to 9 MeV, the detection efficiency was estimated to be about 0.1%. Furthermore, the measured mean free path of fast neutron in BN is about 7.6 cm. Together with the capability of BN for thermal neutron detection, the present results indicate that by incorporating BN with a large thickness, BN neutron detectors are expected to possess the unique capability of directly detecting thermal to fast neutrons as well as outstanding features resulting from the ultrawide bandgap of BN. The identification of a single material that is sensitive to both thermal and fast neutrons is valuable for the development of novel neutron detection technologies.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Toroidal polar topology in strained ferroelectric polymer

Polar topological texture has become an emerging research field for exotic phenomena and potential applications in reconfigurable electronic devices. Here, we report toroidal topological texture self-organized in a ferroelectric polymer, poly(vinylidene fluoride-ran-trifluoroethylene) [P(VDF-TrFE)], that exhibits concentric topology with anticoupled chiral domains. The interplay among the elastic, electric, and gradient energies results in continuous rotation and toroidal assembly of the polarization perpendicular to polymer chains, whereas relaxor behavior is induced along polymer chains. Such toroidal polar topology gives rise to periodic absorption of polarized far-infrared (FIR) waves, enabling the manipulation of the terahertz wave on a mesoscopic scale. Our observations should inform design principles for flexible ferroic materials toward complex topologies and provide opportunities for multistimuli conversions in flexible electronics.

36 MATERIALS SCIENCE↗

Ion Transport in Concentrated Crosslinked Solid Polymer Electrolytes

Crosslinking polymers is a common approach to create mechanically stable solid materials such as polymer electrolytes for lithium batteries. In conventional liquid electrolytes, the solvent molecules move freely to accommodate the field-induced motion of ions. However, in crosslinked polymer electrolytes, the rearrangement of polymer segments is constrained by the deformation limits of the network. Herein, we develop a new transport model that accounts for both the formation of concentration gradients and the elasticity of the electrolyte. The elasticity is incorporated by adding an additional term related to the entropy of crosslinked strands to the electrochemical potential of the salt. The resulting Crosslink Model contains two adjustable parameters: $\mathcal{N}$, the average number of monomers in a strand, and λ crit , the maximum strain the network can sustain. These solid-like constraints produce singularities in the governing transport equations, fundamentally altering the concentration profiles. Plateaus in salt concentrations emerge near the electrodes, and network elasticity introduces a strain overpotential. When compared to a Baseline Model ($\mathcal{N}$ → ∞, equivalent to concentrated solution theory), which predicts steepest gradients near the electrodes, both models yield similar current–voltage relationships. Model predictions are compared to electrochemical data for a poly(ethylene oxide)-based crosslinked polymer electrolyte.

Patel, Vivaan [University of California, Berkeley,↗