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 505 records · Page 28

A method to test the coupling strength of the linear and nonlinear contributions to higher-order flow harmonics via event shape engineering

A Multi-Phase Transport (AMPT) model is used to study the efficacy of shape-engineered events to delineate the degree of coupling between the linear and nonlinear contributions to the higher-order flow harmonics v 4 and v 5 . Here, the study shows that the nonlinear contributions are strongly shape-dependent while the linear contributions are shape-independent, indicating little if any, coupling between the linear and nonlinear flow coefficients. The experimental verification of such patterns could be an invaluable tool for robust extraction of the linear and mode-coupled flow coefficients, especially for beam energies where the charged particle multiplicity and the event statistics precludes the use of current methods to establish the coupling strength.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Poisson-Boltzmann Theory with Non-Linear Ion Correlations

The Poisson-Boltzmann (PB) theory is widely used to depict ionic systems in a mean-field manner without considering the correlations between ions. A modified theory including ion correlation effect while retaining the simplicity of the original PB theory is desired. Although recent theoretical efforts successfully incorporate ion correlations into the PB equation as the self-energy, the results deviate significantly from our molecular dynamics (MD) simulation results due to the linear form of the self-energy. We herein reinterpret the self-energy term from a physics point of view, on the basis of which we present a new way for incorporating ion correlations into the original PB equation by utilizing the Green’s function with a non-linear form of the self-energy. Justified by our MD simulation results, our modified PB theory with non-linear ion correlations surpasses linearly formed modified PB theories.

Su, Mao↗

Structure formation with massive neutrinos: going beyond linear theory

We compute non-linear corrections to the matter power spectrum taking the time- and scale-dependent free-streaming length of neutrinos into account. We adopt a hybrid scheme that matches the full Boltzmann hierarchy to an effective two-fluid description at an intermediate redshift. The non-linearities in the neutrino component are taken into account by using an extension of the time-flow framework. We point out that this remedies a spurious behaviour that occurs when neglecting non-linear terms for neutrinos. This behaviour is related to how efficiently short modes decouple from long modes and can be traced back to the violation of momentum conservation if neutrinos are treated linearly. Furthermore, we compare our results at next to leading order to various other methods and quantify the accuracy of the fluid description. Due to the correct decoupling behaviour of short modes, the two-fluid scheme is a suitable starting point to compute higher orders in perturbations or for resummation methods.

79 ASTRONOMY AND ASTROPHYSICS↗

f(R) gravity on non-linear scales: the post-Friedmann expansion and the vector potential

Many modified gravity theories are under consideration in cosmology as the source of the accelerated expansion of the universe and linear perturbation theory, valid on the largest scales, has been examined in many of these models. However, smaller non-linear scales offer a richer phenomenology with which to constrain modified gravity theories. Here, we consider the Hu-Sawicki form of f(R) gravity and apply the post-Friedmann approach to derive the leading order equations for non-linear scales, i.e. the equations valid in the Newtonian-like regime. We reproduce the standard equations for the scalar field, gravitational slip and the modified Poisson equation in a coherent framework. In addition, we derive the equation for the leading order correction to the Newtonian regime, the vector potential. We measure this vector potential from f(R) N-body simulations at redshift zero and one, for two values of the f{sub R{sub 0}} parameter. We find that the vector potential at redshift zero in f(R) gravity can be close to 50% larger than in GR on small scales for |f{sub R{sub 0}}|=1.289 × 10{sup −5}, although this is less for larger scales, earlier times and smaller values of the f{sub R{sub 0}} parameter. Similarly to in GR, the small amplitude of this vector potential suggests that the Newtonian approximation is highly accurate for f(R) gravity, and also that the non-linear cosmological behaviour of f(R) gravity can be completely described by just the scalar potentials and the f(R) field.

79 ASTRONOMY AND ASTROPHYSICS↗

Non-linearities in the Lyman-α forest and in its cross-correlation with dark matter halos

Here three-dimensional correlations of the Lyman-α (Lyα) forest and cross correlations between the Lyα forest and quasars have been measured on large scales, allowing a precise measurement of the baryon acoustic oscillation (BAO) feature at redshifts z > 2. These 3D correlations are often modelled using linear perturbation theory, but full-shape analyses to extract cosmological information beyond BAO will require more realistic models capable of describing non-linearities present at smaller scales. We present a measurement of the Lyα forest flux power spectrum from large hydrodynamic simulations — the Sherwood simulations — and compare it to different models describing the small-scale deviations from linear theory. We confirm that the model presented in Arinyo-i-Prats et al. (2015) fits the measured 3D power up to k = 10 h Mpc -1 with an accuracy better than 5%, and show that the same model can also describe the 1D correlations with similar precision. We also present, for the first time, an equivalent study for the cross-power spectrum of halos with the Lyα forest, and we discuss different challenges we face when modelling the cross-power spectrum beyond linear scales. We make all our measured power spectra public in https://github.com/andreufont/sherwoo_p3d. This study is a step towards joint analyses of 1D and 3D flux correlations, and towards using the quasar-Lyα cross-correlation beyond BAO analyses.

79 ASTRONOMY AND ASTROPHYSICS↗

A direct detection method of galaxy intrinsic ellipticity-gravitational shear correlation in non-linear regimes using self-calibration

Intrinsic alignment (IA) of galaxies is a challenging source of contamination in the Cosmic shear (GG) signals. The galaxy intrinsic ellipticity-gravitational shear (IG) correlation is generally the most dominant component of such contamination for cross-correlating redshift bins. One of the most effective techniques to mitigate such contamination is the self-calibration (SC) method which extracts the IG correlation and allows for its removal from the GG signal. In a photometric survey, the SC method first extracts the galaxy number density-galaxy intrinsic ellipticity (gI) correlation from the observed galaxy-galaxy lensing correlation using the redshift dependence of lens-source pairs. The IG correlation is computed through a scaling relation using the gI correlation and other lensing observables. The applicability of the SC method has so far been focused on the linear IA scales and the linear galaxy bias. We extend the SC method beyond the linear regime by modifying its scaling relation which can account for the non-linear galaxy bias model and various IA models. In this study, we provide a framework to detect the IG correlation for the redshift bins for source galaxies for the proposed year 1 survey of the Rubin Legacy Survey of Space and Time (LSST Y1). We tested the method for the tidal alignment and tidal torquing (TATT) model of IA and we found that the scaling relation is accurate within 10% and 20% for cross-correlating and auto-correlating redshift bins, respectively. Hence the suppression of IG contamination in observed GG correlation can be accomplished with a factor of 10 and 5, for cross-correlating and auto-correlating redshift bins, respectively. We tested the method's robustness and found that the suppression of IG contamination by a factor of 5 is still achievable for all combinations of cross-correlating bins even with the inclusion of a moderate amount of uncertainties on IA and bias parameters, respectively. We also make available, a branch of the code FAST-PT to provide gI correlations up to 1-loop order term used by the new SC method.

gravitational lensing↗

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↗

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↗

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↗

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)↗

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↗

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↗

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↗

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↗

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↗

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 ↗

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↗

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↗