Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 1,009 records · Page 56

Adaptive Power Flow Approximations With Second-Order Sensitivity Insights

The power flow equations are fundamental to power system planning, analysis, and control. However, the inherent non-linearity and non-convexity of these equations present formidable obstacles in problem-solving processes. To mitigate these challenges, recent research has proposed adaptive power flow linearizations that aim to achieve accuracy over wide operating ranges. The accuracy of these approximations inherently depends on the curvature of the power flow equations within these ranges, which necessitates considering second-order sensitivities. In this paper, we leverage second-order sensitivities to both analyze and improve power flow approximations. We evaluate the curvature across broad operational ranges and subsequently utilize this information to inform the computation of various sample-based power flow approximation techniques. Additionally, we leverage second-order sensitivities to guide the development of rational approximations that yield linear constraints in optimization problems. In conclusion, this approach is extended to enhance accuracy beyond the limitations of linear functions across varied operational scenarios.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Estimation and Visualization of Isosurface Uncertainty from Linear and High-Order Interpolation Methods

Isosurface visualization is fundamental for exploring and analyzing 3D volumetric data. Marching cubes (MC) algorithms with linear interpolation are commonly used for isosurface extraction and visualization. Although linear interpolation is easy to implement, it has limitations when the underlying data is complex and high-order, which is the case for most real-world data. Linear interpolation can output vertices at the wrong location. Its inability to deal with sharp features and features smaller than grid cells can lead to an incorrect isosurface with holes and broken pieces. Despite these limitations, isosurface visualizations typically do not include insight into the spatial location and the magnitude of these errors. We utilize high-order interpolation methods with MC algorithms and interactive visualization to highlight these uncertainties. Our visualization tool helps identify the regions of high interpolation errors. It also allows users to query local areas for details and compare the differences between isosurfaces from different interpolation methods. In addition, we employ high-order methods to identify and reconstruct possible features that linear methods cannot detect. We showcase how our visualization tool helps explore and understand the extracted isosurface errors through synthetic and real-world data.

Ouermi, Timbwaoga↗

Thermodynamic phase transitions of nematic order in magnetic helices

A nematic phase lacks translation order but has orientational order. Nematic phases have been discovered in a variety of systems, including liquid crystals, correlated materials, and superconductors. Here, we report on a magnetic nematic phase, where the basis components are composed of magnetic helices. We directly probed the order parameters associated with the magnetic helices using resonant soft x-ray scattering and find two distinct nematic phases with complex spatiotemporal signatures. Using x-ray correlation spectroscopy, we find that near the phase boundary between the two nematic phases, fluctuations coexist on multiple disparate timescales. Our micromagnetic simulations and density functional theory calculations show that the fluctuations occur concomitantly with a reorientation of the magnetic helices, indicating spontaneous symmetry breaking and the emergence of additional degrees of freedom. Our results provide a framework for characterizing exotic phases that can be extended to a broad class of physical systems.

36 MATERIALS SCIENCE↗

High-Order Mesh r-Adaptivity with Tangential Relaxation and Guaranteed Mesh Validity

High-order meshes are crucial for achieving optimal convergence rates in curvilinear domains, preserving symmetry, and aligning with key flow features in moving mesh simulations [1], but their quality is challenging to control. In prior work, we have developed techniques based on Target-Matrix Optimization Paradigm (TMOP) to adapt a given high-order mesh to the geometry and solution of the partial differential equation (PDE) [2, 3]. Here, we extend this framework to address two key gaps in the literature for highorder mesh 𝑟-adaptivity. First, we introduce tangential relaxation on curved surfaces using solely the discrete mesh representation, eliminating the need for access to underlying geometry (e.g., CAD model). Second, we ensure a continuously positive Jacobian determinant throughout the domain. This determinant positivity is essential for using the high-order mesh resulting from 𝑟-adaptivity with arbitrary quadrature schemes in simulations. The proposed approach is demonstrated to be robust using a variety of numerical experiments.

Mathematics and Computing↗

Multiscale Control of Generic Second Order Traffic Models by Driver-Assist Vehicles

We study the derivation of generic high order macroscopic traffic models from a follow-the-leader particle description via a kinetic approach. First, we recover a third order traffic model as the hydrodynamic limit of an Enskog-type kinetic equation. Next, we introduce in the vehicle interactions a binary control modeling the automatic feedback provided by driver-assist vehicles and we upscale such a new particle description by means of another Enskog-based hydrodynamic limit. The resulting macroscopic model is now a generic second order model (GSOM), which contains in turn a control term inherited from the microscopic interactions. We show that such a control may be chosen so as to optimize global traffic trends, such as the vehicle flux or the road congestion, constrained by the GSOM dynamics. By means of numerical simulations, we investigate the effect of this control hierarchy in some specific case studies, which exemplify the multiscale path from the vehiclewise implementation of a driver-assist control to its optimal hydrodynamic design.

GSOM↗

Lowest-order QED radiative corrections in unpolarized elastic electron–deuteron scattering beyond the ultra-relativistic limit for the proposed deuteron charge radius measurement at Jefferson laboratory

Analogous to the well-known proton charge radius puzzle, a similar puzzle exists for the deuteron charge radius, r d . There are discrepancies observed in the results of r d , measured from electron-deuteron (e - d) scattering experiments, as well as from atomic spectroscopy. In order to help resolve the charge radius puzzle of the deuteron, the PRad collaboration at Jefferson Lab has proposed an experiment for measuring r d , named DRad. This experiment is designed to measure the unpolarized elastic e - d scattering cross section in a low-Q 2 region. To extract the cross section with a high precision, having reliable knowledge of QED radiative corrections is important. Here, in this paper, we present complete numerical calculations of the lowest-order radiative corrections in e - d scattering for the DRad kinematics. The calculations have been performed within a covariant formalism and beyond the ultra-relativistic approximation ($m^2_e$ $\ll$ Q 2 ). Besides, we present a systematic uncertainty on r d arising from higher-order radiative corrections, estimated based on our cross-section results.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Domain-decomposition nonlinear manifold reduced order model

This software combines nonlinear-manifold reduced order models (NM-ROMs) with domain decomposition (DD) techniques. NM-ROMs, which utilize a shallow, sparse autoencoder trained with full order model (FOM) snapshot data, approximate the FOM state on a nonlinear manifold. These models offer advantages over linear-subspace ROMs (LS-ROMs) particularly in scenarios with slowly decaying Kolmogorov n-width. However, the training of NM-ROMs involves a number of parameters that scale with the size of the FOM, and storing high-dimensional FOM snapshots can significantly increase the cost of ROM training for extreme-scale problems. To mitigate these costs, the software employs DD to partition the FOM into smaller subdomains, computes NM-ROMs for each, and then integrates these to form a global NM-ROM. This strategy offers multiple benefits: it enables parallel training of subdomain NM-ROMs, reduces the number of parameters needed, decreases the dimensional requirements of subdomain FOM training data, and allows for customization to the unique characteristics of each FOM subdomain. The use of a shallow, sparse autoencoder architecture in each subdomain NM-ROM facilitates the application of hyper-reduction (HR), simplifying the nonlinear complexities and enhancing computational speed. This software marks the inaugural application of NM-ROM combined with HR to a DD problem. It features an algebraic DD reformulation of the FOM, training of NM-ROMs with HR for each subdomain, and employs a sequential quadratic programming (SQP) solver for the evaluation of the coupled global NMROM. The effectiveness of the DD NM-ROM with HR is numerically demonstrated on the 2D steady-state Burgers' equation, showing an order of magnitude improvement in accuracy over the DD LS-ROM with HR.

Diaz, AlejandroN↗

Large-scale cascading of first-order FBG array in a highly multimode coreless fiber using femtosecond laser for distributed thermal sensing

This research focuses on the performance analysis and characterization of a fiber Bragg gratings (FBGs) array, consisting of 10 first-order FBGs inscribed by a femtosecond (FS) laser in a highly multimode coreless fiber. The study evaluates the FBG array's ability to function as a distributed thermal sensing (DTS) platform, with the coreless fiber chosen as the sensing element due to its immunity to dopant migration at high temperatures. The design of a large cascaded first-order FBG array effectively eliminates unwanted harmonic peaks across a wide spectrum range. In contrast, higher-order FBGs introduce limitations due to the overlapping of Bragg peaks with harmonics. The FBG array's performance is evaluated by measuring the reflection spectrum of each grating at different temperatures, showing a high temperature sensitivity of 15.05 pm/°C at a Bragg wavelength of 1606.3 nm, with a linear response in the temperature range of 24 - 1100 °C. The FBG array was designed for a spatial resolution of 5 mm. A mode scrambler in the sensing network is employed, which suppresses multimodal interference, characterizes FBG peak visibility, and stabilizes the interference spectrum. The stability of the FBG array is also assessed over 24 hrs at 1100 °C, and it is observed to be stable during thermal treatment. Heat treatment at 1100°C improves the signal to noise ratio of the FBG array, demonstrating the robustness and suitability of the proposed FBG array on highly multimode coreless fiber as a potential sensing platform for DTS applications in harsh environmental conditions, overcoming the issues of dopant migration presented by dopes silica optical fibers at high temperatures.

47 OTHER INSTRUMENTATION↗

Third-order photon correlations extract single-nanocrystal multiexciton properties in solution

Colloidal semiconductor nanocrystals are considered promising materials for high-flux optical applications, including lasing, light-emitting diodes, biological imaging, and quantum optics. In high-flux applications, multiexcitons can significantly contribute to emission, influencing its brightness, spectral purity, and kinetics. As a result, understanding and controlling multiexciton emission in colloidal nanocrystal materials is of the utmost importance. In the past, single-nanocrystal photon correlation methods have been applied to understand biexciton and triexciton efficiencies, lifetimes, and spectra. While powerful, such methods suffer from user selection bias and require stable emission from single nanocrystals. To compensate for this shortcoming, second-order correlation methods were developed to extract sample-averaged biexciton properties from a solution of nanocrystals. Until now, however, the analogous third-order solution photon correlation methods remained unexplored. In this work, we present a pair of third-order photon correlation techniques to obtain the sample-averaged single-nanocrystal triexciton quantum yield and lifetime in a solution-phase experiment. These techniques derive from the relationship between the Poisson probability of nanocrystal photon absorption and the intrinsic probability of nanocrystal photon emission. We validate the theoretical background of these techniques by creating a numerical model to simulate the diffusion and emission of many nanocrystals in solution. Our simulations confirm that the average triexciton quantum yield and triexciton lifetime can be extracted from a solution of nanocrystals. These techniques will enable researchers to gain a better understanding of the fundamental multiexciton properties of colloidal nanocrystals.

Horowitz, Jonah R. [Massachusetts Institute of Tec↗

Reconsidering short-range order in complex concentrated alloys

Abstract The seemingly contradictory state of research on short-range order in many-component alloys is addressed through a critical review of the characterization of face-centered-cubic 3 d systems. Despite the paucity of direct observations, the ordering of many widely studied alloys is argued to be primarily interesting for its potential ubiquity. To clarify this situation, future research directions are proposed with reference to historical results, including a review of the fundamental principles of ordering and clustering. Graphical abstract

36 MATERIALS SCIENCE↗

Simultaneous Compensation of Third-Order Resonances at the FNAL Recycler Ring

Third-order resonance lines will have a detrimental effect on the high-intensity operation of the Recycler Ring (RR), under the current Proton Improvement Plan (PIP-II) for the Fermilab Accelerator Complex. Increasing intensity will increase space charge effects, leading to the excitation of normal and skew sextupole lines. Dedicated normal and skew sextupoles have been installed in order to mitigate the effect of these resonance lines. By measuring the response matrix of the third-order Resonance Driving Terms (RDTs) to the currents of these dedicated elements, this study shows how several resonance lines can be compensated simultaneously. Resonance compensation is experimentally verified through loss maps and emittance growth measurements using the Ion Profile Monitor (IPM) system in the Recycler.

43 PARTICLE ACCELERATORS↗

Predicted third-order sweet spots for phi-junction Josephson parametric amplifiers

Hybrid superconductor-semiconductor nanowire Josephson junctions exhibit skewed and phi-shifted current phase relations when an in-plane magnetic field is applied along the weak link's spin-orbit effective field direction. These junctions can have an asymmetric Josephson potential with odd-order nonlinearities. A dominant third-order nonlinearity can be achieved by tuning the magnetic field to a sweet spot. Sweet spots persist when higher order Josephson harmonics are included. This makes it possible to have a single Josephson junction dipole element with three-wave mixing capability, which is favorable for pump-efficient amplification. Electrostatic gate tunability of the semiconductor weak link can make it operable within an extended range of working frequencies, and the inclusion of micromagnets can facilitate near-zero magnetic field operation.

Reza, Tasnum [Univ. of Pittsburgh, PA (United Stat↗

Spectral Equivalence Properties of Higher-Order Tensor Product Finite Elements and Applications to Preconditioning

In this study, we present spectral equivalence results for higher-order tensor product edge-, face- and interior-based finite elements. Specifically, we show for certain choices of shape functions that the mass and stiffness matrices of the higher-order elements are spectrally equivalent to those for an assembly of lowest-order elements on the associated Gauss-Lobatto-Legendre mesh. Based on this equivalence, efficient preconditioners can be designed with favorable computational complexity. Numerical results are presented which confirm the theory and demonstrate the benefits of the equivalence results for overlapping Schwarz preconditioners.

97 MATHEMATICS AND COMPUTING↗

A 6th Order Mehrstellen Finite Volume Discretization of Poisson's Equation in Three Dimensions

We discuss the derivation of a new, sixth-order finite volume scheme for Poisson’s equation on 3D Cartesian equispaced grids. The scheme is based on a discretization of the Laplace operator with a compact (Mehrstellen) 27-point stencil. To achieve sixth order convergence the right hand side of the equation is replaced with a discrete operator that involves the discrete Laplace and Biharmonic operators and the sum of discrete fourth-order cross derivatives applied to the charge function. Numerical tests demonstrate the superiority of the proposed method compared to the well known schemes associated with the 7-point and 19-point discretizations of the Laplacian.

97 MATHEMATICS AND COMPUTING↗

Topological Materials with Complex Long-Range Order

The overarching vision of this proposal is to understand exotic metallic behavior in frustrated quantum magnets. This problem is at the core of understanding a number of mysteries in the study of itinerant magnetism, and of high temperature superconductivity. The central challenge is that magnetic order is generally understood by localized magnetic moments. In these systems the mobile electrons form an ordered magnetic state, sometime breaking sharp symmetry, sometimes forming incommensurate order, and other times forming a spin liquid. The present research focuses on finding complex magnets that are in close proximity to a sea of metallic, itinerant electrons.

36 MATERIALS SCIENCE↗

High-order algorithmic developments and optimizations for large-scale GPU-accelerated simulations (Milestone CEED-MS36)

The goal of this milestone was to improve the high-order software ecosystem for CEED-enabled ECP applications by making progress on efficient matrix-free kernels targeting forthcoming ECP architectures. These kernels included matrix-free preconditioning and the development of new set of CEED solver bake-off problems. As part of this milestone, we also released the next version of the CEED software stack, CEED-4.0, reported on results from several application collaborations, and documented the efforts of porting to AMD GPUs for Frontier and other modern architectures, such as Fugaku. The specific tasks addressed in this milestone were: (1) Port and run CEED benchmarks/miniapps on Frontier EA systems; (2) Demonstrate performant libCEED integration in MFEM, Nek and applications; (3) Matrix-free preconditioning of high-order operators; (4) Benchmark problems for fast high-order solvers on GPU platforms; and (5) Public release of CEED-4.0. The artifacts delivered include the next version of the CEED software stack, CEED-4.0, the next libCEED release, libCEED-0.8, and a number of developments integrated within applications to improve their GPU and CPU performance and capabilities. See the CEED website, https://ceed.exascaleproject.org and the CEED GitHub organization, https://github.com/ceed for more details.

97 MATHEMATICS AND COMPUTING↗

(U) Calculating First-Order Sensitivities to Material Density and Composition in Fixed-Source Problems Using the MCNP6 Perturbation Capability

The Taylor series (differential operator) perturbation method, as implemented in the PERT capability in MCNP6 (Ref. 1), can be used for first-order sensitivity analyses in fixed-source problems.2 Sensitivities of a response with respect to material and nuclide densities can be used for efficient first-order uncertainty quantification. This report discusses how to compute sensitivities of material densities and compositions in order to propagate uncertainties using the standard uncertainty propagation formula.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A Micromorphic Length-Scale Coupling Framework for the Determination of Higher-Order Constitutive Models and the Multi-Scale Simulation of Heterogeneous Materials [Thesis]

Heterogeneous materials and materials with complex microstructures pose a unique challenge in the development of accurate models of their response to external stimuli. These difficulties principally arise due to the difficulty in characterizing and modeling the constituents and their interactions. It is usually possible, though non-trivial, to construct an explicit representation of the microstructure (a direct numeric simulation or DNS) but the method by which the complex modes of deformation and other processes can be homogenized to a reduced-order approximation is frequently unclear. Many homogenization approaches are ad-hoc and lack a strong justification beyond ease of computation. Furthermore, the homogenization approach can, in some cases, not utilize the full breadth of information available in the computation of the macro-scale stresses and deformation measures. It is also noteworthy that homogenization, by its very nature, will tend to obfuscate details of processes occurring at the lower length-scale. It is therefore of interest to include as much information as possible in the construction of the reduced-order model so as to be predictive in a variety of loading environments. We here present a length-scale bridging technique based upon the micromorphic continuum mechanics of Eringen which incorporates volume and surface area averages as a part of its construction. This approach enforces the balance equations at the micro-scale and then studies the effect of the spatially varying nature on the macro-scale. This leads, naturally, to further balance equations which are solved at the macro-scale. This work details a homogenization framework which arises naturally from the micromorphic construction of Eringen attempting to introduce no definitions beyond which are justifiable from micro-structural considerations. One of the results of this effort is the so-called “micromorphic filter” which has been developed and applied to several DNS to demonstrate its effectively. In order to determine the macroscopic degrees of freedom we utilize the special case of an overlap coupling technique where the macro-scale is fully constrained to the micro-scale. This enables us to study the resulting material properties as expected but also allows us to further study the boundary conditions on the additional degrees of freedom at the macro-scale.

42 ENGINEERING↗