Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Physics Regularization”

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 163 records · Page 9

Domain Aware Deep-learning Algorithms Integrated with Scientific-computing Technologies (DADAIST)

This technical report summarized the contribution of the DADAIST project funded by the Data Model Convergence Initiative via the Laboratory Directed Research and Development (LDRD) investments at Pacific Northwest National Laboratory (PNNL). Specifically, we report the development of the NeuroMANCER (Neural Modules with Adaptive Nonlinear Constraints and Efficient Regularizations), a new open-source Scientific Machine Learning library for formulating and solving parametric constrained optimization problems, physics-informed system identification, and parametric optimal control problems. NeuroMANCER is using differentiable programming to combine modern data-driven models and optimization modeling language into a coherent algorithmic and software framework. NeuroMANCER is a Pytorch-based framework and adopts much of its philosophy focused on research and development, rapid prototyping, and streamlined deployment. Strong emphasis is given to extensibility, interoperability with the PyTorch ecosystem, and quick adaptability to custom domain problems. Neuromancer repository contains a comprehensive library of differentiable modules, including custom activation functions, matrix factorizations, deep learning architectures, neural differential equations, differential equation solvers, implicit layers such as iterative solvers, high-level API for symbolic expressions, API for modeling and control of dynamical systems, and extensive set of tutorial code examples in the form of python scripts and jupyter notebooks.

97 MATHEMATICS AND COMPUTING↗

An interpretation and solution of ill-conditioned linear equations

Data insufficiency, poorly conditioned matrices and singularities in equations occur regularly in complex optimization, correlation, and interdisciplinary model studies. This work concerns itself with two methods of obtaining certain physically realistic solutions to ill-conditioned or singular algebraic systems of linear equations arising from such studies. Two efficient computational solution procedures that generally lead to locally unique solutions are presented when there is insufficient data to completely define the model, or a least-squares error formulation of this system results in an ill-conditioned system of equations. If it is assumed that a reasonable estimate of the uncertain data is available in both cases cited above, then we shall show how to obtain realistic solutions efficiently, in spite of the insufficiency of independent data. The proposed methods of solution are more efficient than singular-value decomposition for dealing with such systems, since they do not require solutions for all the non-zero eigenvalues of the coefficient matrix.

Ojalvo, I. U.↗

Modifications to Axially Symmetric Simulations Using New DSMC (2007) Algorithms

Several modifications aimed at improving physical accuracy are proposed for solving axially symmetric problems building on the DSMC (2007) algorithms introduced by Bird. Originally developed to solve nonequilibrium, rarefied flows, the DSMC method is now regularly used to solve complex problems over a wide range of Knudsen numbers. These new algorithms include features such as nearest neighbor collisions excluding the previous collision partners, separate collision and sampling cells, automatically adaptive variable time steps, a modified no-time counter procedure for collisions, and discontinuous and event-driven physical processes. Axially symmetric solutions require radial weighting for the simulated molecules since the molecules near the axis represent fewer real molecules than those farther away from the axis due to the difference in volume of the cells. In the present methodology, these radial weighting factors are continuous, linear functions that vary with the radial position of each simulated molecule. It is shown that how one defines the number of tentative collisions greatly influences the mean collision time near the axis. The method by which the grid is treated for axially symmetric problems also plays an important role near the axis, especially for scalar pressure. A new method to treat how the molecules are traced through the grid is proposed to alleviate the decrease in scalar pressure at the axis near the surface. Also, a modification to the duplication buffer is proposed to vary the duplicated molecular velocities while retaining the molecular kinetic energy and axially symmetric nature of the problem.

Liechty, Derek S.↗

Physics-Infused AI/ML Based Digital-Twin Framework for Flow-Induced-Vibration Damage Prediction in a Nuclear Reactor Heat Exchanger

This report summarizes some of the ongoing work related to the development of an expert-elicitation-digital-twin framework for real time damage state prediction in heat exchanger components of a nuclear reactor. The framework is targeted towards predicting damage associated with coupled low cycle fatigue (associated with regular heat-up, cool-down and power operation transients) and high cycle fatigue (associated with flow induced vibration transients). The overall framework will be based on a NoSQL based database, physics-infused-geometry-dependent virtual-sensor data, different AI/ML techniques-based data-driven-predictive-model applications (Apps) and real-time plant sensor measurements available through few existing sensors. Towards this overall goal, this report updates some of the ongoing work, such as on implementation of a NoSQL Database (such as MongoDB), FE based heat transfer analysis of a heat exchanger (e.g. of a PWR steam generator) for generating geometry-dependent virtual sensor data and evaluation of various AI/ML models such as based on multivariate linear regression, ensembled decision-tree based Random-Forest and Gradient-Boosting regression and high-dimensional-kernel-function-transformation based Support-Vector-Machine regression models. The AI/ML models were evaluated for predicting multi-time-series thermal states at thousands of 3D point-clouds

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Interacting fractons in 2+1-dimensional quantum field theory

We analyze, in perturbation theory, a theory of weakly interacting fractons and non-relativistic fermions in a 2+1 dimensional Quantum Field Theory. In particular we compute the 1-loop corrections to the self energies and interaction vertex, and calculate the associated 1-loop Renormalization Group flows of the coupling constants. Surprisingly, we find that the fracton-fermion coupling does not flow due to an emergent coordinate-dependent symmetry of the effective Lagrangian, making this model a well-defined quantum field theory. We provide additional discussions on the regularization and renormalization of interacting fractonic theories, as well as both qualitative and quantitative remarks regarding the theory at finite temperature and finite chemical potential.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Quantum calculations of the cavity shift in electron magnetic moment measurements

The measurement of the anomalous electron magnetic moment g – 2 through quantum transitions of a single trapped electron is the most stringent test of quantum field theory. These experiments are now so precise that they must account for the effects of the cavity containing the electron. Classical calculations of this “cavity shift” must subtract the electron’s divergent self-field, and thus require knowledge of the exact Green’s function for the cavity’s electromagnetic field. We perform the first fully quantum calculation of the cavity shift in a closed cavity, which instead involves subtracting linearly divergent cavity mode sums and integrals. Using contour integration methods, we find perfect agreement with existing classical results for both spherical and cylindrical cavities, justifying their current use. Moreover, our mode-based results can be naturally generalized to account for systematic effects, necessary to push future measurements to the next order of magnitude in precision.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Flow patterns near a conical sonic line

A study is made in the physical plane of the flow properties near a conical cross flow sonic line. The shapes of the Mach lines at regular and singular sonic points are presented. Through a study of the jump conditions at a singular sonic point, the possibility of a double sonic line is discussed. The properties of bubbles of supersonic cross flow are examined. Numerical results are presented to support some of the theoretical observations.

Salas, M. D.↗

Evaluation of the DAO Retrospective Data Assimilation System

We have developed and implemented a retrospective data assimilation system (RDAS) as an upgrade to the operational DAO/Terra data assimilation system. This formulation aims at improving analysis over filter analysis by the dynamically consistent incorporation of observation information past a given analysis time. The current implementation of the RDAS uses the adjoint of the tangent linear model of a simplified version of the Terra general circulation model and extensions to the physical-space statistical analysis system to propagate observation information back in time. The RDAS adopts the same assumptions of the regular data assimilation system, particularly, no explicit propagation of error covariances are involved therefore rendering a procedure that is computationally affordable. In this study, we show results of experiments conducted to investigate the performance of the 6-hour (lag-1) RDAS. Statistical results obtained over one month during a winter season indicate that the RDAS represents considerable improvement over the regular assimilation. Plans for implementation of the RDAS capability in our new finite-volume data assimilation system will also be presented at the time of the conference.

Zhu, Yanqiu↗

Understanding Evolution of Lithium Trivanadate Cathodes During Cycling via Reformulated Physics-Based Models and Experiments

Degradation of lithium trivanadate ( Li x V 3 O 8 ) cathodes has been widely reported in the literature, but studies have offered little insight towards developing a detailed understanding of the evolution of the active material, and have been inconclusive as to the root cause of degradation. Here, we refit a validated physics-based model to discharge curves over the course of cycling at C/5, and use the evolution of the model parameters to track evolution of the cathode. A regularization penalty for adjusting model parameters from the validated model is introduced as a framework to identify which model parameters can explain a significant portion of the observed change in the voltage profile over the course of cycling. SEM reveals that lithium trivandate particles fracture upon cycling at C/5, consistent with the results of the parameter estimation, deactivation of lithium trivanadate and faster diffusion of lithium within the active particles. The physics-based model is then used to design modified cycling protocols which identify the phase transformation during discharge of lithium trivanadate as the driver of the particle fracture and capacity fade.

25 ENERGY STORAGE↗

Physics of beam-driven ion cyclotron emission in the large plasma device

Abstract Ion cyclotron emission (ICE) is widely observed from spatially localised minority energetic ion populations in toroidal magnetically confined fusion (MCF) plasmas, both tokamaks and stellarators. Its spectral structure is typically regular with narrow suprathermal peaks, whose frequency separation matches a local energetic ion cyclotron frequency. Here we report the first computational (fully nonlinear self-consistent kinetic particle-in-cell code) and analytical (linear magnetoacoustic cyclotron instability (MCI)) studies of ICE observations from cylindrical plasmas contained in the Large Plasma Device (LAPD). Because LAPD is cylindrical, the plasma physics giving rise to the observed ICE spectrum necessarily excludes toroidal effects. Our approach, previously successful for toroidal plasma ICE, assumes slab geometry and hence is well adapted to LAPD. ICE from LAPD is strongly electrostatic, as distinct from electromagnetic, and is driven by 15 keV beam ions for which the ratio of perpendicular speed to the local Alfven speed, v ⊥ / v A , is 0.15, lower than in MCF plasmas from which beam-driven ICE has previously been observed. Our results are in good agreement with these observations. There is congruence between simulated ICE spectra, obtained in the saturated nonlinear regime of our computations, and observed ICE spectra; the underlying physics is essentially the same as in toroidal plasmas; and there is alignment with linear analytical theory where appropriate. The present work establishes a mapping from the cylindrical LAPD ICE observations to toroidal MCF ICE observations. The LAPD spectra are instances of beam-driven MCI-type ICE in its sub-Alfvenic, predominantly electrostatic manifestation, which has precedents in MCF stretching back to the 1990s. An interesting corollary is that, for many purposes, ICE in toroidal MCF plasmas ‘might as well’ be occurring in a cylinder.

Samant, O. (ORCID:0000000226055363)↗

Quantile regression-enriched event modeling framework for dropout analysis in high-temperature superconductor manufacturing

High-temperature superconductor (HTS) tapes have shown promising characteristics of high critical current, which are prerequisites for applications in high-field magnets. Due to the unstable growth conditions in the HTS manufacturing process, however, the frequent occurrences of dropouts in the critical current impede the consistent performance of HTS tapes. To manufacture HTS tapes with large scale, high yield, and uniform performance, it is essential to develop novel data analysis approaches for modeling the dropouts and identifying the related important process parameters. Conventional methods for modeling recurrent events, such as the point process, require the extraction of events from quality measurements. As the critical current is a continuous process, it may not comprehensively represent the drop patterns by transforming the time-series measurements into a set of events. Here, to solve this issue, we develop a novel quantile regression-enriched event modeling (QREM) framework that integrates the non-homogeneous Poisson process for modeling the occurrence of dropouts and the quantile regression for capturing the drop patterns. By incorporating the feature selection and regularization, the proposed framework identifies a set of significant process parameters that can potentially cause the dropouts of HTS tapes. The proposed method is tested on real HTS tapes produced using an advanced manufacturing process, successfully identifying important parameters that influence dropout events including the substrate temperature and voltage. The results demonstrate that the proposed QREM method outperforms the standard point process in predicting the occurrence of dropouts.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

G2Aero (Separable shape tensors for aerodynamic design) [SWR-22-44]

G2Aero is a flexible and practical tool for design and deformation of 2D airfoils and 3D blades using data-driven approaches. G2Aero utilizes the geometry of matrix manifolds—specifically the Grassmannian—to build a novel framework for representing physics-based separable deformations of shapes. G2Aero offers the flexibility to generate perturbations in a customizable way over any portion of the blade. The G2Aero framework utilizes data-driven methods based on a curated database of physically relevant airfoils. Specific tools include: - principal geodesic analysis over normal coordinate neighborhoods of matrix manifolds; - a variety of data-regularized deformations to nominal 2D airfoil shapes; - Riemannian interpolation connecting a sequence of airfoil cross-sections to build 3D blades from 2D data; - consistent perturbations over the span of interpolated 3D blades based on dominant modes from the data-driven analysis.

Doronina, Olga↗

The fundamental constants of orthotropic affine plate/slab equations

The global constants associated with orthotropic slab/plate equations are discussed, and the rotational behavior of the modulus/compliance components associated with orthotropic slabs/plates are addressed. It is concluded that one cluster constant is less than or equal to unity for all physically possible materials. Rotationally anomalous behavior is found in two materials, and a simple inequality which can be used to identify regular or anomalous behavior is presented and discussed in detail.

Brunelle, E. J.↗

Loop corrections in Minkowski spacetime away from equilibrium. Part II. Finite-time results

Loop corrections to finite-time correlation functions in quantum field theories away from equilibrium can be calculated using the in-in path integral approach. In this paper, we calculate the unequal-time two-point correlator for different massless self-interacting scalar quantum field theories on a Minkowski background, starting the field evolution at an arbitrary initial time. We find the counterterms that need to be added to UV-renormalize the result, including usual in-out counterterms in the dynamics and additional initial state counterterms that are required to cancel all UV divergences. We find that the late-time limit of the renormalized correlation function exhibits a linear or logarithmic growth in time, depending on whether the interaction strength is dimension-one or dimensionless, respectively. The late-time correlations match those obtained in our companion paper and, as shown there, the divergences do not indicate a real IR issue, consistent with what one would expect in Minkowski.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Polyakov’s confinement mechanism for generalized Maxwell theory

We study fractional-derivative Maxwell theory, as appears in effective descriptions of, for example, large N f QED 3 , graphene, and some types of surface defects. We argue that when the theory is realized on a lattice, monopole condensation leads to a confining phase via the Polyakov confinement mechanism.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

On ambiguities and divergences in perturbative renormalization group functions

There is an ambiguity in choosing field-strength renormalization factors in the $ \overline{\mathrm{MS}} $ scheme starting from the 3-loop order in perturbation theory. More concerning, trivially choosing Hermitian factors has been shown to produce divergent renormalization group functions, which are commonly understood to be finite quantities. We demonstrate that the divergences of the RG functions are such that they vanish in the RG equation due to the Ward identity associated with the flavor symmetry. It turns out that any such divergences can be removed using the renormalization ambiguity and that the use of the flavor-improved β-function is preferred. We show how our observations resolve the issue of divergences appearing in previous calculations of the 3-loop SM Yukawa β-functions and provide the first calculation of the flavor-improved 3-loop SM β-functions in the gaugeless limit.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Online regularization of Poincaré map of storage rings with Shannon entropy

A measurable chaos indicator is used as the online optimization objective in tuning a complicated nonlinear system—the National Synchrotron Light Source-II storage ring. Through analyzing the Shannon entropy in measured Poincaré maps, not only can the commonly used nonlinear characterizations be extracted, but more importantly, the chaos can be quantified and then used for an online regularization of these maps. The method itself is general and applicable to other tunable nonlinear systems as well. Published by the American Physical Society 2025

36 MATERIALS SCIENCE↗