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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Linear gyrokinetic simulations of microinstabilities within the pedestal region of H-mode NSTX discharges in a highly shaped geometry

Linear (local) gyrokinetic predictions of edge microinstabilities in highly shaped, lithiated and non-lihiated NSTX discharges are reported using the gyrokinetic code GS2. Microtearing modes dominate the non-lithiated pedestal top. The stabilization of these modes at the lithiated pedestal top enables the electron temperature pedestal to extend further inwards, as observed experimentally. Kinetic ballooning modes are found to be unstable mainly at the mid-pedestal of both types of discharges, with unstable trapped electron modes nearer the separatrix region. At electron wavelengths, ETG modes are found to be unstable from mid-pedestal outwards for eta(e,exp)~2.2, with higher growth rates for the lithiated discharge. Near the separatrix, the critical temperature gradient for driving ETG modes is reduced in the presence of lithium, reflecting the reduction of the lithiated density gradients observed experimentally. A preliminary linear study in the edge of non-lithiated discharges shows that the equilibrium shaping alters the electrostatic modes stability, found more unstable at high plasma shaping.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Improving Mathematical Exposition of an Industrial-Scale Linear Program

Industrial-scale models require considerable setup time; hence, once built, they are used in myriad ways to consider closely related cases. In practice, the code for these models frequently evolves without appropriate notational choices, largely as a result of the lengthy development time of, and the number of individuals contributing to, their formulation. This leads to inefficiencies and obfuscates model structures that might be leveraged to expedite solutions. In this paper, we advocate for an emerging literature on model formulation “best practices” and present the reformulation of a widely used industrial-scale linear program. The efficient mathematical expression of this linear program, used to plan capacity expansion in the energy sector, allows for greater transparency of model structures and enhanced ability to identify computational performance improvements, as well as a lucid interpretation of its solutions. This type of formulation is employed in several mathematical programming courses at our university as an example of the advantages of best practices; the model more broadly is used widely to inform policy in the U.S. energy sector.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

ALFA Non-linear Ocean Waves and PTO Control Strategy

Data from Advanced Laboratory and Field Arrays (ALFA) Non-linear Ocean Waves and Power Take-Off (PTO). Control Strategy project conducted at the O.H. Hinsdale Wave Research Laboratory (HWRL) at Oregon State University in 2019/2020. Contains two zip files (ALFANL.zip, ALFANL2.zip) from two phases of testing with data in engineering units and a report detailing the testing. This data collected by HWRL. A readme file in the docs folder explains the data collection and format. Six zip files (foswec-1 to foswec-6) data recorded by the Floating Oscillating Surge Wave Energy Converter (FOSWEC). Data organized by date and time and report describes data. A test report (ALFA Non-linear Ocean Waves Test Report.docx) details the experiments and data recording.

16 TIDAL AND WAVE POWER↗
DOE OSTI2018

Time and Frequency Domain Methods for Basis Selection in Random Linear Dynamical Systems

Polynomial chaos methods have been extensively used to analyze systems in uncertainty quantification. Furthermore, several approaches exist to determine a low-dimensional approximation (or sparse approximation) for some quantity of interest in a model, where just a few orthogonal basis polynomials are required. In this work, we consider linear dynamical systems consisting of ordinary differential equations with random variables. The aim of this paper is to explore methods for producing low-dimensional approximations of the quantity of interest further. We investigate two numerical techniques to compute a low-dimensional representation, which both fit the approximation to a set of samples in the time domain. On the one hand, a frequency domain analysis of a stochastic Galerkin system yields the selection of the basis polynomials. It follows a linear least squares problem. On the other hand, a sparse minimization yields the choice of the basis polynomials by information from the time domain only. An orthogonal matching pursuit produces an approximate solution of the minimization problem. Finally, we compare the two approaches using a test example from a mechanical application.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2013

Optimal Transmission Switching Using the IV-ACOPF Linearization

In this paper we seek to investigate the performance of transmission switching using the iterative linear program approximation to the Current Voltage AC Optimal Power Flow (IV-ACOPF). Several different methods of using this switching are investigated to find a method that addresses the MIP challenges and is both fast and accurate. We consider opening only one line, opening up to five lines, and progressively opening one line at a time. The linear method with switching is much faster than the nonlinear ACOPF and generally finds solutions within 1% of the nonlinear ACOPF.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2024

Computational study of tungsten and depleted uranium photoneutron targets for a 20 MeV electron linear accelerator

Neutron production can be realized with a high energy electron linear accelerator by using Bremsstrahlung and photoneutron converters. In this study, Monte Carlo N-Particle Code (MCNP) was used to evaluate potential photonuclear target designs for a high energy electron linear accelerator for applications such as neutron radiography and neutron resonance spectroscopy. A computational model was developed to inform a target design that would yield a high number of neutrons. It consists of a 20 MeV electron beam incident on a Bremsstrahlung target and a photonuclear target to generate neutrons. This computational model showed that a thickness of 0.75 inches for both tungsten and depleted uranium yields the most neutrons from photoneutron reactions. Saturation in the total number of generated neutrons was observed at over 0.75-inch thickness for both evaluated materials. Depleted uranium yielded approximately twice the number of neutrons overall compared to tungsten. The highest neutron surface flux for Depleted Uranium was 1.06 × 10-4 neutrons/cm2/source electron, and for Tungsten it was 5.12 × 10-5 neutrons/cm2/source electron. The optimal target design for this study’s application would consist of a 0.75 inch-thick block of depleted uranium with the length, width, and/or diameter varying dependent on application.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Novel Solver Algorithms for Nearly Singular Linear Systems Arising in Combustion Modelling

Direct Numerical Simulations of realistic combustion devices are extremely challenging due to the wide separation of scales in the simulation, for example an internal combustion (IC) engine chamber, and the flame thickness of a high-pressure flame. The PeleLMeX solver uses adaptive mesh refinement (AMR) to evolve multi-species reacting flows in the low Mach number limit at the Exascale and relies on an embedded boundary (EB) approach to represent complex geometries. In that framework, the EB geometries often give rise to very small cut-cells along the boundary, which translate into extreme ill-conditioning of the pressure-projection, with eigenvalues that span 15-16 orders of magnitude. In this talk, we focus on the case of a typical IC piston bowl geometry for which we present on a novel approach towards solving these nearly singular linear systems with ILU-based, C-AMG smoothers on massively parallel architectures. In particular, we use scaling and equilibration algorithms to handle the non-normality of the upper triangular factors. This enables us to approximate the highly sequential triangular solve algorithm, embedded in the AMG smoothing-solve phase, with Jacobi iterations. This approximation can be written as a convergent Neumann series whose terms are composed of highly parallel sparse matrix vector multiplications. The result is an algorithm that substantially decreases setup and solve time, compared to state-of-the-art, for these challenging linear systems.

combustion modelling↗
DOE OSTI2023

Sensor Anomaly Detection for Nuclear Reactor Systems Utilizing Linear Regression and K-Means Unsupervised Machine Learning

Nuclear reactors and related systems are becoming increasingly complex due to advancing technologies in next-generation power reactors. This increased complexity necessitates enhanced automation and data management capabilities. To successfully realize autonomous systems, methods must be developed to handle vast volumes of data and effectively distinguish anomalous data from noise and expected data. While impressive models utilizing digital twins and similar approaches are under development, here we propose a simplified model for analyzing fundamental methods and techniques. Initially, we created a general dataset by using initial data from PCTRAN in order to represent ideal steady-state conditions. We then inserted anomalies based on prevalent sensor anomaly types (e.g., point anomalies, linear drift, and downward deviations), along with unusual anomalies such as exponential drift and upward deviations. To detect anomalies, we developed a program that employs data partitioning and linear regression to preprocess and filter the anomalous data. A K-Means machine learning (ML) method was then applied to separate and count the data within the anomalous partition. The results from all datasets—apart from exponential growth—demonstrated positive outcomes, with each returning multiple instances of greaterthan-95% accuracy. We conducted further investigations using Idaho National Laboratory’s RAVEN software to perform a sensitivity analysis on the input variables (R 2 Tolerance, Slope Tolerance, and Window Size) and found that the output variables (Accuracy and Time) were most sensitive to the Window Size. Despite the promising results published, further development is required to effectively apply these methods to nuclear systems. Nevertheless, the strengths of this approach are evident and hold promise for future applications in the field.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2015

Design of Linear Accelerator (LINAC) tanks for proton therapy via Particle Swarm Optimization (PSO) and Genetic Algorithm (GA) approaches

A homemade computer code for designing a Side- Coupled Linear Accelerator (SCL) is written. It integrates a simplified model of SCL tanks with the Particle Swarm Optimization (PSO) algorithm. The computer code main aim is to obtain useful guidelines for the design of Linear Accelerator (LINAC) resonant cavities. The design procedure, assisted via the aforesaid approach seems very promising, allowing future improvements towards the optimization of actual accelerating geometries. (authors)

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Improved Evaluation of Large Network Matrices for Linear Power Flow Within Optimization Problems: Preprint

This work discusses methods for evaluating the Power Transfer Distribution Factor (PTDF) and Line Outage Distribution Factor (LODF) matrices by employing sparse linear algebra for large-scale computing applications. These matrices are critical in many power systems applications, such as the Unit Commitment Problem (UC), pre- and post-contingency power flow analysis, and transmission expansion. These matrices are typically dense, which means they require a significant amount of time and memory to be computed for large networks. However, by analyzing the structure of the matrices and their computation method, it is possible to use reduced memory methods based on sparse matrix operations. This paper shows that sparse linear algebra algorithms are faster and require less memory and time than traditional dense approaches. Additionally, we explore the effect of matrix sparsification by eliminating trailing digits on power flow calculations.

ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2026

Adaptive Control for Load-Following of Boiling Water Reactors Part I: Linear Systems and Fully-Observable Dynamics

Automation control is a key strategy to improve the economic competitiveness of nuclear power plants. Not only does it help reduce operational costs, but it also extends the value proposition of these plants to nontraditional markets, including unattended operations in remote villages and space. However, the dynamics of the operating environments of nuclear reactors are subject to changes over time, and there are no widely adopted methods to ensure that the automation strategy will remain effective over the extended durations required for these applications. Adaptive control is a discipline that offers the possibility to accommodate such changes online. However, it relies on mathematical assumptions that must be respected to ensure robustness and reliability. In this work, we derive an adaptive control formulation for linear systems in which all states are observable and apply it to an instance of load-follow operation for Boiling Water Reactors. We assumed uncertainty in two factors: the temperature coefficient and the control rod worth, both of which are affected over time by the evolution of the nuclear reactor core environment. With an arbitrary penalty factor of 5, we found that the mean absolute and integral time absolute errors can be reduced by more than 90%, underscoring the strength of adaptive control. To extend the application to more challenges, different uncertainties and load-follow trajectories, as well as new formulations that include non-linearity and partial observability, are currently being developed.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2015

Fracture propagation in Indiana Limestone interpreted via linear softening cohesive fracture model

In this study, we examine the use of a linear softening cohesive fracture model (LCFM) to predict single-trace fracture growth in short-rod (SR) and notched 3-point-bend (N3PB) test configurations in Indiana Limestone. The broad goal of this work is to (a) understand the underlying assumptions of LCFM and (b) use experimental similarities and deviations from the LCFM to understand the role of loading paths of tensile fracture propagation. Cohesive fracture models are being applied in prediction of structural and subsurface fracture propagation in geomaterials. They lump the inelastic processes occurring during fracture propagation into a thin zone between elastic subdomains. LCFM assumes that the cohesive zone initially deforms elastically to a maximum tensile stress ( σ max ) and then softens linearly from the crack opening width at σ max to zero stress at a critical crack opening width w 1 . Using commercial finite element software, we developed LCFMs for the SR and N3PB configurations. After fixing σ max with results from cylinder splitting tests and finding an initial Young's modulus ( E ) with unconfined compressive strength tests, we manually calibrate E and w 1 in the SR model against an envelope of experimental data. We apply the calibrated LCFM parameters in the N3PB geometry and compare the model against an envelope of N3PB experiments. For accurate simulation of fracture propagation, simulated off-crack stresses are high enough to require inclusion of damage. Different elastic moduli are needed in tension and compression. We hypothesize that the timing and location of shear versus extensional micromechanical failures control the qualitative macroscopic force-versus-displacement response in different tests. For accurate prediction, the LCFM requires a constant style of failure, which the SR configuration maintains until very late in deformation. The N3PB configuration does not maintain this constancy. To be broadly applicable between geometries and failure styles, the LCFM would require additional physics, possibly including elastoplastic damage in the bulk material and more complicated cohesive softening models.

58 GEOSCIENCES↗
DOE OSTI2018

Site-Selective Catalysis of a Multifunctional Linear Molecule: The Steric Hindrance of Metal–Organic Framework Channels

The site-selective reaction of a multifunctional linear molecule requires a suitable catalyst possessing both uniform narrow channel to limit the molecule rotation and a designed active site in the channel. Recently, nanoparticles (NPs) were incorporated in metal–organic frameworks (MOFs) with the tailorable porosity and ordered nanochannel, which makes these materials (NPs/MOFs) highly promising candidates as catalytic nanoreactors in the field of heterogeneous catalysis. Inspired by a “Gondola” sailing in narrow “Venetian Canal” without sufficient space for a U-turn, a simple heterogeneous catalyst based on NPs/MOFs is developed here that exhibits site-selectivity for the oxidation of diols by restricting the random rotation of the molecule (the “Gondola”) in the limited space of the MOF channel (the narrow “Venetian Canal”), thereby protecting the middle functional group via steric hindrance. This strategy is not limited to the oxidation of diols, but can be extended to the site-selective reaction of many similar multifunctional linear molecules, such as the reduction of alkadienes.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Birefringent Glass‐Engraved Quasi‐Linear Nanograting Metasurface Based on Self‐Organizing Process for Large Aperture High Power Laser Applications

All-glass metasurface “nanograting” structures that exhibit birefringence in the formed layer are reported. The key enabler of this work is ion beam processing at an angle sufficiently off-normal incidence, inducing self-assembly of a deposited metal layer into quasi-linear metallic features that can function as an etching mask. As a result, a fused silica metasurface, monolithic to the underlying substrate, is demonstrated at 375 nm wavelength to exhibit a phase delay angle of 30° between the principal axes. The capability of an angled etch mask replenishment process is also demonstrated for achieving deeper etch depth and for increasing the grating period, another first – to the best of the knowledge. This is the first display of a technology capable of fabricating glass-engraved near-linear grating structure with a feature-to-feature period as small as 118.6 nm. Furthermore, this technology has the potential to generate grating-like structures with periods as small as 12.4 nm, as demonstrated here with reactive ion beam processing assisted mask assembly. Furthermore, these structures are shown to have reflectivity < 0.4% across the wavelength band 350 nm – 1000 nm. Such a technology can enable laser-durable grating structures for the deep-UV and even down to soft X-ray wavelengths.

Ray, Nathan J. [Lawrence Livermore National Labora↗
DOE OSTI2024

Ab initio exploration of low–lying electronic states of linear and bent MNX + (M = Ca, Sr, Ba, Ra; X = O, S, Se, Te, Po) and their origins

High-level multireference and coupled cluster quantum calculations were employed to analyze low-lying electronic states of linear-MNX + and side-bonded-M[NX] + (M = Ca, Sr, Ba, Ra; X = O, S, Se, Te, Po) species. Their full potential energy curves (PECs), dissociation energies (D e s), geometric parameters, excitation energies (T e s), and harmonic vibrational frequencies (ω e s) are reported. The first three chemically bound electronic states of MNX + and M[NX] + are 3 Σ – , 1 Δ, 1 Σ + and 3 A", 1 A', 1 A", respectively. The 3 Σ – , 1 Δ, 1 Σ + of MNX + originate from the M + ( 2 D) + NX( 2 Π) fragments, whereas the 3 A", 1 A', 1A" states of M[NX] + dissociate to M + ( 2 S) + NX( 2 Π) as a result of avoided crossings. The MNX + and M[NX] + are real minima on the potential energy surface and their interconversions are possible. The M 2+ NX – /M 2+ [NX] – ionic structure is an accurate representation for their low-lying electronic states. The D e s of MNX + species were found to depend on the dipole moment (μ) of the corresponding NX ligands and a linear relationship between these two parameters was observed.

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

Topology property analysis and application of stable time-delay regions for linear multiple time-delay systems

This study examines the relationship between the topology of stable time-delay regions and the stability analysis of an linear multiple time delay system. To analyze the topology of stable time-delay regions, we construct a function with a value equal to zero for the time-delay points on the boundaries of stable time-delay regions. The function is continuous and differentiable in the whole defining field with a global minimum of zero so that we can locate the boundaries by minimizing the value of the function. Based on the topology analysis, we proposed a performance validation approach for controllers that are designed to stabilize the system using feedback signals with time delays. The method based on the topology analysis is simple and reliable so that can deal with a linear time-invariant system with high order and multiple time delays (more than 3). The example case study shows that the above method is reliable and enables ultra-low latency coordination and control for the future power grid with ubiquitous power electronics.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2017

Learning from Non-Linear Ecosystem Dynamics Is Vital for Achieving Land Degradation Neutrality

Land Degradation Neutrality is one of the Sustainable Development Goal targets, requiring on-going degradation to be balanced by restoration and sustainable land management. However, restoration and efforts to prevent degradation have often failed to deliver expected benefits,despite enormous investments. Better acknowledging the close relationships between climate, land management and non-linear ecosystem dynamics can help restoration activities to meet their intended goals, while supporting climate change adaptation and mitigation. This paper is the first to link ecological theory of non-linear ecosystem dynamics to Land Degradation Neutrality offering essential insights into appropriate timings, climate-induced windows of opportunities and risks and the financial viability of investments. These novel insights are pre-requisites for meaningful o and monitoring of progress towards Land Degradation Neutrality

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Highly Anisotropic Drude-Weight-Reduction and Enhanced Linear-Dichroism in van der Waals Weyl Semimetal T d -MoTe 2 with Coherent Interlayer Electronic Transport

Weyl semimetal (WSM) states can be achieved by breaking spatial-inversion symmetry or time reversal symmetry. However, the anisotropy of the energy reduction contributing to the emergence of WSM states has seldom been investigated by experiments. A van der Waals metal MoTe 2 exhibits a type-II WSM phase below the monoclinic-to-orthorhombic-phase-transition temperature T c ~ 250 K. Here, in this work, a combined linearly polarized optical-spectroscopy and electrical-transport study of MoTe 2 is reported at different temperatures. The Drude components in the a-axis, b-axis, and c-axis optical conductivity spectra, together with the metallic out-of-plane and in-plane electrical resistivities, indicate the coherent inter-layer and in-plane charge transports. Moreover, the Drude weight in σ 1a (ω), rather than the Drude weights in σ 1b (ω) and σ 1c (ω), decreases dramatically below T c , which exhibits a highly anisotropic decrease in its Drude weight and thus suggests a strongly anisotropic reduction of the electronic kinetic energy in the WSM phase. Furthermore, below T c , due to the in-plane anisotropic spectral-weight transfer from Drude component to high-energy region, the in-plane inter-band-absorption anisotropy increases remarkably ≈770 meV, and has the largest value (≈0.68) of normalized linear dichroism among the reported type-II WSMs. The work sheds light on seeking new WSMs and developing novel photonic devices based on WSMs.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗