Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlinear function”

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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Enhanced extreme ultraviolet high-harmonic generation from chromium-doped magnesium oxide

High-order harmonic generation (HHG) from crystals is emerging as a new ultrashort source of coherent extreme ultraviolet (XUV) light. Doping the crystal structure can offer a new way to control the source properties. Here, we present a study of HHG enhancement in the XUV spectral region from an ionic crystal, using dopant-induced vacancy defects, driven by a laser centered at a wavelength of 1.55 μm. Our numerical simulations based on solutions of the semiconductor Bloch equations and density-functional theory are supported by our experimental observations and demonstrate an increase in the XUV high harmonic yield from doped bulk magnesium oxide (MgO) compared to undoped MgO, even at a low defect concentration. The anisotropy of the harmonic emission as a function of the laser polarization shows that the pristine crystal's symmetry is preserved. Our study paves the way toward the control of HHG in solids with complex defects caused by transition-metal doping.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Dark Energy Survey Year 3 results: Exploiting small-scale information with lensing shear ratios

Using the first three years of data from the Dark Energy Survey (DES), we use ratios of small-scale galaxy-galaxy lensing measurements around the same lens sample to constrain source redshift uncertainties, intrinsic alignments and other systematics or nuisance parameters of our model. Instead of using a simple geometric approach for the ratios as has been done in the past, we use the full modeling of the galaxy-galaxy lensing measurements, including the corresponding integration over the power spectrum and the contributions from intrinsic alignments and lens magnification. We perform extensive testing of the small-scale shear-ratio (SR) modeling by studying the impact of different effects such as the inclusion of baryonic physics, nonlinear biasing, halo occupation distribution descriptions and lens magnification, among others, and using realistic N -body simulations of the DES data. We validate the robustness of our constraints in the data by using two independent lens samples with different galaxy properties, and by deriving constraints using the corresponding large-scale ratios for which the modeling is simpler. The results applied to the DES Y3 data demonstrate how the ratios provide significant improvements in constraining power for several nuisance parameters in our model, especially on source redshift calibration and intrinsic alignments. For source redshifts, SR improves the constraints from the prior by up to 38% in some redshift bins. Such improvements, and especially the constraints it provides on intrinsic alignments, translate to tighter cosmological constraints when shear ratios are combined with cosmic shear and other 2pt functions. In particular, for the DES Y3 data, SR improves S 8 constraints from cosmic shear by up to 31%, and for the full combination of probes ( 3 × 2 pt ) by up to 10%. The shear ratios presented in this work are used as an additional likelihood for cosmic shear, 2 × 2 pt and the full 3 × 2 pt in the fiducial DES Y3 cosmological analysis.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Bayesian Optimization Framework for Imperfect Data or Models

Conventional Bayesian optimization methods implicitly assume that the data and model being optimized are “perfect.” This assumption leads to inaccurate posterior probability distribution functions (PDFs) when applied to “imperfect” data or models. The new Bayesian optimization framework presented in this report provides a way to parameterize the effect of imperfections usually encountered in a prior PDF of generalized data or a model on the posterior PDF. The effects of imperfections are parameterized by a set of constraints imposed on the posterior expectation values of deviations between the data and the model and on their covariance matrix elements. A particular set of values for these constraints conveys an evaluator’s best estimate of the effect of imperfections on the corresponding posterior expectation values. When a prior PDF of generalized data is assumed to be normal, an expression for a posterior PDF satisfying an arbitrary set of constraints is derived analytically for linear models. An analogous iterative algorithm is given for nonlinear models. The corresponding posterior PDF should be used to estimate any posterior expectation values in the presence of imperfections parameterized by that set of constraints. A posterior PDF of a conventional Bayesian optimization method is recovered analytically when all evaluator-specified constraints are set to zero (i.e., in the absence of any imperfections). The analytical expressions derived in this report for normal PDFs and linear models were verified numerically by a Metropolis–Hastings Monte Carlo method. The methods presented herein could be applied to any kind of data or models, including differential cross-section data or integral benchmark experiments.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Heterobilayer moiré magnets: Moiré skyrmions and commensurate-incommensurate transitions

In this work, we study untwisted heterobilayers of ferromagnetic and antiferromagnetic van der Waals materials, with in particular a Dzyaloshinskii-Moriya interaction in the ferromagnetic layer. A continuum low-energy field theory is utilized to study such systems. We develop a phase diagram as a function of the strength of interlayer exchange and Dzyaloshinskii-Moriya interactions, combining perturbative and strong coupling analyses with numerical simulations using Landau-Lifshitz-Gilbert equations. Various moiré-periodic commensurate phases are found, and the commensurate-incommensurate transition is discussed. Among the commensurate phases, we observe an interesting skyrmion lattice phase wherein each moiré unit-cell hosts one skyrmion.

2-dimensional systems↗
DOE OSTI2023

Learning continuous models for continuous physics

Abstract Dynamical systems that evolve continuously over time are ubiquitous throughout science and engineering. Machine learning (ML) provides data-driven approaches to model and predict the dynamics of such systems. A core issue with this approach is that ML models are typically trained on discrete data, using ML methodologies that are not aware of underlying continuity properties. This results in models that often do not capture any underlying continuous dynamics—either of the system of interest, or indeed of any related system. To address this challenge, we develop a convergence test based on numerical analysis theory. Our test verifies whether a model has learned a function that accurately approximates an underlying continuous dynamics. Models that fail this test fail to capture relevant dynamics, rendering them of limited utility for many scientific prediction tasks; while models that pass this test enable both better interpolation and better extrapolation in multiple ways. Our results illustrate how principled numerical analysis methods can be coupled with existing ML training/testing methodologies to validate models for science and engineering applications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Study on Intelligent Compaction-Equipment Logistics Scheduling and Propagation Characteristics of Vibration Wave in Nonlinear Systems with Multistability Based on Field Test

The equipment scheduling and propagation characteristics of vibration wave from vibratory roller ⟶ filling material nonlinear systems with multistability are the core problems of subgrade intelligent construction technology, and the logistics scheduling of the equipment is directly related to the construction efficiency. Aiming at the shortages, one typical subgrade located at the Gu’an station of Beijing-Xiong’an city railway is selected to research and finish the field tests; some findings are shown as follows: first, some valuable suggestions about the logistics scheduling of intelligent equipment are proposed, which can break the barriers between the organizations and improve construction efficiency; second, when the vibration wave propagates from the vibratory roller ⟶ surface of filling material ⟶ different buried depths of filling material, the peak acceleration of vibration wave gradually decreases and is hyperbolic distribution approximately. At the same time, the sensitive of attenuation is shown as follows: Z R 2 >0.9), and the concrete functional relationship among different amplitudes of harmonics can be summarized as y = Ae −BX ; fourth, the vibration energy is mainly concentrated near 10–30 Hz in the vibratory roller, but when the vibration wave propagates from vibratory roller⟶filling material, the vibration energy gradually decreases with the increase in depth, and the marginal spectrum gradually changes from one peak to two peaks, that is, 30–50 Hz and 50–100 Hz; fifth, the vibration energy in the vibrational wheel is distributed averagely in the compaction process, and the effective compaction time is two seconds, which will be helpful for revealing the propagation characteristics of vibration wave, optimizing the compaction quality control models and providing some support for the development of intelligent compaction theory of railway subgrade.

Han, Yixuan↗
DOE OSTI2020

Kinetic effects of self-organization in low-temperature magnetized plasmas (Final Report for DE-SC0019045: 9/1/2018-8/31/2019)

The research goal of this Early Career Research Program project is to advance the science of low temperature magnetized plasmas using first-principles kinetic, e.g., solving the Boltzmann equation, models. A better understanding of the plasma transport properties, such as electrical resistivity and thermal conductivity must be obtained to better control and design magnetized plasma sources. A key physical phenomenon in low-temperature magnetized plasmas is the self-organization of plasma flows, e.g., coherent plasma structures and transition between different discharge modes. Understanding the physics of coherent plasma structures has been identified as one of the four important frontier topics in a recent workshop by the Office of Fusion Energy Science. These self-emerging plasma patterns results from the nonlinear coupling between plasma constituents (ions, electrons, and neutral atoms), electromagnetic fields, and materials. The physics of low-temperature magnetized plasmas is particularly complex because (i) the collisionless instabilities affect collisional phenomena (ionization, transport, etc.) and vice versa and (ii) the presence of materials (sheath, electron injection, etc.) influences the bulk plasma properties in low-temperature plasma (LTP) devices. Computational models of these plasma flows remain challenging because the plasma density can vary several orders of magnitudes, a wide range of temporal and spatial scales must be resolved, the phenomena are inherently three-dimensional, and various physical processes, such as wall interaction, instabilities, inelastic collisions, etc., must be simultaneously and self-consistently taken into account. Understanding the fundamental transport mechanisms in cross-field configurations enables additional controllability of plasma properties and electron energy distribution functions (EEDFs).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Towards interpretable Cryo-EM: disentangling latent spaces of molecular conformations

Molecules are essential building blocks of life and their different conformations (i.e., shapes) crucially determine the functional role that they play in living organisms. Cryogenic Electron Microscopy (cryo-EM) allows for acquisition of large image datasets of individual molecules. Recent advances in computational cryo-EM have made it possible to learn latent variable models of conformation landscapes. However, interpreting these latent spaces remains a challenge as their individual dimensions are often arbitrary. The key message of our work is that this interpretation challenge can be viewed as an Independent Component Analysis (ICA) problem where we seek models that have the property of identifiability. That means, they have an essentially unique solution, representing a conformational latent space that separates the different degrees of freedom a molecule is equipped with in nature. Thus, we aim to advance the computational field of cryo-EM beyond visualizations as we connect it with the theoretical framework of (nonlinear) ICA and discuss the need for identifiable models, improved metrics, and benchmarks. Moving forward, we propose future directions for enhancing the disentanglement of latent spaces in cryo-EM, refining evaluation metrics and exploring techniques that leverage physics-based decoders of biomolecular systems. Moreover, we discuss how future technological developments in time-resolved single particle imaging may enable the application of nonlinear ICA models that can discover the true conformation changes of molecules in nature. The pursuit of interpretable conformational latent spaces will empower researchers to unravel complex biological processes and facilitate targeted interventions. This has significant implications for drug discovery and structural biology more broadly. More generally, latent variable models are deployed widely across many scientific disciplines. Thus, the argument we present in this work has much broader applications in AI for science if we want to move from impressive nonlinear neural network models to mathematically grounded methods that can help us learn something new about nature.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

SIF and Vegetation Indices in the US Midwestern Agroecosystems, 2016-2021, Carbon Monitoring System (CMS)

This dataset contains half-hourly ground solar-induced chlorophyll fluorescence (SIF) and vegetation indices including NDVI, EVI, Red edge chlorophyll index, green chlorophyll index, and photochemical reflectance index at seven crop sites in Nebraska and Illinois for the period 2016-2021. Four sites were located at Eddy Covariance (EC) tower sites (sites US-Ne2, US-Ne3, US-UiB, and US-UiC), and three sites were located on private farms (sites Reifsteck, Rund, and Reinhart). The sites were either miscanthus, corn-soybean rotation or corn-corn-soybean rotation. The spectral data for SIF retrieval and hyperspectral reflectance for vegetation index calculation were collected by the FluoSpec2 system, installed near planting, and uninstalled after harvest to collect whole growing-season data. Raw nadir SIF at 760 nm from different algorithms (sFLD, 3FLD, iFLD, SFM) are included. SFM_nonlinear and SFM_linear represent the Spectral fitting method (SFM) with the assumption that fluorescence and reflectance change with wavelength non-linearly and linearly, respectively. Additional data include two SIF correction factors including calibration coefficient adjustment factor (f_cal_corr_QEPRO) and upscaling nadir SIF to eddy covariance footprint factor (ratio_EC footprint, SIF pixel), and measured FPAR from quantum sensors and Rededge NDVI calculated FPAR. The data are provided in comma-separated values (CSV) format.

Ameriflux↗
DOE OSTI2024

Nonperturbative Zou-Wang-Mandel effect

The Zou-Wang-Mandel (ZWM) effect is a remarkable consequence of photon indistinguishability and continuous-variable entanglement in which an optical phase shift is imprinted on photonic modes associated with optical paths that that do not pass through the phase shift source. By bringing the canonical formalism of continuous-variable Gaussian states to bear on the mode-structure of the ZWM experiment, we show that the physical consequence of implementing optical path identity is a renormalization of quadrature squeezing which governs the entanglement of four effective optical modes. Nonperturbative expressions for the ZWM interference patterns and normalized first-order coherence function are derived. Generalizations to $\mathscr{H}$-graph states with more than four modes directly follow from the general method used to analyze the minimal example. Here, we show that a ZWM interferometer with a laser-seeded signal mode, which estimates an idler phase shift by detecting photons that did not propagate through the phase shift, exhibits an optimal sensitivity comparable to that of a laser-seeded SU(1,1) interferometer if path identity is implemented with high fidelity.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Ferroelectricity in a semiconducting all-inorganic halide perovskite

Ferroelectric semiconductors are rare materials with both spontaneous polarizations and visible light absorptions that are promising for designing functional photoferroelectrics, such as optical switches and ferroelectric photovoltaics. The emerging halide perovskites with remarkable semiconducting properties also have the potential of being ferroelectric, yet the evidence of robust ferroelectricity in the typical three-dimensional hybrid halide perovskites has been elusive. Here, we report on the investigation of ferroelectricity in all-inorganic halide perovskites, CsGeX 3 , with bandgaps of 1.6 to 3.3 eV. Their ferroelectricity originates from the lone pair stereochemical activity in Ge (II) that promotes the ion displacement. This gives rise to their spontaneous polarizations of ~10 to 20 μC/cm 2 , evidenced by both ab initio calculations and key experiments including atomic-level ionic displacement vector mapping and ferroelectric hysteresis loop measurement. Furthermore, characteristic ferroelectric domain patterns on the well-defined CsGeBr 3 nanoplates are imaged with both piezo-response force microscopy and nonlinear optical microscopic method.

36 MATERIALS SCIENCE↗
DOE OSTI2024

On hydromagnetic wave interactions in collisionless, high- β plasmas

We describe the interaction of parallel-propagating Alfvén waves with ion-acoustic waves and other Alfvén waves, in magnetized, high- β collisionless plasmas. This is accomplished through a combination of analytical theory and numerical fluid simulations of the Chew–Goldberger–Low (CGL) magnetohydrodynamic (MHD) equations closed by Landau-fluid heat fluxes. An asymptotic ordering is employed to simplify the CGL-MHD equations and derive solutions for the deformation of an Alfvén wave that results from its interaction with the pressure anisotropy generated either by an ion-acoustic wave or another, larger-amplitude Alfvén wave. The difference in time scales of acoustic and Alfvénic fluctuations at high- β means that interactions that are local in wavenumber space yield little modification to either mode within the time it takes the acoustic wave to Landau damp away. Instead, order-unity changes in the amplitude of Alfvénic fluctuations can result after interacting with frequency-matched acoustic waves. Additionally, we show that the propagation speed of an Alfvén-wave packet in an otherwise homogeneous background is a function of its self-generated pressure anisotropy. This allows for the eventual interaction of separate co-propagating Alfvén-wave packets of differing amplitudes. The results of the CGL-MHD simulations agree well with these predictions, suggesting that theoretical models relying on the interaction of these modes should be reconsidered in certain astrophysical environments. Applications of these results to weak Alfvénic turbulence and to the interaction between the compressive and Alfvénic cascades in strong, collisionless turbulence are also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Effect of Soil Properties and Input Motion on Site Amplification Using Validated Nonlinear Soil Model

Here, three-dimensional (3-D) nonlinear site response analyses are conducted using finite element models of actual soil profiles from ten nuclear power plant (NPP) sites in the United States to investigate the effects of soil properties and input motions on site amplification. The modeling approach developed in this study combines several novel elements, such as 3-D analysis (including vertical motions), nonlinear inelastic behavior of soil (strain-dependent shear modulus reduction and hysteretic damping), formulation of nonreflecting boundary conditions at the base, and generation of realistic outcrop ground motions for specific sites. All these elements of the modeling approach are first validated using actual data from five earthquakes at three downhole array stations recorded in the Kiban-Kyoshin network (KiK-net), Japan. The same approach is then used to develop site models of ten NPP sites in the United States and corresponding ground motions that are spectrally matched to the site hazard spectra. Eight sets of three-component input motions are used in the study and are categorized on the basis of presence or absence of a near-field pulse in the seed ground motions used for spectral matching. It is found that all sites retain a definite site amplification function regardless of the input motion, provided that the seed motion is spectrally matched to the site hazard spectra. The magnitude of site amplification and frequencies at which they occur depend upon soil properties, particularly the shear wave velocity profile and the constitutive relationship (strain-dependent shear modulus reduction and hysteretic damping) of soil. Amplification of spectral acceleration in the vertical direction (up-down motion) is found to be just as much as, if not more than, the amplification in the horizontal direction. Peak shear strain is found to be about 20% larger for near-field motions compared to far-field motions whereas maximum horizontal site amplification for far-field motions is found to be consistently larger than that of near-field motions, even though the differences between the two remain within the scatter resulting from individual ground motions.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2021

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 nonlinear. The literature describes the numerical methods that are used to resolve this classical system, including modern tools, such as the Runge–Kutta fourth–order method. Finding the derivate of closed-form solutions for the EOM is more difficult and, as mentioned in the literature, discovering the solution is not always possible for all the EOM. Fortunately, a few examples are available that serve as a guide to move further in this topic. The purpose of this paper is to find a methodology that will produce the solutions for a given subset of EOMs that fulfill certain requisites. This paper summarizes the literature available on this topic and then follows with the derivation of the EOM using the Euler–Lagrange method. The Routhian method will be used to reduce the size of the expression, and it continues with the formulation of the classical cubic function, ƒ(u), through a novel process. The roots of ƒ(u) are of the utmost importance in finding 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.

Laos, Hector↗
DOE OSTI2020

Angle-resolved nonresonant two-photon single ionization of argon using 9.3-eV photons produced via high-order harmonic generation

Here we present an experimental study on the photoionization dynamics of nonresonant one-color two-photon single valence ionization of neutral argon atoms. Using 9.3-eV photons produced via high harmonic generation and a three-dimensional momentum imaging spectrometer, we detect the photoelectrons and ions produced from nonresonant two-photon ionization in coincidence. Photoionization from the $3p$ orbital produces a photoelectron scattering wave function with $p$- and $f$-partial-wave components, which interfere and result in a photoelectron angular distribution with peak amplitude perpendicular to the vacuum ultraviolet polarization. The comparison between the present results and the two previous sets of theoretical calculations [Phys. Rev. A 44, 324 (1991), and J. Phys. B: At. Mol. Phys. 16, 2737 (1983)] indicates that electron-electron correlation contributes appreciably to the two-photon ionization dynamics.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Hydrokinetic Power Conversion Using Vortex-Induced Oscillation with Cubic Restoring Force

A cubic-spring restoring function with high-deformation stiffening is introduced to passively improve the harnessed marine hydrokinetic power by using flow-induced oscillations/vibrations (FIO/V) of a cylinder. In these FIO/V experiments, a smooth, rigid, single-cylinder on elastic end-supports is tested at Reynolds numbers ranging from 24,000 < Re < 120,000. The parameters of the tested current energy converter (CEC) are cubic stiffness and linear damping. Using the second generation of digital virtual spring-damping (Vck) controller developed by the Marine Renewable Energy Laboratory (MRELab), the cubic modeling of the oscillator stiffness is tested. Experimental results show the influence of the parameter variation on the amplitude, frequency, energy conversion, energy efficiency, and power of the converter. All experiments are conducted in the low turbulence-free surface water (LTFSW) channel of the MRELab of the University of Michigan. The main conclusions are: (1) The nonlinearity in the cubic oscillator is an effective way to extend the vortex-induced vibration (VIV) upper branch, which results in higher harnessing power and efficiency compared to the linear stiffness cylinder converter. (2) Compared to the linear converter, the overall power increase is substantial. The nonlinear power optimum, occurring at the end of the VIV upper branch, is 63% higher than its linear counterpart. (3) The cubic stiffness converter with low harnessing damping achieves consistently good performance in all the VIV regions because of the hardening restoring force, especially at higher flow velocity.

13 HYDRO ENERGY↗
DOE OSTI2020

Variational nonlinear WKB in the Eulerian frame

Nonlinear WKB is a multiscale technique for studying locally plane-wave solutions of nonlinear partial differential equations (PDEs). Its application comprises two steps: (1) replacement of the original PDE with an extended system separating the large scales from the small and (2) reduction of the extended system to its slow manifold. In the context of variational fluid theories with particle relabeling symmetry, nonlinear WKB in the mean Eulerian frame is known to possess a variational structure. This much has been demonstrated using, for instance, the theoretical apparatus known as the generalized Lagrangian mean. On the other hand, the variational structure of nonlinear WKB in the conventional Eulerian frame remains mysterious. By exhibiting a variational principle for the extended equations from step (1) above, here we demonstrate that nonlinear WKB in the Eulerian frame is in fact variational. Remarkably, the variational principle for the extended system admits loops of relabeling transformations as a symmetry group. Noether’s theorem therefore implies that the extended Eulerian equations possess a family of circulation invariants parameterized by S 1 . As an illustrative example, we use our results to systematically deduce a variational model of high-frequency acoustic waves interacting with a larger-scale compressible isothermal flow.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Multilevel Hierarchical Decomposition of Finite Element White Noise with Application to Multilevel Markov Chain Monte Carlo

In this work we develop a new hierarchical multilevel approach to generate Gaussian random field realizations in an algorithmically scalable manner that is well suited to incorporating into multilevel Markov chain Monte Carlo (MCMC) algorithms. This approach builds off of other partial differential equation (PDE) approaches for generating Gaussian random field realizations; in particular, a single field realization may be formed by solving a reaction-diffusion PDE with a spatial white noise source function as the right-hand side. While these approaches have been explored to accelerate forward uncertainty quantification tasks, e.g., multilevel Monte Carlo, the previous constructions are not directly applicable to multilevel MCMC frameworks which build fine-scale random fields in a hierarchical fashion from coarse-scale random fields. Our new hierarchical multilevel method relies on a hierarchical decomposition of the white noise source function in $L^2$ which allows us to form Gaussian random field realizations across multiple levels of discretization in a way that fits into multilevel MCMC algorithmic frameworks. After presenting our main theoretical results and numerical scaling results to showcase the utility of this new hierarchical PDE method for generating Gaussian random field realizations, this method is tested on a four-level MCMC algorithm to explore its feasibility.

algebraic multigrid↗