Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite differences”

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,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Role of tilt grain boundaries on the structural integrity of WSe 2 monolayers

Transition metal dichalcogenides (TMDCs) are potential materials for future optoelectronic devices. Grain boundaries (GBs) can significantly influence the optoelectronic properties of TMDC materials. Here, we have investigated the mechanical characteristics of tungsten diselenide (WSe 2 ) monolayers and failure process with symmetric tilt GBs using ReaxFF molecular dynamics simulations. In particular, the effects of topological defects, loading rates, and temperatures are investigated. We considered nine different grain boundary structures of monolayer WSe 2 , of which six are armchair (AC) tilt structures, and the remaining three are zigzag (ZZ) tilt structures. Our results indicate that both tensile strength and fracture strain of WSe 2 with symmetric tilt GBs decrease as the temperature increases. We revealed an interfacial phase transition for high-angle GBs reduces the elastic strain energy within the interface at finite temperatures. Furthermore, brittle cracking is the dominant failure mode in the WSe 2 monolayer with tilted GBs. WSe 2 GB structures showed more strain rate sensitivity at high temperatures than at low temperatures.

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

Analysis of Dynamical Peculiarities in Nanoalloys at Subsystems Level: Dynamical Degrees of Freedom, Temperature Differences, and the Chameleon Effect

Abstract A novel analysis of the dynamical behavior of nanoalloy systems, as represented by model Ni/Al 13‐atom clusters, over a broad range of energies that cover the stage‐wise transition of the systems from their solid‐like to liquid‐like state is presented. Conceptually, the analysis is rooted in partitioning the systems into judiciously chosen subsystems and characterizing the latter in terms of subsystem‐specific dynamical descriptors that include dynamical degrees of freedom , root‐mean‐square bond‐length fluctuation , and element‐specific subsystem temperature . The analysis reveals a host of intriguing new peculiarities in the dynamical behavior of the Ni/Al 13‐mers, among which are what we call the chameleon effect and the difference in the temperatures of the Ni and Al subsystems at high energies, a difference that strongly depends on the cluster composition and also changes with energy. These do not have an analog in pure Ni 13 and Al 13 and are explained in terms of the coupled effects of the difference between the masses of the Ni and Al atoms (the mass effect) and of the difference in the anharmonicity of the overall interaction potential as experienced by the Ni and Al subsystems of the clusters (the potential effect).

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

Energy-relaxation pathways in finite-sized artificial square-ice structures

The energy-relaxation pathway in artificial spin ice structures refers to the series of individual moment reorientations the system undergoes as it relaxes from an initial state to a final low-energy state. The probability of the artificial spin ice structure reaching a ground state depends on the energy of the intermediate states as approximated by an Arrhenius-type relation. In this work, we report that by modifying the geometry of artificial square-ice structures, we can influence the probability of observing the ground state in different configurations. Furthermore, our findings are supported by dipolar energy calculations for various energy-relaxation pathways and the corresponding experimental data obtained by magnetic force microscopy.

Analog computation↗
DOE OSTI2022

Density functional thermodynamic description of spin, phonon and displacement degrees of freedom in antiferromagnetic-to-paramagnetic phase transition in YNiO 3

This work herein demonstrates a direct density functional description of the finite-temperature thermodynamic properties of solids exhibiting phase transitions through positional and spin symmetry breaking degrees of freedom. A classic example addressed here is the rare-earth (R) nickelates RNiO 3 where the ground state is characterized by crystallographic and magnetic (e.g., antiferromagnetic) long-range order (LRO), whereas the higher temperature paramagnetic phase manifests a range of local spin and positional symmetry breaking motifs with short-range order (SRO). Unlike time-dependent simulations of spin and positional degrees of freedom, in the present work, phases are described via a superposition of static configurations constructed by populating a periodic base lattice supercell allowing for the formation of energy lowing distribution of positional and spin local motifs. The thermal populations of the configurations in such a superposition phase are obtained from the energy-minimized Density Functional Theory (DFT)-calculated partition functions at different temperatures. This approach offers flexible inclusion of different physical contributions to the free energy, such as elastic, electronic and phonon free energies, all obtained from the same underlying DFT total energy calculations of periodic structures. The thermodynamic and magnetic properties of both LRO and SRO crystallographic and spin phases, including antiferromagnetic (AFM) to paramagnetic (PM) Néel phase transition in YNiO3 are studied. Including spin and phonon contributions, we find a DFT-calculated Néel temperature to be 144 K in satisfactory agreement with the experimental value of 145 K; whereas omitting the phonon contribution, one obtains a Néel temperature of 81 K. We present phonon contributions to the DFT-calculated temperature-dependent SRO, heat capacities, and the polymorphous distribution of nonzero local magnetic moments in the PM phase. This approach thus extends to finite temperatures the symmetry-broken DFT description of both the AFM and PM phases, demonstrating that a thermodynamic superposition approach based on symmetry broken configurations evaluated by a mean-field like DFT is sufficient to obtain a consistent description of the thermal physics of the AFM, PM phases and their interconversion in 3d oxides illustrated by YNiO 3 .

36 MATERIALS SCIENCE↗
DOE OSTI2021

Plastic work constrained elastoplastic topology optimization

An elastoplastic topology optimization framework for limiting plastic work generation while maximizing stiffness is presented. The kinematics and constitutive model are based on finite strain linear isotropic hardening plasticity, and the balance laws are solved using a total Lagrangian finite element formulation. Aggregation of the specific plastic work combined with an adaptive normalization scheme efficiently constrains the maximum specific plastic work. The optimization problem is regularized using an augmented partial differential equation filter, and is solved by the method of moving asymptotes where path-dependent sensitivities are derived using the adjoint method. The numerical examples show a clear dependence on the optimized maximum stiffness structures for different levels of constrained specific plastic work. It is also shown that due to the history dependency of the plasticity, the load path significantly influences the structural performance and optimized topology.

42 ENGINEERING↗
DOE OSTI2023

A ductility metric for refractory-based multi-principal-element alloys

We propose a quantum-mechanical dimensionless metric, the local-lattice distortion (LLD), as a reliable predictor of ductility in refractory multi-principal-element alloys (RMPEAs). The LLD metric is based on electronegativity differences in localized chemical environments and combines atomic-scale displacements due to local lattice distortions with a weighted average of valence-electron count. To evaluate the effectiveness of this metric, we examined body-centered cubic (bcc) refractory alloys that exhibit ductile-to-brittle behavior. Our findings demonstrate that local-charge behavior can be tuned via composition to enhance ductility in RMPEAs. With finite-sized cell effects eliminated, the LLD metric accurately predicted the ductility of arbitrary alloys, which compares well with existing tensile-elongation experiments. To validate further, we qualitatively evaluated the ductility of two refractory RMPEAs, i.e., NbTaMoW and Mo 72 W 13 Ta 10 Ti 2.5 Zr 2.5 , through the observation of crack formation under indentation, again showing excellent agreement with LLD predictions. Additionally, a comparative study of three refractory alloys provides further insights into the electronic-structure origin of ductility in refractory RMPEAs. This proposed metric enables rapid and accurate assessment of ductility behavior in the vast RMPEA composition space.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Ensemble variational Fokker-Planck methods for data assimilation

Particle flow filters solve Bayesian inference problems by smoothly transforming a set of particles into samples from the posterior distribution. Particles move in state space under the flow of an McKean-Vlasov-Itˆo process. This work introduces the Variational Fokker-Planck (VFP) framework for data assimilation, a general approach that includes previously known particle flow filters as special cases. The McKean-Vlasov-Itˆo process that transforms particles is defined via an optimal drift that depends on the selected diffusion term. It is established that the underlying probability density - sampled by the ensemble of particles - converges to the Bayesian posterior probability density. For a finite number of particles the optimal drift contains a regularization term that nudges particles toward becoming independent random variables. Based on this analysis, we derive computationally-feasible approximate regularization approaches that penalize the mutual information between pairs of particles, and avoid particle collapse. Moreover, the diffusion plays a role akin to a particle rejuvenation approach that aims to alleviate particle collapse. The VFP framework is very flexible. Different assumptions on prior and intermediate probability distributions can be used to implement the optimal drift, and localization and covariance shrinkage can be applied to alleviate the curse of dimensionality. A robust implicit-explicit method is discussed for the efficient integration of stiff McKean- Vlasov-Itˆo processes. Here, the effectiveness of the VFP framework is demonstrated on three progressively more challenging test problems, namely the Lorenz ’63, Lorenz ’96 and the quasi-geostrophic equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Half-closed discontinuous Galerkin discretisations

Here we introduce the concept of half-closed nodes for nodal discontinuous Galerkin (DG) discretisations. Unlike more commonly used closed nodes in DG, where on every element nodes are placed on all of its boundaries, half-closed nodes only require nodes to be placed on a subset of the element's boundaries. The effect of using different nodes on DG operator sparsity is studied and we find in particular for there to be no difference in the sparsity pattern of the Laplace operator whether closed or half-closed nodes are used. On quadrilateral/hexahedral elements we use the Gauss-Radau points as the half-closed nodes of choice, which we demonstrate is able to speed up DG operator assembly in addition to leverage previously known superconvergence results. We also discuss in this work some linear solver techniques commonly used for Finite Element or discontinuous Galerkin methods such as static condensation and block-based methods, and how they can be applied to half-closed DG discretisations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Engineering Thermal Transport across Layered Graphene–MoS 2 Superlattices

Layering two-dimensional van der Waals materials provides a high degree of control over atomic placement, which could enable tailoring of vibrational spectra and heat flow at the sub-nanometer scale. In this work, using spatially resolved ultrafast thermoreflectance and spectroscopy, we uncover the design rules governing cross-plane heat transport in superlattices assembled from monolayers of graphene (G) and MoS 2 (M). Using a combinatorial experimental approach, we probe nine different stacking sequences, G, GG, MG, GGG, GMG, GGMG, GMGG, GMMG, and GMGMG, and identify the effects of vibrational mismatch, interlayer adhesion, and junction asymmetry on thermal transport. Pure G sequences display evidence of quasi-ballistic transport, whereas adding even a single M layer strongly disrupts heat conduction. The experimental data are described well by molecular dynamics simulations, which include thermal expansion, accounting for the effect of finite temperature on the interlayer spacing. The simulations show that an increase of ~2.4% in the layer separation of GMGMG, relative to its value at 300 K, can lead to a doubling of the thermal resistance. Using these design rules, we experimentally demonstrate a five-layer GMGMG superlattice “thermal metamaterial” with an ultralow effective cross-plane thermal conductivity comparable to that of air.

2D materials↗
DOE OSTI2025

Cabibbo-Kobayashi-Maskawa unitarity deficit reduction via finite nuclear size

We revisit the extraction of the |𝑉 𝑢⁢𝑑 | Cabibbo-Kobayashi-Maskawa (CKM) matrix element from the superallowed transition decay rate of 26⁢𝑚 Al → 26 Mg, focusing on finite nuclear size effects. The decay rate dependence on the 26⁢𝑚 Al charge radius is found to be four times higher than previously believed, necessitating precise determination. However, for a short-lived isotope of an odd 𝑍 element such as 26⁢𝑚 Al , radius extraction relies on challenging many-body atomic calculations. We performed the needed calculations, finding an excellent agreement with previous ones, which used a different methodology. This sets a new standard for the reliability of isotope shift factor calculations in many-electron systems. The ℱ⁡𝑡 value obtained from our analysis is lower by 2.2𝜎 than the corresponding value in the previous critical survey, resulting in an increase in |𝑉 𝑢⁢𝑑 | 2 by 0.9𝜎. Adopting |𝑉 𝑢⁢𝑑 | from this decay alone reduces the CKM unitarity deficit by one standard deviation, irrespective of the choice of |𝑉 𝑢⁢𝑠 |.

beta-decay↗
DOE OSTI2020

Near, far, wherever you are: simulations on the dose efficiency of holographic and ptychographic coherent imaging

Different studies in X-ray microscopy have arrived at conflicting conclusions about the dose efficiency of imaging modes involving the recording of intensity distributions in the near (Fresnel regime) or far (Fraunhofer regime) field downstream of a specimen. A numerical study is presented on the dose efficiency of near-field holography, near-field ptychography and far-field ptychography, where ptychography involves multiple overlapping finite-sized illumination positions. Unlike what has been reported for coherent diffraction imaging, which involves recording a single far-field diffraction pattern, it is found that all three methods offer similar image quality when using the same fluence on the specimen, with far-field ptychography offering slightly better spatial resolution and a lower mean error. These results support the concept that (if the experiment and image reconstruction are done properly) the sample can be near or far; wherever you are, photon fluence on the specimen sets one limit to spatial resolution.

42 ENGINEERING↗
DOE OSTI2023

BARC Element Classifier (barc)

SAND2024-02120O This program enumerates and classifies all possible functional elements in several different categories within the Ballistic Asynchronous Reversible Computing (barc) model of computation. barc software is a research tool that is being used to help document the possible digital behaviors of primitive functional elements in the ABRC a.k.a. BARC model of computation. The method of operation of this software leverages elementary concepts and methods from discrete mathematics: combinatorics (permutations), partitions and equivalence classes, symmetry transformations, and the theory of finite groups. Implementation is in Python programming language. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Frank, Michael↗
DOE OSTI2024

Numerical modeling and experimental validation of low velocity impact of woven GFRP/CFRP composites

Low-velocity impact of 2D woven glass fiber reinforced polymer (GFRP) and carbon fiber reinforced polymer (CFRP) composite laminates was studied experimentally and numerically. Hybrid laminates containing blocked layers of GFRP/CFRP/GFRP with all plies oriented at 0° were investigated. Relatively high impact energies were used to obtain full perforation of the laminate in a low-velocity impact setup. Numerical simulations were carried out using the in-house transient dynamics finite element code, Sierra/SM, developed at Sandia National Laboratories. A three-dimensional continuum damage model was used to describe the response of a woven composite ply. Two methods for handling delamination were considered and compared: (1) cohesive zone modeling and (2) continuum damage mechanics. The reduced model size achieved by omission of the cohesive zone elements produced acceptable results at reduced computational cost. Further, the comparison between different modeling techniques can be used to inform modeling decisions relevant to low velocity impact scenarios. The modeling was validated by comparing with the experimental results and showed good agreement in terms of predicted damage mechanisms and impactor velocity and force histories.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Statistical methods for design and testing of 3D-printed polymers

Different statistical methods are used in various fields to qualify processes and products, especially in emerging technologies like Additive Manufacturing (AM) or 3D printing. Since several statistical methods are being employed to ensure quality production of the 3D-printed parts, an overview of these methods used in 3D printing for different purposes is presented in this paper. The advantages and challenges, to understanding the importance it brings for design and testing optimization of 3D-printed parts are also discussed. The application of different metrology methods is also summarized to guide future researchers in producing dimensionally-accurate and good-quality 3D-printed parts. Here, this review paper shows that the Taguchi Methodology is the commonly-used statistical tool in optimizing mechanical properties of the 3D-printed parts, followed by Weibull Analysis and Factorial Design. In addition, key areas such as Artificial Intelligence (AI), Machine Learning (ML), Finite Element Analysis (FEA), and Simulation require more research for improved 3D-printed part qualities for specific purposes. Future perspectives are also discussed, including other methods that can help further improve the overall quality of the 3D printing process from designing to manufacturing.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Parametric estimation of Poisson's ratio for thin hinged-hinged plates

Cornu's method is an elegant calculation besieged by an impractical approach to obtain accurate estimates for Poisson's ratio. Conventionally, Cornu's method requires several components, and each component can adversely affect the accuracy of the measurement. Furthermore, Cornu's conventional method requires a long beam because beams with short length-to-width ratios cause the estimate of Poisson's ratio to diverge from the true value of Poisson's ratio. We believe that, with the right modifications, Cornu's method can become an attractive approach to obtaining precise estimates for Poisson's ratio from mode shapes. Here we use finite element simulations to show how to use Cornu's method to estimate Poisson's ratio from a mode shape. Our modified Cornu's method removes knife-edges and loading components for a hinged-hinged plate under steady-state excitation. Given the true value of Poisson's ratio, as a performance specification, we show that simple parametric expressions can fit estimates for Poisson's ratio for different length-to-width ratios of thin hinged-hinged plates. Additionally, we show that estimates for Poisson's ratio from higher modes align with the results from the first mode and explain our expectation for this outcome. Furthermore, our results challenge the idea that anticlastic, monoclastic, and synclastic deformation uniquely correspond to positive, zero, and negative estimates of Poisson's ratio, respectively. With the emergence of materials-by-design, we expect that this parametric technique will be able to assist in experimental qualification of thin beam and plate structures with respect to the desired value of Poisson's ratio.

42 ENGINEERING↗
DOE OSTI2020

A linearity preserving nodal variation limiting algorithm for continuous Galerkin discretization of ideal MHD equations

A stabilized continuous Galerkin (CG) method for magnetohydrodynamics (MHD) is presented herein. Ideal, compressible inviscid MHD equations are discretized in space on unstructured meshes using piecewise linear or bilinear finite element bases to get a semi-discrete scheme. Stabilization is then introduced to the semi-discrete method in a strategy that follows the algebraic flux correction paradigm. This involves adding some artificial diffusion to the high order, semi-discrete method and mass lumping in the time derivative term. The result is a low order method that provides local extremum diminishing properties for hyperbolic systems. The difference between the low order method and the high order method is scaled element-wise using a limiter and added to the low order scheme. The limiter is solution dependent and computed via an iterative linearity preserving nodal variation limiting strategy. The stabilization also involves an optional consistent background high order dissipation that reduces phase errors. The resulting stabilized scheme is a semi-discrete method that can be applied to inviscid shock MHD problems and may be even extended to resistive and viscous MHD problems. To satisfy the divergence free constraint of the MHD equations, we add parabolic divergence cleaning to the system. Various time integration methods can be used to discretize the scheme in time. We demonstrate the robustness of the scheme by solving several shock MHD problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Embedding PbS Quantum Dots (QDs) in Pb-Halide Perovskite Matrices: QD Surface Chemistry and Antisolvent Effects on QD Dispersion and Confinement Properties

Hybrid materials of metal chalcogenide colloidal quantum dots (QDs) embedded in metal halide perovskites (MHPs) have led to composites with synergistic properties. Here, we investigate how QD size, surface chemistry, and MHP film formation methods affect the resulting optoelectronic properties of QD/MHP “dot-in-matrix” systems. We monitor the QD absorption and photoluminescence throughout synthesis, ligand exchange, and transfer into the MHP ink, and we characterize the final QD/MHP films via electron microscopy and transient absorption. In addition, we are the first to globally map how PbS QDs are distributed on the micrometer scale within these dot-in-matrix systems, using three-dimensional (3D) tomography time-of-flight secondary ion mass spectrometry. The surface chemistry imparted during synthesis directly affects the optical properties of the dot-in-matrix composites. Pb-halide passivation leads to QD/MHP dot-in-matrix samples with optical properties that are well-described by a theoretical model, based on a Type I finite-barrier heterostructure between the PbS QD and the MHP matrix. Samples without Pb-halide passivation show complicated size-dependent behavior, indicating a transition from a Type I heterostructure between the PbS QD wells and MHP barriers for small-sized QDs to PbS QDs that are electronically decoupled from the MHP matrix for larger QDs. Furthermore, the choice in perovskite antisolvent crystallization method leads to a difference in the spatial QD distribution within the perovskite matrix, differences in carrier lifetime, and photoluminescence shifts of up to 180 meV for PbS in methylammonium lead iodide. This work establishes an understanding of such emerging synergistic systems relevant for technologies such as photovoltaics, infrared emitters and detectors, and other unexplored technological applications.

92 SC - Basic Energy Sciences (BES) (SC-22)↗
DOE OSTI2022

Probing the axion–nucleon coupling with the next generation of axion helioscopes

Abstract A finite axion–nucleon coupling, nearly unavoidable for QCD axions, leads to the production of axions via the thermal excitation and subsequent de-excitation of $$^{57}$$ 57 Fe isotopes in the sun. We revise the solar bound on this flux adopting the up to date emission rate, and investigate the sensitivity of the proposed International Axion Observatory IAXO and its intermediate stage BabyIAXO to detect these axions. We compare different realistic experimental options and discuss the model dependence of the signal. Already BabyIAXO has sensitivity far beyond previous solar axion searches via the nucleon coupling and IAXO can improve on this by more than an order of magnitude.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗