Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linearity”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
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 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 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 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↗
DOE OSTI2020

Application of radiochromic gel dosimetry to commissioning of a megavoltage research linear accelerator for small–field animal irradiation studies

To develop and implement an efficient and accurate commissioning procedure for small-field static beam animal irradiation studies on an MV research linear accelerator (Linatron-M9) using radiochromic gel dosimetry. The research linear accelerator (Linatron-M9) is a 9 MV linac with a static fixed collimator opening of 5.08 cm diameter. Lead collimators were manually placed to create smaller fields of 2 × 2 cm 2 , 1 × 1 cm 2 , and 0.5 × 0.5 cm 2 . Relative dosimetry measurements were performed, including profiles, percent depth dose (PDD) curves, beam divergence, and relative output factors using various dosimetry tools, including a small volume ionization chamber (A14), GAFCHROMIC™ EBT3 film, and Clearview gel dosimeters. The gel dosimeter was used to provide a 3D volumetric reference of the irradiated fields. The Linatron profiles and relative output factors were extracted at a reference depth of 2 cm with the output factor measured relative to the 2 × 2 cm 2 reference field. Absolute dosimetry was performed using A14 ionization chamber measurements, which were verified using a national standards laboratory remote dosimetry service. Absolute dosimetry measurements were confirmed within 1.4% (k = 2, 95% confidence = 5%). The relative output factor of the small fields measured with films and gels agreed with a maximum relative percent error difference between the two methods of 1.1 % for the 1 × 1 cm 2 field and 4.3 % for the 0.5 × 0.5 cm 2 field. These relative errors were primarily due to the variability in the collimator positioning. The measured beam profiles demonstrated excellent agreement for beam size (measured as FWHM), within approximately 0.8 mm (or less). Film measurements were more accurate in the penumbra region due to the film’s finer resolution compared with the gel dosimeter. Following the van Dyk criteria, the PDD values of the film and gel measurements agree within 11% in the buildup region starting from 0.5 cm depth and within 2.6 % beyond maximum dose and into the fall-off region for depths up to 5 cm. The 2D beam profile isodose lines agree within 0.5 mm in all regions for the 0.5 × 0.5 cm 2 and the 1 × 1 cm 2 fields and within 1 mm for the larger field of 2 × 2 cm 2 . The 2D PDD curves agree within approximately 2% of the maximum in the typical therapy region (1–4 cm) for the 1 × 1 cm 2 and 2 × 2 cm 2 and within 5% for the 0.5 × 0.5 cm 2 field. This work provides a commissioning process to measure the beam characteristics of a fixed beam MV accelerator with detailed dosimetric evaluation for its implementation in megavoltage small animal irradiation studies. Radiochromic gel dosimeters are efficient small-field relative dosimetry tools providing 3D dose measurements allowing for full representation of dose, dosimeter misalignment corrections and high reproducibility with low inter-dosimeter variability. Altogether, radiochromic gels are valuable for fast, full relative dosimetry commissioning in comparison to films for application in high-energy small-field animal irradiation studies.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2024

Renormalon cancellation and linear power correction to threshold-like asymptotics of space-like parton correlators

Abstract In this paper, we show that the common hard kernel of double-log-type or threshold-type factorization for certain space-like parton correlators that arise in the context of lattice parton distributions, theheavy-light Sudakov hard kernel, has linear infrared (IR) renormalon. We explicitly demonstrate how this IR renormalon correlates with ultraviolet (UV) renormalons of next-to-leading power operators in two explicit examples: threshold asymptotics of space-like quark-bilinear coefficient functions and transverse momentum dependent (TMD) factorization of quasi wave function amplitude. Theoretically, the pattern of renormalon cancellation complies with general expectations to marginal asymptotics in the UV limit. Practically, this linear renormalon explains the slow convergence of imaginary parts observed in lattice extraction of the Collins-Soper kernel and signals the relevance of next-to-leading power contributions. Fully factorized, fully controlled threshold asymptotic expansion for space-like quark-bilinear coefficient functions in coordinate and moment space has also been proposed.

Physics↗
DOE OSTI2023

An analytic approach to quasinormal modes for coupled linear systems

Quasinormal modes describe the ringdown of compact objects deformed by small perturbations. In generic theories of gravity that extend General Relativity, the linearized dynamics of these perturbations is described by a system of coupled linear differential equations of second order. We first show, under general assumptions, that such a system can be brought to a Schrödinger-like form. We then devise an analytic approximation scheme to compute the spectrum of quasinormal modes. We validate our approach using a toy model with a controllable mixing parameter ε and showing that the analytic approximation for the fundamental mode agrees with the numerical computation when the approximation is justified. The accuracy of the analytic approximation is at the (sub-) percent level for the real part and at the level of a few percent for the imaginary part, even when ε is of order one. Our approximation scheme can be seen as an extension of the approach of Schutz and Will [1] to the case of coupled systems of equations, although our approach is not phrased in terms of a WKB analysis, and offers a new viewpoint even in the case of a single equation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Linearizing the BPS equations with vector and tensor multiplets

We analyse the BPS equations of N = (1, 0) supergravity theory in six dimensions coupled to a vector and tensor multiplet. We show how these BPS equations can be reduced to a set of linear differential equations. This system is triangular in that each layer of equations, while linear, is quadratically sourced by the solutions of the previous layers. We examine several explicit examples and discuss the construction of new families of microstate geometries. We expect that the result presented here will open up new branches of superstrata in which the momentum is encoded in a new class of charge carriers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Non-linear dynamics of jet quenching

We develop a comprehensive analytic framework for jet quenching in QCD media, based on a medium-induced parton cascade sourced by collinear virtual splittings. We show that the energy flow out of the jet cone, driven by turbulent gluon cascades, is governed by a non-linear rate equation that resums gluon splittings at arbitrary angles and is enhanced by the medium length, L. The solution of this equation sets the initial condition for a non-linear DGLAP-like evolution equation, which describes the collinear early vacuum cascade resolved by the medium at angles exceeding the medium resolution angle, θ c . For asymptotic jet energies, the medium-induced cascade displays an exponential behavior that generalizes the Poisson-like distribution of parton energy loss. This formulation enables the resummation of leading contributions in α s ln(1/R), and α s ln(R/θ c ), and powers of α s L. We briefly explore the limit of strong quenching, where analytic treatments are feasible, offering insights into the impact of parton cascades on jet quenching. These results provide guidance for future numerical simulations and analytical investigations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Higher harmonic non-linear flow modes of charged hadrons in Pb-Pb collisions at $\sqrt{s_{\mathrm{NN}}}$ = 5.02 TeV

Anisotropic flow coefficients, v n , non-linear flow mode coefficients, χ n,mk , and correlations among different symmetry planes, ρ n,mk are measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV. Results obtained with multi-particle correlations are reported for the transverse momentum interval 0.2 < p T < 5.0 GeV/c within the pseudorapidity interval 0.4 < |η| < 0.8 as a function of collision centrality. The v n coefficients and χ n,mk and ρ n,mk are presented up to the ninth and seventh harmonic order, respectively. Calculations suggest that the correlations measured in different symmetry planes and the non-linear flow mode coefficients are dependent on the shear and bulk viscosity to entropy ratios of the medium created in heavy-ion collisions. The comparison between these measurements and those at lower energies and calculations from hydrodynamic models places strong constraints on the initial conditions and transport properties of the system.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Development of Transition-Edge Sensor X-ray Microcalorimeter Linear Array for Compton Scattering and Energy Dispersive Diffraction Imaging

In this paper we present a strip transition-edge sensor microcalorimeter linear array detector developed for energy dispersive X-ray diffraction imaging and Compton scattering applications. The prototype detector is an array of 20 transition-edge sensors with absorbers in strip geometry arranged in a linear array. We discuss the fabrication steps needed to develop this array including Mo/Cu bilayer, Au electroplating, and proofof-principle fabrication of long strips of SiNx membranes. We demonstrate minimal unwanted effect of strip geometry on X-ray pulse response and showlinear relationship of 1/pulse height and pulse decay times with absorber length. For the absorber lengths studied, our preliminary measurements show energy resolutions of 40-180eV near 17 keV. Furthermore, we show that the heat flow to the cold bath is nearly independent of the absorber area and depends on the SiNx membrane geometry.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Efficient Trainability of Linear Optical Modules in Quantum Optical Neural Networks

The existence of “barren plateau landscapes” for generic discrete-variable quantum neural networks, which obstructs efficient gradient-based optimization of cost functions defined by global measurements, would be surprising in the case of generic linear optical modules in quantum optical neural networks due to the tunability of the intensity of continuous variable states and the relevant unitary group having exponentially smaller dimension. We demonstrate that coherent light in m modes can be generically compiled efficiently if the total intensity scales sublinearly with m, and extend this result to cost functions based on homodyne, heterodyne, or photon detection measurement statistics, and to noisy cost functions in the presence of attenuation. Here, we further demonstrate efficient trainability of m mode linear optical quantum circuits for variational mean field energy estimation of positive quadratic Hamiltonians for input states that do not have energy exponentially vanishing with m.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

The Cauchy Problem for Boltzmann Bi-linear Systems: The Mixing of Monatomic and Polyatomic Gases

Abstract From a unified vision of vector valued solutions in weighted Banach spaces, this paper establishes the existence and uniqueness for space homogeneous Boltzmann bi-linear systems with conservative collisional forms arising in complex gas dynamical structures. This broader vision is directly applied to dilute multi-component gas mixtures composed of both monatomic and polyatomic gases. Such models can be viewed as extensions of scalar Boltzmann binary elastic flows, as much as monatomic gas mixtures with disparate masses and single polyatomic gases, providing a unified approach for vector valued solutions in weighted Banach spaces. Novel aspects of this work include developing the extension of a general ODE theory in vector valued weighted Banach spaces, precise lower bounds for the collision frequency in terms of the weighted Banach norm, energy identities, angular or compact manifold averaging lemmas which provide coerciveness resulting into global in time stability, a new combinatorics estimate for p -binomial forms producing sharper estimates for the k -moments of bi-linear collisional forms. These techniques enable the Cauchy problem improvement that resolves the model with initial data corresponding to strictly positive and bounded initial vector valued mass and total energy, in addition to only a $$2^+$$ 2 + moment determined by the hard potential rates discrepancy, a result comparable in generality to the classical Cauchy theory of the scalar homogeneous Boltzmann equation.

Physics↗