SEARCH · Engineering Papers

Results for “fixed point method”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Closed-Form Solutions for the Equations of Motion of the Heavy Symmetrical Top with One Point Fixed

The Equations of Motion (EOM) for the Heavy Symmetrical Top with One Point Fixed are highly non-linear. The literature describes the numerical methods that are used to resolve this Classical system including modern tools i.e. the Runge-Kutta Fourth Order method. It is more difficult to derivate closed-form solutions for the EOM and as mentioned in the literature it is not always possible to find the close-form solution for all the EOM. Fortunately, there are a few examples available that will serve as a guide to move further on in this topic. It is the purpose of this paper to find a methodology that will produce the solutions for a given subset of EOM that fulfill certain requisites. The report is organized as follows: it starts with a very short summary of the literature available on this topic and quickly follows into the derivation of the EOM using the Euler-Lagrange method. The Routhian will be used to reduce the size of the expression. It continues with the formulation of the Classical cubic function (f(u)) through a novel process. The roots of f(u) are of the outmost importance to be able to find the EOM closed-form solution, and once the final roots are selected the general method that will produce the closed-form solutions is presented. Two sets of examples are included to show the validity of the process and comparisons of the results from the closed-form solutions vs. the numerical results for these examples are shown.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Highly accelerated life testing (HALT): A review from a statistical perspective

Despite its use in one form or another for at least four decades, HALT and related techniques [e.g., highly accelerated-stress screening (HASS) and stress audits (HASA)] are not well understood within the statistical community and remain controversial. This largely reflects a conflict in motivation between engineers, testing under harsh conditions to discover and eliminate failure modes, and statisticians, taking a more cautious approach to develop quantitative estimates of parameters such as mean time between failures (MTBF). Here, this review article will clarify HALT concepts and methods and explain where it fits within the universe of methods that involve the application of accelerating factors to compress the time required to evaluate or enhance product reliability. A major distinction is between methods such as HALT, a high-stress test-analyze-fix-test iterative process directed at improving reliability by discovering and fixing weak points in a design, and quantitative accelerated life testing (QALT), whose goal is the estimation of product life for a fixed design. We discuss methods such as physics of failure that offer some hope of bridging the gap between the qualitative nature of HALT, and purely quantitative statistical methods. We present a variety of engineering applications of HALT including metal fatigue, piping and pressure vessels, structural damage, radiation damage, and rotating machinery. We also discuss potential synergies between HALT and QALT, such as rapid identification, through HALT, of failure modes requiring quantitative analysis. For further study, extensive references to the applicable literature are provided as well as an appendix that describes related methods.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

On the Development of Order and Interfaces during the Growth of Ultrathin La 2 CuO 4 Films by Molecular Beam Epitaxy

While the atomic structure of interfaces in complex oxide heterostructures created by epitaxial growth has been investigated extensively, few studies have been conducted on how interfaces form and restructure at the initial stage of film growth. The dynamic aspects of the growth behavior can strongly influence the final interfacial atomic structure, which may lead to the emergence of interface-specific properties, such as the rise of interfacial superconductivity between certain Ruddlesden-Popper oxide materials. In this report, the structural development of La2CuO 4 thin films grown by molecular beam epitaxy on a LaSrAlO 4 substrate is investigated by X-ray diffraction measurements with rapid scans over a volume of reciprocal space. This method provides far more detailed information on interface formation than traditional fixed-point measurements. The results show that the atomic structure of the interface becomes fully established after just a single unit cell of growth. Interestingly, restructuring continues to occur within the topmost half to one unit cell of the film during the deposition process. However, diffraction intensity oscillations from both reflection high-energy electron and X-ray measurements stabilize only after the growth of two unit cells, indicating that the growth front morphology continues to evolve until the start of the third unit cell. This multimodal investigation provides insights into the atomic processes taking place during layered oxide interface formation, including the dynamical rearrangement of LaO and CuO 2 layers. Such information is not only relevant to the engineering and optimization of functional layer structures but can also be critical for ultrathin films.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Calculation of permanent magnet arrangements for stellarators: a linear least-squares method

A problem arising in several engineering areas is to design magnets outside a volume that produce a desired magnetic field inside it. One instance of this problem is stellarator design, where it has recently been shown that permanent magnets can provide the required shaping of the magnetic field. Here we demonstrate a robust and efficient algorithm REGCOIL_PM to calculate the spatial distribution of these permanent magnets. The procedure involves a small number of fixed-point iterations, with a linear least-squares problem solved at each step. Furthermore, the method exploits the Biot–Savart Law's exact linearity in magnetization density and approximate linearity in magnet size, for magnets far from the target region. No constraint is placed on the direction of magnetization, so Halbach solutions are found naturally, and the magnitude of the magnetization can be made uniformly equal to a target value.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Fast explicit solutions for neutrino-electron scattering: Explicit asymptotic methods

Here, we present results of explicit asymptotic approximations applied to neutrino-electron scattering in a representative model of neutrino population evolution under conditions characteristic of core-collapse supernova explosions or binary neutron star mergers. It is shown that this approach provides stable solutions of these stiff systems of equations, with accuracy and time stepping comparable to that for standard implicit treatments such as backward Euler, fixed point iteration, and Anderson-accelerated fixed point iteration. Because each time step can be computed more rapidly with the explicit asymptotic approximation than with implicit methods, this suggests that algebraically stabilized explicit integration methods could be used to compute neutrino evolution coupled to hydrodynamics more efficiently in stellar explosions and mergers than the methods currently in use.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Hydrogen Detection Strategies to Support H2@SCALE - The NREL Sensor Laboratory

Hydrogen represents a major pathway to decarbonize and stabilize the national and international energy industry and select manufacturing markets. To facilitate the development of hydrogen markets, the US Department of Energy initiated H2@Scale to bring together stakeholders to advance affordable hydrogen production, transport, storage, and utilization to increase revenue opportunities in multiple energy sectors. One major impediment to hydrogen implementation is cost. To expedite the use of hydrogen in energy and other markets, the United States announced in 2021 the Hydrogen Shot, which seeks to reduce the cost of clean hydrogen by 80% to $1 per 1 kilogram in 1 decade ("1 1 1"). As the cost of hydrogen drops, new applications will emerge that will require unique configurations of existing equipment and infrastructure, and eventually lead to advances in the generation and utilization of hydrogen. As the hydrogen economy expands, sensors and detection methods will need to adapt to changing infrastructure demands to address the primary targets of health & safety, emissions monitoring, and process control. The NREL Sensor Laboratory is playing a pivotal role in advancing the use of hydrogen sensors and detection methodologies in each of these categories to support DOE's mission for safe and efficient utilization in emerging markets. Health & safety monitors are required to ensure that operators and facilities can react to unintended hydrogen releases, either as GH2, LH2, or as a constituent of blends (e.g., natural gas or ammonia). Current detection methodologies focus on safety applications to detect near its lower flammable limit (4 vol %), and typically include point sensors in applications such as fixed or mobile detectors (e.g., personal gas monitors). Methodologies amenable for area detection include acoustic, emerging optical imaging methods, and flame detectors. Comparable detection strategies can be utilized for emissions monitoring and quantization, however few methods can simultaneously cover both low (emissions) and high (health & safety) levels. Deployment of emission level detectors will be required to 1) reduce product loss through small but potentially significant leaks from an environmental or cost perspective, 2) reduce downtime of high demand systems by early identification of eminent system failures (leaks through pump or compressor seals indicative of impending failure), and 3) address potential emission monitoring requirements that may be set by regulating bodies. The first two points should be adopted by industry to reduce the cost-of-goods-sold. The third main category for hydrogen detection relates to process control and may be advantageous for many existing applications. Two main applications are emerging. For example, the purity requirements for hydrogen that is dispensed from refueling systems for hydrogen fuel cell electric vehicles (FCEV) is rigorously regulated by the Standard SAE J2719, which prescribes maximum allowable levels of multiple impurities in the hydrogen fuel and must be verified by a regulatory body. Hydrogen contaminant detectors (HCD) integrated to the fueling station can assure this compliance. HCDs must be able operate in 100% H2 backgrounds and be able to distinguish between multiple contaminants at low ppm to low ppb levels. Secondly, as a strategy to decarbonize the natural gas grid, there are proposals to blend hydrogen with natural gas. This blending will affect transport applications (pipeline infrastructure), stationary combustion systems (turbines), and consumer and commercial appliances. In the short-term, hydrogen levels up to 20% are proposed. Variations in the hydrogen level can have dramatic impact on the combustion process and on the potential response of safety sensors. These mixtures may be regulated so that the concentration at a delivery point must be monitored with high precision. However, routine maintenance may introduce background gases such as ambient air (with water) or maintenance gases (introduced with welding processes or adhesive outgassing.) Therefore, the detection methodology must be robust enough to recover or respond to various contaminants. Several reviews can be found in literature addressing sensing and detection technologies, including their limitations and applications. However, for most applications, limitations can be alleviated by combining various detection techniques either through system integration or implementation of machine learning methods (artificial intelligence). In this presentation, we will discuss several applications, highlight their current approach for hydrogen detection, and suggest detection strategies to supplement their limitations.

ENERGY STORAGE,HYDROGEN↗

DOE OSTI2025

A Novel Approach for Computing Rigid Body Motion Using Linear Accelerations

Here, a novel approach is presented for computing general rigid body motion based on a few known linear accelerations. This method utilizes linear acceleration data obtained from three distinct points on the body, all within a body-fixed reference frame. The only requirement is that the three chosen points must not be collinear. A system of differential-algebraic equations is derived, combining principles of rigid body kinematics with theory of the rotation group SO(3). These equations provide a framework for numerically computing various motion parameters, including angular velocity, angular acceleration, body orientation, velocity field, acceleration field, and displacement field. By numerically solving this system of equations, we can fully characterize rigid body motion in three-dimensional space. A numerical example is provided to demonstrate the practical implementation and efficacy of the proposed technique, illustrating its potential for accurate motion computation in various applications.

42 ENGINEERING↗

DOE OSTI2024

Neutron transport methods for multiphysics heterogeneous reactor core simulation in Griffin

Griffin is a reactor physics application based on the Multiphysics Object-Oriented Simulation Environment (MOOSE). This work discloses the methods, algorithms, and implementation for simulating heterogeneous reactor dynamics models. Griffin utilizes a discontinuous finite-element method with discrete ordinates (DFEM-S ) to discretize the field variable of the multigroup neutron transport equation. Multiphysics feedback is handled using two-step tabulated cross-section methodology. Feedback quantities are evaluated using the MOOSE-MultiApp system to couple various engineering phenomena, such as heat conduction and thermal fluids. The multiphysics DFEM-S system is solved using fixed-point iteration with a fully asynchronous parallel sweeper, unstructured coarse-mesh finite difference acceleration, and a multi-timescale improved quasi-static method scheme. The implementation is applied to a multiphysics microreactor model, with two transients: one initiated by a single heat-pipe failure and another by control drum rotation. Importantly, these examples demonstrate the ability of Griffin to tractably solve the neutron transport equation considering seven independent variables and feedback.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

High precision tests of QCD without scale or scheme ambiguities: The 40th anniversary of the Brodsky–Lepage–Mackenzie method

A key issue in making precise predictions in QCD is the uncertainty in setting the renormalization scale μ r and thus determining the correct values of the QCD running coupling α s (μ r ) at each order in the perturbative expansion of a QCD observable. It has often been conventional to simply set the renormalization scale to the typical scale of the process Q and vary it in the range μ r $\in$ [Q/2, 2Q] in order to estimate the theoretical error. This is the practice of Conventional Scale Setting (CSS). The resulting CSS prediction will however depend on the theorist’s choice of renormalization scheme and the resulting pQCD series will diverge factorially. It will also disagree with renormalization scale setting used in QED and electroweak theory thus precluding grand unification. A solution to the renormalization scale-setting problem is offered by the Principle of Maximum Conformality (PMC), which provides a systematic way to eliminate the renormalization scale-and-scheme dependence in perturbative calculations. The PMC method has rigorous theoretical foundations, it satisfies Renormalization Group Invariance (RGI) and preserves all self-consistency conditions derived from the renormalization group. The PMC cancels the renormalon growth, reduces to the Gell-Mann–Low scheme in the N c → 0 Abelian limit and leads to scale- and scheme-invariant results. The PMC has now been successfully applied to many high-energy processes. In this article we summarize recent developments and results in solving the renormalization scale and scheme ambiguities in perturbative QCD. In particular, we present a recently developed method the PMC ∞ and its applications, comparing the results with CSS. The method preserves the property of renormalizable SU(N)/U(1) gauge theories defined as Intrinsic Conformality (iCF). This property underlies the scale invariance of physical observables and leads to a remarkably efficient method to solve the conventional renormalization scale ambiguity at every order in pQCD. This new method reflects the underlying conformal properties displayed by pQCD at NNLO, eliminates the scheme dependence of pQCD predictions and is consistent with the general properties of the PMC. A new method to identify conformal and β-terms, which can be applied either to numerical or to theoretical calculations is also shown. We present results for the thrust and C-parameter distributions in e + e - annihilation showing errors and comparison with the CSS. We also show results for a recent innovative comparison between the CSS and the PMC ∞ applied to the thrust distribution investigating both the QCD conformal window and the QED N c → 0 limit. In order to determine the thrust distribution along the entire renormalization group flow from the highest energies to zero energy, we consider the number of flavors near the upper boundary of the conformal window. In this flavor-number regime the theory develops a perturbative infrared interacting fixed point. These results show that PMC ∞ leads to higher precision and introduces new interesting features in the PMC. In fact, this method preserves with continuity the position of the peak, showing perfect agreement with the experimental data already at NNLO. We also show a detailed comparison of the PMC ∞ with the other PMC approaches: the multi-scale-setting approach (PMCm) and the single-scale-setting approach (PMCs) by comparing their predictions for three important fully integrated quantities R e+e- , R$_Τ$ and Γ (H→$b\bar{b}$) up to the four-loop accuracy.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2025

The topology of non-resonant stellarator divertors

We apply topological methods to better understand how the magnetic field in the stellarator edge can be diverted away from the confined region. Our primary method is calculating the winding numbers of closed contours, which gives information on the number and nature of fixed points within a bounded region. We first apply this to the non-resonant divertor Hamiltonian system, and present a simple explanation for the system’s diversion: trajectories are guided away from the confined region by X-points which are ‘unpaired’ i.e. do not have corresponding O-points and therefore do not resemble an island chain. We show how similar phenomena can occur in a similar, axisymmetric Hamiltonian system. Secondly, we find examples of neoclassically optimised stellarators in the quasi-symmetric stellarator repository database which divert the magnetic field via unpaired X-points. We present and discuss three examples, each containing novel phenomena which might be desirable for stellarator divertors. These examples serve as an illustration of new divertor possibilities which exist in realistic stellarators, which may ultimately have application for future experiments and reactors.

X-points, topology↗

DOE OSTI2025

Designing Yield Grids for Airborne Nuclear Explosions to Maintain Decrements in Rayleigh Waveform Amplitudes

This document details a method for selecting grid points that sample the yields of airborne, low yield nuclear explosions that constrains the relative size of seismic waveforms sourced by explosions that are parameterized by this grid. In detail, we impose this constraint so that the vertical component of ground displacement that is sourced by explosions with a fixed height of burst (HoB) at grid point k change by a fixed factor of a when compared to the waveform amplitudes sourced by explosions at grid point k + 1, as observed at a fixed standoff distance. This gridding thereby preserves amplitude ratios along the grid so that they are agnostic to k. We assert that this method is useful in tests that apply algorithms against synthetic or augmented datasets for which the performance of the algorithm must be measured against source size, and waveform amplitudes must change uniformly (e.g., Berg et al, 2025).

58 GEOSCIENCES↗

DOE OSTI2020

Fixed-point quantum circuits for quantum field theories

Renormalization group ideas and effective operators are introduced to efficiently prepare ground states of massive lattice field theories on digital quantum devices. This is accomplished with a systematic approximation through localized unitaries that removes an exponentially costly barrier in the spatial volume of the quantum simulation. With these methods, classically computed ground states in a spatial volume L, containing a few Compton wavelengths, can be used to determine operators for preparing the ground state toward the thermodynamic limit with a precision improving as e –mL on beyond-classical quantum registers. Here, due to the exponential spatial decay of correlations in massive theories and the double exponential suppression of digitization artifacts in the number of qubits representing the scalar field, the derived fixed-point quantum circuits are expected to be relevant for simulations of quantum field theories throughout the evolution from small-scale near-term quantum devices to large-scale fault-tolerant quantum computers.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

Analysis of a New Implicit Solver for a Semiconductor Model

Here, we present and analyze a new iterative solver for implicit discretizations of a simplified Boltzmann--Poisson system. The algorithm builds on recent work that incorporated a sweeping algorithm for the Vlasov--Poisson equations as part of nested inner-outer iterative solvers for the Boltzmann--Poisson equations. The new method eliminates the need for nesting and requires only one transport sweep per iteration. It arises as a new fixed-point formulation of the discretized system which we prove to be contractive for a given electric potential. We also derive an accelerator to improve the convergence rate for systems in the drift-diffusion regime. We numerically compare the efficiency of the new solver, with and without acceleration, against a recently developed nested iterative solver.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Qubit Regularization and Qubit Embedding Algebras

Qubit regularization is a procedure to regularize the infinite dimensional local Hilbert space of bosonic fields to a finite dimensional one, which is a crucial step when trying to simulate lattice quantum field theories on a quantum computer. When the qubit-regularized lattice quantum fields preserve important symmetries of the original theory, qubit regularization naturally enforces certain algebraic structures on these quantum fields. We introduce the concept of qubit embedding algebras (QEAs) to characterize this algebraic structure associated with a qubit regularization scheme. We show a systematic procedure to derive QEAs for the O(N) lattice spin models and the SU(N) lattice gauge theories. While some of the QEAs we find were discovered earlier in the context of the D-theory approach, our method shows that QEAs are far richer. A more complete understanding of the QEAs could be helpful in recovering the fixed points of the desired quantum field theories.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Effect of topology changes on the breakup of a periodic liquid jet

Here, the breakup of a periodic jet is examined computationally, using a front-tracking/finite-volume method, where the interface is represented by connected marker points moving with the fluid, while the governing equations are solved on a fixed grid. Tracking the interface allows control of whether topology changes take place or not. The Reynolds and Capillary numbers are kept relatively low ($Re = 150$ and $Ca = 2$) so most of the flow is well resolved. The effect of topology changes is examined by following the jet until it has mostly disintegrated, for different “coalescence criterion,” based on the thickness of thin films and threads. The evolution of both two-dimensional and fully three-dimensional flows is examined. It is found that although there is a significant difference between the evolution when no breakup takes place and when it does, once breakup takes place the evolution is relatively insensitive to exactly how it is triggered for a range of coalescence criterion, and any differences are mostly confined to the smallest scales.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

Delayed Neutron Precursor Group Parameter and Spectra Generation from Fast Fission of 235U in SCALE

Delayed neutron precursor (DNP) group data are important for modeling reactor dynamics. Although the data for individual DNPs have been developed over time, the DNP group data present in the Evaluated Nuclear Data Files (ENDF) have not been updated in the past 20 years. In this work, we use SCALE to recreate the Godiva experiment that was used to generate the original DNP group structure for fast fission of 235U. However, each DNP is modeled using up-to-date data, and the results are then converted into a newly updated group structure. This conversion uses an iterative linear least squares solver to minimize chi-squared. The approaches used in this work also enable energy spectrum generation and uncertainty tracking. The method used in this paper for fast 235U fission DNP group structure updating can be applied to different energies and fissile nuclides. Demonstration of the uncertainty tracking in reactor kinetics and dynamics simulations is shown using point reactor kinetics simulations. Results show that there are data discrepancies between the International Atomic Energy Agency database and data used in ORIGEN, which are currently being fixed. Results also show that the proposed method for group spectra generation performs well.

Seifert, Luke [University of Illinois]↗

DOE OSTI2025

Realizability-preserving discontinuous Galerkin method for spectral two-moment radiation transport in special relativity

Here we present a realizability-preserving numerical method for solving a spectral two-moment model to simulate the transport of massless, neutral particles interacting with a steady background material moving with relativistic velocities. The model is obtained as the special relativistic limit of a four-momentum-conservative general relativistic two-moment model. Using a maximum-entropy closure, we solve for the Eulerian-frame energy and momentum. The proposed numerical method is designed to preserve moment realizability, which corresponds to moments defined by a nonnegative phase-space density. The realizability-preserving method is achieved with the following key components: (i) a discontinuous Galerkin phase-space discretization with specially constructed numerical fluxes in the spatial and energy dimensions; (ii) a strong stability-preserving implicit-explicit time-integration method; (iii) a realizability-preserving conserved to primitive moment solver; (iv) a realizability-preserving implicit collision solver; and (v) a realizability-enforcing limiter. Component (iii) is necessitated by the closure procedure, which closes higher order moments nonlinearly in terms of primitive moments. The nonlinear conserved to primitive and the implicit collision solves are formulated as fixed-point problems, which are solved with custom iterative solvers designed to preserve the realizability of each iterate. With a series of numerical tests, we demonstrate the accuracy and robustness of this discontinuous-Galerkin-implicit-explicit method.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Entanglement properties of disordered quantum spin chains with long-range antiferromagnetic interactions

Entanglement measures are useful tools in characterizing otherwise unknown quantum phases and indicating transitions between them. Here we examine the concurrence and entanglement entropy in quantum spin chains with random long-range couplings, spatially decaying with a power-law exponent $\textit{α}$. Using the strong disorder renormalization group (SDRG) technique, we find by analytical solution of the master equation a strong disorder fixed point, characterized by a fixed point distribution of the couplings with a finite dynamical exponent, which describes the system consistently in the regime $α > \frac{1}{2}$ . A numerical implementation of the SDRG method yields a power-law spatial decay of the average concurrence, which is also confirmed by exact numerical diagonalization. However, we find that the lowest-order SDRG approach is not sufficient to obtain the typical value of the concurrence. We therefore implement a correction scheme which allows us to obtain the leading-order corrections to the random singlet state. This approach yields a power-law spatial decay of the typical value of the concurrence, which we derive both by a numerical implementation of the corrections and by analytics. Next, using numerical SDRG, the entanglement entropy (EE) is found to be logarithmically enhanced for all $\textit{α}$, corresponding to a critical behavior with an effective central charge $\textit{c}$ = ln(2), independent of $\textit{α}$. This is confirmed by an analytical derivation. Using numerical exact diagonalization (ED), we confirm the logarithmic enhancement of the EE and a weak dependence on $\textit{α}$. For a wide range of partition size $\textit{l}$, the EE fits a critical behavior with a central charge close to $\textit{c}$ = 1, which is the same as for the clean Haldane-Shastry model with a power-law-decaying interaction with $\textit{α}$ = 2. Only for small $l \ll L$, in a range which increases with the number of spins N, we find deviations which are rather consistent with the strong disorder fixed point central charge $\textit{c}$ = ln(2). Furthermore, we find using ED that the concurrence shows power-law decay, albeit with smaller power exponents than obtained by SDRG. Finally, we also present results obtained with DMRG and find agreement with ED for sufficiently small $\textit{α}$ < 2, whereas for larger $\textit{α}$ DMRG tends to underestimate the entanglement entropy and finds a faster decaying concurrence.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗