Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “method of multipliers”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Data for Monitoring Fracture Saturation with Internal Transportable Seismic Sources and Twin Neural Networks

Seismic coda-wave analysis is a well-developed method for detecting subtle physical changes in complex media by measuring arrival times in the late-arriving energy from multiply-scattered or reflected waves. However, a challenge arises when multiply-scattered waves are not sufficiently separated in time from the direct arrivals to provide a clear coda wave train. Additional complications for monitoring changes in fracture systems arise when the signals originate from unsynchronized internal sources, such as natural induced seismicity, from acoustic emission, or from transportable intra-fracture sources (chattering dust), that generate uncontrolled signals that vary in arrival time, amplitude and frequency content. Here, we use a twin neural network (also known as a Siamese neural network) for dimensionality reduction to analyze signals from chattering dust to classify the fluid saturation state of a synthetic fracture system. The twin neural network with shared weights generates a low-dimensional representation of the data input by minimizing contrastive loss, serving as the input to a multiclass classifier that accurately classifies whether multiple fractures in a fracture system are fully saturated or partially saturated, or whether a change in saturation has occurred in different fractures in the system. These results show that information buried in unresolved codas from uncontrolled sources can be extracted using machine learning to monitor the evolution of fracture systems caused by physical and chemical processes even when the scattered and direct wave fields overlap.

58 GEOSCIENCES↗
DOE OSTI2021

Derivation of the Future Time Equation for Analog, Non-Multiplying Monte Carlo Simulation

The expected computational time required to simulate a particle from a point in phase space through a Monte Carlo history, termed the expected future time, is found by solving the Future Time Equation (FTE). The expected future time may be useful when generating variance reduction parameters for a Monte Carlo simulation with a method such as Consistent Adjoint Driven Importance Sampling (CADIS). This report presents a detailed derivation of the Future Time Probability Density Function (FTPDF) and FTE for neutral particle Monte Carlo transport to aid future researchers. For simplicity, this derivation only considers analog transport in non-multiplying media.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Discovering strongly lensed quasar candidates with catalogue-based methods from DESI Legacy Surveys

The Hubble tension, revealed by a ~5σ discrepancy between measurements of the Hubble-Lemaitre constant among observations of the early and local Universe, is one of the most significant problems in modern cosmology. In order to better understand the origin of this mismatch, independent techniques to measure H 0 , such as strong lensing time delays, are required. Notably, the sample size of such systems is key to minimising the statistical uncertainties and cosmic variance, which can be improved by exploring the datasets of large-scale sky surveys such as Dark Energy Spectroscopic Instrument (DESI). We identify possible strong lensing time-delay systems within DESI by selecting candidate multiply imaged lensed quasars from a catalogue of 24 440 816 candidate QSOs contained in the ninth data release of the DESI Legacy Imaging Surveys (DESI-LS). Using a friend-of-friends-like algorithm on spatial co-ordinates, our method generates an initial list of compact quasar groups. This list is subsequently filtered using a measure of the similarity of colours among a group’s members and the likelihood that they are quasars. A visual inspection finally selects candidate strong lensing systems based on the spatial configuration of the group members. We identified 620 new candidate multiply imaged lensed quasars (101 grade-A, 214 grade-B, 305 grade-C). This number excludes 53 known spectroscopically confirmed systems and existing candidate systems identified in other similar catalogues. When available, these new candidates will be further checked by combining the spectroscopic and photometric data from DESI.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

A priori error analysis of high-order LL* (FOSLL*) finite element methods

A number of non-standard finite element methods have been proposed in recent years, each of which derives from a specific class of PDE-constrained norm minimization problems. The most notable examples are LL* methods. In this work, we argue that all high-order methods in this class should be expected to deliver substandard uniform h-refinement convergence rates. In fact, one may not even see rates proportional to the polynomial order p > 1 when the exact solution is a constant function. Here, we show that the convergence rate is limited by the regularity of an extraneous Lagrange multiplier variable which naturally appears via a saddle-point analysis. In turn, limited convergence rates appear because the regularity of this Lagrange multiplier is determined, in part, by the geometry of the domain. Numerical experiments support our conclusions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Timing resolution and rate capability of Photonis miniPlanacon XPM85212/A1-S MCP-PMT

In this paper we report on a set of experiments that focus on measuring the timing performance of the Photonis miniPlanacon XPM85212/A1-S microchannel plate photomultiplier tubes along with the detector response at high event rates. The detector has a single photoelectron timing resolution of 30 ps at low rates ( 10 kHz ). We show that both the gain and the timing start to deteriorate around 1 MHz with 20 photoelectrons per channel at 1 0 4 target gain, giving the current per unit of area limits of 1.38 μ A / cm 2 and 2.38 μ A / cm 2 for the two tested devices. Photo-multiplier tubes with lower micro-channel plate resistivity exhibit a better rate capability due to the faster replenishment of the charge inside the microchannels. As we demonstrate, another method to improve the rate capability is to operate the photon detector tubes at a lower gain by reducing the supply voltage and compensating for the loss of amplitude by an additional amplification stage at the expense of the achieved timing resolution. The tube active area over which the beam is spread also plays a role since the same amount of light yields a better rate capability when spread over a larger area. The tubes exposed to high event rates showed an uncharacteristically slow recovery back to its initial gain with 80 % of the initial gain recovered only after several tens of minutes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Variational, stable, and self-consistent coupling of 3D electromagnetics to 1D transmission lines in the time domain

This work presents a new multiscale method for coupling the 3D Maxwell's equations to the 1D telegrapher's equations. While Maxwell's equations are appropriate for modeling complex electromagnetics in arbitrary-geometry domains, simulation cost for many applications (e.g. pulsed power) can be dramatically reduced by representing less complex transmission line regions of the domain with a 1D model. By assuming a transverse electromagnetic (TEM) ansatz for the solution in a transmission line region, we reduce the Maxwell's equations to the telegrapher's equations. Here, we propose a self-consistent finite element formulation of the fully coupled system that uses boundary integrals to couple between the 3D and 1D domains and supports arbitrary unstructured 3D meshes. Additionally, by using a Lagrange multiplier to enforce continuity at the coupling interface, we allow for an absorbing boundary condition to also be applied to non-TEM modes on this boundary. We demonstrate that this feature reduces non-physical reflection and ringing of non-TEM modes off of the coupling boundary. By employing implicit time integration, we ensure a stable coupling, and we introduce an efficient method for solving the resulting linear systems. We demonstrate the accuracy of the new method on two verification problems, a transient O-wave in a rectilinear prism and a steady-state problem in a coaxial geometry, and show the efficiency and weak scalability of our implementation on a cold test of the Z-machine MITL and post-hole convolute.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

A sharp interface Lagrangian-Eulerian method for flexible-body fluid-structure interaction

This paper introduces a sharp-interface approach to simulating fluid-structure interaction (FSI) involving flexible bodies described by general nonlinear material models and across a broad range of mass density ratios. This new flexible-body immersed Lagrangian-Eulerian (ILE) scheme extends our prior work on integrating partitioned and immersed approaches to rigid-body FSI. Our numerical approach incorporates the geometrical and domain solution flexibility of the immersed boundary (IB) method with an accuracy comparable to body-fitted approaches that sharply resolve flows and stresses up to the fluid-structure interface. Unlike many IB methods, our ILE formulation uses distinct momentum equations for the fluid and solid subregions with a Dirichlet-Neumann coupling strategy that connects fluid and solid subproblems through simple interface conditions. As in earlier work, we use approximate Lagrange multiplier forces to treat the kinematic interface conditions along the fluid-structure interface. This penalty approach simplifies the linear solvers needed by our formulation by introducing two representations of the fluid-structure interface, one that moves with the fluid and another that moves with the structure, that are connected by stiff springs. This approach also enables the use of multi-rate time stepping, which allows us to use different time step sizes for the fluid and structure subproblems. Our fluid solver relies on an immersed interface method (IIM) for discrete surfaces to impose stress jump conditions along complex interfaces while enabling the use of fast structured-grid solvers for the incompressible Navier-Stokes equations. The dynamics of the volumetric structural mesh are determined using a standard finite element approach to large-deformation nonlinear elasticity via a nearly incompressible solid mechanics formulation. This formulation also readily accommodates compressible structures with a constant total volume, and it can handle fully compressible solid structures for cases in which at least part of the solid boundary does not contact the incompressible fluid. Selected grid convergence studies demonstrate second-order convergence in volume conservation and in the pointwise discrepancies between corresponding positions of the two interface representations as well as between first and second-order convergence in the structural displacements. The time stepping scheme is also demonstrated to yield second-order convergence. To assess and validate the robustness and accuracy of the new algorithm, comparisons are made with computational and experimental FSI benchmarks. Test cases include both smooth and sharp geometries in various flow conditions. Furthermore, we also demonstrate the capabilities of this methodology by applying it to model the transport and capture of a geometrically realistic, deformable blood clot in an inferior vena cava filter.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Analysis of Lobe Power Calculator and Indication System with Physics and Cycle Based Models

The Advanced Test Reactor (ATR) at INL measures reactor power through two methods, thermal and Nitrogen-16 (N-16) activity. Water power calculator (WPC) is a thermal power system that measures flow and temperature to determine the thermal quadrant powers. The N-16 system utilizes a beta chamber detector that outputs reactivity levels to calculate lobe power through an algorithm called lobe power calculation and indication system (LPCIS). The LPCIS utilizes the N-16 system and the WPC system to determine reactor core power levels. The WPC provides accurate calculations of quadrant and total reactor thermal power. With the use of WPC measurements, thermal-to-N-16 (T2N) power ratios are produced to determine if the two indication systems agree on core power. Relative magnitude equations are used to utilize N-16 coefficients and multipliers to improve the indications of the LPCIS. This is crucial for maintaining safety limits. Currently, the LPCIS system uses linear equations and matrices to calculate lobe power through multipliers and coefficients. Advancements in technology and system upgrades have assisted system engineers at INL with the objective to reach a more dynamic system. In return the system demonstrates an increase in the accuracy of power reading while maintaining safety margins. The new proposed coefficient and multiple method implements a cycle specific, and a physics-based model intended for changing coefficients during operation. This calibration experiment focuses on power splitting, outer shim and neck shim, as well as fuel burning into the RDAS weighting factors. Results show that the Physics learning method yields a smaller error margin inside of the desired power range for the data set 166-A. This proved true for both constrained and unconstrained testing. This is most likely due to the physics data fitting approach that resulted in favorable coefficients and multipliers for 166-A. Continuing to improve the physics-based model will help improve the power accuracy of the LPCIS system.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

Recent Improvements to the NWChem COSMO Module

This work presents recent developments in the COSMO solvation model implementation in NWChem. A new cavity construction approach, based on the solvent-excluding surface (SES) and utilizing the well-established GEPOL algorithm, has been introduced. Additionally, a straightforward procedure to merge surface segments that are too close─often a source of numerical artifacts─has been implemented. The available methods for correcting outlying surface charges have also been reviewed and improved. To validate the new implementation, we computed dielectric solvation energies for a chemically diverse set of approximately 100 molecules, including neutral species, small ions, and common ionic liquid components. Results were compared to those from GAMESS using the double-cavity method as a reference. Although the double-cavity approach can be regarded as more accurate, the simpler correction schemes available in NWChem─based on scaling factors or Lagrange multipliers─can achieve excellent agreement if the potential is also properly corrected, with mean unsigned deviations of around 0.15 kcal/mol. Predictions of typical vapor–liquid and liquid–liquid equilibria using a COSMO-SAC variant based on NWChem also showed very promising results.

Cavities↗
DOE OSTI2021

In crystallo organometallic chemistry

X-ray crystallography is an invaluable tool in design and development of organometallic catalysis, but application typically requires species to display sufficiently high solution concentrations and lifetimes for single crystalline samples to be obtained. In crystallo organometallic chemistry relies on chemical reactions that proceed within the single-crystal environment to access crystalline samples of reactive organometallic fragments that are unavailable by alternate means. Furthermore, this highlight describes approaches to in crystallo organometallic chemistry including (a) solid–gas reactions between transition metal complexes in molecular crystals and diffusing small molecules, (b) reactions of organometallic complexes within the extended lattices of metal–organic frameworks (MOFs), and (c) intracrystalline photochemical transformations to generate reactive organometallic fragments. Application of these methods has enabled characterization of catalytically important transient species, including σ-alkane adducts of transition metals, metal alkyl intermediates implicated in metal-catalyzed carbonylations, and reactive M–L multiply bonded species involved in C–H functionalization chemistry. Opportunities and challenges for in crystallo organometallic chemistry are discussed.

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

Efficient Scalable Contact Network Generation from Population Data

Modeling the contacts among a population is critical to understanding the dynamics of a disease outbreak. Contact networks, where nodes are individuals and edges are contacts among them, are used to represent these complex individual-level interactions. In this work, we are given the daily activity schedules of an urban population that represent the activity location and time of individuals in a population during a single twenty four hour period over multiple days. Using collocation to determine contact between individuals, our goal is to extract hourly contact networks from large-scale activity data. We improve upon the existing adjacency matrix-based method by implementing our custom sparse matrix multiplication algorithm. Starting with a Python implementation, we achieve a 1600x speed up in the computation with a fast custom designed sparse matrix multiplier algorithm implemented in the C++ language. This work is central to future parallel designs of the problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Unconventional compute methods and future challenges for superconducting digital computing

Superconducting digital computing (SDC) based on Josephson junctions (JJs) offers significant potential for enhancing compute throughput and reducing energy consumption compared to conventional room-temperature CMOS-based approaches. Current superconducting logic families exhibit diverse characteristics in clocking strategies, power management, and information encoding techniques. This paper reviews recent advancements in unconventional computing methods specifically designed for superconducting digital circuits, emphasizing temporal computing and pulse-train representations. Notable techniques include race logic (RL), temporal pulse train computing (U-SFQ), and temporal multipliers, each offering unique performance and area advantages suited to superconducting implementations. Additionally, this paper reviews innovations in superconducting coarse-grain reconfigurable architectures (CGRA), superconducting-specific on-chip communication architectures, cryogenic sensor interfaces, and quantum computing control electronics. Finally, we highlight research challenges that should be addressed to facilitate the widespread adoption of superconducting digital computing.

EDA tools↗
DOE OSTI2020

A High Precision Measurement of the Proton Charge Radius at JLab

The elastic electron-proton (e-p) scattering and the spectroscopy of hydrogen atoms are the two traditional methods to determine the proton charge radius (r p ). In 2010, a new method using muonic hydrogen (µH) 1 spectroscopy reported a r p result that was nearly ten times more precise but significantly smaller than the values from the compilation of all previous r p measurements. This discrepancy is often referred to as the "proton charge radius puzzle". In order to investigate the puzzle, the PRad experiment (E12-11-106) was first proposed in 2011 and performed in 2016 in Hall B at the Thomas Jefferson National Accelerator Facility, with both 1.1 and 2.2 GeV electron beams. The experiment measured the e-p elastic scattering cross sections in an unprecedented low values of momentum transfer squared region (Q 2 = 2.1× 10 -4 -0.06 (GeV/c) 2 ), with a sub-percent precision. The PRad experiment utilized a calorimetric method that was magnetic-spectrometer-free. Its detector setup included a large acceptance and high resolution calorimeter (HyCal), and two large-area, high spatial-resolution Gas Electron Multiplier (GEM) detectors. To have a better control over the systematic uncertainties, the absolute e-p elastic scattering cross section was normalized to that of the well-known Møller scattering process, which was measured simultaneously during the experiment. For each beam energy, all data with different Q 2 were collected simultaneously with the same detector setup, therefore sharing the same integrated luminosity. The windowless H 2 gas-flow target utilized in the experiment largely removed a typical background source, the target cell windows. The proton charge radius was determined as r p = 0.831±0.007 stat. ±0.012 p fm, which is smaller than the average r p from previous e-p elastic scattering experiments, but in agreement with the µH spectroscopic results within the experimental uncertainties.

Source record↗
DOE OSTI2021

Distributed control of electric power grids

A method for determining power output levels of a plurality of nodes in an electric power system includes receiving, at a first node of the plurality of nodes, voltage information and multipliers of all neighboring nodes of the first node within the electric power system, determining, by the first node, a local power generation and a local voltage using the voltage information and the multipliers of the neighboring nodes and distributing the local power generation and the local voltage to the neighboring nodes, determining, by the first node, an estimated voltage of each of the neighboring nodes and distributing the estimated voltage to each of the neighboring nodes, and updating, by the first node, a local multiplier using the voltage information received from the neighboring nodes and the estimated voltage of each of the neighboring nodes determined by the node.

Ghosh, Soumyadip↗
DOE OSTI2021

Improved motional Stark effect signal processing using fast Fourier transform spectral analysis

A Fast Fourier Transform (FFT) based method has been developed, which improves the frequency response of the Motional Stark Effect (MSE) system by about a factor of 10 over the conventional analog lock-in method. The method uses fits to rigorously derived analytic expressions for the FFT spectral components of the MSE signal to accurately obtain the amplitudes and phases of the 2f1 and 2f2 photo-elastic modulator (PEM) frequencies that encode the polarization angle. Since no frequency filtering is used in the FFT method, the frequency response is limited by fundamental measurement properties: the frequency response of the detector, photon statistics, sample rate, and the ability to resolve the spectral components. In contrast, the frequency response of the analog lock-in is limited by a low pass filter with a cutoff of around 500 Hz. In the case of the DIII-D MSE system, the output of the photo-multiplier tube detector was sampled at 500 kHz and FFTs with as few as 100 points were used to obtain the amplitudes of the 2f1 and 2f2 PEM frequency components. This corresponds to a frequency response of 5 kHz, about ten times faster than the analog lock-in amplifier system. Details of the FFT method will be presented and compared to those of the analog lock-in system.

Makowski, M. A.↗
DOE OSTI2022

Numerical algorithms for water waves with background flow over obstacles and topography

Abstract We present two accurate and efficient algorithms for solving the incompressible, irrotational Euler equations with a free surface in two dimensions with background flow over a periodic, multiply connected fluid domain that includes stationary obstacles and variable bottom topography. One approach is formulated in terms of the surface velocity potential while the other evolves the vortex sheet strength. Both methods employ layer potentials in the form of periodized Cauchy integrals to compute the normal velocity of the free surface, are compatible with arbitrary parameterizations of the free surface and boundaries, and allow for circulation around each obstacle, which leads to multiple-valued velocity potentials but single-valued stream functions. We prove that the resulting second-kind Fredholm integral equations are invertible, possibly after a physically motivated finite-rank correction. In an angle-arclength setting, we show how to avoid curve reconstruction errors that are incompatible with spatial periodicity. We use the proposed methods to study gravity-capillary waves generated by flow around several elliptical obstacles above a flat or variable bottom boundary. In each case, the free surface eventually self-intersects in a splash singularity or collides with a boundary. We also show how to evaluate the velocity and pressure with spectral accuracy throughout the fluid, including near the free surface and solid boundaries. To assess the accuracy of the time evolution, we monitor energy conservation and the decay of Fourier modes and compare the numerical results of the two methods to each other. We implement several solvers for the discretized linear systems and compare their performance. The fastest approach employs a graphics processing unit (GPU) to construct the matrices and carry out iterations of the generalized minimal residual method (GMRES).

Ambrose, David M.↗
DOE OSTI2022

Multielectron Organic Redoxmers for Energy-Dense Redox Flow Batteries

Redox flow battery is a highly promising stationary energy storage method but the limited energy density and high chemical cost are among the main barriers for commercialization. Multi-electron organic redoxmers represent a family of structurally tailorable candidates that can achieve multiplied energy density with decreased materials consumption, potentially resulting in a viable solution to address these challenges. In this work, the recent development of organic molecules with reversible multi-redox activities in both aqueous and nonaqueous electrolytes is reviewed. The major focus is on the fundamental correlation between the chemical structures and the functional properties of reported multi-electron organic molecules. Valuable insights are offered on rational structural design strategies for improving the relevant physico- and electro-chemical properties. Finally, the current challenges are discussed to suggest future research needs along the avenue of using the multi-electron approach to achieve energy-dense, stable, cost-effective redox flow batteries.

25 ENERGY STORAGE↗
DOE OSTI2022

Apparatus and method for generating medical isotopes

An apparatus for generating medical isotopes provides an annular fissile solution vessel surrounding a neutron generator. The annular fissile solution vessel provides for good capture of the emitted neutrons and a geometry that provides enhanced stability in an aqueous reactor. A neutron multiplier and/or a neutron moderator may be used to improve the efficiency and control the criticality of the reaction in the annular fissile solution vessel.

Piefer, Gregory Richard↗