Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear projection”

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 37 records · Page 2

Gradient-based constrained optimization using a database of linear reduced-order models

A methodology grounded in model reduction is presented for accelerating the gradient-based solution of a family of linear or nonlinear constrained optimization problems where the constraints include at least one linear Partial Differential Equation (PDE). A key component of this methodology is the construction, during an offline phase, of a database of pointwise, linear, Projection-based Reduced-Order Models (PROM)s associated with a design parameter space and the linear PDE(s). A parameter sampling procedure based on an appropriate saturation assumption is proposed to maximize the efficiency of such a database of PROMs. A real-time method is also presented for interpolating at any queried but unsampled parameter vector in the design parameter space the relevant sensitivities of a PROM. The practical feasibility, computational advantages, and performance of the proposed methodology are demonstrated for several realistic, nonlinear, aerodynamic shape optimization problems governed by linear aeroelastic constraints.

97 MATHEMATICS AND COMPUTING↗

Nucleon-pair truncation of the shell model for medium-heavy nuclei

Background: Establishing computationally tractable models of atomic nuclei is a long-time goal of nuclear structure physics. A flexible framework which easily includes excited states and many-body correlations is the configuration-interaction shell model (SM), but the exponential growth of the basis means one needs an efficient truncation scheme, ideally one that includes both deformation and pairing correlations. Purpose: Here, we propose an efficient truncation scheme of the SM: starting from a pair condensate variationally defined by Hartree-Fock single-particle states and the particle-number conserved Bardeen-Cooper-Schrieffer (NBCS) approximation, we carry out projection of states with good angular momentum. Methods: After generating Hartree-Fock single-particle states with Kramers degeneracy in a SM space, we optimize the pair amplitudes in the NBCS by minimizing the energy, and then use linear algebra projection (LAP) of states with good angular momentum. Both NBCS and LAP are computationally fast. Results: Our calculations yield good agreement with full configuration-interaction SM calculations for low-lying states of transitional and rotational nuclei with axially symmetric and triaxial deformation in medium- and heavy-mass regions: 44,46,48 Ti, 48,50 Cr, 52 Fe, 60,62,64 Zn, 66,68 Ge, 68 Se, and 108,110 Xe. We predict low-lying states of 112–114 Ba and 116–120 Ce, nuclei difficult to reach by large-scale SM calculations. Conclusions: Both pair correlation and the configuration mixing between different intrinsic states play a key role in reproducing collectivity and shape coexistence, demonstrating the utility of this truncation scheme of the SM to study transitional and deformed nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Manifold projection image segmentation for nano-XANES imaging

As spectral imaging techniques are becoming more prominent in science, advanced image segmentation algorithms are required to identify appropriate domains in these images. We present a version of image segmentation called manifold projection image segmentation (MPIS) that is generally applicable to a broad range of systems without the need for training because MPIS uses unsupervised machine learning with a few physically motivated hyperparameters. We apply MPIS to nanoscale x-ray absorption near edge structure (XANES) imaging, where XANES spectra are collected with nanometer spatial resolution. We show the superiority of manifold projection over linear transformations, such as the commonly used principal component analysis (PCA). Moreover, MPIS maintains accuracy while reducing computation time and sensitivity to noise compared to the standard nano-XANES imaging analysis procedure. Finally, we demonstrate how multimodal information, such as x-ray fluorescence data and spatial location of pixels, can be incorporated into the MPIS framework. We propose that MPIS is adaptable for any spectral imaging technique, including scanning transmission x-ray microscopy, where the length scale of domains is larger than the resolution of the experiment.

36 MATERIALS SCIENCE↗

Multifidelity Ensemble Kalman Filtering Using Surrogate Models Defined by Theory-Guided Autoencoders

Data assimilation is a Bayesian inference process that obtains an enhanced understanding of a physical system of interest by fusing information from an inexact physics-based model, and from noisy sparse observations of reality. The multifidelity ensemble Kalman filter (MFEnKF) recently developed by the authors combines a full-order physical model and a hierarchy of reduced order surrogate models in order to increase the computational efficiency of data assimilation. The standard MFEnKF uses linear couplings between models, and is statistically optimal in case of Gaussian probability densities. This work extends the MFEnKF into to make use of a broader class of surrogate model such as those based on machine learning methods such as autoencoders non-linear couplings in between the model hierarchies. We identify the right-invertibility property for autoencoders as being a key predictor of success in the forecasting power of autoencoder-based reduced order models. We propose a methodology that allows us to construct reduced order surrogate models that are more accurate than the ones obtained via conventional linear methods. Numerical experiments with the canonical Lorenz'96 model illustrate that nonlinear surrogates perform better than linear projection-based ones in the context of multifidelity ensemble Kalman filtering. We additionality show a large-scale proof-of-concept result with the quasi-geostrophic equations, showing the competitiveness of the method with a traditional reduced order model-based MFEnKF.

97 MATHEMATICS AND COMPUTING↗

Long-Term Corrosion Studies Of Pine Derived Bio-Oil And Blends With Heavy Fuel Oil

Biomass derived liquid fuels offer a means to reduce greenhouse gas emissions compared to those produced by combustion of petroleum derived liquid fuels. However, the corrosivity of bio-oils toward the less expensive structural materials creates a material selection problem for designers of storage tanks and combustion systems. Samples of candidate structural materials are being exposed for thousands of hours in multiple fast pyrolysis bio-oils and conditions to evaluate the corrosion resistance of these materials. One method to mitigate the corrosivity of bio-oils and speed their adoption, while also decreasing the pollution issues associated with low quality petroleum derived fuels, is to utilize blends of bio-oil and heavy fuel oil in engines currently solely burning a petroleum-based fuel. In addition to the corrosion studies conducted in 100% bio-oil, studies were also conducted with blends of the fast pyrolysis bio-oil produced from pine tree components with a heavy fuel oil that is used in ocean-going ships. This bio-oil had a very high carboxylic acid content which made it very corrosive to carbon and 2¼ Cr-1 Mo steel and even 409 stainless steel. The heavy fuel oil was not corrosive to carbon steel, but its sulfur content makes it a significant pollution producer and particularly undesirable for use near coastlines. Corrosion tests were conducted with the individual components and with various blends of the two liquid fuels. Studies showed a significantly lower corrosivity of the blends than would be projected assuming linear mixing behavior. Adoption of such blends holds the potential to reduce production of sulfur-containing exhaust gases as well as carbon dioxide from non-renewable fuels.

Keiser, Jim↗

Long-Term Corrosion Studies of Pine Derived Bio-Oil and Blends with Heavy Fuel Oil

Biomass derived liquid fuels offer a means to reduce greenhouse gas emissions compared to those produced by combustion of petroleum derived liquid fuels. However, the corrosivity of bio-oils toward the less expensive structural materials creates a material selection problem for designers of storage tanks and combustion systems. Samples of candidate structural materials are being exposed for thousands of hours in multiple fast pyrolysis bio-oils and conditions to evaluate the corrosion resistance of these materials. One method to mitigate the corrosivity of bio-oils and speed their adoption, while also decreasing the pollution issues associated with low quality petroleum derived fuels, is to utilize blends of bio-oil and heavy fuel oil in engines currently solely burning a petroleum-based fuel. In addition to the corrosion studies conducted in 100% bio-oil, studies were also conducted with blends of the fast pyrolysis bio-oil produced from pine tree components with a heavy fuel oil that is used in ocean-going ships. This bio-oil had a very high carboxylic acid content which made it very corrosive to carbon and 2¼ Cr-1 Mo steel and even 409 stainless steel. The heavy fuel oil was not corrosive to carbon steel, but its sulfur content makes it a significant pollution producer and particularly undesirable for use near coastlines. Corrosion tests were conducted with the individual components and with various blends of the two liquid fuels. Studies showed a significantly lower corrosivity of the blends than would be projected assuming linear mixing behavior. Adoption of such blends holds the potential to reduce production of sulfur-containing exhaust gases as well as carbon dioxide from non-renewable fuels.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

RPI LINAC refurbishment and upgrade project

The Electron Linear Accelerator at Rensselaer Polytechnic Institute (RPI) was constructed in the late 1950s and started operation in 1961 [1]. The facility was designed and used primarily for measurements of Nuclear Data (ND) with emphasis on the resonance energy region. As such it was designed to produce short pulses of neutrons (5-5000 ns) and equipped with detection stations at different neutron flight path distances from the source (15 m-250 m). When the LINAC started operation it was one of the most powerful accelerators in the world for high energy-resolution measurements in the resonance region. The original LINAC with some upgrades is still in operation but no longer provides world class capabilities for measurements with short pulses. In 2011 discussions between the RPI ND group, Naval Reactors, and the Nuclear Criticality Safety Program about the future of the LINAC concluded that maintaining the capability to perform world-class nuclear data measurements at RPI was desirable. To that end, it was decided to refurbish the facility and upgrade the accelerator capabilities to enable high accuracy measurement in the resonance region with narrow neutron pulse and high electron beam power. Located at a major US technical university, this project provides a needed US nuclear data capability for both research and education, and also includes external users and industrial applications.

43 PARTICLE ACCELERATORS↗

Status of the International Linear Collider

This paper is not a proposal for a CERN future project but provides information on the International Linear Collider (ILC) considered for Japan in order to facilitate the European Strategy discussion in a global context. It describes progress to date, ongoing engineering studies, updated cost estimate for the machine at $\sqrt{s}=250~\rm GeV$ and the situation in Japan. The physics of the ILC is not presented here, but jointly for all Linear Collider projects in a separate document ``A Linear Collider Vision for the Future of Particle Physics'' submitted for the forthcoming European Strategy deliberations.

Abe, Y. [KEK, Tsukuba]↗

Gluon unpolarized, polarized, and transversity GPDs from lattice QCD: Lorentz-covariant parametrization

We identify the matrix elements necessary to determine the leading-twist gluon generalized parton distributions (GPDs) H g , E g , H ˜ g , E ˜ g , H g T , E g T , H ˜ g T , E ˜ g T in lattice QCD calculations. We present a method to achieve a Lorentz-covariant parameterization of the matrix elements in terms of a linearly independent basis of tensor structures. This parameterization is crucial for projecting lattice QCD matrix elements onto light cone distributions. For the first time, we determine the corresponding components that project onto the linear combinations of invariant amplitudes, which reduce to the different gluon GPDs in the light cone limit and enable their separation in a lattice QCD calculation for spin-0 and spin- 1 2 hadrons. Hence, this work lays the foundation for the numerical determination of the gluon GPDs from first-principle lattice QCD calculations, directly advancing our understanding of the mass and spin structures and mechanical properties of the nucleon, as well as the physics underlying deeply virtual Compton scattering and deeply virtual meson production in a range of experimental processes. Published by the American Physical Society 2025

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Enabling Efficient Sparse Computations using Linear Algebra Aware Compilers

This project developed the LAPIS compiler framework, built on the Multilevel Intermediate Representation (MLIR), to optimize sparse linear algebra operations and support performance portability across diverse architectures. The main innovation of LAPIS is the Kokkos dialect, which allows for lowering codes from a high productivity language to different architectures in an elegant way. The dialect also allows the conversion of lower-level MLIR code to C++ Kokkos code, facilitating the integration of scientific machine learning (SciML) models into applications. To extend LAPIS for distributed memory architectures, a new partition dialect was created to manage the distribution of sparse tensors and express communication patterns for sparse linear algebra operations. This dialect also supports the distributed execution of operators and includes algorithmic optimizations to minimize communication to improve performance. The project also demonstrates that MLIR can enable effective linear algebra-level optimizations, improving performance on different GPUs for both sparse and dense linear algebra kernels. Key applications of LAPIS include sparse linear algebra and graph kernels, TenSQL, a relational database management solution built on GraphBLAS, and the development of subgraph isomorphism and monomorphism kernels, showcasing performance portability. In summary, the LAPIS framework supports productivity, performance, portability, and distributed memory execution, while also enabling linear algebra-level optimizations that are challenging in traditional programming languages, with successful applications ranging from simple sparse linear algebra to complex graph kernels.

97 MATHEMATICS AND COMPUTING↗

Randomized algorithms for accelerating linear algebraic computations

The project supported the development of new methodologies for performing matrix computations that form key building blocks in modern scientific computing, such as low rank approximation of matrices, and efficient representations of global operators that arise in simulations of physical phenomena.

97 MATHEMATICS AND COMPUTING↗

LCLS-II Cryomodules Production Experience and Lessons Learned Towards LCLS-II-HE Project

LCLS-II is an upgrade project for the linear coherent light source (LCLS) at SLAC. The LCLS-II linac consists of thirty-five 1.3 GHz and two 3.9 GHz superconducting RF (SRF) continuous wave (CW) cryomodules with high quality factor cavities. Cryomodules were produced at Fermilab and at Jefferson Lab in collaboration with SLAC. Fermilab successfully completed the assembly, testing and delivery of seventeen 1.3 GHz and three 3.9 GHz cryomodules. LCLS-II-HE is a planned upgrade project to LCLS-II. The LCLS-II-HE linac will consist of twenty-three 1.3 GHz cryomodules with high gradient and high quality factor cavities. This paper presents LCLS-II-HE cryomodule production plans, emphasizing the improvements done based on the challenges, mitigations, and lessons learned from LCLS-II.

43 PARTICLE ACCELERATORS↗

Verification of an improved equation-free projective integration method for neoclassical plasma-profile evolution in tokamak geometry

A brute-force, long-time gyrokinetic simulation of plasma profile evolution in magnetic fusion devices is not desirable due to large computational resource requirements and a possible accumulation of numerical error. The equation-free projective integration method of Keverekidis et al. [Commun. Math. Sci. 1(4), 715–762 (2003)] is one of the outstanding candidates in projecting micro-scale simulations to a longer timescale. However, its application to tokamak plasma has not been fruitful due to the appearance of spurious transient oscillations in the lifting process, which are present when the kinetic simulations are initialized with a simplified model distribution function and which make the kinetic simulations to deviate from the desired paths. In this work, a kinetically informed lifting algorithm is added to the equation-free projective integration method, which is then verified in the electrostatic gyrokinetic particle-in-cell code XGCa [R. Hager and C. S. Chang, Phys. Plasmas 23, 042503 (2016)] for a neoclassical ion heat transport problem with adiabatic electrons. This new lifting operator is demonstrated to control spurious transients, enabling an over four-times reduction in the overall computing time in the time-evolution of the ion temperature profile in an axisymmetric toroidal plasma. Further reduction in the computing time is found to be limited due to the stability properties of the linear least squares projective integrator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Planar slicing for nonextrusion AM processes

Extrusion-based AM processes, including material extrusion and directed energy deposition, construct objects by continuously depositing a relatively small amount of molten feedstock to a specific location. Nonextrusion AM processes, on the other hand, tend to form each layer by linearly or areally projecting either energy onto a vat of photocurable liquid or a binding agent into a bed of loose powder. Such systems vary significantly from extrusion-based AM processes but do share similarities. The general approach to path planning is the same, but nonextrusion processes typically favor the use of bitmap representations, instead of polygons, for pathing computation and do not use a g-code file to command the printer. In this chapter, a high-level discussion of these differences will be covered.

Macdonald, Eric↗

Latent map Gaussian processes for mixed variable metamodeling

Gaussian processes (GPs) are ubiquitously used in sciences and engineering as metamodels. Standard GPs, however, can only handle numerical or quantitative variables. Here we introduce latent map Gaussian processes (LMGPs) that inherit the attractive properties of GPs and are also applicable to mixed data which have both quantitative and qualitative inputs. The core idea behind LMGPs is to learn a continuous, low-dimensional latent space or manifold which encodes all qualitative inputs. To learn this manifold, we first assign a unique prior vector representation to each combination of qualitative inputs. We then use a low-rank linear map to project these priors on a manifold that characterizes the posterior representations. As the posteriors are quantitative, they can be directly used in any standard correlation function such as the Gaussian or Matern. Hence, the optimal map and the corresponding manifold, along with other hyperparameters of the correlation function, can be systematically learned via maximum likelihood estimation. Through a wide range of analytic and real-world examples, we demonstrate the advantages of LMGPs over state-of-the-art methods in terms of accuracy and versatility. In particular, we show that LMGPs can handle variable-length inputs, have an explainable neural network interpretation, and provide insights into how qualitative inputs affect the response or interact with each other. We also employ LMGPs in Bayesian optimization and illustrate that they can discover optimal compound compositions more efficiently than conventional methods that convert compositions to qualitative variables via manual featurization.

42 ENGINEERING↗