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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Optimization to Generate Equations of State for Hydrogen Production

On a high level, the larger project in question, HydroGEN, aims to develop software used for finding equations of state (EOS) to optimize catalyst configuration for H 2 production through water splitting. In particular, this summer project focused on solving the nonlinear equations used in fitting the equations. This problem involved using Python to solve a linear system with nonlinear constraints. In order for this to be achieved, Pyomo was used to build a model and the solver Ipopt, interior point optimizer, was used. Pyomo is a Python-based language developed at Sandia; it is an optimization modeling language. Rather than solving the entire problem at once, a toy problem was created, simplifying the problem down to the most important focus. This problem had a known solution, comparable to the calculated solution to assess accuracy and as progress was made towards finding solutions, complexity was gradually added to the problem. After building and solving the toy problem, it was found that it gave reasonably accurate solutions, better compared to the two existing solvers previously used with this project in terms of functionality. The solver is now ready for implementation into the project’s main software.

08 HYDROGEN↗
DOE OSTI2025

Chemical Control of Symmetry and Bandgap in Tungsten Oxyhalide van der Waals Semiconductors

Tunability in solid-state materials is essential for testing theory, discovering quantum phases, and enabling functionality. Layered van der Waals (vdW) semiconductors offer a unique platform, providing new degrees of freedom at the two-dimensional (2D) limit through exfoliation and external controls. Here, in this study, we demonstrate tunability of symmetry and electronic structure via halogen substitution in a family of layered vdW tungsten oxyhalides. Substituting the halogens in WO 2 X 2 (X = I, Br, Cl) tunes the bandgap across a broad energy range and modifies the structural symmetry from centrosymmetric to noncentrosymmetric. By alloying WO 2 I 2– y Br y , we continuously tune the polar distortion and optical gap across the visible range. These insights into halogen substitution effects on symmetry and electronic structure lay the foundation for new tunable vdW semiconductors for optoelectronics and nonlinear optics.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Shifts in valence states in bimetallic MXenes revealed by electron energy-loss spectroscopy (EELS)

MXenes are an emergent class of two-dimensional materials with a very wide spectrum of promising applications. The synthesis of multiple MXenes, specifically solid-solution MXenes, allows fine tuning of their properties, expands their range of applications, and leads to enhanced performance. The functionality of solid-solution MXenes is closely related to the valence state of their constituents: transition metals, oxygen, carbon, and nitrogen. However, the impact of changes in the oxidation state of elements in MXenes is not well understood. In this work, three interrelated solid-solution MXene systems (Ti$_{2–y}$Nb$_y$CT$_x$, Nb$_{2–y}$V$_y$CT$_x$, and Ti$_{2–y}$V$_y$CT$_x$) were investigated with scanning transmission electron microscopy and electron energy-loss spectroscopy to determine the localized valence states of metals at the nanoscale. The analysis demonstrates changes in the electronic configuration of V upon modification of the overall composition and within individual MXene flakes. These shifts of oxidation state can explain the nonlinear optical and electronic features of solid-solution MXenes. Vanadium appears to be particularly sensitive to modification of the valence state, while titanium maintains the same oxidation state in Ti–Nb and Ti–V MXenes, regardless of stoichiometry. Here, the study also explains Nb's influential role in the previously observed electronic properties in the Nb–V and Nb–Ti systems.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Multifunctional Fiber-Reinforced Polymer Composites for Damage Detection and Memory

Self-structural health monitoring (SHM) functionalities for fiber-reinforced polymer composites have become highly sought after to ensure the structural safety of newly advancing components in the automotive, civil, mechanical, and aerospace industries. This paper introduces a self-damage detection and memory (SDDM) hybrid composite material, where the structural carbon fiber tow is transformed into a piezoresistive sensor network, and the structural glass fiber operates as electrical insulation. In this study, SDDM specimens were fabricated, and tensile and impact tests were performed. The tensile tests of SDDM specimens find two distinct loading peaks: first where the carbon fiber fails, and second where the glass fiber fails. A linear correlation was observed between the carbon fiber resistance and composite strain up to a threshold, beyond which a sharp nonlinear increase in resistance occurred. The resistance then approached infinity, coinciding with the first loading peak and failure of the carbon fiber elements. This demonstrates the potential for a damage early warning threshold. Additionally, the effect of stitching the sensor tow in a zig-zag pattern over a large area was investigated using tailored fiber placement (TFP) of 1-loop, 3-loop, and 5-loop specimens. Tensile testing found that increasing the number of loops improved the sensor’s accuracy for strain sensing. Furthermore, impact tests were conducted, and as the impact energy progressively increased, the sensor resistance permanently increased. This illustrates a capability for self-memory of microdamage throughout the life cycle of the structure, potentially useful for predicting the remaining life of the composite.

Demo, Luke B. (ORCID:0000000208806522)↗
DOE OSTI2023

New Developments and Capabilities Within WEC-Sim: Preprint

WEC-Sim is an open-source software for simulating wave energy converters, which has been actively developed and applied since its initial release in 2014 to simulate a wide variety of device archetypes. WEC-Sim is developed jointly by the National Renewable Energy Laboratory (NREL) and Sandia National Laboratories (Sandia) within the MATLAB/SIMULINK environment. A general wave-to-wire model begins with a deployment site resource characterization, which is used to complete the hydrodynamic simulation of wave energy converters (WEC), with the power generation profile imported to a grid simulator to understand the influence on the local electrical network. While modeling the entire wave-to-wire is difficult and encompasses multiple time scales and physics, WEC-Sim is focused on the hydrodynamics simulation to predict, analyze, and optimize WEC dynamics and power performance. WEC-Sim simulations are performed in the time domain based on the radiation and diffraction method using hydrodynamics coefficients derived from boundary element method (BEM)-based frequency-domain potential flow solvers (e.g., WAMIT, NEMOH, Capytaine, or ANSYS-AQWA). With this level of modeling fidelity, WEC-Sim can handle floating body hydrodynamics, mechanical and electrical power generation methods, advanced control implementation, mooring systems, and other unique applications such as desalination. Additional WEC-Sim functionalities include pre-built Simulink blocks and MATLAB scripts that can simulate a wide range of floating systems and the corresponding auxiliary subsystems. The developers of WEC-Sim continue to release new versions of the software, at least annually, with our latest release in September 2022. These releases include bug fixes, updates to software documentation, as well as new features to expand WEC-Sim's capabilities to model a wide range of WEC concepts. This publication will highlight the new features added to WEC-Sim between versions 4.1.0 to 5.0.1 which spans over a two year period from June 2020 to September 2022. New features to be described will include topics such as continuous integration checks, revised Morison Element and nonlinear hydro implementations, run directly from Simulink (required for hardware-in-the-loop execution), BEMIO updates to import Capytaine BEM hydrodynamics, addition of cable blocks, and new wave visualization features.

TIDAL AND WAVE POWER↗
DOE OSTI2021

Validation and Demonstration of Control System Functional Capabilities within the IES Plug-and-Play Simulation Environment

The concept of an integrated energy system (IES) is meant to combine different energy technologies in synergistic ways to achieve a more secure and economical energy supply. The RAVEN-based HYBRID framework is used to find the optimal installed capacity and the optimal economical dispatch of each component of the IES. The new RAVEN plugin for grid and capacity optimization (HERON) only addresses the limits that affect the production variables and the corresponding rates of variation (explicit constraints). However, other variables are subject to constraints, and the associated limits should be accounted for (implicit constraints). In particular, for the power dispatch problem, the optimization algorithm takes into account the limits on the electrical power output and the corresponding hourly power variations but does not consider other constraints on process variables whose response affects the service life of the IES. This report describes a scheme that allows accounting for implicit constraints without increasing the size of the optimization problem. To obtain a more accurate approximation of the nonlinear dynamic behavior, a parametric version of the dynamic mode decomposition with control (DMDc) algorithm was developed to derive the state-space representation matrices of the IES components at different scheduling parameter. Thanks to this approach, a more accurate approximation of the system response can be obtained, the limits imposed by thermal mechanical implicit constraints can be translated into power dispatch limits, and the feedbacks to HERON power dispatcher can be provided. To assess the developed methodology, a power dispatching test case composed of three power generating and storage units (Balance of Plant, Secondary Energy Source, Thermal Energy Storage) was developed. The power output of each one of the three units was optimized to meet the imposed time-dependent load demand trajectory and to maximize the IES profitability by meeting both the explicit and implicit constraints.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

Deep Learning enabled spectral energy conversion for in situ exposure measurements

A detector-specific deep learning (DL) approach is presented for spectra-to-exposure conversion using large-format sodium iodide (NaI(Tl)) detectors deployed for in situ environmental radiation measurements in emergency response scenarios. Accurate determination of exposure from NaI spectra is challenging due to poor energy resolution, partial energy absorption, and the strong sensitivity of traditionally deployed analytical conversion methods to calibrated source geometry and pre-deployment assumptions. Here, to address these limitations, a multi-layer perceptron model was trained on a hybrid in situ /Monte Carlo dataset constructed to span a broad range of photon energies, spatial extents, and realistic deployment variability, representative of general in situ emergency response conditions. The DL model was evaluated against commonly fielded analytical approaches under matched simulation conditions, including a single-factor method, a G-function method, and a modeled pressurized ion chamber (PIC) baseline. This study was intentionally computational in scope to enable controlled, like-for-like comparisons between conversion techniques while minimizing confounding real-world variability. Comparison to the modeled PIC provides contextual benchmarking and is not intended as a field inter-comparison with deployed instruments. Across the evaluated 20 keV to 3 MeV energy range, the DL approach consistently exhibited higher accuracy and reduced variance relative to the analytical methods against a deterministically calculated exposure. This may indicate improved robustness to spectral complexity without reliance on source-, geometric-, or spectral region-specific optimization. While results do not represent real-world validation, the presented work demonstrates that deep learning may effectively learn the nonlinear detector response-to-exposure relationship for asymmetric NaI(Tl) detectors and offers a promising pathway for improving in situ exposure estimation using spectroscopic systems already integrated into initial real-time emergency response operations.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2020

A Physics-Driven Deep-Learning Network for Solving Nonlinear Inverse Problems

Geosteering inversion, which can be viewed as a nonlinear inverse problem, is an important technique used by directional drilling. Traditional methods that rely on iterative procedures and regularization are sensitive to the selection of initial values and can be slow due to convergence issues. In industrial applications, a lookup table is used to produce fast predictions. However, performance is not guaranteed by this approach due to the limitation of the hardware. In this paper, we propose a novel physics-driven deep-learning framework for providing a fast, accurate surrogate to solve the inverse problem. Particularly, leveraged by the forward physical model and 1D convolutional neural network (1D-CNN), the proposed method provides more reliable solutions to the inverse problem with improved performance. In addition, a new physics-driven loss function is introduced to accommodate both the model misfit and the data misfit. In conclusion, our experiments demonstrate the effectiveness of our method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

High‐order cubature rules for tetrahedra

Summary In this paper, we construct new high‐order numerical integration schemes for tetrahedra, with positive weights and integration points that are in the interior of the domain. The construction of cubature rules is a challenging problem, which requires the solution of strongly nonlinear algebraic (moment) equations with side conditions given by affine inequality constraints. We present a robust algorithm based on a sequence of three modified Newton procedures to solve the constrained minimization problem. In the literature, numerical integration rules for the tetrahedron are available up to order p =15. We obtain integration rules for the tetrahedron from p =2 to p =20, which are computed using multiprecision arithmetic. For p ≤15, our approach provides integration rules that have the same or fewer number of integration points than existing rules; for p =16 to p =20, our rules are new. Numerical tests are presented that verify the polynomial‐precision of the cubature rules. Convergence studies are performed for the integration of exponential, rational, weakly singular and trigonometric test functions over tetrahedra with flat and curved faces. In all tests, improvements in accuracy is realized as p is increased, though in some cases nonmonotonic convergence is observed.

Jaśkowiec, Jan↗
DOE OSTI2023

Ideal MHD induced temperature flattening in spherical tokamaks

This paper extends the analysis first presented in Jardin et al. [Phys. Rev. Lett. 128, 245001 (2022)] to more thoroughly examine the stability of spherical torus equilibrium to ideal magnetohydrodynamic (MHD) infernal modes and their nonlinear consequences. Here we demonstrate that in a 3D resistive magnetohydrodynamic (MHD) simulation of a NSTX discharge, anomalous transport can occur due to these instabilities. Here we generate a family of equilibrium of differing β and use this to show that these instabilities could explain the experimentally observed flattening of the electron temperature profile at modest β. The modes studied in this paper are found to occur with poloidal mode number m and toroidal mode number n when the ratio m/ n is in the range of 1.2–1.5, when the central safety factor is in this range or slightly lower, and when the central region has very low magnetic shear. Our analysis gives some insight as to why the unstable linear growth rates are oscillatory functions of the toroidal mode number n. We present a simulation of an initially stable configuration that passes through a stability boundary at a critical β as it is heated. We also show that a particular NSTX discharge is unstable to these modes over a timescale of several hundred ms. We conclude that these modes must be taken into account when performing predictive modeling. An appendix shows that similar modes can be found in R / a = 4 tokamaks for certain q-profiles and β values.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Nonlinear electrical transport near the metal–insulator transition in V 4 O 7 thin films

Nonlinear electrical transport associated with correlated electronic states has been widely investigated in transition-metal oxides near metal–insulator transitions. Here, in this study, we investigate nonlinear transport and threshold switching in sputter-deposited V 4 O 7 thin films grown on fused silica substrates. Temperature-dependent transport measurements show a metal–insulator transition near 240 K with negligible thermal hysteresis, defining the temperature scale that governs the nonlinear electrical response. Current–voltage measurements reveal reproducible bipolar threshold switching over a broad temperature range, with the threshold voltage decreasing systematically as the transition temperature is approached. Electrothermal finite-element simulations reproduce the measured switching characteristics and show that switching occurs when localized Joule heating drives a confined region of the device toward the transition temperature. The simulated peak local temperature at threshold lies close to the intrinsic metal–insulator transition temperature of V 4 O 7 , establishing that the nonlinear electrical response originates from electrothermal feedback acting on the strongly temperature-dependent conductivity near the transition. These results identify V 4 O 7 thin films as a model correlated-oxide system in which volatile threshold switching emerges from proximity to an extended metal–insulator transition.

25 ENERGY STORAGE↗
DOE OSTI2021

Mitigation of the internal p-n junction in CoS 2 -contacted FeS 2 single crystals: Accessing bulk semiconducting transport

Pyrite FeS 2 is an outstanding candidate for a low-cost, nontoxic, sustainable photovoltaic material, but efficient pyrite-based solar cells are yet to materialize. Recent studies of single crystals have shed much light on this by uncovering a p-type surface inversion layer on n-type (S-vacancy doped) crystals, and the resulting internal p-n junction. This leaky internal junction likely plays a key role in limiting efficiency in pyrite-based photovoltaic devices, also obscuring the true bulk semiconducting transport properties of pyrite crystals. Here, we demonstrate complete mitigation of the internal p-n junction in FeS 2 crystals by fabricating metallic CoS 2 contacts via a process that simultaneously diffuses Co (a shallow donor) into the crystal, the resulting heavy n doping yielding direct Ohmic contact to the interior. Low-temperature bulk transport studies of controllably Co- and S-vacancy doped semiconducting crystals then enable a host of previously inaccessible observations and measurements, including determination of donor activation energies (which are as low as 5 meV for Co), observation of an unexpected second activated transport regime, realization of electron mobility up to 2100 cm 2 V –1 s –1 , elucidation of very different mobilities in Co- and S-vacancy-doped cases, and observation of an abrupt temperaturedependent crossover to bulk Efros-Shklovskii variable-range hopping, accompanied by an unusual form of nonlinear Hall effect. Aspects of the results are interpreted with the aid of first-principles electronic structure calculations on both Co- and S-vacancy-doped FeS 2 . Furthermore, this work thus demonstrates unequivocal mitigation of the internal p-n junction in pyrite single crystals, with important implications for both future fundamental studies and photovoltaic devices

36 MATERIALS SCIENCE↗
DOE OSTI2022

A modified model parametrization algorithm for solving a special type of heat and mass transfer systems

A new method for solving nonlinear heat and mass transfer design tasks was considered. Systems using the Number of Transfer Units (NTU) method are a special type of mathematical model of heat and mass exchangers. It was observed, that the NTU models in a form of differential-algebraic equations (DAEs) cannot be directly solved with higher values of NTU. The requirements for consistent initial conditions, as well as numerical limitations of DAEs solvers, result, that the solution to the considered design problems that cannot be obtained by a classical direct shooting procedure. To overcome the presented difficulties, the αDAE model optimization algorithm was adjusted for solving NTU-based models. The new approach consists of 3 main steps: 1) task discretization by a multiple-shooting approach, 2) design an appropriate function $f_{NTU}$(α) to effectively influence the variability of the state variables described by dynamical relations, 3) the iterative numerical optimization algorithm for the new parametrized system. Moreover, computations can be performed by a chosen numerical optimization approach, which can be communicated with an available outer procedure for solving differential-algebraic equations. The presented algorithm was implemented and applied to solve the design task with the NTU model of a counter-flow exchanger. Here, the new approach was used to modify the system dynamics to influence the difficulty of the considered problem. Finally, the presented method enabled failure-free numerical computations for the higher values of the NTU parameter.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

The Aemulus Project. V. Cosmological Constraint from Small-scale Clustering of BOSS Galaxies

We analyze clustering measurements of BOSS galaxies using a simulation-based emulator of two-point statistics. We focus on the monopole and quadrupole of the redshift-space correlation function, and the projected correlation function, at scales of 0.1 ~ 60 h -1 Mpc. Although our simulations are based on wCDM with general relativity (GR), we include a scaling parameter of the halo velocity field, γ f , defined as the amplitude of the halo velocity field relative to the GR prediction. We divide the BOSS data into three redshift bins. After marginalizing over other cosmological parameters, galaxy bias parameters, and the velocity scaling parameter, we find fσ 8 (z = 0.25) = 0.413 ± 0.031, fσ 8 (z = 0.4) = 0.470 ± 0.026, and fσ 8 (z = 0.55) = 0.396 ± 0.022. Compared with Planck observations using a flat Lambda cold dark matter model, our results are lower by 1.9σ, 0.3σ, and 3.4σ, respectively. These results are consistent with other recent simulation-based results at nonlinear scales, including weak lensing measurements of BOSS LOWZ galaxies, two-point clustering of eBOSS LRGs, and an independent clustering analysis of BOSS LOWZ. All these results are generally consistent with a combination of ${\gamma }_{f}^{1/2}{\sigma }_{8}\approx 0.75$. We note, however, that the BOSS data is well fit assuming GR, i.e., γf = 1. We cannot rule out an unknown systematic error in the galaxy bias model at nonlinear scales, but near-future data and modeling will enhance our understanding of the galaxy–halo connection, and provide a strong test of new physics beyond the standard model.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Active Site Engineering via Optimizing the Heterogeneous Support Structure for Single-Atom Catalysis

Supported single-atom catalysts show a large range of activities and selectivities that depend on the local environment of the catalytic sites. Here a theory-based optimization strategy is presented that is based on a density functional theory determination of the transition states and intermediates for a low-dimensional coordinate representation of the heterogeneity of the active sites. The approach is applied to a vanadium catalyst on an amorphous SiO 2 support that involves a large kinetic network described using a full chemistry model. Without assuming a priori scaling relations or mechanism reduction, the optimal state of heterogeneity is found to lie at atomic configurations where the activation energies for two distinct key chemical processes are equal. It is found a posteriori that the behavior of the system is consistent with linear free energy scaling relations in the randomness parameters. The energetic span theory proves quite useful in reducing the full chemistry model to a small number of key reactions. The use of a nonlinear optimization algorithm in combination with energetic span theory provides significant simplification in treating disordered systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Continuous Integration, In-Code Documentation, and Automation for Nuclear Quality Assurance Conformance

The Multiphysics Object Oriented Simulation Environment (MOOSE) is an open-source, finite element framework for solving highly coupled sets of nonlinear equations. The development of the framework and applications occurs concurrently using an agile, continuous-integration software package. Included in the framework is an in-code, extensible documentation system. Using these two tools in union with the repository management tools GitHub and GitLab, a software quality plan was created and followed such that MOOSE and a MOOSE-based application (BISON) have been shown to meet the American Society of Mechanical Engineers’ Nuclear Quality Assurance-1 standard. The approach relies heavily on automation for both testing and documentation. The resulting effort demonstrates that a rigorous software quality plan may be implemented that incurs a minimal impact on day-to-day development of the software, satisfying the stringent guidelines necessary to operate the software in a safety function within a nuclear facility.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Soft Switching Over the Entire Line Cycle for a Quadruple Active Bridge DCX in a DC to Three-Phase AC Module

This paper is focused on a transformer-isolated quadruple active bridge (QAB) dc-dc converter loaded by three full-bridge dc-ac inverters. The QAB functions as a "DC trans-former" (DCX) in the dc-to-three-phase ac module. The QAB outputs provide time-varying power at twice the line frequency, which presents challenges in maintaining zero voltage switching (ZVS) on the secondary sides during low-power portions of the line cycle. It is shown how ZVS can be achieved in a phase, even at zero-power transfer, using a relatively small circulating current provided by the magnetizing inductance of the high-frequency transformer. The approach is particularly effective in high-voltage applications using SiC MOSFETs, where reductions in switching loss outweigh additional conduction losses due to the circulating currents. A detailed analysis of ZVS switching waveforms at the zero power transfer is presented, including effects of nonlinear device capacitances. Analytical expressions are given for the optimal values of the magnetizing inductance and dead times of the QAB primary and secondary bridges. The approach is verified by experimental results on a 600 V, 4 kW prototype, demonstrating greater than 98.5% efficiency from 1 kW to 4 kW, with a peak efficiency of 99.0%. Compared to conventionally operated prototype, a 54% reduction in total loss is achieved at rated power.

41 EE - Solar Energy Technologies Office (EE-4S)↗
DOE OSTI2023

Nonlinear second order electromagnetic gyrokinetic theory for a tokamak plasma

The steep plasma pressure gradient that forms at the edge of the high confinement, H-mode regime of tokamak operation provides free energy to drive electromagnetic micro-instabilities that are widely believed to influence the transport processes in this so-called pedestal region. This high pressure gradient also provides a high current density (bootstrap current), known to influence ballooning mode stability and to be important for driving kink modes in the ideal magneto-hydrodynamic plasma model (so-called peeling-ballooning modes). Furthermore, efficient, steady state future tokamak power plants must operate with a large bootstrap current in the core and especially concerning spherical tokamaks, confinement will be influenced by electromagnetic turbulence. To accommodate these important situations, conventional electromagnetic gyrokinetic theory is extended to incorporate neoclassical effects in the equilibrium drives, allowing $B_{\vartheta} ~ B_0$ (B0 is the confining magnetic field, and $B_{\vartheta}$ is its poloidal component). This provides a global gyrokinetic model that self-consistently captures the consequences of large bootstrap current fractions on the equilibrium distribution functions.

kink modes↗