Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “orientational order”

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 253 records · Page 14

Fast and accurate reduced-order modeling of a MOOSE-based additive manufacturing model with operator learning

One predominant challenge in additive manufacturing (AM) is to achieve specific material properties by manipulating manufacturing process parameters during the runtime. Such manipulation tends to increase the computational load imposed on existing simulation tools employed in AM. The goal of the present work is to construct a fast and accurate reduced-order model (ROM) for an AM model developed within the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework, ultimately reducing the time/cost of AM control and optimization processes. Our adoption of the operator learning (OL) approach enabled us to learn a family of differential equations produced by altering process variables in the laser’s Gaussian point heat source. More specifically, we used the Fourier neural operator (FNO) and deep operator network (DeepONet) to develop ROMs for time-dependent responses. Furthermore, we benchmarked the performance of these OL methods against a conventional deep neural network (DNN)-based ROM. Ultimately, we found that OL methods offer comparable performance and, in terms of accuracy and generalizability, even outperform DNN at predicting scalar model responses. The DNN-based ROM afforded the fastest training time. Furthermore, all the ROMs were faster than the original MOOSE model yet still provided accurate predictions. FNO had a smaller mean prediction error than DeepONet, with a larger variance for time-dependent responses. Unlike DNN, both FNO and DeepONet were able to simulate time series data without the need for dimensionality reduction techniques. Finally, the present work can help facilitate the AM optimization process by enabling faster execution of simulation tools while still preserving evaluation accuracy.

36 MATERIALS SCIENCE↗

The nano- and meso-scale structure of amorphous calcium carbonate

Abstract Understanding the underlying processes of biomineralization is crucial to a range of disciplines allowing us to quantify the effects of climate change on marine organisms, decipher the details of paleoclimate records and advance the development of biomimetic materials. Many biological minerals form via intermediate amorphous phases, which are hard to characterize due to their transient nature and a lack of long-range order. Here, using Monte Carlo simulations constrained by X-ray and neutron scattering data together with model building, we demonstrate a method for determining the structure of these intermediates with a study of amorphous calcium carbonate (ACC) which is a precursor in the bio-formation of crystalline calcium carbonates. We find that ACC consists of highly ordered anhydrous nano-domains of approx. 2 nm that can be described as nanocrystalline. These nano-domains are held together by an interstitial net-like matrix of water molecules which generate, on the mesoscale, a heterogeneous and gel-like structure of ACC. We probed the structural stability and dynamics of our model on the nanosecond timescale by molecular dynamics simulations. These simulations revealed a gel-like and glassy nature of ACC due to the water molecules and carbonate ions in the interstitial matrix featuring pronounced orientational and translational flexibility. This allows for viscous mobility with diffusion constants four to five orders of magnitude lower than those observed in solutions. Small and ultra-small angle neutron scattering indicates a hierarchically-ordered organization of ACC across length scales that allow us, based on our nano-domain model, to build a comprehensive picture of ACC formation by cluster assembly from solution. This contribution provides a new atomic-scale understanding of ACC and provides a framework for the general exploration of biomineralization and biomimetic processes.

59 BASIC BIOLOGICAL SCIENCES↗

Multimode two-dimensional vibronic spectroscopy. I. Orientational response and polarization-selectivity

Two-dimensional Electronic–Vibrational (2D EV) spectroscopy and two-dimensional Vibrational–Electronic (2D VE) spectroscopy are among the newest additions to the coherent multidimensional spectroscopy toolbox, and they are directly sensitive to vibronic couplings. In this first of two papers, the complete orientational response functions are developed for a model system consisting of two coupled anharmonic oscillators and two electronic states in order to simulate polarization-selective 2D EV and 2D VE spectra with arbitrary combinations of linearly polarized electric fields. Here, we propose analytical methods to isolate desired signals within complicated spectra and to extract the relative orientation between vibrational and vibronic dipole moments of the model system using combinations of polarization-selective 2D EV and 2D VE spectral features. Time-dependent peak amplitudes of coherence peaks are also discussed as means for isolating desired signals within the time-domain. This paper serves as a field guide for using polarization-selective 2D EV and 2D VE spectroscopies to map coupled vibronic coordinates on the molecular frame.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Learning interpretable surface elasticity properties from bulk properties via neural network equation learners

Surface elasticity is central to understanding the mechanics and stability of surfaces and interfaces. It is characterized by quantities such as surface tension, residual surface stress, and surface stiffness. However their analytical expressions are typically difficult to derive from atomistic data, and depend strongly on modeling choices. This work presents a neural network-based equation learner which combines customized activation functions and connection-based pruning to discover parsimonious, closed-form equations for surface elasticity from atomistic simulations. Applying the method to seven face-centered cubic (FCC) metals, our equation learner uncovers interpretable equations that describe both low-Miller index and high-Miller index surface properties, capturing long-tail property distributions accurately. The discovered expressions are decoupled into two components: a universal, geometry-driven orientation function, and material-specific baseline coefficients. We find that lower-order properties such as surface tension are fundamentally geometry dependent, while higher-order properties such as surface stress and elasticity show more complex geometry and material dependence. We also relate material dependent coefficients to bulk properties, forming a clear map from bulk material properties to surface elasticity. Overall, this approach demonstrates that interpretable neurosymbolic machine learning can bridge the gap between atomistic simulations and physical laws, enabling the discovery of generalizable structure–property relationships for materials science phenomena such as surface elasticity.

Equation learning↗

An implicit, sharp numerical treatment of viscous terms at arbitrarily shaped liquid-gas interfaces in evaporative flows

We report a novel methodology for accurate treatment of viscous terms in evaporation problems. The proposed scheme is an extension of the sharp viscous treatment of Kang et al. (2000) [7] to 3D phase change problems. To ensure accuracy and grid converging solutions, a new implicit approach to computing viscous fluxes across the phase interfaces is proposed, which was previously unavailable in fixed grid numerical schemes. Analytical relations were derived for the jump in velocity gradients across the 2D and 3D phase interfaces, an important constituent of the proposed scheme. The relations show a non-vanishing jump in the tangential gradients across the phase interface that are associated with evaporative flux and interfacial curvature. The proposed methodology demonstrated first order accuracy in canonical test cases. It is general and applicable to arbitrarily oriented interfaces, and can be readily implemented in existing evaporation flow solvers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Coupling of organic cation and inorganic lattice in methylammonium lead halide perovskites: Insights into a pressure-induced isostructural phase transition

Organic-inorganic hybrid metal-halide perovskite (OIHP) materials provide a tunable platform for engineering their optoelectronic properties. Although several high-pressure studies have been conducted from the OIHP family of single crystals and films, the exact nature of the dynamic coupling of the CH 3 NH 3 (MA) cation with the octahedral lattice framework and the mechanisms responsible for the structural phase transformation under pressure are not well captured. By combined photoluminescence (PL), synchrotron-based x-ray diffraction, and Raman-scattering studies as a function of pressure from methylammonium lead bromide (MAPbBr 3 ), we shed light on an isostructural phase transition due to the coupling of the MA cation and the PbBr 6 lattice through hydrogen bonding. The sharp discontinuities at ~ 1 GPa and ~ 3 GPa in the PL peak positions correlate with the structural changes observed in high-pressure XRD and Raman-scattering studies. The electronic band edge as a function of pressure is calculated within density-functional theory. The PL peak position, intensity and width of the excitonic peak show significant changes at 2 GPa, which corroborate the changes observed in high-pressure Raman-scattering studies. The frequencies of the lattice modes and the C-H/N-H bending and stretching modes of the MA cation show anomalous changes and other nuances at 2 GPa. The suppression of rotational and orientational disorder of the organic moiety is initiated at 2 GPa and the ordering is completed by 3.0 GPa, leading to an order-disorder type cubic II (Im$\overline{3}$) to orthorhombic ( Pnma ) phase transition. Along with the revelation of an isostructural transformation at 2 GPa, this paper highlights the impact of molecular vibrations on the electronic properties of MAPbBr 3 under pressure.

36 MATERIALS SCIENCE↗

Idiomatic Correctness-Checking via Julienne in Fortran 2023

This paper presents a unified approach to unit testing and runtime assertion checking using Fortran 2023. The paper describes the support for our approach in the Julienne framework. Julienne leverages recent Fortran standards to implement object-oriented design patterns, support testing parallel programs, and implement functional programming patterns in order to craft idioms inspired by natural-language expressions. The presented idioms employ novel operators to write expressions that evaluate to a test-diagnosis object encapsulating two components: (1) the test outcome or assertion outcome and (2) an automatically generated diagnostic string. Two other novel aspects of the approach include (1) the ability to enforce assertions inside pure procedures and (2) the ability to output rich diagnostic information inside pure procedures during error termination when assertions fail. The latter capability mitigates against a reason that Fortran programmers commonly cite for not writing pure procedures: difficulty obtaining useful program output inside pure procedures when debugging code. This paper demonstrates how the adoption of the proposed idioms leads naturally to a unifying theme across two otherwise disparate technologies: unit testing and runtime assertion checking. Finally, this paper describes the usage of the Julienne testing framework for writing unit tests and assertions in the Matcha high-performance computing application and the Fiats deep learning library.

Rouson, Damian↗

Imaging phonon dynamics with ultrafast electron microscopy: Kinematical and dynamical simulations

Ultrafast x-ray and electron scattering techniques have proven to be useful for probing the transient elastic lattice deformations associated with photoexcited coherent acoustic phonons. Indeed, femtosecond electron imaging using an ultrafast electron microscope (UEM) has been used to directly image the influence of nanoscale structural and morphological discontinuities on the emergence, propagation, dispersion, and decay behaviors in a variety of materials. Here, we describe our progress toward the development of methods ultimately aimed at quantifying acoustic-phonon properties from real-space UEM images via conventional image simulation methods extended to the associated strain-wave lattice deformation symmetries and extents. Using a model system consisting of pristine single-crystal Ge and a single, symmetric Lamb-type guided-wave mode, we calculate the transient strain profiles excited in a wedge specimen and then apply both kinematical- and dynamical-scattering methods to simulate the resulting UEM bright-field images. While measurable contrast strengths arising from the phonon wavetrains are found for optimally oriented specimens using both approaches, incorporation of dynamical scattering effects via a multi-slice method returns better qualitative agreement with experimental observations. Contrast strengths arising solely from phonon-induced local lattice deformations are increased by nearly an order of magnitude when incorporating multiple electron scattering effects. We also explicitly demonstrate the effects of changes in global specimen orientation on the observed contrast strength, and we discuss the implications for increasing the sophistication of the model with respect to quantification of phonon properties from UEM images.

36 MATERIALS SCIENCE↗

Control-Oriented Modeling of Lithium-Ion Batteries

Battery management systems (BMSs), which monitor and optimize performance while ensuring safety, require control-oriented models, i.e., models tailored to the design and implementation of estimation and control algorithms. Physics-based electrochemical models describe detailed battery phenomena, but are too computationally intensive for use in estimation and control applications. Single particle models (SPMs), which retain some of the physics of electrochemical models, are often used for control-oriented battery modeling since they are computationally efficient; however, they are only valid over very low frequency ranges and C-rates. Empirical equivalent circuit models (ECMs) are also used for control-oriented battery modeling since they are computationally efficient and can describe battery behavior over wide frequency ranges; however, they provide no physical understanding of the battery and, therefore, have limited applicability. Further, fractional order terms (e.g., Warburg impedances) are often employed, making the models unwieldy for use in the time domain. This work provides a control-oriented battery model that combines the benefits of SPM and ECM models, while overcoming their limitations. The proposed model incorporates some of the battery physics found in electrochemical models, can easily be used in both the time and frequency domains, and describes battery behavior over its entire frequency range. A linearized single particle model, which incorporates key electrochemical parameters, is used for modeling battery physics at very low frequencies. For low frequencies, integer-order linear systems are used to approximate diffusion physics described by Warburg impedances, and high frequency behavior is modeled by the double layer capacitance effect. The proposed battery model is more computationally efficient than full electrochemical models since it does not require the solution of PDEs, is accurate for a wider frequency range than the SPM considered in this paper, and does not suffer from the unwieldiness and limited applicability of empirical ECMs. Finally, the model is validated in the time and frequency domains via a comparison to pseudo-two-dimensional (P2D) model simulations and experimental data.

25 ENERGY STORAGE↗

Pore-Scale Modeling of Electrokinetics in Geomaterials

Pore-scale finite-volume continuum models of electrokinetic processes are used to predict the Debye lengths, velocity, and potential profiles for two-dimensional arrays of circles, ellipses and squares with different orientations. The pore-scale continuum model solves the coupled Navier–Stokes, Poisson, and Nernst–Planck equations to characterize the electro-osmotic pressure and streaming potentials developed on the application of an external voltage and pressure difference, respectively. Here, this model is used to predict the macroscale permeabilities of geomaterials via the widely used Carmen–Kozeny equation and through the electrokinetic coupling coefficients. The permeability results for a two-dimensional X-ray tomography-derived sand microstructure are within the same order of magnitude as the experimentally calculated values. The effect of the particle aspect ratio and orientation on the electrokinetic coupling coefficients and subsequently the electrical and hydraulic tortuosity of the porous media has been determined. These calculations suggest a highly tortuous geomaterial can be efficient for applications like decontamination and desalination.

36 MATERIALS SCIENCE↗

Hybridized Discontinuous Galerkin Methods for Computational Fluid Dynamics

Hybridizable Discontinuous Galerkin (HDG) methods hold promise for any applications with significant advection character, including thermal hydraulics in light water reactors and advanced reactor concepts and fluid models of plasmas in magnetic confinement fusion. Its features include natural upwinding, local element conservation, and extensibility to arbitrarily high order accuracy. In the last fiscal year we have implemented HDG in the Multiphysics Object-Oriented Simulation Environment (MOOSE). We developed a first-of-its-kind automatic static condensation system in MOOSE’s underlying finite element library libMesh which can condense out arbitrarily many internal variables. Finally, we developed the first preconditioner for HDG discretizations of the Navier-Stokes equations which shows robust performance across a wide range of problem sizes and Reynolds numbers. This preconditioner yields solution times that are equivalent to the fastest developed for industry standard finite volume methods. Moreover, the arbitrarily high-order nature of HDG makes it a prime candidate for acceleration via graphical processing units (GPUs). We believe these developments will hold significant importance for future DOE Nuclear Energy (NE) and Fusion Energy Science (FES) programs.

97 MATHEMATICS AND COMPUTING↗

Quantifying the resistance to dislocation glide in single phase FeCrAl alloy

The resistance to dislocation glide associated with slip systems {110}<111> and {112}<111> in single phase FeCrAl alloy is measured via micromechanical testing in a scanning electronic microscopy at room temperature. Two important factors, the shape and orientation of a pillar, are discussed with respect to the glide resistance and stress-strain response. Maximizing Schmid factor of one specific slip system while minimizing others is recommended in order to diminish obvious dislocations-induced hardening during in-situ testing. Apparent Schmid factor analysis is conducted to select grains with preferred orientations. Two types of pillars with conventional cylindrical shape or dog-bone shape are tested under compression to estimate the resistance to dislocation glide and evaluate the effect of pillar shape on the compression stress-strain response. One dog-bone pillar is tested under tension to check the tension-compression isotropy of dislocation slip. We discover that the shape of a pillar to a smaller extent affects the measured resistance but strongly influences the stress strain response. Cylindrical pillars exhibit apparent hardening associated with early yielding due to stress or strain concentration at contact region, while dog-bone pillars show an obvious yielding and continuous shearing without hardening. The resistance is 220 MPa for slip system {110}<111> and 230 MPa for slip system {112}<111>. Finite element analysis is performed to account for the influence of pillar shape and contact condition on mechanical response.

36 MATERIALS SCIENCE↗

Observation of Giant Optical Linear Dichroism in a Zigzag Antiferromagnet FePS 3

Direct optical probing of the antiferromagnetic order parameter in atomically thin samples is challenging, for example, via magneto-optical spectroscopy, due to the lack of net magnetization. Here, we report zigzag-antiferromagnetism (AFM) induced optical linear dichroism (LD) in layered transition-metal thiophosphate FePS3 down to the monolayer limit. The observed LD is giant despite having the optical wave vector parallel to the Neel vector. The LD is at least one order of magnitude larger than those reported in other antiferromagnetic systems, where the optical wave vector is orthogonal to the Neel vector. The large LD enables the probe of 60 degrees orientated zigzag-AFM domains. The optical anisotropy in FePS 3 originates from an electronic anisotropy associated with the zigzag direction of the AFM order and is independent of the spin-pointing direction. Furthermore, our findings point to a new optical approach for the investigation and control of zigzag or stripe magnetic order in strongly correlated systems.

2D magnet↗

Magnetic anisotropy and spin dynamics in the kagome magnet Fe 4 Si 2 Sn 7 O 16 : NMR and magnetic susceptibility study on oriented powder

Fe 4 Si 2 Sn 7 O 16 hosts an undistorted kagome lattice of Fe 2+ (3d 6 , S = 2) ions. We present results of bulk magnetization and Sn nuclear magnetic resonance (NMR) measurements on an oriented Fe 4 Si 2 Sn 7 O 16 powder sample oriented in geometries parallel (∥) and perpendicular (⊥) to the external applied magnetic field used for orienting the powder (Bori). The bulk susceptibility χ shows a broad peak at T N ~ 3 K associated with antiferromagnetic ordering. NMR spectra indicate the presence of planar anisotropy in the kagome planes. From an analysis of the static NMR shift (K) and dynamic spin-lattice relaxation rate (1/T 1 ) we conclude the presence of dominant magnetic fluctuations in the kagome planes. For the ∥ orientation, K scales linearly with the bulk susceptibility for temperatures down to ~ 4 K, while in the ⊥ orientation K starts to deviate strongly below T ~ 30 K. We associate this deviation with the onset of spin-tilting towards the kagome planes. These correlations are also reflected in the 1/T 1 data for the ∥ orientation, which starts to decrease below T ~ 30 K. Here, in this correlated regime, T N < T < ~ 30 K, we discuss the formation of positive chiral spin correlations in the kagome planes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Dynamic Structural Change of Plant Epidermal Cell Walls under Strain

Abstract The molecular foundations of epidermal cell wall mechanics are critical for understanding structure–function relationships of primary cell walls in plants and facilitating the design of bioinspired materials. To uncover the molecular mechanisms regulating the high extensibility and strength of the cell wall, the onion epidermal wall is stretched uniaxially to various strains and cell wall structures from mesoscale to atomic scale are characterized. Upon longitudinal stretching to high strain, epidermal walls contract in the transverse direction, resulting in a reduced area. Atomic force microscopy shows that cellulose microfibrils exhibit orientation‐dependent rearrangements at high strains: longitudinal microfibrils are straightened out and become highly ordered, while transverse microfibrils curve and kink. Small‐angle X‐ray scattering detects a 7.4 nm spacing aligned along the stretch direction at high strain, which is attributed to distances between individual cellulose microfibrils. Furthermore, wide‐angle X‐ray scattering reveals a widening of (004) lattice spacing and contraction of (200) lattice spacing in longitudinally aligned cellulose microfibrils at high strain, which implies longitudinal stretching of the cellulose crystal. These findings provide molecular insights into the ability of the wall to bear additional load after yielding: the aggregation of longitudinal microfibrils impedes sliding and enables further stretching of the cellulose to bear increased loads.

59 BASIC BIOLOGICAL SCIENCES↗

Interplay of pressure and temperature in modulating the magnetic properties of Cr 2 O 3

The spin-lattice interaction is pivotal in tailoring materials properties and is crucial for developing novel spintronic devices. Here, in the current work, we aim to explore the effects of pressure and temperature on the magnetic properties of Cr 2 O 3 to gain insights into its spin-lattice interaction. Through high-pressure neutron diffraction experiments, we observed an enhancement of the magnetic Bragg intensity with increasing pressure and it becomes more pronounced as the temperature increases. Results from first-principles calculations reveal a strengthening of the easy-axis magnetic anisotropy, doubling from 0 to 20 GPa. The exchange parameters, calculated based on this spin orientation, show an enhancement of the dominant magnetic interactions and an increase of magnetic ordering temperature when pressure increases. These findings suggest that the contributions from these two mechanisms are responsible for the observed increase in magnetic Bragg intensity.

Cr2O3↗

Large-scale 2D heterostructures from hydrogen-bonded organic frameworks and graphene with distinct Dirac and flat bands

The current strategies for building 2D organic-inorganic heterojunctions involve mostly wet-chemistry processes or exfoliation and transfer, leading to interface contaminations, poor crystallizing, or limited size. Here we show a bottom-up procedure to fabricate 2D large-scale heterostructure with clean interface and highly-crystalline sheets. As a prototypical example, a well-ordered hydrogen-bonded organic framework is self-assembled on the highly-oriented-pyrolytic-graphite substrate. The organic framework adopts a honeycomb lattice with faulted/unfaulted halves in a unit cell, resemble to molecular “graphene”. Interestingly, the topmost layer of substrate is self-lifted by organic framework via strong interlayer coupling, to form effectively a floating organic framework/graphene heterostructure. The individual layer of heterostructure inherits its intrinsic property, exhibiting distinct Dirac bands of graphene and narrow bands of organic framework. Our results demonstrate a promising approach to fabricate 2D organic-inorganic heterostructure with large-scale uniformity and highly-crystalline via the self-lifting effect, which is generally applicable to most of van der Waals materials.

36 MATERIALS SCIENCE↗