Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SERIES EXPANSION”

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

Interaction models and configurational entropies of binary MoTa and the MoNbTaW high entropy alloy

We introduce a simplified method to model the interatomic interactions of high entropy alloys based on a lookup table of cluster energies. Furthermore, these interactions are employed in replica exchange Monte Carlo simulations with histogram analysis to obtain thermodynamic properties across a broad temperature range. Kikuchi's cluster variation method entropy formalism is applied to directly calculate entropy from statistics on short- and long-range chemical order, and we discuss the convergence of the entropy as clusters of differing size are included. A high temperature series expansion aids in our understanding of the convergence. Computer codes implementing these methods, and supporting data, are freely available on the internet.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Almost strong zero modes at finite temperature

Interacting fermionic chains exhibit extended regions of topological degeneracy of their ground states as a result of the presence of Majorana or parafermionic zero modes localized at the edges. In the opposite limit of infinite temperature, the corresponding nonintegrable spin chains, obtained via generalized Jordan-Wigner mapping, are known to host so-called almost strong zero modes, which are long-lived with respect to any bulk excitations. Here we study the fairly unexplored territory that bridges these two extreme cases of zero and infinite temperature. We blend two established techniques for states, the Lanczos series expansion and a tensor network ansatz, uplifting them to the level of operator algebra. This allows us to efficiently simulate large system sizes for arbitrarily long timescales and to extract the temperature-dependent decay rates. We observe that for the Kitaev-Hubbard model, the decay rate of the edge mode depends exponentially on the inverse temperature 𝛽, and on an effective energy scale Δ eff that is greater than the thermodynamic gap of the system Δ.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

An iterative reconstruction algorithm for pole figure inversion using total variation regularization

This paper describes a new discrete method for inverting X-ray pole figures to estimate the orientation distribution function (ODF). The method employs the equal-volume `cubochoric' representation for uniform discretization of orientation space, SO(3). The forward-projection model is combined with an anisotropic total variation term to iteratively determine the ODF. Here, the efficacy of the new method is evaluated with both model and experimental data and compared with existing discrete and series expansion methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Off-axis representation of hyperbolic mirror shapes for X-ray beamlines

Mirror-centered, closed-form expressions for hyperbolic surfaces used in X-ray beamlines have been derived. Hyperbolic mirrors create a virtual focus or source point and can be used to lengthen or shorten the effective focal distance of a compound optical system. The derivations here express off-axis segments of a hyperbolic surface in terms of the real and virtual focal distances and the incident glancing angle at the center of the mirror. Conventional mathematical expressions of hyperbolic shapes describe the surfaces in Cartesian or polar coordinates centered on an axis of symmetry, necessitating cumbersome rotation and translation to mirror-centered coordinates. The representation presented here, with zero slope and the origin at the central point, is most convenient for modeling, metrology, aberration correction, and general surface analysis of off-axis configurations. The direct derivation avoids the need for nested coordinate transforms. A series expansion provides a helpful approximation; the coefficients of the implicit equation are also provided.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Dynamic Phasor Modeling of Various Multipulse Rectifiers and a VSI Fed by 18-Pulse Asymmetrical Autotransformer Rectifier Unit for Fast Transient Analysis

To study the fast modelling of dynamic characteristics caused by power electronic switching, the dynamic phasor (DP) theory based on the time-varying Fourier decomposition and frequency shifting is firstly applied to the modelling of various multipulse rectifiers. The DP model for symmetrical 12-pulse phase-shifting reactor rectifier unit (PSR-RU) is derived with model order reduction, relations between ac and dc terminals, and Taylor series expansion. The DP model of asymmetrical 18-pulse autotransformer rectifier unit (AT-RU) is proposed based on the switching functions expressed in the DP domain. Meanwhile, the DP model of voltage source inverter (VSI) fed by asymmetrical 18-pulse AT-RU is built with the harmonic state-space (HSS) equations. Under both balanced and unbalanced conditions, the good calculation accuracy and rapid simulation speed of the developed DP models are validated by the detailed time-domain (TD) simulation.

42 ENGINEERING↗
DOE OSTI2017

Optimal Bandwidth for Multitaper Spectrum Estimation

A systematic method for bandwidth parameter selection is desired for Thomson multitaper spectrum estimation. We give a method for determining the optimal bandwidth based on a mean squared error (MSE) criterion. When the true spectrum has a second-order Taylor series expansion, one can express quadratic local bias as a function of the curvature of the spectrum, which can be estimated by using a simple spline approximation. This is combined with a variance estimate, obtained by jackknifing over individual spectrum estimates, to produce an estimated MSE for the log spectrum estimate for each choice of time-bandwidth product. The bandwidth that minimizes the estimated MSE then gives the desired spectrum estimate. Additionally, the bandwidth obtained using our method is also optimal for cepstrum estimates. We give an example of a damped oscillatory (Lorentzian) process in which the approximate optimal bandwidth can be written as a function of the damping parameter. Furthermore, the true optimal bandwidth agrees well with that given by minimizing estimated the MSE in these examples.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Uncertainty reduction in residual stress measurements by an optimised inverse solution using nonconsecutive polynomials

Many destructive methods for measuring residual stresses such as the slitting method require an inverse analysis to solve the problem. The accuracy of the result as well as an uncertainty component (the model uncertainty) depends on the basis functions used in the inverse solution. The use of a series expansion as the basis functions for the inverse solution was analysed in a previous work for the particular case where functions orders grew consecutively. The present work presents a new estimation of the model uncertainty and a new improved methodology to select the final basis functions for the case where the basis is composed of polynomials. Including nonconsecutive polynomial orders in the basis generates a larger space of possible solutions to be evaluated and allows the possibility to include higher-order polynomials. The paper includes a comparison with two other inverse analyses methodologies applied to synthetically generated data. With the new methodology, the final error is reduced and the uncertainty estimation improved.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Dual energy CT: How well can pseudo-monochromatic imaging reduce metal artifacts?

Purpose: Dual Energy CT (DECT) provides so-called monoenergetic images based on a linear combination of the original polychromatic images. At certain patient-specific energy levels, corresponding to certain patient- and slice-dependent linear combination weights, e.g., E = 160 keV corresponds to α = 1.57, a significant reduction of metal artifacts may be observed. The authors aimed at analyzing the method for its artifact reduction capabilities to identify its limitations. The results are compared with raw data-based processing. Methods: Clinical DECT uses a simplified version of monochromatic imaging by linearly combining the low and the high kV images and by assigning an energy to that linear combination. Those pseudo-monochromatic images can be used by radiologists to obtain images with reduced metal artifacts. The authors analyzed the underlying physics and carried out a series expansion of the polychromatic attenuation equations. The resulting nonlinear terms are responsible for the artifacts, but they are not linearly related between the low and the high kV scan: A linear combination of both images cannot eliminate the nonlinearities, it can only reduce their impact. Scattered radiation yields additional noncanceling nonlinearities. This method is compared to raw data-based artifact correction methods. To quantify the artifact reduction potential of pseudo-monochromatic images, they simulated the FORBILD abdomen phantom with metal implants, and they assessed patient data sets of a clinical dual source CT system (100, 140 kV Sn) containing artifacts induced by a highly concentrated contrast agent bolus and by metal. In each case, they manually selected an optimal α and compared it to a raw data-based material decomposition in case of simulation, to raw data-based material decomposition of inconsistent rays in case of the patient data set containing contrast agent, and to the frequency split normalized metal artifact reduction in case of the metal implant. For each case, the contrast-to-noise ratio (CNR) was assessed. Results: In the simulation, the pseudo-monochromatic images yielded acceptable artifact reduction results. However, the CNR in the artifact-reduced images was more than 60% lower than in the original polychromatic images. In contrast, the raw data-based material decomposition did not significantly reduce the CNR in the virtual monochromatic images. Regarding the patient data with beam hardening artifacts and with metal artifacts from small implants the pseudo-monochromatic method was able to reduce the artifacts, again with the downside of a significant CNR reduction. More intense metal artifacts, e.g., as those caused by an artificial hip joint, could not be suppressed. Conclusions: Pseudo-monochromatic imaging is able to reduce beam hardening, scatter, and metal artifacts in some cases but it cannot remove them. In all cases, the CNR is significantly reduced, thereby rendering the method questionable, unless special post processing algorithms are implemented to restore the high CNR from the original images (e.g., by using a frequency split technique). Raw data-based dual energy decomposition methods should be preferred, in particular, because the CNR penalty is almost negligible.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Toward extracting $$\gamma $$ from $$B\rightarrow DK$$ without binning

Abstract $$B^\pm \rightarrow DK^\pm $$ B ± → D K ± transitions are known to provide theoretically clean information about the CKM angle $$\gamma $$ γ , with the most precise available methods exploiting the cascade decay of the neutral D into CP self-conjugate states. Such analyses currently require binning in the D decay Dalitz plot, while a recently proposed method replaces this binning with the truncation of a Fourier series expansion. In this paper, we present a proof of principle of a novel alternative to these two methods, in which no approximations at the level of the data representation are required. In particular, our new strategy makes no assumptions about the amplitude and strong phase variation over the Dalitz plot. This comes at the cost of a degree of ambiguity in the choice of test statistic quantifying the compatibility of the data with a given value of $$\gamma $$ γ , with improved choices of test statistic yielding higher sensitivity. While our current proof-of-principle implementation does not demonstrate optimal sensitivity to $$\gamma $$ γ , its conceptually novel approach opens the door to new strategies for $$\gamma $$ γ extraction. More studies are required to see if these can be competitive with the existing methods.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

A New Analytical Expression for Estimating the Adiabatic Temperature Rise in Lithium-Ion Batteries During High-Power Pulses

A coupled, thermal-electrochemical model is used to explain why the Joule heating assumption (i.e., I 2 R) does not provide a good representation of the temperature rise during high power pulses in lithium-ion batteries, even in cases where the reversible heat generation can be neglected. The poor agreement occurs because the internal resistance changes during the pulse due to the opposing effects of mass transport in the electrolyte, which raises the resistance, and heat generation (temperature rise), which lowers the resistance. These insights are used to propose a new analytical expression for predicting the temperature rise during adiabatic pulses with limited experimental and physical data. The expression accounts for thermal effects using a Taylor series expansion of an Arrhenius-type equation. It accounts for transport effects using an approximate solution to the one-dimensional diffusion equation. The new expression is shown to accurately estimate the simulated, adiabatic temperature rise across a range of loadings (1 to 4 mAh cm -2 ) and C-rates (1C to 10C) for cells containing LiNi 0.5 Mn 0.3 Co 0.2 O 2 positive electrodes and graphite negative electrodes. In conclusion, it is also shown to accurately estimate the experimental, adiabatic temperature rise measured for cells with LiFePO 4 positive electrodes and mesocarbon negative electrodes with minimal changes to the fitted parameters.

25 ENERGY STORAGE↗
DOE OSTI2016

Incorporation of the Joule Heating of highly conducting materials into the Truchas code via an asymptotic approach

The EM module of the Truchas code currently lacks the capability to model the Joule (Ohmic) heating of highly conducting materials that are inserted into induction furnaces from time to time to change the heating profile. This effect is difficult to simulate directly because of the requirement to resolve the extremely thin skin depth of good conductors, which is computationally costly. For example, copper has a skin depth, δ ~ 1 mm, for an oscillation frequency of tens of kHz. The industry is interested in determining what fraction of the heating power is lost to the Joule heating of these good conductors inserted inside the furnaces. The approach presented in this document is one of asymptotics where the leading order (unperturbed) solution is taken as that which emerges from solving the EM problem for a perfectly conducting insert. The conductor is treated as a boundary of the domain. The perturbative correction enters as a series expansion in terms of the dimensionless skin depth δ/L, where L is the characteristic size of the EM system. The correction at each order depends on the previous. This means that the leading order correction only depends on the unperturbed solution, in other words, it does not require Truchas to perform an additional EM field solve. Thus, the Joule heating can be captured by a clever leveraging of the existing tools in Truchas with only slight modifications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Investigation of an Ultrafast Harmonic Resonant RF Kicker

An Energy Recovery Linac (ERL) based multi-turn electron Circulator Cooler Ring (CCR) is envisaged in the proposed Jefferson Lab Electron Ion Collider (JLEIC) to cool the ion bunches with high energy (55 MeV), high current (1.5 A), high repetition frequency (476.3 MHz), high quality magnetized electron bunches. A critical component in this scheme is a pair of ultrafast kickers for the exchange of electron bunches between the ERL and the CCR. The ultrafast kicker should operate with the rise and fall time in less than 2.1 ns, at the repetition rate of ~10s MHz, and should be able to run continuously during the whole period of cooling. These -and-fall time being combined together, are well beyond the state-of-art of traditional pulsed power supplies and magnet kickers. To solve this technical challenge, an alternative method is to generate this high repetition rate, fast rise-and-fall time short pulse continuous waveform by summing several finite number of (co)sine waves at harmonic frequencies of the kicking repetition frequency, and these harmonic modes can be generated by the Quarter Wave Resonater (QWR) based multifrequency cavities. Assuming the recirculator factor is 10, 10 harmonic modes (from 47.63 MHz to 476.3 MHz) with proper amplitudes and phases, plus a DC offset are combined together, a continuous short pulse waveform with the rise-and-fall time in less than 2.1 ns, repetition rate of 47.63 MHz waveform can be generated. With the compact and matured technology of QWR cavities, the total cost of both hardware development and operation can be reduced to a modest level. Focuse on the technical scheme, three main topics will be discussed in this thesis: the synthetization of the kicking pulse, the design and optimization of the deflecting QWR multi-integer harmonic frequency resonator and the fabrication and bench measurements of a half scale copper prototype. In the kicking pulse synthetization part, we begin with the Fourier Series expansion of an ideal square pulse, and get a Flat-Top waveform which will give a uniform kick over the bunch length of the kicked electron bunches, thus the transverse emittance of these kicked electron bunches can be maintained. By using two identical kickers with the betatron phase advance of 180 degree or its odd multiples, the residual kick voltage wave slopes at the unkicked bunch position will be totally cancelled out. Flat-Top waveform combined with two kicker scheme, the transverse emittance of the cooling electron bunches will be conserved during the whole injection, recirculation, and ejection processes. In the cavity design part, firstly, the cavity geometry is optimized to get high transverse shunt impedance thus less than 100 W of RF losses on the cavity wall can be achieved for all these 10 harmonic modes. To support all these 10 harmonic modes, group of four QWRs are adopted with the mode distribution of 5:3:1:1. In the multi-frequency cavities such as the five-mode-cavity and the three-mode-cavity, tunings are required to achieve the design frequencies for each mode. Slight segments of taper design on the inner conductor help to get the frequencies to be exactly on the odd harmonic modes. Stub tuners equal to the number of resonant modes are inserted to the outer conductor wall to compensate the frequency shifts due manufacturing errors and other perturbations during the operation such as the change of the cavity temperature. Single loop couple is designed for all harmonic modes in each cavity. By adjusting its loop size, position and rotation, it is possible to get the fundamental mode critical coupled and other higher harmonic modes slightly over coupled. A broadband circulator will be considered for absorbing the reflected power. Finally in this part, multipole field components due to the asymmetric cylindrical structure around the beam axis of the cavity as well as the beam-induced higher order mode (HOM) issues will be analyzed and discussed in this thesis. A half-scale copper prototype cavity (resonant frequencies from 95.26 MHz to 857.34 MHz) was fabricated to validate the electromagnetic characteristics. With this half scale prototype, the tuning processes of multiple harmonic frequencies, unloaded quality factor measurements of each mode, and bead-pull measurements are performed. The bench measurement results matched well with the simulation results, which have validated our cavity design and construction methods. Finally, a simple mode combining experiment with five separate signal generators was performed on this prototype cavity and the desired fast rise/fall time (1.2 ns), high repetition rate (95.26 MHz) waveform was captured, which finally proved our design of this ultrafast harmonic kicker.

43 PARTICLE ACCELERATORS↗
DOE OSTI2017

Development of an ejecta particle size measurement diagnostic based on Mie scattering

The goal of this work is to determine the feasibility of extracting the size of particles ejected from shocked metal surfaces (ejecta) from the angular distribution of light scattered by a cloud of such particles. The basis of the technique is the Mie theory of scattering, and implicit in this approach are the assumptions that the scattering particles are spherical and that single scattering conditions prevail. The meaning of this latter assumption, as far as experimental conditions are concerned, will become clear later. The solution to Maxwell’s equations for spherical particles illuminated by a plane electromagnetic wave was derived by Gustav Mie more than 100 years ago, but several modern treatises discuss this solution in great detail. The solution is a complicated series expansion of the scattered electric field, as well as the field within the particle, from which the total scattering and absorption cross sections as well as the angular distribution of scattered intensity can be calculated numerically. The detailed nature of the scattering is determined by the complex index of refraction of the particle material as well as the particle size parameter, x, which is the product of the wavenumber of the incident light and the particle radius, i.e. x = 2rπ= λ. Figure 1 shows the angular distribution of scattered light for different particle size parameters and two orthogonal incident light polarizations as calculated using the Mie solution. It is obvious that the scattering pattern is strongly dependent on the particle size parameter, becoming more forward-directed and less polarizationdependent as the particle size parameter increases. This trend forms the basis for the diagnostic design.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Linear and Nonlinear Constitutive Model for Piezoelectricity in ALEGRA-FE

This report develops and documents linear and nonlinear constitutive relations implemented in ALEGRA-FE (ferroelectric). A thermodynamic framework is created to describe the electromechanical system in the form of a free energy functional. Constitutive relations are derived by taking series expansions of the free energy functional with respect to the independent fields. First order expansion terms yield linear constitutive relations and higher order expansion terms yield non-linear constitutive relations. This document serves as supplement to Section 4 of Sandia Report SAND2013-7363, Rev 3. Methods for implementation of kinematic relations of piezoelectric models and rotation of material principal axes are covered in the supplemented report. Additional discussion on phase velocity calculation is also presented.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2016

Analysis of Electromagnetic Scattering From a Dihedral Cut Into a Sphere Using Vector Spherical Harmonics

Proper calibration or verification of the calibration of fully - polarimetric radar requires more than a simple specular target, such as a conducting sphere or corner reflector. The dihedral reflector is a useful radar - scattering target for calibrating the polarization response of a radar system, because it possesses the ability to change the direction of the scattered electric - field vector in a special way, while a sphere or trihedral target does not. However, the scattered electromagnetic field of a sphere can be computed with arbitrary precision using the well - known Mie series, which is obtained by expanding the scattered field as a series of vector spherical harmonics. In contrast, the computation of the scattered electromagnetic field of dihedral reflectors typically relies on numerical methods or various approximate solutions of Maxwell's equations. Precision calibration of wide - bandwidth polarimetric radar has not been adequately achieved for some applications using this approach. This report describes a method for computing the scattered electromagnetic field of a specially shaped dihedral, derived through the application of a series expansion in vector spherical harmonics, with the goal of achieving the same precision that is available for the sphere.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Theory and Implementation of the Spectacular Nonlinear Viscoelastic Constitutive Model

This report is a comprehensive guide to the nonlinear viscoelastic Spectacular model, which is an isotropic, thermo-rheologically simple constitutive model for glass-forming materials, such as amorphous polymers. Spectacular is intermediate in complexity to the previous PEC and SPEC models (Potential Energy Clock and Simplified Potential Energy Clock models, respectively). The model form consists of two parts: a Helmholtz free energy functional and a nonlinear material clock that controls the rate of viscoelastic relaxation. The Helmholtz free energy is derived from a series expansion about a reference state. Expressions for the stress and entropy functionals are derived from the Helmholtz free energy following the Rational Mechanics approach. The material clock depends on a simplified expression for the potential energy, which itself is a functional of the temperature and strain histories. This report describes the thermo-mechanical theory of Spectacular, the numerical methods for time-integrating the model, model verification for its implementation in LAMÉ, a user guide for its implementation in LAMÉ, and ideas for future work. A number of appendices provide supplementary mathematical details and a description of the procedure used to derive the simplified potential energy from the full expression for the potential energy. The goal of this report is create a convenient point-of-entry for engineers who wish to learn more about Spectacular, but also to serve as a reference manual for advanced users of the model.

36 MATERIALS SCIENCE↗