Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “relaxation”

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,549 records · Page 86
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Atomistic modeling of energy band alignment in CdTe(1 0 0) and CdTe(1 1 1) surfaces

An atomic-scale perspective of energy band alignment in CdTe surfaces has not been spatially studied despite the major role surfaces play in forming interfaces within CdTe-based thin film photovoltaic devices. Atomistic modeling based on density functional theory coupled with surface Green’s function is used for calculating energy band alignment of CdTe surfaces. The CdTe(1 0 0) ((1×1) and c(2×2) reconstruction) and CdTe(1 1 1) ((1×1) and (2×2) reconstruction) facets without and with surface relaxation provide insightful band bending characteristics that influence charge carrier transport. Results show that unrelaxed (1×1) CdTe(1 0 0) and CdTe(1 1 1) surfaces bend the valence band downward with surface polarity dictating the surface potential magnitude. The reconstructed CdTe(1 0 0) c(2×2) and CdTe(1 1 1) (2×2) surfaces result in favorable surface electronic features in relation to their unreconstructed variants. In addition, the structurally relaxed CdTe(1 1 1) surfaces develop an internal energy cusp potential that may enhance hole charge transport toward the back of CdTe solar cell devices. In conclusion, energy band alignments calculated within the study lead to a detailed understanding of how CdTe surfaces may affect CdTe-based thin film photovoltaic applications.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Probing the meta-stability of oxide core/shell nanoparticle systems at atomic resolution

Hybrid nanoparticles allow exploiting the interplay of confinement, proximity between different materials and interfacial effects. However, to harness their properties an in-depth understanding of their (meta)stability and interfacial characteristics is crucial. This is especially the case of nanosystems based on functional oxides working under reducing conditions, which may severely impact their properties. Here, the in-situ electron-induced selective reduction of Mn 3 O 4 to MnO is studied in magnetic Fe 3 O 4 /Mn 3 O 4 and Mn 3 O 4 /Fe 3 O 4 core/shell nanoparticles by means of high-resolution scanning transmission electron microscopy combined with electron energy-loss spectroscopy. Such in-situ transformation allows mimicking the actual processes in operando environments. A multi-stage image analysis using geometric phase analysis combined with particle image velocity enables direct monitoring of the relationship between structure, chemical composition and strain relaxation during the Mn 3 O 4 reduction. In the case of Fe 3 O 4 /Mn 3 O 4 core/shell the transformation occurs smoothly without the formation of defects. However, for the inverse Mn 3 O 4 /Fe 3 O 4 core/shell configuration the electron beam-induced transformation occurs in different stages that include redox reactions and void formation followed by strain field relaxation via formation of defects. This study highlights the relevance of understanding the local dynamics responsible for changes in the particle composition in order to control stability and, ultimately, macroscopic functionality.

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

The synergistic effects of central core size and end group engineering on performance of narrow bandgap nonfullerene acceptors

Understanding the relationship between the molecular structure and photoelectric properties of fused-ring non-fullerene acceptors (NFAs) is of far-reaching significance to the development of organic solar cells (OSCs). Herein, six NFAs based on multiple thiophenes (4 T, 6 T and 8 T) are employed to systematically probe the synergistic effects of extending central core size and terminal fluorination. The absorption results manifest that the molecular absorption is comprehensively affected by molecular crystallinity, planarity, conjugation length and intramolecular electron push–pull effect, simultaneously. The intensity of electron push–pull effect is not only related to the electron-donating ability of central core and the electron-withdrawing ability of end group, but also may be related to the distance between the positive and negative centers. The extension of central core leads to the more planar backbone and stronger crystallinity of NFAs, and less energy loss (E loss ) in its based OSC. Compared with the extension of central core, terminal fluorination has a greater impact on molecular photoelectric properties. The terminal fluorination significantly enhances the push–pull effect, lowers the energy levels, and slightly increases the vibrational relaxation. As a result, the strongest crystallinity and coplanarity of 8TIC-4F lead to a low vibrational relaxation of 0.18 eV, which makes PTB7-Th:8TIC-4F device exhibit a small E loss of 0.51 eV and a high efficiency of 10.4%. In addition, the fluorinated 6TIC-4F with suitable core size exhibits suitable energy level, absorption, crystallization, and phase separation morphology, making its as-cast device up to 11.61% efficiency. Finally, to the best of the authors’ knowledge, the PCE of 11.61% for PTB7-Th:6TIC-4F based device without any treatment is one of the highest values reported for the NAFs with over 1000 nm absorption.

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

Curing process and pore structure of metakaolin-based geopolymers: Liquid-state 1H NMR investigation

Geopolymers are emerging construction materials with lower carbon dioxide emissions compared to the conventional cementitious materials. The knowledge of the curing process and the related pore structures are important for optimizing the properties of these materials for different applications. The curing process and final pore structure are sensitive to the amount of used water, however the specifics are unclear. The curing process and pore structures of metakaolin-based geopolymers with a narrow water-to-solid (w/s) ratio (0.59–0.66) were monitored by nuclear magnetic resonance (NMR) relaxometry and cryoporometry. The 14-day curing process was investigated by monitoring the change of T{sub 2} and T{sub 1} relaxation times and water signal intensity. After the curing, the pore structures were characterized by 2D T{sub 1}-T{sub 2} correlation and T{sub 2}-T{sub 2} exchange measurements of absorbed water. The pore size distributions (PSDs) were measured with NMR cryoporometry and compared to nitrogen physisorption and mercury intrusion porosimetry (MIP) results. We found that the relaxation times decreased as the pore structure of the geopolymers matured during the curing while the dissolution and the condensation periods of the curing were distinguished by the changes in signal amplitude reflecting the proton density. After the curing, three distinct pore sizes and connectivity between pores were identified from T{sub 1}-T{sub 2} and T{sub 2}-T{sub 2} spectra. Their PSDs were measured, and they were found to correspond to two different pore sizes originating from the arrangement of clusters and defective pores. In the narrow w/s ratio (0.59–0.66), the curing times were the same for all samples when cured at 24 °C while the pore sizes were observed to increase as a function of the w/s ratio.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Plastic deformation mechanism of calcium-silicate hydrates determined by deviatoric-stress Raman spectroscopy

Highlights: • C-(A-)S-H samples under stress were characterized at the bond and grain level. • Small angle scattering showed that 720 MPa reduced grain thickness by 30%. • In-situ Raman showed that cross-linked C-A-S-H slips along its intralayer. • However, in-situ Raman showed that C-S-H slips along the interlayer. • In-situ Raman also showed a growth in CaCO{sub 3}, indicating fracture in both samples. Creep of the cement matrix affects the structural stability of concrete. In Portland cements, the creep is largely controlled by the binding phase calcium-(aluminum-)silicate-hydrate, or C-(A-)S-H. This phase has a lamellar structure and under deviatoric stress aligns its c-axis with the principal stress. However, the limiting resistance to this reorientation is unknown at the nanocrystalline level. Small-angle X-ray scattering shows that the lamellae thickness decreases under 100's MPa deviatoric stress. Deviatoric stress Raman spectroscopy shows that there are two ways that this break-up can occur. If the material's silicate chains are cross-linked, then strain in SiO bonds does not increase above certain stresses, indicating a relaxation adjacent to the SiO bond. If the chains are not cross-linked, then the silicate chains are broken up by rastering against each other, introducing defects. These results show that the plastic deformation of C-(A-)S-H is relevant for Portland cement creep.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Reducing the computational burden of a microgrid energy management system

As renewable technology advances and decreases in cost, microgrids are becoming an appealing means of distributed generation both for isolated communities and integrated with existing electrical grid systems. Due to their small size, however, microgrids may have financial limitations which preclude them from using commercial software to optimize control of their assets. Open-source optimization solvers are a viable alternative, but increase computation time. This work expands on a rolling horizon optimization framework for economic dispatch within an existing residential microgrid located in Hoover, Alabama. The microgrid has an open-source solver requirement and a need for quick solution time on a rolling horizon as opposed to a day-ahead commitment. We present a method of reducing integer variables by relaxation which completes two goals: reduction in computation time for real-time operations, and reduction in daily operational cost for the microgrid. Seasonal data for load and photovoltaic (PV) power was also collected from the microgrid to facilitate simulation testing. Computation time was successfully reduced using multiple variations of the relaxation method, while obtaining solution quality with operational cost similar to or better than the original model.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

An optimization-based approach to tailor the mechanical response of soft metamaterials undergoing rate-dependent instabilities

An optimization-based design framework is proposed to tune the response of soft metamaterials involving both geometric instabilities and nonlinear viscoelastic material behavior. Designing the response of soft metamaterials to harness instabilities and undergo large, tailored configuration changes will enable advancements in soft robotics, shock and vibration mitigation, and flexible electronics. In line with the metamaterial concept, the response of these materials is governed to a large extent by the geometric and topological makeup of their small-scale features. However, the link between structure and response is less intuitive for soft metamaterials due to their reliance upon highly nonlinear responses triggered by geometric instabilities. This is further complicated by the effects of viscoelastic relaxation, which recent studies have shown to alter the emergence of instabilities in non-intuitive ways. Here, these effects are accounted for in our framework to achieve various design objectives, including tailored force–displacement response and maximized energy absorption from both geometric and material effects. To fully automate this process, it is essential to have a completely robust equation solver for forward problems involving instabilities and viscoelastic relaxation. We achieve this by casting the search for stable mechanical equilibrium — i.e. the forward problem — as a minimization problem and utilize a trust region algorithm to robustly handle instabilities and follow energetically-favorable equilibrium paths through critical points.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Fast meta-solvers for 3D complex-shape scatterers using neural operators trained on a non-scattering problem

Three-dimensional target identification using scattering techniques requires high accuracy solutions and very fast computations for real-time predictions in some critical applications. We first train a deep neural operator (DeepONet) to solve wave propagation problems described by the Helmholtz equation in a domain without scatterers but at different wavenumbers and with a complex absorbing boundary condition. We then design two classes of fast meta-solvers by combining DeepONet with either relaxation methods, such as Jacobi and Gauss-Seidel, or with Krylov methods, such as GMRES and BiCGStab, using the trunk basis of DeepONet as a coarse-scale preconditioner. We leverage the spectral bias of neural networks to account for the lower part of the spectrum in the error distribution while the upper part is handled inexpensively using relaxation methods or fine-scale preconditioners. The meta-solvers are then applied to solve scattering problems with different shape of scatterers, at no extra training cost. We first demonstrate that the resulting meta-solvers are shape-agnostic, fast, and robust, whereas the standard standalone solvers may even fail to converge without the DeepONet. We then apply both classes of meta-solvers to scattering from a submarine, a complex three-dimensional problem. We achieve very fast solutions, especially with the DeepONet-Krylov methods, which require orders of magnitude fewer iterations than any of the standalone solvers.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Monte-Carlo modeling of phonon thermal transport using DFT-based anisotropic dispersion relations over the full Brillouin zone

In this work, we present a Monte Carlo (MC) approach to solve the phonon Boltzmann transport equation (BTE) in which the anisotropic phonon dispersion relations over the full Brillouin zone (BZ) are used. In this approach, the discretization of the BZ used to compute the phonon relaxation time places constraints on the direction of scattered phonons in the real-space simulation domain. The phonon dispersion and phonon relaxation times are calculated using the density functional theory (DFT) approach. The modified MC approach is validated by a close examination of its ability to simulate phonon transport in both the ballistic and diffusive regimes for multiple materials including GaAs, InAs, ThO 2 , and α-U. In doing so, the phonon thermal conductivities from 100 K to 1000 K are calculated and compared with traditional non-transport solution of the phonon BTE. It is found that the phonon thermal conductivities of α-U and ThO 2 obtained from MC simulations using isotropic dispersion are larger than the values obtained using anisotropic phonon dispersion relations over the full BZ. The effect of phonon-defect scattering on the thermal conductivity of ThO 2 is also studied as an application of the current MC approach and found to agree with previously computed values in the literature. The MC solver developed here has been parallelized as a step to demonstrate its potential to solving computationally intensive phonon thermal transport problems at the mesoscale.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Graph neural networks predict energetic and mechanical properties for models of solid solution metal alloy phases

Here, we developed a PyTorch-based architecture called HydraGNN that implements graph convolutional neural networks (GCNNs) to predict the formation energy and the bulk modulus for models of solid solution alloys for various atomic crystal structures and relaxed volumes. We trained the GCNN surrogate model on a dataset for nickel–niobium (NiNb) generated by the embedded atom model (EAM) empirical interatomic potential for demonstration purposes. The dataset was generated by calculating the formation energy and the bulk modulus as a prototypical elastic property for optimized geometries starting from initial body-centered cubic (BCC), face-centered cubic (FCC), and hexagonal compact packed (HCP) crystal structures, with configurations spanning the possible compositional range for each of the three types of initial crystal structures. Numerical results show that the GCNN model effectively predicts both the formation energy and the bulk modulus as function of the optimized crystal structure, relaxed volume, and configurational entropy of the model structures for solid solution alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Multi-material hydrodynamics with algebraic sharp interface capturing

A finite volume method for Eulerian multi-material hydrodynamics with sharp interface capturing is presented here. The pressure-temperature non-equilibrium multi-material system with finite-rate pressure relaxation in mixed-cells is considered here. This pressure closure facilitates material-property-dependent pressure relaxation, rather than instantaneous pressure equilibration, which in turn allows the use of unsplit high-order time-integrators. A modified tangent of hyperbola for interface capturing (THINC) method is used to reconstruct multi-material (> 2) interfaces, on three-dimensional unstructured meshes. A simple modification which extends the THINC reconstruction to interfaces between more than two materials is proposed. It is demonstrated that the modified THINC can capture multi-material interfaces within 2–4 tetrahedral cells. Since no geometric reconstructions are required by the THINC method, the presented multi-material method is algorithmically simple, and computationally efficient. Consistent reconstructions of conserved quantities at material interfaces ensure that conservation and closure laws are satisfied at the discrete level. Through a suite of test problems solved on unstructured meshes, it is demonstrated that the presented method is a promising candidate for accurate and efficient multi-material hydrodynamics computations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Particle collisionality in scaled kinetic plasma simulations

Kinetic plasma processes, such as magnetic reconnection, collisionless shocks, and turbulence, are fundamental to the dynamics of astrophysical and laboratory plasmas. Simulating these processes often requires particle-in-cell (PIC) methods, but the computational cost of fully kinetic simulations can necessitate the use of artificial parameters, such as a reduced speed of light and ion-to-electron mass ratio, to decrease expense. While these approximations can preserve overall dynamics under specific conditions, they introduce nontrivial impacts on particle collisionality that are not yet well understood. In this work, we develop a method to scale particle collisionality in simulations employing an artificial speed of light and/or an artificial ion-to-electron mass ratio. By introducing species-dependent scaling factors, we independently adjust inter- and intra-species collision rates to better replicate the collisional properties of the physical system. Our approach maintains the fidelity of electron and ion transport properties while preserving critical relaxation rates, such as energy exchange timescales, within the limits of weakly collisional plasma theory. Furthermore, we demonstrate the accuracy of this scaling method through benchmarking tests against theoretical relaxation rates and connecting to fluid theory, highlighting its ability to retain key transport properties. Existing collisional PIC implementations can be easily modified to include this scaling, which will enable deeper insights into the behavior of marginally collisional plasmas across various contexts.

Totorica, S. R. [Princeton Univ., NJ (United State↗
DOE OSTI2020

Toward a scalable robust security-constrained optimal power flow using a proximal projection bundle method

Robust security-constrained optimal power flow (rSCOPF) aims to find the worst-case contingencies of alternating current optimal power flow (ACOPF) in power systems. With the rise of GPU architectures on the upcoming supercomputer architectures, optimization algorithms that rely on sparse linear algebra and indefinite linear systems are becoming increasingly hard to solve efficiently (e.g. interior-point method). To address this we revisit a maximin optimization formulation of the rSCOPF and the single-level mixed-integer semidefinite programming (MISDP) reformulation, which is obtained by taking the Lagrangian relaxation of the inner minimization ACOPF problem. In this paper, we focus on the development of a proximal projection bundle method (PPBM) for solving continuous relaxation node subproblems of the MISDP problem, based primarily on the well-known alternating direction method of multipliers. Cutting planes reminiscent of bundle method ideas are also applied in coordination with updates of the proximal parameter. The cutting-plane method can generate a large number of linear inequalities, leading to a large scale but decomposable quadratic programming (QP) subproblem that is amenable to GPUs. We present the numerical results on the IEEE 30, 57, 118, and 300-bus systems by using our PBMM method. We discuss the main computational bottleneck of our method, which is the time taken to solve each iteration of a QP subproblem instance of the PPBM, and how GPU architectures can accelerate this solution process.

bundle method↗
DOE OSTI2022

Call attention to using DRT and EIS to quantify the contributions of solid oxide cell components to the total impedance

The contributions of anode and cathode processes to solid oxide fuel cell (SOFC) impedance were studied using electrochemical impedance spectroscopy (EIS) and distribution of relaxation time (DRT). Specifically, the role of gas composition at both anode and cathode was explored in a systematic study aimed at deconvoluting, identifying and quantifying the contributions of different anode and cathode processes. Further, five distinct relaxation time distributions were observed in the DRT spectra; although the peaks were ascribed to diffusion and charge transfer processes at the electrodes, it was found that both electrodes contributed to various peaks at the same time. Moreover, it was found that two distinctive operating conditions could return equivalent DRT spectra. While DRT analysis allows to obtain useful information regarding cell performance, extra consideration is needed when assessing and quantifying anodic and cathodic resistances within the SOFC.

08 HYDROGEN↗
DOE OSTI2020

Towards understanding the development and characteristics of under-expanded flash boiling jets

Here in this research, the n-hexane jets discharged from a straight channel were studied numerically to reveal the development and characteristics of under-expanded flashing jets. Firstly, the development of flash boiling jets was investigated and compared to that of two-phase jets without phase transition. It was shown that the two-phase jets without phase transition can be under-expanded only if the flow is choked. In contrast, flash boiling jets can be under-expanded whether or not the choked conditions are verified. Secondly, the factors affecting the flashing jet behaviors were discussed. The relaxation time in phase transition, the expansion ratio, and the injection pressure have been identified as three key factors determining the jet characteristics. With high enough expansion ratio and low enough relaxation time, the expansion was so strong that the flow became transonic and shock structures similar to those of under-expanded gaseous jets appeared near the orifice exit. Finally, two types of choking have been identified for flash boiling jets, i.e., vaporization-induced choking and expansion-induced choking. The former one is mainly ascribed to the vapor formation in the orifice which reduces the local speed of sound, while the latter one is due to the expansion of under-expanded two-phase flow and the resultant flow acceleration.

42 ENGINEERING↗
DOE OSTI2021

A statistical mechanics approach to macroscopic limits of car-following traffic dynamics

Herein we study the derivation of macroscopic traffic models from car-following vehicle dynamics by means of hydrodynamic limits of an Enskog-type kinetic description. We consider the superposition of Follow-the-Leader (FTL) interactions and relaxation towards a traffic-dependent Optimal Velocity (OV) and we show that the resulting macroscopic models depend on the relative frequency between these two microscopic processes. If FTL interactions dominate then one gets an inhomogeneous Aw–Rascle–Zhang model, whose (pseudo) pressure and stability of the uniform flow are precisely defined by some features of the microscopic FTL and OV dynamics. Conversely, if the rate of OV relaxation is comparable to that of FTL interactions then one gets a Lighthill–Whitham–Richards model ruled only by the OV function. We further confirm these findings by means of numerical simulations of the particle system and the macroscopic models. Unlike other formally analogous results, our approach builds the macroscopic models as physical limits of particle dynamics rather than assessing the convergence of microscopic to macroscopic solutions under suitable numerical discretisations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Viscoelastic necking dynamics between attractive microgels

Hypothesis: Microgels can deform and interpenetrate and display colloid/polymer duality. The effective interaction of microgels in the collapsed state is governed by the interplay of polymer–solvent interfacial tension and bulk elasticity. A connecting neck is shown to mediate microgel interaction, but its temporal evolution has not been addressed. Here we hypothesize that the necking dynamics of attractive microgels exhibits liquid-like or solid-like behavior over different time and length scales. Experiments: We simulate the merging and pinching of attractive microgels with different crosslinking densities in explicit solvent using dissipative particle dynamics. The temporal coalescence dynamics of microgels is investigated and compared with simple liquid and polymeric droplets. We model the neck growth on long time scales using Maxwell model of polymer relaxation and compare the theoretical prediction with simulation data. The mechanical strength of the neck is characterized systematically via simulated pinch-off of microgels by steered molecular dynamics. Findings: We evidence a crossover in the coalescence dynamics reflecting the viscoelastic signature of microgels. In contrast to the common knowledge that viscoelastic materials respond elastically on short time scales, the early expansion of the microgel neck exhibits a linear behavior, similar to the viscous coalescence of liquid droplets. However, the late regime with arrested dynamics resembles sintering of solid particles. Through an analytical model relating microgel dynamics to neck growth, we show that the long-term behavior is governed by stress relaxation of the polymers in the neck region and predict an exponential decay in the rate of growth, which agrees favorably with the simulation. Different from coalescence, the thread thinning in microgel breakup primarily highlights its polymeric characteristics.

42 ENGINEERING↗
DOE OSTI2022

Magnetoelastic properties of multiferroic hexagonal ErMnO 3

Here, the strength and dynamics of magnetoelastic coupling through the paramagnetic (PM) – antiferromagnetic (AFM) – ferrimagnetic (FIM) transitions in multiferroic hexagonal ErMnO 3 have been investigated by Resonant Ultrasound Spectroscopy. Elastic stiffening by up to 2% below the PM – AFM transition at 80 K arises from biquadratic coupling between strain and the magnetic order parameter with relaxation times longer than ~10 -6 s for the response of spins to changes in strain. In contrast with YMnO 3 , the PM – AFM transition in ErMnO 3 is accompanied by a peak in acoustic loss immediately below the Néel point which is interpreted in terms of strain relaxation accompanying ordering of spins of Er 3+ at 4b sites. Changes in the magnetic ordering scheme at the AFM – FIM transition near 3 K are accompanied by elastic softening of ~0.03 %. During poling of the low temperature ferrimagnetic structure round magnetic hysteresis loops, small changes in elastic stiffness which arise due to the contribution of piezomagnetic and/or piezoelectric moduli are detected. Contributions of piezoelectric moduli to acoustic resonance frequencies also permit changes in the configuration of ferroelectric domains to be detected in response both to cycling through this transition and to application of a magnetic field. A peak in acoustic loss in the vicinity of 250 K is attributed to strain-mediated pinning/freezing of some aspect of the domain microstructure with an activation energy of ~0.25–0.3 eV. A return to the original elastic properties on heating to temperatures above ~250 K is interpreted in terms of backswitching of domains to the configuration they had at the start. These observations confirm the existence of subtle variations in magnetoelastic coupling behaviour relating to both the magnetic order parameters and magnetic domain structures.

36 MATERIALS SCIENCE↗