Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear systems”

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,171 records · Page 65

Design and Optimization of a Fast-Switching Electromagnetic Valve Actuator for Off-highway Vehicle Hydraulic Systems

Electrification of off-highway vehicles offers improved efficiency and reduction in greenhouse gas emissions. However, the extremely low power densities of electric machines relative to the hydraulics used in conventional off-highway vehicles makes electrification challenging. Although the recently proposed Hybrid Hydraulic Electric Architecture (HHEA) can overcome this limitation, it requires frequent switching between the pressure rails using switching valves resulting in losses. Fast switching minimizes the switching losses in these valves, which otherwise significantly reduce the overall efficiency of the system. The primary contribution of this paper is the design and optimization of a linear electromagnetic actuator for a fast-switching valve. This paper identifies the actuator requirements to enable fast switching and presents a design and optimization technique to realize actuators that meet the design objectives. In addition, a multi-physics approach is proposed to model the performance of the complete switching valve by considering the electromagnetic, fluidic, and mechanical physics. An optimal actuator design is obtained using the developed optimization framework, and the multi-physics model demonstrates that the selected design achieves a 30% reduction in throttling-related switching losses by switching 70 % faster than comparable commercial valve designs.

Linear actuators, off-highway vehicle electrificat↗

Analog In-Memory Computing for the Synthetic Aperture Radar Polar Format Algorithm

As the utility of synthetic aperture radar (SAR) systems increases in autonomous vehicles, satellites, and other power- and space-constrained edge applications, there is a growing need for processors that can form SAR images at low power. In recent years, analog in-memory compute (AIMC) has shown immense promise for accelerating neural networks and other matrix-vector multiplication (MVM) heavy workloads at the edge. Here, in this work, we examine how the polar format algorithm (PFA), a popular SAR image formation algorithm, can be mapped to these AIMC systems. The PFA maps readily onto analog MVMs because it primarily consists of two linear operations: interpolation of frequency-domain data to a Cartesian grid, followed by a 2-D Fourier transform. This work presents two approaches to map the interpolation operation onto MVMs in analog hardware: a chirp transform and a modified form of sinc interpolation. These mappings introduce algorithmic errors, and their effect on the quality of SAR image formation is examined, both quantitatively and qualitatively. In addition, the impact of errors introduced by the analog hardware is explored to determine which approach is optimal under varying assumptions about the underlying analog memory devices and circuits.

Analog computing↗

A Trajectory-Driven Algorithm for Differentiating SRB Measures on Unstable Manifolds

Sinai-Ruelle-Bowen (SRB) measures are limiting stationary distributions describing the statistical behavior of chaotic dynamical systems. Directional derivatives of SRB measure densities conditioned on unstable manifolds are critical in the sensitivity analysis of hyperbolic chaos. These derivatives, known as the SRB density gradients, are by-products of the regularization of Lebesgue integrals appearing in the original linear response expression. In this paper, we propose a novel trajectory- driven algorithm for computing the SRB density gradient defined for systems with high-dimensional unstable manifolds. We apply the concept of measure preservation together with the chain rule on smooth manifolds. Due to the recursive one-step nature of our derivations, the proposed procedure is memory-efficient and can be naturally integrated with existing Monte Carlo schemes widely used in computational chaotic dynamics. Here, we numerically show the exponential convergence of our scheme, analyze the computational cost, and present its use in the context of Monte Carlo integration.

97 MATHEMATICS AND COMPUTING↗

CMLM (Co-Optimized Machine-Learned Manifolds) [SWR-23-41]

Co-optimized Machine-Learned Manifolds (CMLM) is a data-driven approach for developing reduced-order manifold models for high-dimensional chemically reacting systems. It involves a specially designed neural network, the training of which simultaneously optimizes linear combinations of species that define the manifold, nonlinear mapping to outputs of interest such as reaction rates, and (optionally) subfilter closure for large eddy simulation. This software package provides an implementation of the CMLM approach in Python using the PyTorch machine learning library. A few example cases are included, showing how the tool can be applied to different types of data from 0D and 1D reacting simulations performed using Cantera. The neural networks can be saved in a format that is readable by the Pele suite of combustion solvers for use in reacting computational fluid dynamics simulations. This software repository contains several python scripts to perform various tasks associated with the Co-optimized Machine Learned Manifolds (CMLM) model, which is described in Perry, Henry de Frahan, and Yellapantula, CNF, 2022 (https://doi.org/10.1016/j.combustflame.2022.112286). This includes not only the code that defines the CMLM model, but also scripts to generate suitable training data, scripts to pre-process the data, scripts to train the CMLM model, and scripts to plot the output, as well as various other helper files. The scripts depend on several commonly used python libraries for data analysis and chemical reaction computations. The trained models that result from this tool are designed to work with the an interface being implemented in the Pele suite of reacting flow solvers (https://github.com/AMReX-Combustion).

Perry, Bruce↗

A Statistical Interpolation Code for Ocean Analysis and Forecasting

Abstract We present a data assimilation package for use with ocean circulation models in analysis, forecasting, and system evaluation applications. The basic functionality of the package is centered on a multivariate linear statistical estimation for a given predicted/background ocean state, observations, and error statistics. Novel features of the package include support for multiple covariance models, and the solution of the least squares normal equations either using the covariance matrix or its inverse—the information matrix. The main focus of this paper, however, is on the solution of the analysis equations using the information matrix, which offers several advantages for solving large problems efficiently. Details of the parameterization of the inverse covariance using Markov random fields are provided and its relationship to finite-difference discretizations of diffusion equations are pointed out. The package can assimilate a variety of observation types from both remote sensing and in situ platforms. The performance of the data assimilation methodology implemented in the package is demonstrated with a yearlong global ocean hindcast with a 1/4° ocean model. The code is implemented in modern Fortran, supports distributed memory, shared memory, multicore architectures, and uses climate and forecasts compliant Network Common Data Form for input/output. The package is freely available with an open source license from www.tendral.com/tsis/ .

Srinivasan, Ashwanth↗

Massively parallel and universal approximation of nonlinear functions using diffractive processors

Nonlinear computation is essential for a wide range of information processing tasks, yet implementing nonlinear functions using optical systems remains a challenge due to the weak and power-intensive nature of optical nonlinearities. Overcoming this limitation without relying on nonlinear optical materials could unlock unprecedented opportunities for ultrafast and parallel optical computing systems. Here, we demonstrate that large-scale nonlinear computation can be performed using linear optics through optimized diffractive processors composed of passive phase-only surfaces. In this framework, the input variables of nonlinear functions are encoded into the phase of an optical wavefront—e.g., via a spatial light modulator (SLM)—and transformed by an optimized diffractive structure with spatially varying point-spread functions to yield output intensities that approximate a large set of unique nonlinear functions–all in parallel. We provide proof establishing that this architecture serves as a universal function approximator for an arbitrary set of bandlimited nonlinear functions, also covering wavelength-multiplexed nonlinear functions as well as multi-variate and complex-valued functions that are all-optically cascadable. Our analysis also indicates the successful approximation of typical nonlinear activation functions commonly used in neural networks, including the sigmoid, tanh, ReLU (rectified linear unit), and softplus. We numerically demonstrate the parallel computation of one million distinct nonlinear functions, accurately executed at wavelength-scale spatial density at the output of a diffractive optical processor. Furthermore, we experimentally validated this framework using in situ optical learning and approximated 35 unique nonlinear functions in a single shot using a compact setup consisting of an SLM and an image sensor. These results establish diffractive optical processors as a scalable platform for massively parallel universal nonlinear function approximation, paving the way for new capabilities in analog optical computing based on linear materials.

Rahman, Md Sadman Sakib [University of California,↗

Closed-Form Solutions for the Equations of Motion of the Heavy Symmetrical Top with One Point Fixed

The Equations of Motion (EOM) for the Heavy Symmetrical Top with One Point Fixed are highly non-linear. The literature describes the numerical methods that are used to resolve this Classical system including modern tools i.e. the Runge-Kutta Fourth Order method. It is more difficult to derivate closed-form solutions for the EOM and as mentioned in the literature it is not always possible to find the close-form solution for all the EOM. Fortunately, there are a few examples available that will serve as a guide to move further on in this topic. It is the purpose of this paper to find a methodology that will produce the solutions for a given subset of EOM that fulfill certain requisites. The report is organized as follows: it starts with a very short summary of the literature available on this topic and quickly follows into the derivation of the EOM using the Euler-Lagrange method. The Routhian will be used to reduce the size of the expression. It continues with the formulation of the Classical cubic function (f(u)) through a novel process. The roots of f(u) are of the outmost importance to be able to find the EOM closed-form solution, and once the final roots are selected the general method that will produce the closed-form solutions is presented. Two sets of examples are included to show the validity of the process and comparisons of the results from the closed-form solutions vs. the numerical results for these examples are shown.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

LDRD23-0730: Invoking Multilayer Networks to Develop a Paradigm for Security Science—Summary Report

Current approaches to securing high consequence facilities (HCF) and critical assets are linear and static and therefore struggle to adapt to emerging threats (e.g., unmanned aerial systems) and changing environmental conditions (e.g., decreasing operational control). The pace of change in technological, organizational, societal, and political dynamics necessitates a move toward codifying underlying scientific principles to better characterize the rich interactions observed between HCF security technology, infrastructure, digital assets, and human or organizational components. The promising results of Laboratory Directed Research and Development (LDRD) 20-0373—“Developing a Resilient, Adaptive, and Systematic Paradigm for Security Analysis”—suggest that when compared to traditional security analysis, invoking multilayer network (MLN) modeling for HCF security system components captures unexpected failure cases and unanticipated interactions.

97 MATHEMATICS AND COMPUTING↗

Interstellar Extinction, Polarization, and Grain Alignment in the Sh 2-185 (IC 59 and IC 63) Region

Optical and infrared continuum polarization from the interstellar medium is driven by radiative processes aligning the grains with the magnetic field. While a quantitative, predictive theory of radiative alignment torques (RATs) exists and has been extensively tested, several parameters of the theory remain to be fully constrained. In a recent paper, Medan & Andersson showed that the polarization efficiency (and therefore grain alignment efficiency) at different locations in the wall of the Local Bubble (LB) could be modeled as proportional to the integrated light intensity from the surrounding stars and OB associations. Here we probe that relationship at high radiation field intensities by studying the extinction and polarization in the two reflection nebulae IC 59 and IC 63 in the Sh 2-185 H ii region, illuminated by the B0 IV star γ Cassiopeia. We combine archival visual polarimetry with new seven-band photometry in the Vilnius system, to derive the polarization efficiency from the material. We find that the same linear relationship seen in the LB wall also applies to the Sh 2-185 region, strengthening the conclusion from the earlier study.

47 OTHER INSTRUMENTATION↗

Design and evaluation of a low-cost sensor node for near-background methane measurement

Abstract. We developed a low-cost methane sensing node incorporating two metal oxide (MOx) sensors (Figaro Engineering TGS2611-E00 and TGS2600), humidity and temperature sensing, data storage, and telemetry. We deployed the prototype sensor alongside a reference methane analyzer at two sites: one outdoors and one indoors. We collected data at each site for several months across a range of environmental conditions (particularly temperature and humidity) and methane levels. We explored calibration models to investigate the performance of our system and its suitability for methane background monitoring and enhancement detection, first selecting a linear regression to fit a sensor baseline response and then fitting methane response by the sensor deviation from baseline. We achieved moderate accuracy in a 2 to 10 ppm methane range compared to data from the reference analyzer (RMSE < 0.6 ppm), but we found that the sensor response varied over time, possibly as a result of changes in non-targeted gas concentrations. We suggest that this cross sensitivity may be responsible for mixed results in previous studies. We discuss the implications of our results for the use of these and similar inexpensive MOx sensors for methane monitoring in the 2 to 10 ppm range.

Furuta, Daniel↗

Tunable plasmonic color device and method of making the same

A plasmonic system is disclosed. The system includes at least one polarizer that is configured to provide at least one linearly polarized broadband light beam, an anisotropic plasmonic metasurface (APM) assembly having a plurality of nanoantennae each having a predetermined orientation with respect to a global axis representing encoded digital data, the APM assembly configured to receive the at least one linearly polarized broadband light beam and by applying localized surface plasmon resonance reflect light with selectable wavelengths associated with the predetermined orientations of the nanoantennae, and at least one analyzer that is configured to receive the reflected light with selectable wavelength, wherein the relative angles between each of the at least one analyzers and each of the at least one polarizers are selectable with respect to the global axis, thereby allowing decoding of the digital data.

Kildishev, Alexander V.↗

Risk-averse optimization for resilience enhancement of complex engineering systems under uncertainties

With the growth of complexity and extent, large scale interconnected network systems, e.g., transportation networks or infrastructure networks, become more vulnerable to external disturbances. Hence, managing potential disruptive events during the design, operating, and recovery phase of an engineered system and therefore improving the system’s resilience is an important yet challenging task. Here, to ensure system resilience after the occurrence of failure events, this study proposes a mixed-integer linear programming (MILP) based restoration framework using heterogeneous dispatchable agents. The scenario-based stochastic optimization (SO) technique is adopted to deal with the inherent uncertainties imposed on the recovery process from nature. Moreover, different from conventional SO using deterministic equivalent formulations, the CVaR risk measure is implemented for this study because of the temporal sparsity of the decision making in applications such as the recovery from extreme events. The resulting restoration framework involves a large-scale MILP problem and thus an adequate decomposition technique i.e. modified Lagrangian dual decomposition, is also employed to achieve tractable computational complexity. Case study results based on the IEEE 37-bus test feeder demonstrate the benefits of using the proposed framework for resilience improvement as well as the advantages of adopting SO formulations.

42 ENGINEERING↗

Multimodal (Non)Linear Optical Nanoimaging and Nanospectroscopy

This article highlights recent advances in linear and nonlinear spectral nanoimaging. The described developments are motivated by the need to characterize molecular and material systems noninvasively with nanometer spatial and femtosecond temporal resolution. Indeed, the ability to image and chemically characterize heterogeneous interfaces with joint nano–femto resolution is a prerequisite to advancing our fundamental understanding of processes as diverse as heterogeneous catalysis, microbial communication, and energy flow in pristine/defect-containing low-dimensional quantum materials, to name a few. We describe pioneering work and recent demonstrations of (non)linear optical nanoimaging and nanospectroscopy, with an emphasis on high spatial resolution measurements conducted under ambient laboratory conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Impact of high-temperature annealing on hafnia-silica composite coatings deposited via ion beam sputtering for high-peak power 1064 nm lasers

The maximum power handling fluence of high-peak and average power laser systems is often limited by the laser damage of the coatings on optical components. Furthermore, these multilayer dielectric coatings are limited in their maximum power handling due to laser-damage-prone defects in the lower optical bandgap, higher optical index material. Some of these defects can be mitigated by thermal annealing to high temperatures, which can greatly reduce the linear absorbative precursors. Typically, hafnia and silica are the materials of choice for high-peak and average power laser systems in the ultraviolet through infrared spectral range; however, hafnia crystalizes readily when annealed at high temperatures. In this study, we prepare composite HfO 2 -SiO 2 coatings by co-sputtering hafnia and silica in an ion beam sputtering system and compare them to pure hafnia-based coatings. We demonstrate that crystallinity in hafnia can be completely suppressed when it is mixed with silica, such as the composite coatings in this study. High reflectors were fabricated and annealed, demonstrating that the multilayer dielectric stacks can survive high-temperature annealing and exhibit an excellent linear absorption of 0.2 +/- 1 ppm at 1064 nm. Short- and long-pulse laser damage was explored, demonstrating the complex relationship between linear absorption and the non-linear absorption which drives pulsed laser damage. These results provide an excellent route to the creation of very low linear absorption optical coatings, which also utilize low scattering materials that are best suited for high-peak and average power applications.

Harthcock, Colin [Lawrence Livermore National Labo↗

Towards Ab-Initio Simulations of Crystalline Defects at the Exascale Using Spectral Quadrature Density Functional Theory

Defects in crystalline solids play a crucial role in determining properties of materials at the nano, meso- and macroscales, such as the coalescence of vacancies at the nanoscale to form voids and prismatic dislocation loops or diffusion and segregation of solutes to nucleate precipitates, phase transitions in magnetic materials via disorder and doping. First principles Density Functional Theory (DFT) simulations can provide a detailed understanding of these phenomena. However, the number of atoms needed to correctly simulate these systems is often beyond the reach of many widely used DFT codes. The aim of this article is to discuss recent advances in first principles modeling of crystal defects using the spectral quadrature method. The spectral quadrature method is linear scaling with respect to the number of atoms, permits spatial coarse-graining, and is capable of simulating non-periodic systems embedded in a bulk environment, which allows the application of appropriate boundary conditions for simulations of crystalline defects. In this article, we discuss the state-of-the-art in ab-initio modeling of large metallic systems of the order of several thousand atoms that are suitable for utilizing exascale computing resourses.

36 MATERIALS SCIENCE↗

Holographic Kolmogorov-Sinai entropy and the quantum Lyapunov spectrum

In classical chaotic systems the entropy, averaged over initial phase space distributions, follows a universal behavior. While approaching thermal equilibrium it passes through a stage where it grows linearly, while the growth rate, the Kolmogorov-Sinai entropy (rate), is given by the sum over all positive Lyapunov exponents. A natural question is whether a similar relation is valid for quantum systems. We argue that the Maldacena-Shenker-Stanford bound on quantum Lyapunov exponents implies that the upper bound on the growth rate of the entropy, averaged over states in Hilbert space that evolve towards a thermal state with temperature T, should be given by πT times the thermal state’s von Neumann entropy. Strongly coupled, large N theories with black hole duals should saturate the bound. To test this we study a large number of isotropization processes of random, spatially homogeneous, far from equilibrium initial states in large N, $\mathcal{N}$ = 4 Super Yang Mills theory at strong coupling and compute the ensemble averaged growth rate of the dual black hole’s apparent horizon area. We find both an analogous behavior as in classical chaotic systems and numerical evidence that the conjectured bound on averaged entropy growth is saturated granted that the Lyapunov exponents are degenerate and given by λ i = ±2πT. This fits to the behavior of classical systems with plus/minus symmetric Lyapunov spectra, a symmetry which implies the validity of Liouville’s theorem.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Midfidelity Model Verification for a Point-Absorbing Wave Energy Converter with Linear Power Takeoff: Preprint

In the preliminary design stage of a wave energy converter (WEC), fast and reliable simulation tools are required. High-fidelity numerical models are usually employed to study the wave-structure interaction, yet the computational cost is demanding. Instead, mid-fidelity models provide simulations in the order of real time. In this study, Uppsala University WEC operates in a relatively mild sea state and is modeled using WEC-Sim. The model is verified based on OpenFOAM simulations. To analyze the ability of mid-fidelity model to capture the WEC dynamics, the system is investigated separately in 1, 2 and 3 DoF. The contribution of viscous phenomena is examined. Moreover, linear and weakly non-linear solutions provided by WEC-Sim are studied. The results obtained indicate that the viscous effects in heave and surge motion can be neglected but not for the pitch. The weakly non-linear WEC-Sim solution successfully agrees with the CFD, whereas the linear solution could drive to misleading results.

CFD↗

Midfidelity Model Verification for a Point-Absorbing Wave Energy Converter with Linear Power Take-Off

In the preliminary design stage of a wave energy converter (WEC), fast and reliable simulation tools are required. High-fidelity numerical models are usually employed to study the wave-structure interaction, yet the computational cost is demanding. Instead, mid-fidelity models provide simulations in the order of real time. In this study, Uppsala University WEC operates in a relatively mild sea state and is modeled using WEC-Sim. The model is verified based on OpenFOAM simulations. To analyze the ability of mid-fidelity model to capture the WEC dynamics, the system is investigated separately in 1, 2 and 3 DoF. The contribution of viscous phenomena is examined. Moreover, linear and weakly non-linear solutions provided by WEC-Sim are studied. The results obtained indicate that the viscous effects in heave and surge motion can be neglected but not for the pitch. The weakly non-linear WEC-Sim solution successfully agrees with the CFD, whereas the linear solution could drive to misleading results.

CFD↗