Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Order parameters”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Recent progress in the phase-field dislocation dynamics method

The phase-field dislocation dynamics (PFDD) method, originated in 2002, is a continuum dislocation model that uses order parameters to describe dislocation slips in crystalline materials. In the past two decades, and especially since it was last reviewed in 2016, PFDD was advanced significantly in terms of the mathematical formulation, numerical implementation, and applicability. The main purpose of this short review is to summarize recent progress made to improve the energy functional formulation and numerical techniques of PFDD as well as its recent applications. Additionally, some recommendations for future work to further extend the PFDD method are presented.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Quantifying the atomistic free-volume morphology of materials with graph theory

Here, we introduce a new computational methodology for the identification and characterization of free volume within/around atomistic configurations. This scheme employs a three-stage workflow, by which spheres are iteratively grown inside of voxels, and ultimately converted to planar graphs, which are then characterized via a graph-based order parameter. Our approach is computationally efficient, physically intuitive, and universally transferable to any material system. Validation of our methodology is performed on several sets of materials problems: (1) classification of unique free volumes in various crystal phases, (2) autonomous detection and classification of complex surface defects during epitaxial growth simulations, (3) characterization of free volume defects in metals/alloys, and (4) quantification of the spatio-temporal behavior of nano-scale free volume morphologies as a function of both temperature and free-volume size. Our method accurately identifies and characterizes unique free volumes over a multitude of systems and length scales, indicating its potential for future use in understanding the relationship between free volume morphology and material properties under both static and dynamic conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A generalized and adaptable tensor-contraction-based cluster expansion formalism for multicomponent solids

Density functional theory (DFT)-based simulations of materials have first-principles accuracy, but are very computationally expensive. For simulating various properties of multi-component alloys, the cluster expansion (CE) technique has served as the standard workaround to improve computational efficiency. However, the standard CE technique is difficult to extend to exotic and/or low-symmetry lattices, often implemented via iteration over particular cluster types, which must be enumerated per lattice structure. In this work, we introduce the tensor cluster expansion (TCE), implemented in the open-source code tce-lib, which maps correlation functions to mixed tensor contractions, eliminating the need to iterate over cluster types and additionally making the calculation of correlation functions well-suited for massively parallel architectures like GPUs. We show that local interaction energies are an immediate consequence of the TCE formalism, yielding nearly $\mathcal{O}$(1) energy difference calculations. We then use this formalism to fit CE models for the TaW and CoNiCrFeMn systems, and use these models to respectively compute the enthalpy of mixing curve and Cowley short-range order parameters, showing excellent agreement with ground truth data.

Cluster expansion↗
DOE OSTI2023

Thermo4PFM: Facilitating Phase-field simulations of alloys with thermodynamic driving forces

Phase-field modeling is a popular front-tracking approach used to model solidification. Its time-evolution equations are often coupled to alloy composition and/or thermal diffusion in high-resolution multiphysics approaches. Materials thermodynamic properties tabulated in CALPHAD databases can be used for phase-field modeling to parameterize bulk energies of alloys. In addition, they can be naturally integrated into models such as the Kim-Kim-Suzuki (KKS) model where driving forces depend on the differences between chemical potentials of co-existing phases. In that case, a small system of coupled nonlinear equations needs to be solved at every point in space where the phase-field order parameter is to be updated and evolved in time. Here we present Thermo4PFM, a solver for the KKS equations for binary and ternary alloys, with two or three phases, and parameterized with CALPHAD models. Thermo4PFM is open source, written in C++, and can take advantage of Graphics Processing Units (GPU) accelerators. Using OpenMP offload capabilities for C++ classes, an excellent performance is demonstrated on GPU using the LLVM compiler. CALPHAD data is read from simple JSON files using an open source parser from the boost library.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Fouling behavior of zwitterionic membranes compared to polyamide membranes

Membrane fouling remains a critical bottleneck for reverse osmosis (RO) desalination, driving energy consumption and reducing membrane lifetime. Here, we employ all-atom molecular dynamics simulations to investigate the antifouling behavior of random zwitterionic amphiphilic copolymer (r-ZAC) membranes composed of sulfobetaine methacrylate (SBMA) and allyl methacrylate (AMA), benchmarked against conventional polyamide (PA) RO membranes. Structural and dynamical analyses—including radial distribution functions, coordination numbers, tetrahedral order parameters, vector orientation, and residence-time correlation functions—reveal that r-ZAC surfaces sustain tightly bound, long-lived hydration layers with preserved tetrahedrality and anisotropic water orientation, in sharp contrast to the weak and disordered hydration of PA. Steered molecular dynamics simulations demonstrate that r-ZAC membranes impose substantial free-energy barriers to foulant approach (alginate ≈ 90 kcal/mol, sucrose ≈ 35 kcal/mol, humic acid ≈ 15 kcal/mol), whereas PA membranes exhibit negligible barriers (< 1 kcal/mol) and thermodynamically favorable adsorption. Detailed foulant–surface interaction analyses show that zwitterionic hydration and electrostatic heterogeneity in r-ZAC suppress adhesion, except in the case of amphiphilic humic acid, which exploits multiple binding modes. Together, these results establish molecular-level design principles for antifouling membranes: the combination of zwitterionic hydration, structured interfacial water, and controlled amphiphilic balance in r-ZAC membranes provides superior resistance to organic fouling relative to PA.

Cross-linked polyamide↗
DOE OSTI2022

Implementing contact angle boundary conditions for second-order Phase-Field models of wall-bounded multiphase flows

In the present work, a general formulation is proposed to implement the contact angle boundary conditions for the second-order Phase-Field models, which is applicable to N-phase (N ≥ 2) moving contact line problems. To remedy the issue of mass change due to the contact angle boundary condition, a source term or Lagrange multiplier is added to the original second-order Phase-Field models, which is determined by the consistent and conservative volume distribution algorithm so that the summation of the order parameters and the consistency of reduction are not influenced. To physically couple the proposed formulation to the hydrodynamics, especially for large-density-ratio problems, the consistent formulation is employed. The reduction-consistent conservative Allen-Cahn models are chosen as examples to illustrate the application of the proposed formulation. The numerical scheme that preserves the consistency and conservation of the proposed formulation is employed to demonstrate its effectiveness. Results produced by the proposed formulation are in good agreement with the exact and/or asymptotic solutions. The proposed method captures complex dynamics of moving contact line problems having large density ratios.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Self-induced spin-orbit torques in metallic ferromagnets

Herein we present a phenomenological theory of spin-orbit torques in a metallic ferromagnet with spin-relaxing boundaries. The model is rooted in the coupled diffusion of charge and spin in the bulk of the ferromagnet, where we account for the anomalous Hall effects as well as the anisotropic magnetoresistance in the corresponding constitutive relations for both charge and spin sectors. The diffusion equations are supplemented with suitable boundary conditions reflecting the spin-sink capacity of the environment. In inversion-asymmetric heterostructures, the uncompensated spin accumulation exerts a dissipative torque on the order parameter, giving rise to a current-dependent linewidth in the ferromagnetic resonance with a characteristic angular dependence. We compare our model to recent spin-torque ferromagnetic resonance measurements, illustrating how rich self-induced spin-torque phenomenology can arise even in simple magnetic structures.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A phase-field fracture formulation for generalized standard materials: The interplay between thermomechanics and damage

Accurately modeling fracture of ductile materials poses open challenges in the field of computational mechanics due to the multiphysics nature of their failure processes. Integrating the interplay between thermodynamics and damage into ductile fracture models is vital for predicting critical failure modes. Here, in this paper, we develop a versatile phase-field (PF) framework for modeling ductile fracture, taking into account finite-strain elasto-plasticity. The framework stems from a variational formulation of constitutive relations for generalized standard materials (GSMs), whose response is described by a Helmholtz free energy and a dissipation pseudo-potential. Its variational structure is based on a minimum principle for a functional that expresses the sum of power densities for reversible and irreversible processes. By minimizing this functional with a constraint on a von Mises yield function, we derive the evolution equation for the equivalent plastic strain and an associative flow rule. This constrained optimization problem is analytically solved for a wide class of thermo-viscoplasticity models. The key innovations of the current work include (i) a cubic plastic degradation function that accounts for a non-vanishing damage-dependent yield stress, (ii) closed-form expressions of the Helmholtz free energy and dissipation pseudo-potential for three thermo-viscoplasticity models, (iii) an extended Johnson–Cook plasticity model with a nonlinear hardening law, and (iv) a plastic work heat source that depends on the plastic degradation function and a variable Taylor–Quinney (TQ) coefficient. The capabilities of the proposed framework are tested with the aid of four ductile fracture problems, including the Sandia Fracture Challenge. In each of these problems, we examine the evolution of relevant field variables such as the PF order parameter, the equivalent plastic strain, the temperature, and the internal power dissipation density, in addition to the overall structural response quantified by the force–displacement curve. These numerical studies demonstrate that the proposed framework effectively represents ductile fracture, yielding computational results that exhibit good agreement with experimental data.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Impact of grain boundary and surface diffusion on predicted fission gas bubble behavior and release in UO 2 fuel

In this work, we quantify the impact of grain boundary (GB) and surface diffusion on fission gas bubble evolution and fission gas release in UO 2 nuclear fuel using simulations with a hybrid phase field/cluster dynamics model. Here, we begin with a comprehensive literature review of uranium vacancy and xenon atom diffusivity in UO 2 through the bulk, along GBs, and along surfaces. In our model we represent fast GB and surface diffusion using a heterogeneous diffusivity that is a function of the order parameters that represent bubbles and grains. We find that the GB diffusivity directly impacts the rate of gas release via GB transport, and that the GB diffusivity is likely below 104 times the lower value from Olander and van Uffelen. We also find that the surface diffusivity impacts bubble coalescence and mobility, and that the bubble surface diffusivity is likely below 10 -4 times the value from Zhou and Olander.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Towards the design of nature-inspired materials: Impact of complex pore morphologies via higher-order homogenization

Even though the development of novel materials that mimic nature is widely used in a variety of engineering and scientific fields, the relationship between effective material properties and underlying, often complex pore morphology is still not fully understood. To address this knowledge gap and accelerate the development of novel nature-inspired materials, this paper adopts a higher-order asymptotic homogenization method to numerically investigate the effect of complex micropore morphology on the effective mechanical properties of a porous system. Specifically, we create unique pore morphologies with varying levels of complexity that serve as a more realistic representation of natural materials. Here, we then use the second-order homogenization method to capture the role of pore size, shape, orientation, and distribution on effective properties. By creating different pore morphologies, we systematically studied the relationship between morphology and effective mechanical properties. The results highlight the necessity of higher-order parameters to fully capture the role of realistic pore morphologies on effective mechanical properties and provide a path forward in the design of nature-inspired materials.

36 MATERIALS SCIENCE↗
DOE OSTI2021

From classical thermodynamics to phase-field method

Phase-field method is a density-based computational method at the mesoscale for modeling and predicting the temporal microstructure and property evolution during materials processes. The focus of this article is on connecting the most common phase-field equations to the very basic first and second laws of classical thermodynamics through rudimentary irreversible thermodynamics. It briefly discusses the relations of the continuum phase-field equations to their counter parts at the microscopic and atomic levels. It attempts to clarify the contributions of long-range elastic, electrostatic, and magnetic interactions to domain structure evolution during structural, ferroelectric, and ferromagnetic phase transformations by separating order parameter changes due to the presence of quasi-static fields and those arising from phase transformations. A few examples are presented to demonstrate the possibility of employing the phase-field method to provide guidance to designing materials for optimum properties or discovering novel mesoscale phenomena or new materials functionalities. Here, the article ends with a brief perspective on a number of potential future directions on the development and applications of phase-field method beyond its traditional applications to structural alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Electrons matter in ferromagnetic criticality

When elemental iron is cooled from high temperatures, magnetic moments carried by its d-electrons begin a continuous transition to a state of parallel aligned moments below its Curie temperature TC of about 1040 K. As established in the mid-1900s, the statistical physics of this second order phase transition is universal, independent of the ordering temperature, material structure and chemical composition as well as the character of moment-bearing electrons. Finally, the thermodynamics of the phase transition is determined solely by symmetries of the transition’s order parameter – magnetization in the case of a ferromagnet.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Fragmentation of the edge of a terminated Cu nanolayer within a Nb matrix upon annealing

Here, we investigate the evolution of the tip of a terminated Cu nano-layer embedded within a Nb matrix upon annealing at 800°C. The tip of the terminated Cu layer fragments progressively into a series of single-crystal Cu particles upon continued annealing. This finding shows that layer fragmentation is the second step, after layer pinchoff, in the thermal coarsening of nanolayered metal composites. Comparison to prior phase-field modeling permits us to infer that the mobility at 800°C of Cu-Nb, as modeled by the Cahn-Hillard equation, is M(c)=(M 0 = 7.3 ± 3.6 nm 5 / s • eV) |1 – c 2 |, where the order parameter, c, is 1 for Cu and -1 for Nb.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Facile alignment estimation in carbon nanotube films using image processing

Whether a macroscopic assembly of carbon nanotubes can exhibit the one-dimensional properties expected from individual nanotubes critically depends on how well the nanotubes are aligned inside the assembly. Therefore, a simple and accurate method for assessing the degree of alignment is desired for the rapid characterization of carbon nanotube films and fibers. Here, we present an end-to-end solution for determining the global and local spatial orientation of carbon nanotubes in films within a short amount of time using a fast, precise, and economical approach based on an image processing method applied to scanning electron microscopy images. Further, we first use Laplacian edge enhancement filtering for improving the appearance of edge regions, which is followed by image partitioning into multiple blocks to capture the nanoscale orientation characteristics and total variation-based image decomposition of these image blocks. We then perform a 2D-fast Fourier transform on the image decomposed textural components of these edge-enhanced image blocks to determine the orientation distribution, which is utilized to estimate the 2D nematic order parameter. To show the effectiveness of our method, we corroborated our results against results obtained with other state-of-the-art image processing and experimental techniques.

2D-FFT↗
DOE OSTI2022

To grind or not to grind? The influence of mechanical and thermal treatments on the $\mathrm{Cu/Zn}$ disorder in $\mathrm{Cu}$ 2 $\mathrm{ZnSn}$($\mathrm{S}$ x $\mathrm{Se}$ 1-x ) 4 monograins

We report that Kesterite-type based thin film solar cell technologies are mainly based on polycrystalline absorber layers. A promising low cost alternative technology uses Cu 2 ZnSn(S x Se 1-x ) 4 (CZTSSe) monograins (single crystals of 20–100 μm size) which are fixed in a polymer matrix to form a flexible solar cell. The Cu/Zn disorder is discussed as a possible reason for band tailing causing voltage losses limiting the efficiency of CZTSSe-based devices. The experimental determination of the order parameter Q which is a quantitative measure of Cu/Zn disorder, requires a differentiation between the isoelectronic cations Cu + and Zn 2+ . An in-depth analysis of neutron diffraction data allows the determination of type and concentration of intrinsic point defects including a distinction between Cu and Zn. Neutron diffraction requires large sample volumes, thus monograins offer the unique possibility to correlate structural disorder in CZTSSe with device performance parameters. In this study we tackle the influence of grinding the monograins on stoichiometry deviations, the Cu/Zn disorder as well as intrinsic point defects and optoelectronic properties of CZTSSe monograins. Moreover, an easy methodology based on Raman scattering spectroscopy is proposed for the assessment of Cu/Zn disorder in the CZTSSe compounds.

14 SOLAR ENERGY↗
DOE OSTI2021

Establishing a unified framework for ion solvation and transport in liquid and solid electrolytes

Electrolytes used in rechargeable batteries must enable rapid translation of the working ion between macroscopically separated electrodes. These electrolytes are, however, usually designed and synthesized using atomic-level insights. Whether the ideal electrolyte for a particular battery is a solid or a liquid remains an important unresolved question, especially as solids with conductivities comparable with liquids are discovered. Here, to help resolve such questions, we present the first steps toward a unified framework for relating atomic and continuum scale phenomena. Solvation shells in liquids are entities that translate with the working ion for a short while before they break up due to Brownian motion. By contrast, solvation cages in classical solids and polymers cannot not translate with the working ion. Mobility of the entities that make up the cages and shells, which is quantified by an order parameter, is shown to influence translation of the working ion on continuum length scales.

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

Identification and correction of temporal and spatial distortions in scanning transmission electron microscopy

Scanning transmission electron microscopy (STEM) has become the technique of choice for quantitative characterization of atomic structure of materials, where the minute displacements of atomic columns from high-symmetry positions can be used to map strain, polarization, octahedra tilts, and other physical and chemical order parameter fields. The latter can be used as inputs into mesoscopic and atomistic models, providing insight into the correlative relationships and generative physics of materials on the atomic level. However, these quantitative applications of STEM necessitate understanding the microscope induced image distortions and developing the pathways to compensate them both as part of a rapid calibration procedure for in situ imaging, and the post-experimental data analysis stage. Here, we explore the spatiotemporal structure of the microscopic distortions in STEM using multivariate analysis of the atomic trajectories in the image stacks. Based on the behavior of principal component analysis (PCA), we develop the Gaussian process (GP)-based regression method for quantification of the distortion function. The limitations of such an approach and possible strategies for implementation as a part of in-line data acquisition in STEM are discussed. Here, the analysis workflow is summarized in a Jupyter notebook that can be used to retrace the analysis and analyze the reader's data.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Predictive Design of Hybrid Improper Ferroelectric Double Perovskite Oxides

The computational design of suitable multiferroic double perovskite oxides requires finding materials that exhibit sizable polarization, magnetization, and coupling between them. Oxides with the chemical formula of AA'BB'O 6 with building blocks of ABO 3 single perovskite oxides in centrosymmetric Pnma symmetry are strong candidates that have been reported to satisfy such criteria. The system lowers to noncentrosymmetric, polar P2 1 symmetry if A/A' layered and B/B' rocksalt cation orderings are imposed. A detailed compositional search over a variety of chemical spaces followed by evaluating their polarization may lead to the identification of more of these compounds with ferroelectric ordering. The standard density functional theory practices to estimate polarization within the Berry phase formalism require the systems to be perfectly insulating. The number of compounds that can be evaluated using this method is therefore limited. In this work, we introduce a predictive learning strategy based on importance sampling to build a series of machine learning models using results from first-principles simulations to predict polarization and the corresponding switching barrier. The geometry-driven features related to charge states and cationic radii play key roles in predicting the switching barrier with complementary contributions from the key structural mode-based order parameters. These modes become important to draw reasonable predictions of polarization components from machine learning models. In conclusion, our predictive models identify candidates with high polarizations and low switching barriers from a pool of double perovskite oxides, suitable for future investigation for their potential applications in spintronic devices.

36 MATERIALS SCIENCE↗