Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linearization”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

A digital twin based on OpenFAST linearizations for real-time load and fatigue estimation of land-based turbines

Monitoring a wind turbine requires intensive instrumentation, which would be too cost-prohibitive to deploy to an entire wind plant. This work presents a technique that uses readily available measurements to estimate signals that would otherwise require additional instrumentation. This study presents a digital twin concept with a focus on estimating wind speed, thrust, torque, tower-top position, and loads in the tower using supervisory control and data acquisition (SCADA) measurements. The model combines a linear state-space model obtained using OpenFAST linearizations, a wind speed estimator, and a Kalman filter algorithm that integrates measurements with the state model to perform state estimations. The measurements are: top acceleration, generator torque, pitch, and rotational speed. The article extends previous work that derived the linear state-space model using a different method. The new implementation, based on OpenFAST linearization capability, allows for a systematic extension of the method to more states, inputs, outputs, and to the offshore environment. Results from the two methods are compared, and the validation is made with additional measurements using the GE 1.5-MW turbine located at the National Renewable Energy Laboratory test site. Real-time damage equivalent loads of the tower bottom moment are estimated with an average accuracy of approximately 10%. Overall, the results from this proof of concept are encouraging, and further application of the model will be considered.

17 WIND ENERGY↗
DOE OSTI2015

SIMULTANEOUS LINEAR AND CIRCULAR OPTICAL POLARIMETRY OF ASTEROID (4) VESTA

From a single 3.8 hr observation of the asteroid (4) Vesta at 13.°7 phase angle with the POlarimeter at Lick for Inclination Studies of Hot jupiters 2 (POLISH2) at the Lick Observatory Shane 3 m telescope, we confirm rotational modulation of linear polarization in the B and V bands. We measure the peak-to-peak modulation in the degree of linear polarization to be ΔP = (294 ± 35) × 10{sup −6} (ppm) and time-averaged ΔP/P = 0.0575 ± 0.0069. After rotating the plane of linear polarization to the scattering plane, asteroidal rotational modulation is detected with 12σ confidence and observed solely in Stokes Q/I. POLISH2 simultaneously measures Stokes I, Q, U (linear polarization), and V (circular polarization), but we detect no significant circular polarization with a 1σ upper limit of 78 ppm in the B band. Circular polarization is expected to arise from multiple scattering of sunlight by rough surfaces, and it has previously been detected in nearly all other classes of solar system bodies except for asteroids. Subsequent observations may be compared with surface albedo maps from the Dawn Mission, which may allow the identification of compositional variation across the asteroidal surface. These results demonstrate the high accuracy achieved by POLISH2 at the Lick 3 m telescope, which is designed to directly detect scattered light from spatially unresolvable exoplanets.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Effects of linear trends on estimation of noise in GNSS position time-series

A thorough understanding of time-dependent noise in Global Navigation Satellite System (GNSS) position time-series is necessary for computing uncertainties in any signals found in the data. However, estimation of time-correlated noise is a challenging task and is complicated by the difficulty in separating noise from signal, the features of greatest interest in the time-series. In this study, we investigate how linear trends affect the estimation of noise in daily GNSS position time-series. We use synthetic time-series to study the relationship between linear trends and estimates of time-correlated noise for the six most commonly cited noise models. We find that the effects of added linear trends, or conversely de-trending, vary depending on the noise model. The commonly adopted model of random walk (RW), flicker noise (FN) and white noise (WN) is the most severely affected by de-trending, with estimates of low-amplitude RW most severely biased. FN plus WN is least affected by adding or removing trends. Non-integer power-law noise estimates are also less affected by de-trending, but are very sensitive to the addition of trend when the spectral index is less than one. We derive an analytical relationship between linear trends and the estimated RW variance for the special case of pure RW noise. Finally, overall, we find that to ascertain the correct noise model for GNSS position time-series and to estimate the correct noise parameters, it is important to have independent constraints on the actual trends in the data.

58 GEOSCIENCES↗
DOE OSTI2023

A non-linear solution to the S 8 tension – II. Analysis of DES Year 3 cosmic shear

ABSTRACT Weak galaxy lensing surveys have consistently reported low values of the S8 parameter compared to the Planck lambda cold dark matter (ΛCDM) cosmology. Amon & Efstathiou used KiDS-1000 cosmic shear measurements to propose that this tension can be reconciled if the matter fluctuation spectrum is suppressed more strongly on non-linear scales than assumed in state-of-the-art hydrodynamical simulations. In this paper, we investigate cosmic shear data from the Dark Energy Survey (DES) Year 3. The non-linear suppression of the matter power spectrum required to resolve the S8 tension between DES and the Planck ΛCDM model is not as strong as inferred using KiDS data, but is still more extreme than predictions from recent numerical simulations. An alternative possibility is that non-standard dark matter contributes to the required suppression. We investigate the redshift and scale dependence of the suppression of the matter power spectrum. If our proposed explanation of the S8 tension is correct, the required suppression must extend into the mildly non-linear regime to wavenumbers $k\sim 0.2 \, h\, {\rm Mpc}^{-1}$. In addition, all measures of S8 using linear scales should agree with the Planck ΛCDM cosmology, an expectation that will be testable to high precision in the near future.

Preston, Calvin (ORCID:0000000281136916)↗
DOE OSTI2015

Non-linear violent disc instability with high Toomre's Q in high-redshift clumpy disc galaxies

We utilize zoom-in cosmological simulations to study the nature of violent disc instability in clumpy galaxies at high redshift, z = 1–5. Our simulated galaxies are not in the ideal state assumed in Toomre instability, of linear fluctuations in an isolated, uniform, rotating disc. There, instability is characterized by a Q parameter below unity, and lower when the disc is thick. Instead, the high-redshift discs are highly perturbed. Over long periods they consist of non-linear perturbations, compact massive clumps and extended structures, with new clumps forming in interclump regions. This is while the galaxy is subject to frequent external perturbances. We compute the local, two-component Q parameter for gas and stars, smoothed on a ~1 kpc scale to capture clumps of 10 8-9 M ⊙ . The Q < 1 regions are confined to collapsed clumps due to the high surface density there, while the interclump regions show Q significantly higher than unity. Tracing the clumps back to their relatively smooth Lagrangian patches, we find that Q prior to clump formation typically ranges from unity to a few. This is unlike the expectations from standard Toomre instability. We discuss possible mechanisms for high- Q clump formation, e.g. rapid turbulence decay leading to small clumps that grow by mergers, non-axisymmetric instability, or clump formation induced by non-linear perturbations in the disc. Alternatively, the high- Q non-linear VDI may be stimulated by the external perturbations such as mergers and counter-rotating streams. In conclusion, the high Q may represent excessive compressive modes of turbulence, possibly induced by tidal interactions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Variational quantum and neural quantum states algorithms for the linear complementarity problem

Variational quantum algorithms (VQAs) are promising hybrid quantum-classical methods designed to leverage the computational advantages of quantum computing while mitigating the limitations of current noisy intermediate-scale quantum (NISQ) hardware. Although VQAs have been demonstrated as proofs of concept, their practical utility in solving real-world problems—and whether quantum-inspired classical algorithms can match their performance—remains an open question. We present a novel application of the variational quantum linear solver (VQLS) and its classical neural quantum states-based counterpart, the variational neural linear solver (VNLS), as key components within a minimum map Newton solver for a complementarity-based rigid-body contact model. We demonstrate using the VNLS that our solver accurately simulates the dynamics of rigid spherical bodies during collision events. These results suggest that quantum and quantum-inspired linear algebra algorithms can serve as viable alternatives to standard linear algebra solvers for modelling certain physical systems.

neural quantum states↗
DOE OSTI2022

Topological alloy engineering and locally linearized gap dependence on concentration

Alloy engineering is a well-established approach to tune various materials’ properties, but its application to topological alloys remains rudimentary. Of special interest is the band gap, the most defining property of topological materials; however, the concentration dependence of energy gaps in topological alloys remains unknown. Here we systematically investigate the band gap evolution of a topological alloy as a function of alloy concentration, using KZnSb 1-x Bi x as a prototype, based on first-principles calculations. In contrast to the well-established smooth bowing curve for a trivial gap in semiconductor alloys, we found that the topological gap evolves generally with a complex fragmented pattern due to topological phase transitions, and most strikingly a linear dependence on concentration locally in each distinct phase. Such gap linearization is fundamentally rooted in the linear dependence on alloy concentration of spin-orbit coupling (SOC) that predominantly determines a topological gap. Furthermore, we demonstrate topological alloy engineering as a general approach to tune the topological order by modulating the band edge composition and degeneracy through the alloying-induced interplay of SOC and atomic orbital on-site energy, while the linear gap dependence on alloy concentration remains independent of the degree of topological order.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Controlled skyrmion ratchet in linear protrusion defects

Using atomistic simulations, we investigate the dynamical behavior of a single skyrmion interacting with an asymmetric linear protrusion array under external ac driving. The linear protrusion array is placed in the xy plane and is composed of magnetic walls near the edge of the nanotrack that have angled linear protrusions into the nanotrack. The structure forms a hard axis for –x direction motion and an easy axis for +x direction motion. When the ac drive is applied along the x direction, the skyrmion moves along the hard direction of the substrate asymmetry in three phases: a pinned phase with localized skyrmion orbits, a constant velocity phase where the orbits become delocalized, and a reentrant pinned phase with larger localized orbits. We measure the dependence of the skyrmion velocity on the frequency and amplitude of the ac drive. All three phases appear for all frequency values simulated here, and in the constant velocity phase the skyrmion velocity depends only on the frequency and not on the amplitude of the ac drive. When ac driving is applied in the y direction, the skyrmion moves along the easy direction of the substrate asymmetry and exhibits the same three phases as for x-direction driving along with a fourth phase which, at high driving frequencies, consists of a series of constant velocity phases, each with different average skyrmion velocities. For low frequencies, the constant velocity phase is lost and the skyrmion speed increases linearly with increasing ac drive amplitude due to a Magnus boost effect. Finally, our findings suggest new ways to create reliable data transport for spintronic devices using skyrmions as information carriers, where the skyrmion direction and speed can be controlled by varying only the ac drive amplitude and frequency.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Temperature scaling behavior of the linear magnetoresistance observed in high-temperature superconductors

An analytical model invoking variations in the charge-carrier density is used to generate magnetoresistance curves that are almost indistinguishable from those produced by sophisticated numerical models. This demonstrates that, though disorder is pivotal in causing linear magnetoresistance, the form of the magnetoresistance thus generated is insensitive to details of the disorder. Additionally, taken in conjunction with the temperature (T) dependence of the zero-field resistivity, realistic levels of disorder are shown to be sufficient to explain the linear magnetoresistance and field-T resistance scaling observed in high-temperature pnictide and cuprate superconductors. Hence, though the T-linear zero-field resistance is a definite signature of the “strange metal” state of high-temperature superconductors, their linear magnetoresistance and its scaling is unlikely to be so.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Adaptive linear predictor FIR filter based on the Cyclone V FPGA with HPS to reduce narrow band RFI in AERA radio detection of cosmic rays

We present the new approach to a filtering of radio frequency interferences (RFI) in the Auger Engineering Radio Array (AERA) which study the electromagnetic part of the Extensive Air Showers. The radio stations can observe radio signals caused by coherent emissions due to geomagnetic radiation and charge excess processes. AERA observes frequency band from 30 to 80 MHz. This range is highly contaminated by human-made RFI. In order to improve the signal to noise ratio RFI filters are used in AERA to suppress this contamination. The first kind of filter used by AERA was the Median one, based on the Fast Fourier Transform (FFT) technique. The second one, which is currently in use, is the infinite impulse response (IIR) notch filter. The proposed new filter is a finite impulse response (FIR) filter based on a linear prediction (LP). A periodic contamination hidden in a registered signal (digitized in the ADC) can be extracted and next subtracted to make signal cleaner. The FIR filter requires a calculation of n=32, 64 or even 128 coefficients (dependent on a required speed or accuracy) by solving of n linear equations with coefficients built from the covariance Toeplitz matrix. This matrix can be solved by the Levinson recursion, which is much faster than the Gauss procedure. The filter has been already tested in the real AERA radio stations on Argentinean pampas with a very successful results. The linear equations were solved either in the virtual soft-core NIOSR processor (implemented in the FPGA chip as a net of logic elements) or in the external Voipac PXA270M ARM processor. The NIOS processor is relatively slow (50 MHz internal clock), calculations performed in an external processor consume a significant amount of time for data exchange between the FPGA and the processor. Test showed a very good efficiency of the RFI suppression for stationary (long-term) contaminations. However, we observed a short-time contaminations, which could not be suppressed either by the IIR-notch filter or by the FIR filter based on the linear predictions. For the LP FIR filter the refreshment time of the filter coefficients was to long and filter did not keep up with the changes of a contamination structure, mainly due to a long calculation time in a slow processors. We propose to use the Cyclone V SE chip with embedded micro-controller operating with 925 MHz internal clock to significantly reduce a refreshment time of the FIR coefficients. The lab results are promising. (authors)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Studying the Peak Force of Magnetic Linear Motion Devices Using Dimensional Neutrality

The merits of different sizing equations for linear motion magnetic force devices are assessed. A peak force scaling analysis for a linear magnetic coupling with respect to the volumetric force density, mass force density, magnetic shear stress and energy density is compared. By utilizing a 3-D analytic based model of a Halbach cylinder magnetic coupling it is shown that the maximizing of volumetric and mass force density is not an appropriate sizing metric for a linear magnetic actuator. Only energy density and magnetic shear stress are dimensionally neutral sizing metrics for linear actuators.

16 TIDAL AND WAVE POWER↗
DOE OSTI2024

Linear Solvers for Collector Systems of Generalized Large-scale Inverter-Based Resources

Collector systems for inverter-based resources (IBRs) are typically represented by equivalent circuits for electromagnetic transient (EMT) simulations. Recent studies have revealed that modeling a detailed collector system is essential to accurately represent the behavior of IBRs, especially when dealing with partial tripping during external disturbances. However, there are several challenges in simulating a detailed EMT model of a collector system due to the time required to simulate such systems. Thus, this paper investigates the modeling of a detailed collector system, taking into account its configuration and components as defined in IEEE standard 2800. The configurations include the collector systems of generalized large-scale IBR plants. The components include the main IBR transformer, collector bus, and feeders with lines and/or cables. The EMT model of the collector system is represented by differential algebraic equations (DAEs) that are discretized to form linear equations that are solved using linear solvers. In this paper, linear solvers are proposed based on the Schur complement method, which are utilized for simulation of the EMT model of collector systems of generalized large-scale IBRs to accelerate simulation speed while maintaining the accuracy of the results. The proposed solvers are verified by comparing the performance to that of linear solvers provided in MATLAB.

Choi, Jongchan↗
DOE OSTI2021

Learning Optimal Power Flow Solutions using Linearized Models in Power Distribution Systems

Solving nonlinear optimal power flow (OPF) problem is computationally expensive, and poses scalability challenges for power distribution networks. An alternative to solving the original nonlinear OPF is the linear approximated OPF models. Although, these linear approximated OPF models are fast, the resulting solutions may result in significant optimality gap. Lately, the application of machine learning (ML) methods in successfully solving the nonlinear OPF has been reported. These methods learn and estimate the nonlinear control policies using a purely data-driven approach. In this paper, we propose an approach to complements the ML based approach to solving OPF using solutions from known linearized OPF model. Specifically, we use supervised learning to map the solutions of linear OPF to nonlinear control variables. Unlike, the traditional ML based methods for OPF that approximate the full distribution feeder model using function approximation, our approach uses a two-node approximation of radial networks. The proposed approach is validated using IEEE 123 bus test system for OPF solutions obtained using the nonlinear OPF models.

optimal power flow, power distribution systems, su↗
DOE OSTI2024

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↗
DOE OSTI2024

Complex-valued universal linear transformations and image encryption using spatially incoherent diffractive networks

As an optical processor, a diffractive deep neural network (D2NN) utilizes engineered diffractive surfaces designed through machine learning to perform all-optical information processing, completing its tasks at the speed of light propagation through thin optical layers. With sufficient degrees of freedom, D2NNs can perform arbitrary complex-valued linear transformations using spatially coherent light. Similarly, D2NNs can also perform arbitrary linear intensity transformations with spatially incoherent illumination; however, under spatially incoherent light, these transformations are nonnegative, acting on diffraction-limited optical intensity patterns at the input field of view. Here, we expand the use of spatially incoherent D2NNs to complex-valued information processing for executing arbitrary complex-valued linear transformations using spatially incoherent light. Through simulations, we show that as the number of optimized diffractive features increases beyond a threshold dictated by the multiplication of the input and output space-bandwidth products, a spatially incoherent diffractive visual processor can approximate any complex-valued linear transformation and be used for all-optical image encryption using incoherent illumination. The findings are important for the all-optical processing of information under natural light using various forms of diffractive surface-based optical processors.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Investigation of Microwave Transducers for Linearity Dependence and Applications in Quantum Networking

Quantum devices have the potential to revolutionize applications in computing, communications, and sensing; however, current state-of-art resources must operate at extremely low temperatures, making the routing of microwave control and readout signals challenging to scale. Interest in microwave photonic solutions to this problem has grown in recent years, in which control signals are delivered to the cold stage via optical fiber, where they are converted to electrical signals through photodetection. Overall link performance depends strongly on the characteristics of the photodiode, yet detailed measurements of many detector properties remain lacking at cold temperatures. In this work, we examine and compare the performance of a modified uni-traveling carrier photodiode (MUTC-PD) at both room (300 K) and liquid nitrogen (80 K) temperatures, focusing in particular on responsivity, bandwidth, and linearity. In line with previous work, we find a sharp reduction in responsivity at 1550 nm as temperature decreases, while RF bandwidth remains steady. Interestingly, our linearity tests reveal that the RF output saturates more quickly at 80 K, suggesting reduced linearity with lower temperature, the cause of which is still under investigation. Our results should help contribute to the understanding and future design of highly linear cryogenic quantum links.

Bolton, Summer↗
DOE OSTI2014

TH-E-BRE-01: A 3D Solver of Linear Boltzmann Transport Equation Based On a New Angular Discretization Method with Positivity for Photon Dose Calculation Benchmarked with Geant4

Purpose: The Linear Boltzmann Transport Equation (LBTE) solved through statistical Monte Carlo (MC) method provides the accurate dose calculation in radiotherapy. This work is to investigate the alternative way for accurately solving LBTE using deterministic numerical method due to its possible advantage in computational speed from MC. Methods: Instead of using traditional spherical harmonics to approximate angular scattering kernel, our deterministic numerical method directly computes angular scattering weights, based on a new angular discretization method that utilizes linear finite element method on the local triangulation of unit angular sphere. As a Result, our angular discretization method has the unique advantage in positivity, i.e., to maintain all scattering weights nonnegative all the time, which is physically correct. Moreover, our method is local in angular space, and therefore handles the anisotropic scattering well, such as the forward-peaking scattering. To be compatible with image-guided radiotherapy, the spatial variables are discretized on the structured grid with the standard diamond scheme. After discretization, the improved sourceiteration method is utilized for solving the linear system without saving the linear system to memory. The accuracy of our 3D solver is validated using analytic solutions and benchmarked with Geant4, a popular MC solver. Results: The differences between Geant4 solutions and our solutions were less than 1.5% for various testing cases that mimic the practical cases. More details are available in the supporting document. Conclusion: We have developed a 3D LBTE solver based on a new angular discretization method that guarantees the positivity of scattering weights for physical correctness, and it has been benchmarked with Geant4 for photon dose calculation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Nonequilibrium molecular dynamics study of ring polymer melts under shear and elongation flows: A comparison with their linear analogs

We present detailed results for the structural and rheological properties of unknotted and unconcatenated ring polyethylene (PE) melts under shear and elongation flows via direct atomistic nonequilibrium molecular dynamics simulations. Short (C{sub 78}H{sub 156}) and long (C{sub 400}H{sub 800}) ring PE melts were subjected to planar Couette flow (PCF) and planar elongational flow (PEF) across a wide range of strain rates from linear to highly nonlinear flow regimes. The results are analyzed in detail through a direct comparison with those of the corresponding linear polymers. We found that, in comparison to their linear analogs, ring melts possess rather compact chain structures at or near the equilibrium state and exhibit a considerably lesser degree of structural deformation with respect to the applied flow strength under both PCF and PEF. The large structural resistance of ring polymers against an external flow field is attributed to the intrinsic closed-loop configuration of the ring and the topological constraint of nonconcatenation between ring chains in the melt. As a result, there appears to be a substantial discrepancy between ring and linear systems in terms of their structural and rheological properties such as chain orientation, the distribution of chain dimensions, viscosity, flow birefringence, hydrostatic pressure, the pair correlation function, and potential interaction energies. The findings and conclusions drawn in this work would be a useful guide in future exploration of the characteristic dynamical and relaxation mechanisms of ring polymers in bulk or confined systems under flowing conditions.

42 ENGINEERING↗