Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Weight function”

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 415 records · Page 23

Higher-order particle representation for particle-in-cell simulations

In this paper we present an alternative approach to the representation of simulation particles for unstructured electrostatic and electromagnetic PIC simulations. In our modified PIC algorithm we represent particles as having a smooth shape function limited by some specified finite radius, r 0 . A unique feature of our approach is the representation of this shape by surrounding simulation particles with a set of virtual particles with delta shape, with fixed offsets and weights derived from Gaussian quadrature rules and the value of r 0 . As the virtual particles are purely computational, they provide the additional benefit of increasing the arithmetic intensity of traditionally memory bound particle kernels. The modified algorithm is implemented within Sandia National Laboratories' unstructured EMPIRE-PIC code, for electrostatic and electromagnetic simulations, using periodic boundary conditions. We show results for a representative set of benchmark problems, including electron orbit, a transverse electromagnetic wave propagating through a plasma, numerical heating, and a plasma slab expansion. In this work, good error reduction across all of the chosen problems is achieved as the particles are made progressively smoother, with the optimal particle radius appearing to be problem-dependent.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Assembly of Close–Packed Ferroelectric Polymer Nanowires via Interface–Epitaxy with ReS 2

The flexible, transparent, and low-weight nature of ferroelectric polymers makes them promising for wearable electronic and optical applications. To reach the full potential of the polarization-enabled device functionalities, large-scale fabrication of polymer thin films with well-controlled polar directions is called for, which remains a central challenge. The widely exploited Langmuir–Blodgett, spin-coating, and electrospinning methods only yield polymorphous or polycrystalline films, where the net polarization is compromised. Here, an easily scalable approach is reported to achieve poly(vinylidene fluoride-trifluoroethylene) P(VDF-TrFE) thin films composed of close-packed crystalline nanowires via interface-epitaxy with 1T'-ReS 2 . Upon controlled thermal treatment, uniform P(VDF-TrFE) films restructure into about 10 and 35 nm-wide (010)-oriented nanowires that are crystallographically aligned with the underlying ReS 2 , as revealed by high-resolution transmission electron microscopy. Piezoresponse force microscopy studies confirm the out-of-plane polar axis of the nanowire films and reveal coercive voltages as low as 0.1 V. Reversing the polarization can induce a conductance switching ratio of >10 8 in bilayer ReS 2 , over six orders of magnitude higher than that achieved by an untreated polymer gate. Lastly, this study points to a cost-effective route to large-scale processing of high-performance ferroelectric polymer thin films for flexible energy-efficient nanoelectronics.

crystalline nanowires↗

A novel ignition model for low velocity impact of heterogeneous explosives based on interacting hot spots

While numerous studies have focused on the ignition of explosives occurring in high velocity impact and the associated shock-to-detonation transition, there has been growing interest in developing computational models focused on low-velocity impact regimes. A predictive low-velocity impact ignition model will be important for analyzing high explosive safety and potential accident scenarios. This work introduces a novel ignition model based on the concept of thermally interacting hot spots to simulate low velocity impacted heterogeneous explosives where observed ignition times are on the order of milliseconds. The model asserts that relevant hot spots are micron-sized, the typical separation between neighboring hot spots is on the order of a hundred microns, and that neighbors interact thermally through heat conduction across the interstitial region between them. To achieve tractable numerical solutions, hot spots are assumed to form a periodic array as opposed to the highly irregular positioning in an actual explosive. This idealization allows a single two hotspot system to characterize the ignition process. Consequently, the model is referred to as the two hot spot Frank-Kamenetskii ignition model. In the present study, hot spots are modeled as constant heat sources terms, but this can be extended to include grain-scale phenomena like frictional heating of micron-sized growing cracks that are confined under high pressure. Because the micron-sized features are below the scale that can be efficiently resolved at a systems level, an efficient subscale scheme based on the Method of Weighted Residuals (MWR) is used to efficiently solve the equations. In conclusion, we carry out numerical examples and analytic predictions illustrating the accuracy and the functioning of the model.

97 MATHEMATICS AND COMPUTING↗

Using ground-based measurements to recover the spectra of radiation escaping from distant light-pollution sources

ABSTRACT A professional dark sky assessment normally requires the numerical modelling of light propagation from ground-based sources into the atmosphere. The spectral power distribution (SPD) of light escaping from an artificial source is one of the key parameters needed in predicting models to simulate the optical signatures of the night sky as a function of distance. SPD at zero distance is due to a non-trivial superposition of light beams, but does not necessarily coincide with the weighted SPD of individual lights. This is because most photons can be removed by obstacles or can interact with many surfaces before being directed into the atmosphere. We have developed an experimental and theoretical method for obtaining the initial SPD from ground-based spectroscopy of the night sky. The method is applicable to distant sources of light and is demonstrated in field measurements conducted on the light dome of Chicago with a new sky glow spectrometer. We show that the method requires analysing the relative impact of local light sources scattered around the measuring site on the radiance of the light dome. Our theoretical method along with ground-based spectroscopy represents a new approach to characterizing light-pollution sources which are otherwise difficult to study due to lack of information on either light-source inventory or a heterogeneous 3D structure of light-emitting and blocking environments.

Kocifaj, M.↗

Emergence of topological and trivial interface states in VSe 2 films coupled to Bi 2 Se 3

Coupling ordinary metals with topological Bi 2 Se 3 can instigate the long-range migration of the spin-polarized Dirac states. Instead, for trivial metallic VSe 2 films on Bi 2 Se 3 , topological and trivial Rashba-type interface states emerge, each strongly localized at the VSe 2 /Bi 2 Se 3 interface. Their rapidly decaying spectral weights are uncovered by thickness-dependent band mappings of VSe 2 films and replicated by a phenomenological first-principles model of the spectral function. In this work, the results are pertinent to spintronic devices requiring spin transmission across metal/topological insulator interfaces.

36 MATERIALS SCIENCE↗

Binary pseudo-random array (BPRA) for inspection and calibration for cylindrical wavefront interferometry

High-accuracy metrology is vitally important in manufacturing ultra-high-quality free-form mirrors designed to manipulate X-ray light with nanometer-scale wavelengths. However, surface topography measurements are instrument dependent, and without the knowledge of how the instrument performs under the practical usage conditions, the measured data contain some degree of uncertainty. Binary Pseudo Random Array (BPRA) “white noise” artifact are effective and useful for characterizing the Instrument Transfer Function (ITF) of surface topography metrology tools and wavefront measurement instrument. BPRA artifact contains features with all spatial frequencies in the instrument bandpass with equal weight. As a result, power spectral density of the patterns has a deterministic white-noise-like character that allows direct determination of the ITF with uniform sensitivity over the entire spatial frequency range. The application examples include electron microscopes, x-ray microscopes, interferometric microscopes, and large field-of-view Fizeau Interferometers. Furthermore, we will introduce the application of BPRA method to characterizing the ITF of Cylindrical Wavefront Interferometry (CWI), by developing the BPRA artifact which matches the radius of curvature of the cylindrical wavefront. The data acquisition and analysis procedures for different applications of the ITF calibration technique developed are also discussed.

Munechika, K↗

A Novel Feature of Valence Quark Distributions in Hadrons

Examining the evolution of the maximum of valence quark distribution, qV, weighted by Bjorken x, h(x,t)≡xqV(x,t), it is observed that h(x,t) at the peak becomes a one-parameter function; h(xp,t)=Φ(xp(t)), where xp is the position of the peak, t=logQ2, and Q2 is the resolution scale. This observation is used to derive a new model-independent relation which connects the partial derivative of the valence parton distribution functions (PDFs) in xp to the quantum chromodynamics (QCD) evolution equation through the xp derivative of the logarithm of the function Φ(xp(t)). A numerical analysis of this relation using empirical PDFs results in an observation of the exponential form of the Φ(xp(t))=h(xp,t)=CeDxp(t) for leading to next-to-next leading order approximations of PDFs for the range of Q2, covering four orders in magnitude. The exponent, D, of the observed “height-position” correlation function converges with the increase in the order of approximation. This result holds for all the PDF sets considered. A similar relation is observed also for the pion valence quark distribution, indicating that the obtained relation may be universal for any non-singlet partonic distribution. The observed “height-position” correlation is used also to indicate that no finite number of exchanges can describe the analytic behavior of the valence quark distribution at the position of the peak at fixed Q2.

Leon, Christopher↗

Origin of light-induced metastability in ZrTe 5

Here we study the nonequilibrium electronic structure of a model Dirac semimetal ZrTe 5 by using time-and-angle resolved photoemission spectroscopy and density functional theory–based electron and phonon calculations. By measuring the electronic dispersion near the Γ point at time delays up to 10 picoseconds, we discovered that the band spectral weight does not recover during the measured temporal window, revealing the existence of a light-induced metastable state in the electronic structure of this material. Our calculations find that the photoexcited A 1⁢g phonon mode leads to a band renormalization that both supports our experimental observations at the zone center and predicts changes to the band structure outside of our experimental window, ultimately showing the evolution from a direct to an indirect gap semimetal; such band renormalization dramatically reduces the electron-hole recombination rate giving rise to the metastability in this system.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Train Like a (Var)Pro: Efficient Training of Neural Networks with Variable Projection

Deep neural networks (DNNs) have achieved state-of-the-art performance across a variety of traditional machine learning tasks, e.g., speech recognition, image classification, and segmentation. The ability of DNNs to efficiently approximate high-dimensional functions has also motivated their use in scientific applications, e.g., to solve partial differential equations and to generate surrogate models. In this paper, we consider the supervised training of DNNs, which arises in many of the above applications. We focus on the central problem of optimizing the weights of the given DNN such that it accurately approximates the relation between observed input and target data. Devising effective solvers for this optimization problem is notoriously challenging due to the large number of weights, nonconvexity, data sparsity, and nontrivial choice of hyperparameters. To solve the optimization problem more efficiently, we propose the use of variable projection (VarPro), a method originally designed for separable nonlinear least-squares problems. Our main contribution is the Gauss--Newton VarPro method (GNvpro) that extends the reach of the VarPro idea to nonquadratic objective functions, most notably cross-entropy loss functions arising in classification. These extensions make GNvpro applicable to all training problems that involve a DNN whose last layer is an affine mapping, which is common in many state-of-the-art architectures. In our four numerical experiments from surrogate modeling, segmentation, and classification, GNvpro solves the optimization problem more efficiently than commonly used stochastic gradient descent (SGD) schemes. Finally, GNvpro finds solutions that generalize well, and in all but one example better than well-tuned SGD methods, to unseen data points.

97 MATHEMATICS AND COMPUTING↗

Robust scalable initialization for Bayesian variational inference with multi-modal Laplace approximations

Predictive modeling typically relies on Bayesian model calibration to provide uncertainty quantification. Variational inference utilizing fully independent (“mean-field”) Gaussian distributions are often used as approximate probability density functions. This simplification is attractive since the number of variational parameters grows only linearly with the number of unknown model parameters. However, the resulting diagonal covariance structure and unimodal behavior can be too restrictive to provide useful approximations of intractable Bayesian posteriors that exhibit highly non-Gaussian behavior, including multimodality. High-fidelity surrogate posteriors for these problems can be obtained by considering the family of Gaussian mixtures. Gaussian mixtures are capable of capturing multiple modes and approximating any distribution to an arbitrary degree of accuracy, while maintaining some analytical tractability. Unfortunately, variational inference using Gaussian mixtures with full-covariance structures suffers from a quadratic growth in variational parameters with the number of model parameters. The existence of multiple local minima due to strong nonconvex trends in the loss functions often associated with variational inference present additional complications, These challenges motivate the need for robust initialization procedures to improve the performance and computational scalability of variational inference with mixture models. In this work, we propose a method for constructing an initial Gaussian mixture model approximation that can be used to warm-start the iterative solvers for variational inference. The procedure begins with a global optimization stage in model parameter space. In this step, local gradient-based optimization, globalized through multistart, is used to determine a set of local maxima, which we take to approximate the mixture component centers. Around each mode, a local Gaussian approximation is constructed via the Laplace approximation. Finally, the mixture weights are determined through constrained least squares regression. The robustness and scalability of the proposed methodology is demonstrated through application to an ensemble of synthetic tests using high-dimensional, multimodal probability density functions. Here, the practical aspects of the approach are demonstrated with inversion problems in structural dynamics.

97 MATHEMATICS AND COMPUTING↗

Effect of Confinement on the Density and Diffusivity of Organic Molecules in Single-Digit Nanopores Relative to Bulk Fluid Behavior

Prediction of fluid properties in microporous and mesoporous confinement is non-trivial due to distinct influences from pore surface functionality and shape and fluid molecular characteristics. The vast majority of literature focuses nearly exclusively on a small number of components that are generally low molecular weight gases and vapors. Missing are studies that utilize a relatively large number of compounds that are liquid at ambient conditions to probe the impact of molecular shape and size on confined fluid properties. Such materials are of great interest in the extraction of oil and gas from shale reservoirs or for the design of selective adsorbents for chemical separation. In this work, we reveal non-intuitive relationships that emerge when comparing bulk fluid properties to those when confined in amorphous silica single-digit nanopores (approx. 5 nm SBA-15), specifically the confined density and self-diffusivity. A counterintuitive result was obtained where larger molecules had relatively faster self-diffusivities; however, this behavior was revealed to be caused by these larger compounds exhibiting a relatively low density when confined. As the prediction of confined fluid properties impacts a diversity of disciplines, the results highlight the need for additional investigations where consistent and large-scale experimental procedures are employed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Vinylsilanes as Chain-Transfer Agents in Ethylene Polymerization: Direct Synthesis of Heterotelechelic Polyolefins

A single-step route to heterotelechelic polyethylene is enabled by palladium-diimine-catalyzed polymerization of ethylene using vinylsilanes as chain-transfer agents. The reaction affords α-alkenyl-ω-silyl end-capped heterotelechelic polymers whose molecular weights are controllable over a wide range by adjusting the [ethylene]/[vinylsilane] molar ratios. Notably, highly efficient end-capping with silanes carrying a variety of functionalities can be achieved under the optimized conditions. Furthermore, the alkenyl- and silyl-terminated telechelics serve as polymer precursors for further reactions and can be converted into additional telechelic functionalized polyolefins in good yields.

Chen, Zhi-Hao [University of Houston, TX (United S↗

Spectral Properties and Coding Transitions of Haar-Random Quantum Codes

A quantum error-correcting code with a nonzero error threshold undergoes a mixed-state phase transition when the error rate reaches that threshold. We explore this phase transition for Haar-random quantum codes, in which the logical information is encoded in a random subspace of the physical Hilbert space. We focus on the spectrum of the encoded system density matrix as a function of the rate of uncorrelated, single-qudit errors. For low error rates, this spectrum consists of well-separated bands, representing errors of different weights. As the error rate increases, the bands for high-weight errors merge. The evolution of these bands with increasing error rate is well described by a simple analytic ansatz. Using this ansatz, as well as an explicit calculation, we show that the threshold for Haar-random quantum codes saturates the hashing bound, and thus coincides with that for random stabilizer codes. For error rates that exceed the hashing bound, typical errors are uncorrectable, but postselected error correction remains possible until a much higher detection threshold. Postselection can in principle be implemented by projecting onto subspaces corresponding to low-weight errors, which remain correctable past the hashing bound.

decoherence↗

Accelerating Scientific Simulations with Bi-Fidelity Weighted Transfer Learning

High-fidelity modeling is an essential design tool for many engineering applications. However, for complex systems, computational cost can be a limiting factor. Analyzing parameter sensitivity, uncertainty quantification, and design optimization require many model evaluations. Surrogate models are often used to develop the relationship between model parameters and quantities of interest. However, in the case of complex systems, surrogate models require several degrees of freedom and, thus, a large number of data points to determine the correct dependencies. For many applications, this may be prohibitively expensive. The reduction of computational requirements can be achieved by leveraging low-fidelity models. Low-fidelity models represent the system at a coarser resolution with the advantage of computational efficiency. Therefore, a bi-fidelity modeling paradigm, which augments the accuracy of a low-fidelity model in a computationally efficient manner by invoking limited runs of a high-fidelity model, can be leveraged to sufficiently balance the accuracy and computational requirements. In this work, a bi-fidelity weighted transfer learning method using neural networks was applied to a computational fluid dynamics heat transfer modeling problem. The transfer learning advantage was investigated as a function of hyperparameters. Our main finding is that the use of a bi-fidelity modeling paradigm achieves accuracy close to that of a high-fidelity Gaussian process model while significantly reducing computational cost. The bi-fidelity model achieves comparable performance with 90 high-fidelity samples-that is, 60% less than the samples needed to achieve similar accuracy without the use of bi-fidelity modeling,

Borowiec, Katarzyna↗

Signatures of cysteine oxidation on muscle structural and contractile proteins are associated with physical performance and muscle function in older adults: Study of Muscle, Mobility and Aging (SOMMA)

Oxidative stress is considered a contributor to declining muscle function and mobility during aging; however, the underlying molecular mechanisms remain poorly described. We hypothesized that greater levels of cysteine (Cys) oxidation on muscle proteins are associated with decreased measures of mobility. Herein, we applied a novel redox proteomics approach to measure reversible protein Cys oxidation in vastus lateralis muscle biopsies collected from 56 subjects in the Study of Muscle, Mobility and Aging (SOMMA), a community-based cohort study of individuals aged 70 years and older. We tested whether levels of Cys oxidation on key muscle proteins involved in muscle structure and contraction were associated with muscle function (leg power and strength), walking speed, and fitness (VO 2 peak on cardiopulmonary exercise testing) using linear regression models adjusted for age, sex, and body weight. Higher oxidation levels of select nebulin Cys sites were associated with lower VO 2 peak, while greater oxidation of myomesin-1, myomesin-2, and nebulin Cys sites was associated with slower walking speed. Higher oxidation of Cys sites in key proteins such as myomesin-2, alpha-actinin-2, and skeletal muscle alpha-actin were associated with lower leg power and strength. We also observed an unexpected correlation (R = 0.48) between a higher oxidation level of eight Cys sites in alpha-actinin-3 and stronger leg power. Despite this observation, the results generally support the hypothesis that Cys oxidation of muscle proteins impairs muscle power and strength, walking speed, and cardiopulmonary fitness with aging.

59 BASIC BIOLOGICAL SCIENCES↗

The Cosmological Bootstrap: Spinning correlators from symmetries and factorization

We extend the cosmological bootstrap to correlators involving massless spinning particles, focusing on spin-1 and spin-2. In de Sitter space, these correlators are constrained both by symmetries and by locality. In particular, the de Sitter isometries become conformal symmetries on the future boundary of the spacetime, which are reflected in a set of Ward identities that the boundary correlators must satisfy. We solve these Ward identities by acting with weight-shifting operators on scalar seed solutions. Using this weight-shifting approach, we derive three- and four-point correlators of massless spin-1 and spin-2 fields with conformally coupled scalars. Four-point functions arising from tree-level exchange are singular in particular kinematic configurations, and the coefficients of these singularities satisfy certain factorization properties. We show that in many cases these factorization limits fix the structure of the correlators uniquely, without having to solve the conformal Ward identities. The additional constraint of locality for massless spinning particles manifests itself as current conservation on the boundary. We find that the four-point functions only satisfy current conservation if the s, t, and u-channels are related to each other, leading to nontrivial constraints on the couplings between the conserved currents and other operators in the theory. For spin-1 currents this implies charge conservation, while for spin-2 currents we recover the equivalence principle from a purely boundary perspective. For multiple spin-1 fields, we recover the structure of Yang--Mills theory. Finally, we apply our methods to slow-roll inflation and derive a few phenomenologically relevant scalar-tensor three-point functions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Stakeholder-guided holistic, Adaptive Framework for enhancing community Energy Resilience (SAFER) (Final Technical Report)

The Stakeholder-guided holistic, Adaptive Framework for enhancing community Energy Resilience (SAFER) project advances resilience science and engineering by addressing challenges in rural Kansas communities where aging infrastructure, extreme weather, and socioeconomic disparities heighten vulnerability to energy disruptions. Traditional approaches often focus on technical performance while overlooking community concerns and priorities. SAFER responds by integrating community perspectives with advanced analytical frameworks to create a holistic model for measuring and improving resilience. Project objectives included developing novel resilience metrics, advancing modeling frameworks that capture interdependencies across infrastructures, and embedding community-centric indicators directly into planning processes for distributed energy resources. The key technical innovations included the creation of self-organizing map (SOM)-based indices for objective resilience quantification, hetero-functional graph theory (HFGT) models linking power, water, transportation, and community assets, and graph neural network (GNN) tools for identifying critical nodes in complex systems. Community-centric energy planning was demonstrated through optimal siting and sizing of (photovoltaic) PV and battery storage, ensuring resilience enhancements also addressed energy burden and energy insecurity. SAFER engaged community partners in Dodge City and Ford County through surveys, focus groups, and workshops, generating more than 600 responses that established baseline measures of energy burden, financial insecurity, and willingness-to-pay to avoid outages. This data, organized in terms of a community capitals framework, informed the development of weighted reliability indices that better reflect community costs than traditional utility metrics. SAFER’s GNN-based critical node identification framework identified expert-labelled critical nodes with over 99% accuracy, while also uncovering additional functionalities essential for proactive resilience planning. The project’s models demonstrated that optimal PV and storage deployment could improve resilience indices by over 11 percent, with dispatch strategies further enhancing outcomes, confirming both the technical effectiveness and economic feasibility of these approaches. Through its combined emphasis on rigorous modeling, community-focused planning, and community engagement, SAFER advances the state of resilience research while delivering direct benefits to rural communities. The project provides tools, guidelines, and resilience heatmaps that help utilities, local governments, and residents better anticipate disruptions, prioritize investments, and strengthen the capacity to withstand and recover from energy-related hazards. Furthermore, the developed HFG and GNN frameworks are designed for transferability, allowing them to be adapted for resilience planning in other communities with minimal retraining. This inductive learning capability provides a scalable pathway to extend the SAFER project’s impact. Thus, creating a foundation for a nationally applicable model of infrastructure resilience. Additionally, the HFG can also be extended to include other FEMA community lifelines.

14 SOLAR ENERGY↗

Adsorption of terbium (III) on DGA and LN resins: Thermodynamics, isotherms, and kinetics

Two commercially available extraction chromatography (EXC) resins containing N,N,N’,N’-tetra-n-octyldiglycolamide (DGA Resin, Normal, 50 – 100 μm) and Bis(2-ethylhexyl) phosphate (LN Resin, 100 – 150 μm) were used as adsorbents to study fundamental adsorption properties such as thermodynamic values, equilibrium isotherms, and kinetic uptake models for terbium(III) adsorption. Weight distribution ratios (D w ) for terbium on DGA and LN resins were measured using a [ 160 Tb]Tb 3+ radiometric tracer in nitric acid as a function of acidity, temperature, initial analyte concentration, and equilibrium time. The D w values showed increasing binding affinity for DGA resin at high nitric acid concentrations and decreasing binding affinity for LN resins. Thermodynamic studies for DGA and LN resins revealed that the Gibbs free energy (ΔG) increased consistently with temperature. To model equilibrium data, increasingly higher parameter equilibrium isotherm models (Henry (1) < Langmuir, Freundlich (2) < Redlich-Peterson (3) < Fritz-Schluender (4)) were compared on their root mean squared errors (RMSE) and adjusted determination coefficients to determine the most applicable model. In all cases, the empirical four-parameter Fritz-Schluender isotherm demonstrated a superior fit. Similar comparisons for reaction-based kinetic models (Pseudo-first-order < Pseudo-second-order < Pseudo-n-order) revealed that the higher-order PNO model yielded a superior fit of kinetic data for both resins. Furthermore, in some cases, adsorption isotherms and kinetic models could also be modeled by a lower-order model with minimal change in error parameters. Weber-Morris plots revealed that two linear sections are observed for each resin, where the first linear segment is attributed to fast (film diffusion) adsorption of terbium, followed by slower intraparticle diffusion of terbium through the pores as the rate-limiting step. Based on the Weber-Morris plot, both film and intraparticle diffusion are involved in controlling the kinetic rate of adsorption for DGA and LN resins.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗