Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical computation”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Design of flexible proton beam imaging energy spectrometers (PROBIES)

A simple repeating step-filtered structure has been designed and deployed for the detection of MeV-energy proton beams produced from high-intensity short-pulse laser-driven experiments. Due to the diagnostic's ability to recover energy-resolved images of proton beams, it is called the proton beam imager and energy spectrometer (or PROBIES). The flexible design has two useful applications. First, it can be fielded in a configuration with a scintillator and imaging setup in order to collect proton beam information in high-repetition rate experiments. Here we show a configuration that collects the equivalent of a nine-film radiochromic film (RCF) stack, and a technique for using a neural network (NN) to rapidly analyze the images on millisecond time-scales. Second, when using RCF stacks as the detector, the energy resolution capability is greatly increased compared to traditional RCF stacks, while simultaneously reducing the total number of films required to make each measurement. We demonstrate a configuration that provides up to 90 spatially- and energy-resolved datapoints using just ten films. Both designs presented here can be easily customized for anticipated experimental outputs and manufactured at low cost. Such designs allow for rapid prototyping using modern high-accuracy 3D printers or on-demand computer numerical control machining to reduce the reliance on time-consuming film media, analysis, or bulky electro-magnetically dispersing spectrometers. The analysis of data is discussed and a methodology for constructing a NN that is capable of rapidly analyzing data on-the-fly with greater than 95% accuracy is presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

The look-elsewhere effect from a unified Bayesian and frequentist perspective

When searching over a large parameter space for anomalies such as events, peaks, objects, or particles, there is a large probability that spurious signals with seemingly high signi ficance will be found. This is known as the look-elsewhere effect and is prevalent throughout cosmology, (astro)particle physics, and beyond. To avoid making false claims of detection, one must account for this effect when assigning the statistical significance of an anomaly. This is typically accomplished by considering the trials factor, which is generally computed numerically via potentially expensive simulations. In this paper we develop a continuous generalization of the Bonferroni and Sidak corrections by applying the Laplace approximation to evaluate the Bayes factor, and in turn relating the trials factor to the prior-to-posterior volume ratio. Here, we use this to define a test statistic whose frequentist properties have a simple interpretation in terms of the global p-value, or statistical significance. We apply this method to various physics-based examples and show it to work well for the full range of p-values, i.e. in both the asymptotic and non-asymptotic regimes. We also show that this method naturally accounts for other model complexities such as additional degrees of freedom, generalizing Wilks' theorem. This provides a fast way to quantify statistical significance in light of the look-elsewhere effect, without resorting to expensive simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Real-time scattering in Ising field theory using matrix product states

We study scattering in Ising field theory (IFT) using matrix product states and the time-dependent variational principle. IFT is a one-parameter family of strongly coupled nonintegrable quantum field theories in 1+1 dimensions, interpolating between massive free fermion theory and Zamolodchikov's integrable massive 𝐸 8 theory. Particles in IFT may scatter either elastically or inelastically. In the postcollision wave function, particle tracks from all final-state channels occur in superposition; processes of interest can be isolated by projecting the wave function onto definite particle sectors, or by evaluating energy density correlation functions. Using numerical simulations we determine the time delay of elastic scattering and the probability of inelastic particle production as a function of collision energy. We also study the mass and width of the lightest resonance near the 𝐸 8 point in detail. Close to both the free fermion and 𝐸 8 theories, our results for both elastic and inelastic scattering are in good agreement with expectations from form-factor perturbation theory. Using numerical computations to go beyond the regime accessible by perturbation theory, we find that the high-energy behavior of the two-to-two particle scattering probability in IFT is consistent with a conjecture of Zamolodchikov. Our results demonstrate the efficacy of tensor-network methods for simulating the real-time dynamics of strongly coupled quantum field theories in 1+1 dimensions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Coulomb corrections for the nonflip and spin-flip electromagnetic 𝑝 ↑⁢ 𝐴 amplitudes

It is demonstrated that, within the eikonal approach, the Coulomb corrections to the elastic electromagnetic nonflip and spin-flip proton-nucleus amplitudes are identical when the two amplitudes share the same exponential form factors. This result allows Coulomb corrections to be computed numerically, and with high precision, for both electromagnetic and hadronic elastic 𝑝 ↑⁢ 𝐴 amplitudes in the massless-photon limit, including the effects of soft magnetic photon exchange. The method relies on analytical expressions and numerical integrations over a finite impact-parameter range with nonsingular integrands, providing a practical and systematically controlled framework for phenomenological applications.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Probing magnetic configuration-mediated topological phases via high harmonic generation in MnBi 2 Te 4

Compared to a space group, a magnetic space group is much more complex, as both magnetic structure and magnetic moment direction can change the original symmetries of materials. The interplay between space symmetry and magnetism can generate versatile novel quantum states. However, detecting these topological phases experimentally, achieved through manipulating the magnetic configuration, has been restricted. It is mainly because the intrinsic link between the theory and the available experimental technique remains elusive. Here, we show that high harmonic generation (HHG) can identify these topological quantum states. In this work, we use rhombohedral MnBi 2 Te 4 film as an example, and analyze the symmetry-dependent harmonic order and signal by combining first-principles calculations and time-dependent Liouville equation numerical computations. Our results provide a fundamental basis for using HHG to study the topological quantum phases mediated by various magnetic configurations that can be easily realized by applying an external magnetic field.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Quasiparton distributions in massive QED2: Toward quantum computation

We analyze the quasiparton distributions of the lightest 𝜂′ meson in massive two-dimensional quantum electrodynamics (QED2) by exact diagonalization. The Hamiltonian and boost operators are mapped onto spin qubits in a spatial lattice with open boundary conditions. The lowest excited state in the exact diagonalization is shown to interpolate continuously between an anomalous 𝜂′ state at strong coupling, and a nonanomalous heavy meson at weak coupling, with a cusp at the critical point. The boosted 𝜂′ state follows relativistic kinematics but with large deviations in the luminal limit. The spatial quasiparton distribution function and amplitude for the 𝜂′ state are computed numerically for increasing rapidity both at strong and weak coupling, and compared to the exact light front results. The numerical results from the boosted form of the spatial parton distributions, compare fairly with the inverse Fourier transformation of the luminal parton distributions, derived in the lowest Fock space approximation. Our analysis points out some of the limitations facing the current lattice program for the parton distributions.

Lattice field theory↗
DOE OSTI2025

Precision string phenomenology

Calabi-Yau compactifications of the E 8 × E 8 heterotic string provide a promising route to recovering the four-dimensional particle physics described by the Standard Model. While the topology of the Calabi-Yau space determines the overall matter content in the low-energy effective field theory, further details of the compactification geometry are needed to calculate the normalized physical couplings and masses of elementary particles. In this work, we present numerical computations of physical Yukawa couplings in a number of heterotic models in the standard embedding and demonstrate the existence of natural hierarchies, a coveted feature in string model building. Published by the American Physical Society 2025

Berglund, Per (ORCID:0000000316754133)↗
DOE OSTI2025

Emergent hydrodynamic mode on SU(2) plaquette chains and quantum simulation

We search for emergent hydrodynamic modes in real-time Hamiltonian dynamics of 2+1-dimensional SU(2) lattice gauge theory on a quasi-one-dimensional plaquette chain, by numerically computing symmetric correlation functions of energy densities on lattice sizes of about 20 with the local Hilbert space truncated at 𝑗 max = $\frac{1}{2}$. Because of the Umklapp processes, we only find a mode for energy diffusion. The symmetric correlator exhibits transport peak near zero frequency with a width approximately proportional to momentum squared at small momentum, when the system is fully quantum ergodic, as indicated by the eigenenergy level statistics. This transport peak leads to a power-law 𝑡 −$\frac{1}{2}$ decay of the symmetric correlator at late time, also known as the long-time tail, as well as diffusionlike spreading in position space. We also introduce a quantum algorithm for computing the symmetric correlator on a quantum computer and find it gives results consistent with exact diagonalization when tested on the IBM emulator. Finally we discuss the future prospect of searching for the sound modes.

Hamiltonian systems↗
DOE OSTI2021

Thermalization in the one-dimensional Salerno model lattice

The Salerno model constitutes an intriguing interpolation between the integrable Ablowitz-Ladik (AL) model and the more standard (nonintegrable) discrete nonlinear Schrödinger (DNLS) one. The competition of local on-site nonlinearity and nonlinear dispersion governs the thermalization of this model. Here, we investigate the statistical mechanics of the Salerno one-dimensional lattice model in the nonintegrable case and illustrate the thermalization in the Gibbs regime. As the parameter interpolating between the two limits (from DNLS toward AL) is varied, the region in the space of initial energy and norm densities leading to thermalization expands. The thermalization in the non-Gibbs regime heavily depends on the finite system size; we explore this feature via direct numerical computations for different parametric regimes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Eikonal spin-dependent odderon and gluon Sivers function of a proton and its small-𝑥 evolution

The matrix element in the proton of the eikonal Odderon operator, with a helicity flip, has been shown to correspond to the dipole gluon Sivers function. We employ a three quark light-front model of the proton to determine the Sivers function at moderately small 𝑥 0 ∼ 0.1 and transverse momentum 𝑘 ⊥ ≲ 1 GeV. The model light-cone (LC) wave function predicts the properties of 𝑥⁢𝑓$^{⊥⁢𝑔}_{1⁢𝑇⁡}$(𝑥,𝑘 ⊥ ) such as its overall magnitude, the position of its peak in 𝑘 ⊥ , and its behavior at small 𝑘 ⊥ . We then compute numerically the Balitsky-Fadin-Kuraev-Lipatov anomalous dimension characterizing the power-law tail at 𝑘 ⊥ ≳ 1.5 GeV of the gluon Sivers function at small (but preasymptotic) LC momentum fractions 𝑥 ∼ 𝑥 0 ⁢𝑒 −1/𝛼 𝑠 : 𝑥⁢𝑓$^{⊥⁢𝑔}_{1⁢𝑇⁡}$⁡(𝑥,𝑘 ⊥ ) ∼ 𝑘$^{−3.3}_{⊥}$.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Open-source electrochemical cell for in situ X-ray absorption spectroscopy in transmission and fluorescence modes

X-ray spectroscopy is a valuable technique for the study of many materials systems. Characterizing reactions in situ and operando can reveal complex reaction kinetics, which is crucial to understanding active site composition and reaction mechanisms. In this project, the design, fabrication and testing of an open-source and easy-to-fabricate electrochemical cell for in situ electrochemistry compatible with X-ray absorption spectroscopy in both transmission and fluorescence modes are accomplished via windows with large opening angles on both the upstream and downstream sides of the cell. Using a hobbyist computer numerical control machine and free 3D CAD software, anyone can make a reliable electrochemical cell using this design. Onion-like carbon nanoparticles, with a 1:3 iron-to-cobalt ratio, were drop-coated onto carbon paper for testing in situ X-ray absorption spectroscopy. Cyclic voltammetry of the carbon paper showed the expected behavior, with no increased ohmic drop, even in sandwiched cells. Chronoamperometry was used to apply 0.4 V versus reversible hydrogen electrode, with and without 15 min of oxygen purging to ensure that the electrochemical cell does not provide any artefacts due to gas purging. The XANES and EXAFS spectra showed no differences with and without oxygen, as expected at 0.4 V, without any artefacts due to gas purging. The development of this open-source electrochemical cell design allows for improved collection of in situ X-ray absorption spectroscopy data and enables researchers to perform both transmission and fluorescence simultaneously. It additionally addresses key practical considerations including gas purging, reduced ionic resistance and leak prevention.

08 HYDROGEN↗
DOE OSTI2025

Securing Smart Manufacturing: Detection of Cyber-Physical Attacks in CNC-Based Systems

As Industry 4.0 advances, the integration of computer numerical control (CNC) machines and advanced manufacturing technologies is transforming production into smart manufacturing systems that blend physical and digital processes as cyber-physical systems. However, this increased cyber-physical connectivity exposes manufacturing systems to cyber threats that can cause severe operational and financial disruptions. This paper presents a comparative study on cyber attacks and anomaly detection techniques in manufacturing, focusing on network traffic from CNC machines. The data extracted from network packets includes machine commands and control signals exchanged between the machine's interface and control system, crucial for maintaining operational integrity. We explore two types of cyber attacks, design modification and command injection, which pose substantial risks to CNC machine productivity and system integrity. Our investigation involves experiments on a real CNC system, highlighting the urgent need for effective detection mechanisms. To address these threats, we evaluate three anomaly detection methods: dynamic time warping (DTW), rolling average, and a deep learning, long short-term memory (LSTM) time-series-based autoencoder. Each is assessed for its effectiveness in identifying anomalous behaviors caused by the attacks. Our findings demonstrate the unique strengths and limitations of each detection technique, providing a deeper understanding of their applicability in realworld manufacturing environments. The comparative analysis indicates that while certain methods are highly effective against specific attack types, others offer broader applicability across different attacks. This study contributes to the accurate detection of anomalies in CNC machining processes, thereby enhancing the reliability and security of smart manufacturing systems against diverse cyber threats.

Williams, Bethanie [Tennessee Technological Univer↗
DOE OSTI2020

A Review on Artificial Intelligence for Grid Stability Assessment

Artificial intelligence provides a convenient route for power grid stability assessment. Compared with simulation-based approaches, artificial intelligence can potentially save time on model development and numerical computation in stability assessment. This paper first reviewed existing literature on using artificial intelligence for power grid stability assessment. Then a machine-leaning-based tool is presented and developed to assess power grid transient stability, frequency stability, and small signals stability. Test results verified the accuracy and effectiveness of the AI tool for power grid stability assessment.

You, Shutang↗
DOE OSTI2021

Comparison of Rotary and Linear Cutting Methodology in Determining Specific Cutting Energy of Granite

Abstract Single polycrystalline diamond compact (PDC) cutting is a practical technique to understand the rock–tool interactions in drag-bit-type geothermal drilling operations. This paper introduces a rotary cutting method to determine specific cutting energy (SCE) and compares it with the conventional linear cutting method. In this work, granite is selected to represent hard rock formations in geothermal drilling. Cutting tests are conducted on a Computer Numerical Control (CNC) machine with a realistic cutting speed of 12.7 m/min and several chip loads ranging from 0.08 to 0.25 mm. The cutting force is measured using a dynamometer and then converted to specific cutting energy (SCE). The results show that the rotary method produces an inverse relationship between SCE and chip load, whereas the linear method shows the opposite. As a result, the produced SCE by the rotary method tends to be lower than that of the linear method at a higher chip load at and over 0.16 mm. The difference may be attributed to the cutting configuration and associated force components.

Engineering↗
DOE OSTI2022

Machine Straightness Error Measurement Based on Optical Fiber Fabry–Pérot Interferometer Monitoring Technique

Abstract In this research, we propose an Error Separation Technique (EST) based on optical fiber sensors for on-machine straightness error measurement. Two fiber optic Fabry–Pérot interferometers have been developed serving as two displacement sensors. The displacement distance is computed according to the reflected spectrum from interferometers, which can achieve a sub-micrometer resolution. The two-point method has been employed to separate the straightness error of the slides and the profile error of a fine-polished standard block. The spacing distance between two interferometers is determined by the diameter of optical fibers so that the EST's resolution has the potential to reach the sub-millimeter scale. In the experiment, the straightness error has been measured on a commercially available computer numerical control machine tool, and the measurement has been conducted on its x-axis. The spacing distance between two optical fiber sensors is 1.5 mm which equals the EST's resolution along the machine tool's x-axis. The separated profile error of the measured standard block is around 30 µm which has been verified by a high precision Coordinate Measurement Machine (CMM). The magnitude of the separated straightness error is around 40 µm. This technique is flexible and simple to be conducted, which can contribute to the micro-machine tool calibration and other straightness error applications.

Engineering↗
DOE OSTI2023

Modeling and Optimal Design of a Desalination System Integrated Between a Glass-Covered Solar Collection Water Chamber and a Heat Dissipating Chimney

Abstract Solar thermal energy is trapped in a glass-covered water chamber/basin to provide the heat for evaporation of brackish water in a novel desalination system. To harvest clean water, a chimney is integrated with the water basin, which draws air into the chamber to be humidified by the vapor and then being ventilated via buoyance force. Uniquely, thermal conductive metal sheet is recommended to make the chimney, which allows the vapor in the humid air condenses easily when flowing up. The condensate at the inner wall of the chimney flows down to be collected as clean water. Mathematical modeling and numerical computation have been carried out to delineate the coupling of the buoyance-force-driven flow with the heat and mass transfer of air and water in the solar collection chamber and the condensation of vapor in the heat-dissipating chimney. The objective of the simulation and optimization of the system is to find the best match of the dimensions of the water chamber with a chimney to maximize the production of clean water and energy efficiency. The model has been used to simulation several cases (of water from 40 °C to 50 °C) with available experimental data from the authors’ previous work, and the agreement was satisfactory. The optimization studies found that there is a maximum air flowrate corresponding to a critical chimney height due to the requirement that the chimney is designated to dissipate heat as much as possible to condense water vapor. With the chimney height greater than the critical height, the airflow rate will have a slight decrease. Nevertheless, higher than a critical height is still needed for a chimney to condense more moisture. Optimized chamber diameters at different chimney heights are provided for reference of optimal system designs.

Engineering↗
DOE OSTI2025

Rapid Commissioning of Large Machine Tools Using Finite Element-Based Correction of Geometric Errors

Large computer numerical control (CNC) machine tools derive their stiffness from monolithic cast iron bases or weldments that are sometimes integral to machine motion systems like box ways or guideways. However, the sheer size of castings and even floor flatness deviations result in dimensional errors in these systems, which manifest as machine motion errors. Typical geometric alignment processes rely on an iterative approach, where measurements are taken to assess alignment (straightness, squareness, and parallelism), followed by adjustment of the machine supports (fixators or leveling pads), which can take weeks even for an experienced operator. Conversely, a novel method is proposed to shorten the correction time by eliminating the trial-and-error process in favor of a more deterministic approach guided by a finite element (FE) method. A feasibility study is conducted on a CNC polymer hybrid machine, with a steel weldment frame, supported by six leveling pads. An FE model of the frame is utilized to obtain recommended leveling pad adjustments, based on measurement of machine errors taken using a laser tracker. After a single adjustment cycle, measurements reveal that geometric errors of the machine tool are reduced from 2.22 mm of flatness deviation to 0.32 mm, achieving an 85.6% reduction. Furthermore, the entire process including measurement, adjustment, and assessment is completed in just 6 h by two operators who are not professional service engineers. In conclusion, this methodology demonstrates feasibility for scaling up, especially to large, high-precision CNC machine tools with bases mounted by fixators, offering the capability for bidirectional adjustment.

42 ENGINEERING↗
DOE OSTI2021

Prediction of Thermal Conditions of DED With FEA Metal Additive Simulation

This paper presents the integration of wire-arc additive manufacturing (WAAM) using Gas Metal Arc Welding (GMAW) into a machine tool to create a retrofit hybrid computer numeric control (CNC) machine tool. GMAW, along with other direct energy deposition systems, has the capacity to deposit material faster than the excess thermal energy can dissipate. This results in the need to allow the part to cool between consecutive layers, which is the most time-consuming part of the additive process. Finite element analysis (FEA) was used in conjunction with monitored build plate surface temperatures during deposition samples to improve adequate dwell time prediction and to develop a cooling system. A deposition was completed where no dwell time was used and the build plate along with the machine table temperatures were monitored. A second deposition was completed where only one bead was deposited and the traverse speed was increased. The GMAW welder was mounted on a 3-axis CNC machine where two square deposition samples were completed. A FEA model was designed and verified using the monitored samples. The model will be used to determine improved depositions speeds and whether forced cooling would allow for an increased deposition rate without structural failure. It was determined the FEA software can be used to accurately model and predict the thermal response of WAAM AM components.

Heinrich, Lauren↗