SEARCH · Engineering Papers

Results for “monotonicity”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 649 records · Page 36

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Variational Quantum Fidelity Estimation

Computing quantum state fidelity will be important to verify and characterize states prepared on a quantum computer. In this work, we propose novel lower and upper bounds for the fidelity F ( ρ , σ ) based on the ``truncated fidelity'' F ( ρ m , σ ) , which is evaluated for a state ρ m obtained by projecting ρ onto its m -largest eigenvalues. Our bounds can be refined, i.e., they tighten monotonically with m . To compute our bounds, we introduce a hybrid quantum-classical algorithm, called Variational Quantum Fidelity Estimation, that involves three steps: (1) variationally diagonalize ρ , (2) compute matrix elements of σ in the eigenbasis of ρ , and (3) combine these matrix elements to compute our bounds. Our algorithm is aimed at the case where σ is arbitrary and ρ is low rank, which we call low-rank fidelity estimation, and we prove that no classical algorithm can efficiently solve this problem under reasonable assumptions. Finally, we demonstrate that our bounds can detect quantum phase transitions and are often tighter than previously known computable bounds for realistic situations.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Grid-Forming Frequency Shaping Control for Low-Inertia Power Systems

As power systems transit to a state of high renewable penetration, little or no presence of synchronous generators makes the prerequisite of well-regulated frequency for grid-following inverters unrealistic. Thus, there is a trend to resort to grid-forming inverters which set frequency directly. We propose a novel grid-forming frequency shaping control that is able to shape the aggregate system frequency dynamics into a first-order one with the desired steady-state frequency deviation and Rate of Change of Frequency (RoCoF) after a sudden power imbalance. The no overshoot property resulting from the first-order dynamics allows the system frequency to monotonically move towards its new steady-state without experiencing frequency Nadir, which largely improves frequency security. We prove that our grid-forming frequency-shaping control renders the system internally stable under mild assumptions. The performance of the proposed control is verified via numerical simulations on a modified Icelandic Power Network test case.

control systems↗

DOE OSTI2021

Ripple-Type Control for Enhancing Resilience of Networked Physical Systems

Distributed control agents have been advocated as an effective means for improving the resiliency of our physical infrastructures under unexpected events. Purely local control has been shown to be insufficient, centralized optimal resource allocation approaches can be slow. In this context, we put forth a hybrid low-communication saturation-driven protocol for the coordination of control agents that are distributed over a physical system and are allowed to communicate with peers over a “hotline” communication network. According to this protocol, agents act on local readings unless their control resources have been depleted, in which case they send a beacon for assistance to peer agents. Our ripple-type scheme triggers communication locally only for the agents with saturated resources and it is proved to converge. Moreover, under a monotonicity assumption on the underlying physical law coupling control outputs to inputs, the devised control is proved to converge to a configuration satisfying safe operational constraints. The assumption is shown to hold for voltage control in electric power systems and pressure control in water distribution networks. Numerical tests corroborate the efficacy of the novel scheme.

distributed control↗

DOE OSTI2025

Physics-Informed Machine Learning Model for Ceramic Matrix Composite Creep

A physics-informed recurrent neural network (RNN) based surrogate model is developed to emulate the nonlinear, time-dependent constitutive behavior of ceramic matrix composites (CMCs) driven by matrix damage and constituent creep at the microscale. Physics-informed constraints are introduced into the surrogate model through regularization to ground the prediction in physics and improve its predictive capabilities. Training data is generated using the high-fidelity generalized method of cells (HFGMC) approach which calls appropriate creep and damage models for each of the constituents. This coupling permits simulating the nonlinear behavior of CMCs based on constituent response at the microscale along with microstructural features such as fiber and porosity volume fraction and fiber radius. The microscale repeating unit cell is loaded under creep fatigue conditions to replicate the material loading experienced in a turbine engine. Therefore, the RNN-based surrogate model is tasked with predicting, as a function of variable input stress sequence, temperature, and microstructural features, the resulting strain history response while satisfying physical constraints related to creep rate, isochoric inelastic deformation, and strain energy density. The trained surrogate model is shown to effectively match the strain history over quantified distributions of microstructural features and relevant loading regimes and temperatures. Neural network based surrogate models can offer efficient alternatives to running computationally intensive multiscale material models to simulate the nonlinear response of large structural models. Therefore, the presented work provides evidence towards the feasibility of developing, training, and running such models for CMCs with complex microstructures, nonlinear time-dependent material response, and under non-monotonic loading conditions.

ceramic matrix composites↗

DOE OSTI2020

How to design a focusing guide: The large moderator case1

As continuously shaped super-mirrors are becoming available, the conceptual design of focusing guides should explore a wider range of possibilities to accomplish an efficient neutron beam extraction. Starting from a desired phase-space volume at the sample position and using an upstream ray-tracing approach, the acceptance diagram of any focusing guide can be calculated at the moderator position. To ensure high brilliance transfer and homogeneous coverage, the acceptance diagram should be fully included in the neutron source emission phase-space volume. Following this idea, the guide system can be scaled into dimensionless geometric figures that convey performance limits for a desired cross-section reduction. Moreover, if we impose a monotonic increase of the reflection angle with divergence angle at the sample position, the shape of the mirror is analytically determined. This approach was applied in the design of a focusing guide for SNAP instrument at SNS, at ORNL, USA. The results of McStas simulations are presented with different options included. Our approach facilitates finding an optimal solution for connecting multiple guide pieces to avoid excessive losses and ensure a homogeneous phase space coverage.

Frost, Matthew↗

DOE OSTI2021

Fracture Network Localization Preceding Catastrophic Failure in Triaxial Compression Experiments on Rocks

We quantify the spatial distribution of fracture networks throughout six in situ X-ray tomography triaxial compression experiments on crystalline rocks at confining stresses of 5–35 MPa in order to quantify how fracture development controls the final macroscopic failure of the rock, a process analogous to those that control geohazards such as earthquakes and landslides. Tracking the proportion of the cumulative volume of fractures with volumes >90th percentile to the total fracture volume, ∑ v 90 / v t o t indicates that the fracture networks tend to increase in localization toward these largest fractures for up to 80% of the applied differential stress. The evolution of this metric also matches the evolution of the Gini coefficient, which measures the deviation of a population from uniformity. These results are consistent with observations of localizing low magnitude seismicity before large earthquakes in southern California. In both this analysis and the present work, phases of delocalization interrupt the general increase in localization preceding catastrophic failure, indicating that delocalization does not necessarily indicate a reduction of seismic hazard. However, the proportion of the maximum fracture volume to the total fracture volume does not increase monotonically. Experiments with higher confining stress tend to experience greater localization. To further quantify localization, we compare the geometry of the largest fractures, with volumes >90th percentile, to the best fit plane through these fractures immediately preceding failure. The r 2 scores and the mean distance of the fractures to the plane indicate greater localization in monzonite than in granite. The smaller mean mineral diameter and lower confining stress in the granite experiments may contribute to this result. Tracking these various metrics of localization reveals a close association between macroscopic yielding and the acceleration of fracture network localization. Near yielding, ∑ v 90 / v t o t and the Gini coefficient increase while the mean distance to the final failure plane decreases. Macroscopic yielding thus occurs when the rate of fracture network localization increases.

58 GEOSCIENCES↗

DOE OSTI2022

Effects of the Order Parameter Anisotropy on the Vortex Lattice in UPt3

We have used small-angle neutron scattering to determine the vortex lattice phase diagram in the topological superconductor UPt 3 for the applied magnetic field along the crystalline c -axis. A triangular vortex lattice is observed throughout the superconducting state, but with an orientation relative to the hexagonal basal plane that changes with field and temperature. At low temperature, in the chiral B phase, the vortex lattice undergoes a non-monotonic rotation with increasing magnetic field. The rotation amplitude decreases with increasing temperature and vanishes before reaching the A phase. Within the A phase an abrupt ±15° vortex lattice rotation was previously reported by Huxley et al. , Nature 406 , 160-164 (2000). The complex phase diagram may be understood from competing effects of the superconducting order parameter, the symmetry breaking field, and the Fermi surface anisotropy. The low-temperature rotated phase, centered around 0.8 T, reported by Avers et al. , Nature Physics 16 , 531-535 (2020), can be attributed directly to the symmetry breaking field.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Electromagnetic response of disordered superconducting cavities

We present the results for the resonant frequency shift and quality factor of disordered niobium (Nb) superconducting radio-frequency cavities driven out of equilibrium by the resonant microwave field. The theory is based on the non-equilibrium theory of superconductivity for the current response to the electromagnetic field at the vacuum–metal interface. We are able to accurately predict the observed frequency shifts with a precision of order fractions of kHz over the full temperature range 0 < T ≤ T c , including the negative frequency shift anomalies that are observed very near T c . The origin of these anomalies is shown to be the competition between the normal metal skin depth and the London penetration depth, which diverges as T→T c ⁻. An analytical approximation to the full current response, valid for |T-Tc| $\ll$ Tc, accounts for the negative frequency shift near T c . The non-monotonic dependence of the quality factor on the quasiparticle scattering rate is related to the pair-breaking effect of disorder on the superfluid fraction and, thus, the London penetration depth.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2021

Fluctuations and Pinning for Individually Manipulated Skyrmions

We numerically examine the dynamics of individually dragged skyrmions interacting simultaneously with an array of other skyrmions and quenched disorder. For drives just above depinning, we observe a broadband noise signal with a 1/ f characteristic, while at higher drives, narrowband or white noise appears. Even in the absence of quenched disorder, the threshold force that must be applied to translate the driven skyrmion is finite due to elastic interactions with other skyrmions. The depinning threshold increases as the strength of the quenched disorder is raised. Above the depinning force, the skyrmion moves faster in the presence of quenched disorder than in a disorder-free system since the pinning sites prevent other skyrmions from being dragged along with the driven skyrmion. For strong pinning, we find a stick-slip motion of the driven skyrmion which produces a telegraph noise signature. The depinning threshold increases monotonically with skyrmion density in the absence of quenched disorder, but when pinning is present, the depinning threshold changes nonmonotonically with skyrmion density, and there are reentrant pinned phases due to a competition between pinning induced by the quenched disorder and that produced by the elastic interactions of the skyrmion lattice.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Protocol-dependent frictional granular jamming simulations: cyclical, compression, and expansion

Granular matter takes many paths to pack in natural and industrial processes. The path influences the packing microstructure, particularly for frictional grains. We perform discrete element modeling simulations of different paths to construct packings of frictional spheres. Specifically, we explore four stress-controlled protocols implementing packing expansions and compressions in various combinations thereof. We characterize the eventual packed states through their dependence of the packing fraction and coordination number on packing pressure, identifying non-monotonicities with pressure that correlate with the fraction of frictional contacts. These stress-controlled, bulk-like particle simulations access very low-pressure packings, namely, the marginally stable limit, and demonstrate the strong protocol dependence of frictional granular matter.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Many-Body Effects in FeN4 Center Embedded in Graphene

We introduce a computational approach to study porphyrin-like transition metal complexes, bridging density functional theory and exact many-body techniques, such as the density matrix renormalization group (DMRG). We first derive a multi-orbital Anderson impurity Hamiltonian starting from first principles considerations that qualitatively reproduce generalized gradient approximation (GGA)+U results when ignoring inter-orbital Coulomb repulsion U ′ and Hund exchange J. An exact canonical transformation is used to reduce the dimensionality of the problem and make it amenable to DMRG calculations, including all many-body terms (both intra- and inter-orbital), which are treated in a numerically exact way. We apply this technique to FeN 4 centers in graphene and show that the inclusion of these terms has dramatic effects: as the iron orbitals become single occupied due to the Coulomb repulsion, the inter-orbital interaction further reduces the occupation, yielding a non-monotonic behavior of the magnetic moment as a function of the interactions, with maximum polarization only in a small window at intermediate values of the parameters. Furthermore, U ′ changes the relative position of the peaks in the density of states, particularly on the iron d z 2 orbital, which is expected to affect the binding of ligands greatly.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Theoretical Analysis of Fractional Viscoelastic Flow in Circular Pipes: General Solutions

Fractional calculus is a relatively old yet emerging field of mathematics with the widest range of engineering and biomedical applications. Despite being an incredibly powerful tool, it, however, requires promotion in the engineering community. Rheology is undoubtedly one of the fields where fractional calculus has become an integral part of cutting-edge research. There exists extensive literature on the theoretical, experimental, and numerical treatment of various fractional viscoelastic flows in constraint geometries. However, the general theoretical approach that unites several most commonly used models is missing. Here we present exact analytical solutions for fractional viscoelastic flow in a circular pipe. We find velocity profiles and shear stresses for fractional Maxwell, Kelvin–Voigt, Zener, Poynting–Thomson, and Burgers models. The dynamics of these quantities are studied with respect to normalized pipe radius, fractional orders, and elastic moduli ratio. Three different types of behavior are identified: monotonic increase, resonant, and aperiodic oscillations. The models developed are applicable in the widest material range and allow for the alteration of the balance between viscous and elastic properties of the materials.

Gritsenko, Dmitry (ORCID:0000000248181512)↗

DOE OSTI2026

Methanol Steam Reforming on Ru/m-ZrO 2 : Sodium Promotion of the CO 2 -Forming Pathway

Sodium (Na) promotion of Ru/m-ZrO 2 was investigated to elucidate how an alkali modification tunes selectivity in methanol steam reforming (MSR). H 2 -TPR/XANES/EXAFS show that Na increases surface basicity and strengthens Ru–O interactions, shifting RuO x reduction and H 2 spillover to a higher temperature. DRIFTS reveals Na-induced red shifts of the formate ν(CH) band and changes in OCO vibrational splitting, consistent with weakening of the formate C–H bond and an altered binding geometry. CO 2 -TPD confirms a monotonic shift toward stronger basic sites with increasing Na concentrations. Under MSR conditions, Na selectively increases CO2 concentration at the expense of CO. At ~80% conversion and 325 °C, CO 2 selectivity increases from 12.0% (unpromoted) to 16.2, 21.0, and 26.5% for 0.5, 1.0, and 1.8% Na, respectively; at ~300 °C and ~66–69% conversion, CO 2 selectivity increases from 8.6% to 23.7% at 1.8% Na. Transient MSR experiments further show earlier and larger H 2 evolution upon Na addition, corroborating the promotion of the dehydrogenation/decarbonylation route to CO 2 + H 2 . We propose that Na increases basicity and modifies the Ru–support interface to favor formate dehydrogenation/decarboxylation, thereby increasing the H 2 yield and lowering CO formation. Ru’s higher-energy, less occupied d-band stabilizes CO and oxygenated intermediates more strongly in the reforming environment, making the CO-forming pathway more resistant to suppression than on Pt.

CO2 selectivity↗

DOE OSTI2020

Effect of Computational Schemes on Coupled Flow and Geo-Mechanical Modeling of CO 2 Leakage through a Compromised Well

Carbon capture, utilization, and storage (CCUS) describes a set of technically viable processes to separate carbon dioxide (CO 2 ) from industrial byproduct streams and inject it into deep geologic formations for long-term storage. Legacy wells located within the spatial domain of new injection and production activities represent potential pathways for fluids (i.e., CO 2 and aqueous phase) to leak through compromised components (e.g., through fractures or micro-annulus pathways). The finite element (FE) method is a well-established numerical approach to simulate the coupling between multi-phase fluid flow and solid phase deformation interactions that occur in a compromised well system. We assumed the spatial domain consists of a three-phases system: a solid, liquid, and gas phase. For flow in the two fluids phases, we considered two sets of primary variables: the first considering capillary pressure and gas pressure (PP) scheme, and the second considering liquid pressure and gas saturation (PS) scheme. Fluid phases were coupled with the solid phase using the full coupling (i.e., monolithic coupling) and iterative coupling (i.e., sequential coupling) approaches. The challenge of achieving numerical stability in the coupled formulation in heterogeneous media was addressed using the mass lumping and the upwinding techniques. Numerical results were compared with three benchmark problems to assess the performance of coupled FE solutions: 1D Terzaghi’s consolidation, Liakopoulos experiments, and the Kueper and Frind experiments. We found good agreement between our results and the three benchmark problems. For the Kueper and Frind test, the PP scheme successfully captured the observed experimental response of the non-aqueous phase infiltration, in contrast to the PS scheme. These exercises demonstrate the importance of fluid phase primary variable selection for heterogeneous porous media. We then applied the developed model to the hypothetical case of leakage along a compromised well representing a heterogeneous media. Considering the mass lumping and the upwinding techniques, both the monotonic and the sequential coupling provided identical results, but mass lumping was needed to avoid numerical instabilities in the sequential coupling. Additionally, in the monolithic coupling, the magnitude of primary variables in the coupled solution without mass lumping and the upwinding is higher, which is essential for the risk-based analyses.

deformation flow↗

DOE OSTI2020

Leveraging Cell Expansion Sensing in State of Charge Estimation: Practical Considerations

Measurements such as current and terminal voltage that are typically used to determine the battery’s state of charge (SOC) are augmented with measured force associated with electrode expansion as the lithium intercalates in its structure. The combination of the sensed behavior is shown to improve SOC estimation even for the lithium ion iron phosphate (LFP) chemistry, where the voltage–SOC relation is flat (low slope) making SOC estimation using measured voltage difficult. For the LFP cells, the measured force has a non-monotonic F–SOC relationship. This presents a challenge for estimation as multiple force values can correspond to the same SOC. The traditional linear quadratic estimator can be driven to an incorrect SOC value. To address these difficulties, a novel switching estimation gain is used based on determining the operating region that corresponds to the actual SOC. Moreover, a drift in the measured force associated with a shift of the cell SOC–expansion behavior over time is addressed with a bias estimator for the force signal. The performance of Voltage-based (V) and Voltage and Force-based (V&F) SOC estimation algorithms are then compared and evaluated against a desired ±5% absolute error bound of the SOC using a dynamic stress test current protocol that tests the proposed estimation scheme across wide range of SOC and current rates.

25 ENERGY STORAGE↗

DOE OSTI2021

Evaluation of Formate Salt PCM’s for Latent Heat Thermal Energy Storage

This work examines formate salts as potential phase change materials (PCMs) for middle-high temperature (≤250 °C) latent heat thermal energy storage applications. The thermophysical properties of three formate salts were characterized: pure sodium formate and binary blends of sodium/potassium formate and sodium/calcium formate. The stability of formate PCM’s was evaluated by thermal cycling using differential scanning calorimetry where sodium formate and sodium/potassium formate appeared stable over 600 cycles, while sodium/calcium formate exhibited a monotonic decrease in heat of fusion over the test period. A longer test with sodium formate led to gas release and decomposition of the salt. FTIR analysis of the PCM showed degradation of formate to oxalate. T-history experiments with 50-g PCM quantities demonstrated a bulk supercooling of only 2–3 °C for these salts. Thermal conductivity enhancement of over 700% was achieved by embedding aluminum in the solid PCM. Finally, mild carbon steel was immersed in molten sodium formate for up to 2000 h. Sodium formate was found to be non-corrosive, as calculated by mass loss and confirmed by cross-sectional high-resolution microscopy. FTIR analysis of the PCM after 2000 h shows oxidation at the free surface, while the bulk PCM remained unchanged, further indicating a need to protect the formate from atmospheric exposure when used as a PCM.phase change materials; formate salts; latent heat thermal energy storage; thermal cycling; supercooling; thermal conductivity enhancement; corrosion

25 ENERGY STORAGE↗

DOE OSTI2021

A Lamination Model for Pressure-Assisted Sintering of Multilayered Porous Structures

This work describes a lamination model for pressure-assisted sintering of thin, multilayered, and porous structures based on the linear viscous constitutive theory of sintering and the classical laminated plate theory of continuum mechanics. A constant out-of-plane normal stress is assumed in the constitutive relation. The lamination relations between the force/moment resultants and the strain/curvature rates are presented. Numerical simulations were performed for a symmetric tri-layer laminate consisting of a 10% gadolinia doped ceria (Ce 0.9 Gd 0.1 O 1.95-δ ) composite structure, where porous layers were adhered to the top and bottom of a denser layer under uniaxially-applied pressures and the sinter forging conditions. The numerical results show that, compared with free sintering, the applied pressure can significantly reduce the sintering time required to achieve given layer thicknesses and porosities. Unlike free sintering, which results in a monotonic decrease of the laminate in-plane dimension, pressure-assisted sintering may produce an in-plane dimension increase or decrease, depending on the applied pressure and sintering time. Finally, the individual layers in the laminate exhibit different stress characteristics under pressure-assisted sintering.

42 ENGINEERING↗

DOE OSTI2023

Effect of Production Bias on Radiation-Induced Segregation in Ni-Cr Alloys

We present an in-depth investigation into the Radiation-Induced Segregation (RIS) phenomenon in Ni-Cr alloys. All the pivotal factors affecting RIS such as surface’s absorption efficiency, grain size, production bias, dose rate, temperature, and sink density were systematically studied. Through comprehensive simulations, the individual and collective impacts of these factors were analyzed, enabling a refined understanding of RIS. A notable finding was the significant influence of production bias on point defects’ interactions with grain boundaries/surfaces, thereby playing a crucial role in RIS processes. Production bias alters the neutrality of these interactions, leading to a preferential absorption of one type of point defect by the boundary and consequent establishment of distinct surface-mediated patterns of point defects. These spatial patterns further result in non-monotonic spatial profiles of solute atoms near surfaces/grain boundaries, corroborated by experimental observations. In particular, a positive production bias, signifying a higher production rate of vacancies over interstitials, drives more Cr depletion at the grain boundary. Moreover, a temperature-dependent production bias must be considered to recover the experimentally reported dependence of RIS on temperature. The severity of radiation damage and RIS becomes more pronounced with increased production bias, dose rate, and grain size, while high temperatures or sink density suppress the RIS severity. Model predictions were validated against experimental data, showcasing robust qualitative and quantitative agreements. The findings pave the way for further exploration of these spatial dependencies in subsequent studies, aiming to augment the comprehension and predictability of RIS processes in alloys.

36 MATERIALS SCIENCE↗