Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “normalization 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.

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

At least 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

SDA: a symbolic differential algebra package in C++

Truncated Power Series Algebra (TPSA), or Differential Algebra (DA), is a well-established tool in accelerator physics, commonly used for generating high-order maps of dynamic systems, as well as in symplectic tracking, normal form analysis, verified integration, optimization, and fast multipole methods. This package is the first to perform symbolic DA computations, enabling traceability of initial condition contributions and runtime reduction for repeated DA calculations, potentially expanding DA’s applications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

A Continuous-Discontinuous Galerkin Method for Electromagnetic Simulations Based on an All-Frequency Stable Formulation

In this paper, a potential-based partial-differential formulation, called the all-frequency stable formulation, is presented for the accurate and robust simulation of electromagnetic problems at all frequencies. Due to its stability from (near) dc to microwave frequencies, this formulation can be applied to simulate wide-band and multiscale problems without encountering the infamous low-frequency breakdown issue or applying basis function decompositions such as the tree-cotree splitting technique. To provide both efficient and flexible numerical solutions to the electromagnetic formulation, a mixed continuous-discontinuous Galerkin (CDG) method is proposed and implemented. In regions with homogeneous media, the continuous Galerkin method is employed to avoid the introduction of duplicated degrees of freedom (DoFs) on the elemental interfaces, while on the interfaces of two different media, the discontinuous Galerkin method is applied to permit the jump of the normal components of the electromagnetic fields. Numerical examples are provided to validate and demonstrate the proposed numerical solver for problems in a wide electromagnetic spectrum.

Yan, Su↗
DOE OSTI2022

The use of subchondral bone topography to approximate the center of rotation of the elbow joint in dogs

Abstract Objectives The aim of this study was to compare the approximate center of rotation in normal and diseased elbows in dogs. Study design Cross‐sectional study Sample population Computer tomography scans of nine dogs with unilateral fragmented medial coronoid process (FMCP). Methods A board certified radiologist confirmed that each dog had unilateral FMCP, and a normal contralateral elbow. Digital 3D models of all elbow joints were uploaded into a surgical planning software package. Four axes approximating the center of rotation (COR) of elbow joints were generated using five geometric shapes based on subchondral topography of the humeral condyle radius and ulna. Images showing the locations where axes exited the medial and lateral cortex of the humeral condyle were captured and imported into a second software package, for measurement of distances between exit points and the origin of a system of axes. Results In normal joints 20/27 (74%) axes exited the medial cortex, and 25/27 (93%) axes exited the lateral cortex cranial and distal to the medial and lateral epicondyles, respectively. In diseased joints 22/27 (81%) axes exited medial cortex and 19/27 (70%) axes exited the lateral cortex, caudal and distal to the medial and lateral epicondyles, respectively. Conclusion Based on CT‐ derived geometry, the COR of elbow affected with FMCP was generally more caudal than normal. Clinical significance External landmarks approximating the location of the elbow COR are provided, and while not validated, may assist in planning, creation, and assessment of procedures for FMCP.

Berger, Chen↗
DOE OSTI2025

A cohesive zone treatment for the material point method involving problems of large deformation and damage

A new algorithm is described that permits the use of cohesive zones in the material point method for problems involving large deformation and fracture. In contrast to previous cohesive zone implementations, this method does not utilize massless surface-element particles. Instead, cohesive tractions are computed using the shape function mappings from a reference grid configuration in combination with explicitly defined particle surface normals and surface positions. These normals and relative surface positions are updated each time step according to particle deformation. The tractions are converted to cohesive forces using the nodal areas and mapped back to particles using the same reference shape function mappings. These forces are then remapped by conventional particle-to-grid interpolation as external forces using the current-configuration shape-function mappings. This allows highly compliant cohesive zones to function over jump displacements larger than a grid cell. Upon damage, these interfaces can revert to conventional multi-field contact surfaces. This approach is general and readily applies to two and three dimensions as well as being compatible with damage-field gradient partitioning offering exceptional computational flexibility. The framework for this method enables other capabilities, such as improved contact precision using explicitly defined surface normals and positions, and a method to mitigate spurious material damage at weak discontinuities between stiff brittle materials and soft or compliant materials.

Cohesive zone↗
DOE OSTI2020

Normal modes in thermal AdS via the Selberg zeta function

The heat kernel and quasinormal mode methods of computing 1-loop partition functions of spin s s fields on hyperbolic quotient spacetimes \mathbb{H}^{3}/\mathbb{Z} ℍ 3 / ℤ are related via the Selberg zeta function. We extend that analysis to thermal \text{AdS}_{2n+1} AdS 2 n + 1 backgrounds, with quotient structure \mathbb{H}^{2n+1}/\mathbb{Z} ℍ 2 n + 1 / ℤ . Specifically, we demonstrate the zeros of the Selberg function encode the normal mode frequencies of spin fields upon removal of non-square-integrable modes. With this information we construct the 1-loop partition functions for symmetric transverse traceless tensors in terms of the Selberg zeta function and find exact agreement with the heat kernel method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Uncertainty quantification of mass models using ensemble Bayesian model averaging

Developments in the description of the masses of atomic nuclei have led to various nuclear mass models that provide predictions for masses across the whole chart of nuclides. These mass models play an important role in understanding the synthesis of heavy elements in the rapid neutron capture ( r ) process. However, it is still a challenging task to estimate the size of uncertainty associated with the predictions of each mass model. In this work, a method called ensemble Bayesian model averaging (EBMA) is introduced to quantify the uncertainty of one-neutron separation energies (S 1 n ) which are directly relevant in the calculations of r -process observables. Here, this Bayesian method provides a natural way to perform model averaging, selection, and uncertainty quantification, by combining the mass models as a mixture of normal distributions whose parameters are optimized against the experimental data, employing the Markov chain Monte Carlo method using the no-u-turn sampler. The EBMA model optimized with all the experimental S 1 n from the AME2003 nuclides are shown to provide reliable uncertainty estimates when tested with the new data in the AME2020.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Efficient perturbation-tracking method for directly probing the spectral phonon properties from molecular dynamics simulations

Existing methods for directly extracting the spectral phonon properties from molecular dynamics (MD) simulations, like the normal mode analysis (NMA) and spectral energy density analysis, all require a very long simulation time to produce reliable results with good convergence. So far, these methods are mainly applied in studies using small systems and with empirical potentials, as the heavy computational load has greatly hindered their further applications. Here we propose a perturbation-tracking (PT) method for directly probing the mode-wise phonon anharmonic frequencies and lifetimes. We show that results obtained from our method are in excellent agreement with those from the conventional NMA approach, using Si as the model material system. Comparing with the NMA approach, the PT method offers a greater accuracy and significant improvement of efficiency. It takes an average of two orders of magnitude and up to three orders of magnitude less simulation time to obtain the same lifetime result of a phonon mode with intermediate to high accuracy. Meanwhile, our method preserves all the dynamics of probed phonon mode from a particular state, which means it is capable of studying the transient thermal transport processes in a nonequilibrium system. Besides the exceptional efficiency, our method also comes with freedom to choose to probe only those modes of interest. This makes it ideal for use with large systems and in computationally demanding applications, such as ab initio MD simulations. Moreover, the PT method we propose here is very straightforward and easy to implement.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

SNS Warm Linac Circulator Breakdown Considerations for the PPU Project

Multipacting in accelerating structures is a complex phenomenon about which there is much to be understood. While multipacting research efforts have primarily been focused on superconducting radio frequency (SRF) systems, normal conducting accelerating structures which have a higher thermal capacity, and a greater vacuum pressure tolerance could benefit from additional investigation. This research details multipacting simulation methods and the results of 3-D electromagnetic simulations of RF vacuum windows used on normal conducting linac (NCL) cavities. Benchmarking of the peak electric fields in these structures, benefits of material processing and possible techniques for reducing or eliminating multipacting activities are discussed.

Toby, George↗
DOE OSTI2021

Multipacting Analysis of Warm Linac RF Vacuum Windows

Multipacting in accelerating structures is a complex phenomenon with which there is much to be understood. While multipacting research efforts have primarily been focused on superconducting radio frequency (SRF) systems, normal conducting accelerating structures that have a higher thermal capacity and a greater vacuum pressure tolerance could benefit from additional investigation. This research details multipacting simulation methods and the results of 3-D electromagnetic simulations of RF vacuum windows used on normal conducting linac (NCL) cavities. Possible techniques for reducing and eliminating multipacting activities in these structures are discussed.

Toby, George↗
DOE OSTI2020

A framework to evaluate IMEX schemes for atmospheric models

We present a new evaluation framework for implicit and explicit (IMEX) Runge–Kutta time-stepping schemes. The new framework uses a linearized nonhydrostatic system of normal modes. We utilize the framework to investigate the stability of IMEX methods and their dispersion and dissipation of gravity, Rossby, and acoustic waves. We test the new framework on a variety of IMEX schemes and use it to develop and analyze a set of second-order low-storage IMEX Runge–Kutta methods with a high Courant–Friedrichs–Lewy (CFL) number. We show that the new framework is more selective than the 2-D acoustic system previously used in the literature. Schemes that are stable for the 2-D acoustic system are not stable for the system of normal modes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Conductive Thin Films over Large Areas by Supramolecular Self-Assembly

We report a “one-step” method for preparing conductive thin films with cylindrical microdomains oriented normal to the surface over large areas using the supramolecular assembly of poly(styrene-block-4-vinylpyridine) (PS19-b-P4VP 5 ) and 5,10,15,20-tetrakis(4-hydroxyphenyl)-21H,23H-porphine (HOTPP). HOTPP interacts with the P4VP block by hydrogen bonding between the hydroxyl group of HOTPP and pyridine ring of PS 19 -b-P4VP 5 , forming cylindrical P4VP(HOTPP) domains having an average diameter of ~17 nm in a PS matrix. Here, dynamic light scattering, contact angle, and in situ grazing incidence small-angle X-ray scattering measurements show a morphological transition from spherical micelles in solution to cylindrical microdomains oriented normal to the substrate surface during the drying process. From the dependence of current on voltage, an average current of ~4.0 nA is found to pass through a single microdomain, pointing to a promising route for organic semiconductor device applications.

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

Statistical White-Line Analysis in High-Throughput TXM-XANES for Chemical State Quantification

The transmission X-ray microscopy (TXM) based X-ray absorption near-edge structure (XANES) technique provides three-dimensional mapping of element-specific chemical states at nanometer-scale spatial resolution and micrometer-scale fields of view. However, compared to conventional volume-averaged XANES (VA-XANES) measurements, the inherently small voxel size in TXM-XANES leads to a lower signal-to-noise ratio, making full-spectrum analysis computationally demanding and less robust. Here, we present the structural and compositional conditions for a statistical white-line analysis framework under which chemical state information can be directly extracted from the white-line peak position in voxel spectra without the need for voxel-wise background subtraction or normalization, under well-defined structural and compositional conditions. The method is validated on layered oxide cathode materials, where low-order polynomial fitting accurately reproduces white-line features, and the extracted energy distributions correlate strongly with VA-XANES results. This statistical approach enables high-throughput, dose-efficient, and noise-robust chemical state quantification in TXM-XANES, offering broad applicability to functional materials requiring nanoscale oxidation-state mapping.

TXM↗
DOE OSTI2019

First-principles electronic structure in second-moment calculation of mode frequencies: Failure of quasiharmonic approximation in silicon

Using first-principles electronic structure calculations (local density approximation), we present second-moment calculations of the temperature dependence of the vibrational normal mode frequencies in Si up to 1500 K . The method is based on simple ensemble averages of displacements and forces and is easy to implement. These efficient calculations incorporate all interactions nonperturbatively, implicitly including all high order force constants, in contrast to typical perturbative approaches. Also in this work, we propose and apply to our classical calculations an estimated quantum correction that could be used for classical approaches other than the moments method. We compare our results to available experiments, including recent inelastic neutron scattering data and conclude that the second-moment gives the frequency shift to within a factor of 2 in contrast to quasiharmonic approaches that fail by an order of magnitude to predict the anharmonic shift. Noting the relative simplicity and speed of the calculations, we conclude that the second-moment approximation reasonably accounts for anharmonicity in a system that has been very closely studied experimentally and may reliably be applied in an exploratory way to more complex systems. Finally, we suggest that further improvement on the results could be obtained by incorporating the fourth moment, which can be implemented with little additional cost.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Quantum Monte Carlo calculations in configuration space with three-nucleon forces

Neutron matter, through its connection to neutron stars as well as systems like cold atom gases, is one of the most interesting yet computationally accessible systems in nuclear physics. The Configuration-Interaction Monte Carlo (CIMC) method is a stochastic many-body technique allowing to tackle strongly coupled systems. In contrast to other Quantum Monte Carlo methods employed in nuclear physics, the CIMC method can be formulated directly in momentum space allowing for an efficient use of non-local interactions. In this work we extend CIMC method to include three-nucleon interactions through the normal-ordered two-body approximation. Here, we present results for the equation of state of neutron matter in line with other many-body calculations that employ low resolution chiral interactions, and provide predictions for the momentum distribution and the static structure factor.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Global quark spin correlations in relativistic heavy ion collisions

The observation of the vector meson’s global spin alignment by the STAR Collaboration reveals that strong spin correlations may exist for quarks and antiquarks in relativistic heavy-ion collisions in the normal direction of the reaction plane. We propose a systematic method to describe such correlations in the quark matter. We classify them as local and long range quark spin correlations in the system. We show in particular that the effective quark spin correlations contain the genuine spin correlations originated directly from the dynamical process as well as those induced by averaging over other degrees of freedom. We also show that such correlations can be studied by measuring the vector meson’s spin density matrix and hyperon-hyperon and hyperon-anti-hyperon spin correlations. We present the relationships between these measurable quantities and spin correlations of quarks and antiquarks. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Bilevel Nodal Behind-the-meter Solar Disaggregation Under Unexpected Extreme Weather Conditions

As the power grid undergoes significant paradigm shift due to the increasing penetration of renewable generation, the ever-growing installation of behind-the-meter (BTM) solar generation in the power grid also has a significant impact on nodal loads, posing challenges on transmission operators. Furthermore, increasing frequent and severe extreme weather events intertwine with ubiquitous BTM solar generations and have amplified the challenges of accurately model nodal load profiles, especially under the lack of ground-truth information for verification. To tackle these challenges, this paper introduces a bilevel model that utilizes year-long data (e.g., proxy solar, zonal load, and individual node load profiles) to disaggregate metered profiles into actual demand and BTM solar generation at each transmission node. The proxy solar not only scales the BTM solar generation of individual nodes but also create a compensation term for enhancing performance on days with unexpected extreme weather events. The proposed algorithm is validated with real-world PJM Interconnection data during unexpected events like the recent Winter Storm Elliott. For quantitative evaluations, a novel Score error is introduced, which is based on mean percentages and load scales and offers a universal assessment method suitable for all nodes and different data formats (e.g., normalized or raw values).

behind-the-meter solar, load disaggregations, load↗
DOE OSTI2021

Baseline [ 18 F]GTP1 tau PET imaging is associated with subsequent cognitive decline in Alzheimer’s disease

Background: The role and implementation of tau PET imaging for predicting subsequent cognitive decline in Alzheimer’s disease (AD) remains uncertain. This study was designed to evaluate the relationship between baseline [ 18 F]GTP1 tau PET and subsequent longitudinal change across multiple cognitive measures over 18 months. Methods: Our analyses incorporated data from 67 participants, including cognitively normal controls (n = 10) and β-amyloid (Aβ)-positive individuals ([ 18 F] florbetapir Aβ PET) with prodromal (n = 26), mild (n = 16), or moderate (n = 15) AD. Baseline measurements included cortical volume (MRI), tau burden ([ 18 F]GTP1 tau PET), and cognitive assessments [Mini-Mental State Examination (MMSE), Clinical Dementia Rating (CDR), 13-item version of the Alzheimer’s Disease Assessment Scale-Cognitive Subscale (ADAS-Cog13), and Repeatable Battery for the Assessment of Neuropsychological Status (RBANS)]. Cognitive assessments were repeated at 6-month intervals over an 18-month period. Associations between baseline [ 18 F]GTP1 tau PET indices and longitudinal cognitive performance were assessed via univariate (Spearman correlations) and multivariate (linear mixed effects models) approaches. The utility of potential prognostic tau PET cut points was assessed with ROC curves. Results: Univariate analyses indicated that greater baseline [ 18 F]GTP1 tau PET signal was associated with faster rates of subsequent decline on the MMSE, CDR, and ADAS-Cog13 across regions of interest (ROIs). In multivariate analyses adjusted for baseline age, cognitive performance, cortical volume, and Aβ PET SUVR, the prognostic performance of [ 18 F]GTP1 SUVR was most robust in the whole cortical gray ROI. When AD participants were dichotomized into low versus high tau subgroups based on baseline [ 18 F]GTP1 PET standardized uptake value ratios (SUVR) in the temporal (cutoff = 1.325) or whole cortical gray (cutoff = 1.245) ROIs, high tau subgroups demonstrated significantly more decline on the MMSE, CDR, and ADAS-Cog13. Conclusions: Our results suggest that [ 18 F]GTP1 tau PET represents a prognostic biomarker in AD and are consistent with data from other tau PET tracers. Tau PET imaging may have utility for identifying AD patients at risk for more rapid cognitive decline and for stratification and/or enrichment of participant selection in AD clinical trials.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Analysis of spent nuclear fuel multipurpose canister dynamics during rail transportation

The primary mode of spent nuclear fuel transportation within the United States will be by railcar. One such system is the Atlas railcar, which is designed to transport 17 different spent nuclear fuel cask systems, including bare fuel systems and canister fuel systems. In the latter configuration, multipurpose canisters containing spent nuclear fuel may be placed within an overpack for storage, or within a cask for transportation. Compared to bare fuel systems, canister fuel systems have additional degrees of freedom for motion during transportation, because clearance between the cask and canister allows for some motion of the canister to occur relative to the cask. This work investigates the effect of canister motion on the shock and vibration imparted to the spent nuclear fuel within. Structural dynamic analyses have been conducted to identify the effects of canister to cask clearance, presence and type of dunnage, and loading direction and frequency. This modeling study calculates anticipated cask motion, canister motion, and spent nuclear fuel structural dynamic response to normal conditions of transportation railcar motion using finite element analysis methods that were developed to model the rail segment of the ENSA/DOE (Equipos Nucleares S.A., U.S. Department of Energy) multimodal transportation test of 2017.

SNF, UNF, LS-DYNA, dynamics, kinematics, shock, vi↗