Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Differentiable physics”

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

An adaptive and stability-promoting layerwise training approach for sparse deep neural network architecture

This work presents a two-stage adaptive framework for progressively developing deep neural network (DNN) architectures that generalize well for a given training data set. In the first stage, a layerwise training approach is adopted where a new layer is added each time and trained independently by freezing parameters in the previous layers. We impose desirable structures on the DNN by employing manifold regularization, sparsity regularization, and physics-informed terms. We introduce a ε – δ – stability-promoting concept as a desirable property for a learning algorithm and show that employing manifold regularization yields a ε – δ stability-promoting algorithm. Further, we also derive the necessary conditions for the trainability of a newly added layer and investigate the training saturation problem. In the second stage of the algorithm (post-processing), a sequence of shallow networks is employed to extract information from the residual produced in the first stage, thereby improving the prediction accuracy. Numerical investigations on prototype regression and classification problems demonstrate that the proposed approach can outperform fully connected DNNs of the same size. Moreover, by equipping the physics-informed neural network (PINN) with the proposed adaptive architecture strategy to solve partial differential equations, we numerically show that adaptive PINNs not only are superior to standard PINNs but also produce interpretable hidden layers with provable stability. As a result, we also apply our architecture design strategy to solve inverse problems governed by elliptic partial differential equations.

42 ENGINEERING↗

Uncovering Backbone Conformation for Rigid DPP-Based Donor–Acceptor Conjugated Polymer Using Deuterium Labeling and Neutron Scattering

The conjugated polymer’s backbone conformation dictates the delocalization of electrons, ultimately affecting its optoelectronic properties. Most conjugated polymers can be viewed as semirigid rods with their backbone embedded among long alkyl side chains. Thus, it is challenging to experimentally quantify the conformation of a conjugated backbone. Here, we performed contrast variation neutron scattering on rigid conjugated donor–acceptor (D–A) diketopyrrolopyrrole (DPP) polymers with selectively deuterated side chains to measure the conjugated backbone conformation. We first synthesized DPPbased polymers with deuterated side chains, confirmed by NMR and FTIR. Using contrast variation neutron scattering, we found that the DPP-based conjugated polymers are much more rigid than poly(3-alkylthiophenes), with persistence length (L p ) at 16–18 nm versus 2–3 nm. More importantly, in contrast to the relatively flexible poly(3-alkylthiophenes) whose backbone is more flexible than the whole polymer, we found that the backbone of DPP-based polymers has the same L p value compared to the whole polymer chain. This indicates that side chain interference on backbone conformation is not present for the semirigid polymer, which is further confirmed by coarse-grained molecular dynamics (CG-MD) simulations. Our work provides a novel protocol to probe polymer’s backbone conformation and paradigm-shifting understanding of the backbone conformation of semirigid conjugated polymers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Subtle Molecular Changes Largely Modulate Chiral Helical Assemblies of Achiral Conjugated Polymers by Tuning Solution-State Aggregation

Understanding the solution-state aggregate structure and the consequent hierarchical assembly of conjugated polymers is crucial for controlling multiscale morphologies during solid thin-film deposition and the resultant electronic properties. However, it remains challenging to comprehend detailed solution aggregate structures of conjugated polymers, let alone their chiral assembly due to the complex aggregation behavior. Herein, we present solution-state aggregate structures and their impact on hierarchical chiral helical assembly using an achiral diketopyrrolopyrrole-quaterthiophene (DPP-T4) copolymer and its two close structural analogues wherein the bithiophene is functionalized with methyl groups (DPP-T2M2) or fluorine atoms (DPP-T2F2). Combining in-depth small-angle X-ray scattering analysis with various microscopic solution imaging techniques, we find distinct aggregate in each DPP solution: (i) semicrystalline 1D fiber aggregates of DPP-T2F2 with a strongly bound internal structure, (ii) semicrystalline 1D fiber aggregates of DPP-T2M2 with a weakly bound internal structure, and (iii) highly crystalline 2D sheet aggregates of DPP-T4. These nanoscopic aggregates develop into lyotropic chiral helical liquid crystal (LC) mesophases at high solution concentrations. Intriguingly, the dimensionality of solution aggregates largely modulates hierarchical chiral helical pitches across nanoscopic to micrometer scales, with the more rigid 2D sheet aggregate of DPP-T4 creating much larger pitch length than the more flexible 1D fiber aggregates. Combining relatively small helical pitch with long-range order, the striped twist-bent mesophase of DPP-T2F2 composed of highly ordered, more rigid 1D fiber aggregate exhibits an anisotropic dissymmetry factor (g-factor) as high as 0.09. This study can be a prominent addition to our knowledge on a solution-state hierarchical assembly of conjugated polymers and, in particular, chiral helical assembly of achiral organic semiconductors that can catalyze an emerging field of chiral (opto)electronics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Resonance cones in cold plasma: Origin, singularities, and power flow

In magnetized tenuous plasma, typical at the plasma edge of fusion devices, a nearly electrostatic wave mode with relatively enhanced electric field can propagate along a specific angle with the magnetic field. For this characteristic, it is known as a “resonance cone.” For instance, these waves can be excited by radio frequency antennas in the ion-cyclotron and lower-hybrid range of frequencies. We consider the resonance cones emitted by idealized spatially extended sources. In 2D, we use a novel geometric construction which generalizes the d'Alembert solution to curved boundaries/moving sources, and show, for the first time, that singular electric fields arise under these conditions, thereby bringing the resonance cones in line with the other resonances of the cold plasma theory. Still in 2D, we give an expression for the amount of power radiated by resonance cones in terms of surface quantities on the source, which is finite despite the singular electric field. In conclusion, we generalize the conclusions regarding the presence and location of singular electric fields to the 3D electromagnetic case.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Theoretical investigation of decoherence channels in athermal phonon sensors

The creation and evolution of nonequilibrium phonons is central in applications ranging from cosmological particle searches to decoherence processes in qubits. However, the fundamental understanding of decoherence pathways for athermal phonon distributions in solid-state systems remains an open question. Using first-principles calculations, we investigate the primary decay channels of athermal phonons in two technologically relevant semiconductors—Si and GaAs. We quantify the contributions of anharmonic, isotopic, and interfacial scattering in these materials. From this, we construct a model to estimate the thermal power in a readout scheme as a function of time. We discuss the implication of our results on noise limitations in current phonon sensor designs and strategies for improving coherence in next-generation phonon sensors.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Neural posterior unfolding

Differential cross section measurements are the currency of scientific exchange in particle and nuclear physics. A key challenge for these analyses is the correction for detector distortions, known as deconvolution or unfolding. Binned unfolding of cross section measurements traditionally rely on the regularized inversion of the response matrix that represents the detector response, mapping pre-detector (`particle level') observables to post-detector (`detector level') observables. In this paper we introduce Neural Posterior Unfolding, a modern, Bayesian approach that leverages normalizing flows for unfolding. By using normalizing flows for neural posterior estimation, NPU offers several key advantages including implicit regularization through the neural network architecture, fast amortized inference that eliminates the need for repeated retraining, and direct access to the full uncertainty in the unfolded result. In addition to introducing NPU, we implement a classical Bayesian unfolding method called Fully Bayesian Unfolding (FBU) in modern Python so it can also be studied. These tools are validated on simple Gaussian examples and then tested on simulated jet substructure examples from the Large Hadron Collider (LHC). We find that the Bayesian methods are effective and worth additional development to be analysis ready for cross section measurements at the LHC and beyond.

Analysis and statistical methods↗

Modelling and simulation of brinicle formation

Below the Arctic sea ice, under the right conditions, a flux of icy brine flows down into the sea. The icy brine has a much lower fusion point and is denser than normal seawater. As a result, it sinks while freezing everything around it, forming an ice channel called a brinicle (also known as ice stalactite). In this paper, we develop a mathematical model for this phenomenon, assuming cylindrical symmetry. The fluid is considered to be viscous and quasi-stationary. The heat and salt transport are weakly coupled to the fluid motion and are modelled with the corresponding conservation equations, accounting for diffusive and convective effects. Finite-element discretization is employed to solve the coupled system of partial differential equations. We find that the model can capture the general behaviour of the physical system and generate brinicle-like structures while also recovering dendrite composition, which is a physically expected feature aligned with previous experimental results. This represents, to our knowledge, the first complete model proposed that captures the global structure of the physical phenomenon even though it has some discrepancies, such as brine accumulation.

97 MATHEMATICS AND COMPUTING↗

Single‐parent expression drives dynamic gene expression complementation in maize hybrids

SUMMARY Single‐parent expression (SPE) is defined as gene expression in only one of the two parents. SPE can arise from differential expression between parental alleles, termed non‐presence/absence (non‐PAV) SPE, or from the physical absence of a gene in one parent, termed PAV SPE. We used transcriptome data of diverse Zea mays (maize) inbreds and hybrids, including 401 samples from five different tissues, to test for differences between these types of SPE genes. Although commonly observed, SPE is highly genotype and tissue specific. A positive correlation was observed between the genetic distance of the two inbred parents and the number of SPE genes identified. Regulatory analysis showed that PAV SPE and non‐PAV SPE genes are mainly regulated by cis effects, with a small fraction under trans regulation. Polymorphic transposable element insertions in promoter sequences contributed to the high level of cis regulation for PAV SPE and non‐PAV SPE genes. PAV SPE genes were more frequently expressed in hybrids than non‐PAV SPE genes. The expression of parentally silent alleles in hybrids of non‐PAV SPE genes was relatively rare but occurred in most hybrids. Non‐PAV SPE genes with expression of the silent allele in hybrids are more likely to exhibit above high parent expression level than hybrids that do not express the silent allele, leading to non‐additive expression. This study provides a comprehensive understanding of the nature of non‐PAV SPE and PAV SPE genes and their roles in gene expression complementation in maize hybrids.

Li, Zhi↗

An Educational Guide for 2D Stellar Structure Calculations of Rapidly Rotating Stars using the ESTER code

The Evolution STEllaire en Rotation (ESTER) code is the first 2D stellar structure code to be made open-source and freely available to the astronomy and astrophysics community. An important and novel advancement of this code is that it can reproduce the distorted shape and observable signatures (e.g., gravity darkening) of rapidly rotating stars. ESTER also calculates the steady-state large-scale flows within the star, namely their differential rotation and associated meridional circulation. In this report, we explore and document the physics implemented within version 1.1.0rc2 of the ESTER code, in a way that complements published descriptions. We illustrate this physics by plotting how stellar structure parameters vary through stellar interiors at a range of latitudes and at different angular velocities. We investigate how the thin convective envelopes of intermediate mass stars vary with latitude when rapidly rotating, becoming deeper and thicker near the equator. Simple comparisons of ESTER model predictions (e.g., central temperature and density, luminosity) with the output from the Modules for Experiments in Stellar Astrophysics (MESA) code [Paxton et al., 2010] shows generally good agreement. Additional comparisons provide important benchmarking and verification for ESTER as a comparatively young code. Finally, we provide a guide for installing and running the code on our local university cluster, aimed at helping students to begin work.

79 ASTRONOMY AND ASTROPHYSICS↗

Predictor-corrector models for lightweight massive machine-type communications in Industry 4.0

Future Industry 4.0 scenarios are characterized by seamless integration between computational and physical processes. To achieve this objective, dense platforms made of small sensing nodes and other resource constraint devices are ubiquitously deployed. All these devices have a limited number of computational resources, just enough to perform the simple operation they are in charge of. The remaining operations are delegated to powerful gateways that manage sensing nodes, but resources are never unlimited, and as more and more devices are deployed on Industry 4.0 platforms, gateways present more problems to handle massive machine-type communications. Although the problems are diverse, those related to security are especially critical. To enable sensing nodes to establish secure communications, several semiconductor companies are currently promoting a new generation of devices based on Physical Unclonable Functions, whose usage grows every year in many real industrial scenarios. Those hardware devices do not consume any computational resource but force the gateway to keep large key-value catalogues for each individual node. In this context, memory usage is not scalable and processing delays increase exponentially with each new node on the platform. In this paper, we address this challenge through predictor-corrector models, representing the key-value catalogues. Models are mathematically complex, but we argue that they consume less computational resources than current approaches. The lightweight models are based on complex functions managed as Laurent series, cubic spline interpolations, and Boolean functions also developed as series. Unknown parameters in these models are predicted, and eventually corrected to calculate the output value for each given key. The initial parameters are based on the Kane Yee formula. An experimental analysis and a performance evaluation are provided in the experimental section, showing that the proposed approach causes a significant reduction in the resource consumption.

97 MATHEMATICS AND COMPUTING↗

Analysis of Material Properties of Aged Rigid Polyurethane Foam

Foams are classified as either rigid or flexible and are commonly used in applications such as pipe insulation, carpet fibers, and chair foam cushions. Polyurethane foams are excellent thermal insulators, have excellent mechanical properties, and are chemically inert. With the addition of fire protective additives, polyurethane foams have the ability to intumescence when exposed to high temperatures. For the future storage of the 9977 nuclear shipping packages, the material properties of the General Plastics 3716 rigid polyurethane foam were evaluated for long term performance by comparing aged samples to unaged samples. General Plastics 3716 rigid polyurethane foam was aged from November 2014 to June 2019 at varying temperatures of: 160 deg. F, 160 deg. F + 50% relative humidity (RH), 185 deg. F, 215 deg. F, and 250 deg. F. The foam was then analyzed via Fourier Transform Infrared Spectroscopy (FTIR), 3D Optical Microscope, Thermal Gravimetric Analysis (TGA), Differential Scanning Calorimetry (DSC). FTIR and Optical Microscope: No significant change to chemical or physical structure was evident by FTIR or optical microscope analysis. DSC: Decrease of heat capacity over time for all oven temperatures except a slight increase for the 185 deg. F conditioned foams. Largest decrease evident with foams conditioned at 250 deg. F. A slight increase in heat capacity for humid samples. Similar trends in heat capacities at 71 deg. C and 102 deg. C. TGA: As conditioning temperature increased, a decline was evident in the average percent weight loss, onset temperature, and first derivative peak temperature. No significant effect due to humidity except for a decrease in first derivative peak temperature. Unaged samples exhibited strong intumescence. Foams conditioned at 160 deg. F exhibited slight intumescence. Results suggest fire protective additives decomposed at higher temperatures, which contributed to a lower average percent mass loss in higher temperature conditioning. Further thermal degradation testing would be beneficial to evaluate the ability of the conditioned foams to intumesce. Evaluate irradiation dose rate effects on conditioned foam. Possible conditioning in an environment with a known humidity or analyzing a new aliquot of sample for each trial to achieve a more representative heat capacity calculation.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

E-PINNs: Epistemic Physics-Informed Neural Networks

Physics-informed neural networks (PINNs) have demonstrated promise as a framework for solving forward and inverse problems involving partial differential equations. Despite recent progress in the field, it remains challenging to quantify uncertainty in these networks. While techniques such as Bayesian PINNs (B-PINNs) provide a principled approach to capturing epistemic uncertainty through Bayesian inference, they can be computationally expensive for large-scale applications. In this work, we propose Epistemic Physics-Informed Neural Networks (E-PINNs), a framework that uses a small network, the epinet, to efficiently quantify epistemic uncertainty in PINNs. The proposed approach works as an add-on to existing, pre-trained PINNs with a small computational overhead. We demonstrate the applicability of the proposed framework in various test cases and compare the results with B-PINNs using Hamiltonian Monte Carlo (HMC) posterior estimation and dropout-equipped PINNs (Dropout-PINNs). In our experiments, E-PINNs achieve calibrated coverage with competitive sharpness at substantially lower cost. We demonstrate that when B-PINNs produce narrower bands, they under-cover in our tests. E-PINNs also show better calibration than Dropout-PINNs in these examples, indicating a favorable accuracy-efficiency trade-off.

AI for Science↗

PINNSTRIPES (Physics-Informed Neural Network SurrogaTe for Rapidly Identifying Parameters in Energy Systems) [SWR-22-12]

Energy systems models typically take the form of complex partial differential equations which make multiple forward calculations prohibitively expensive. Fast and data-efficient construction of surrogate models is of utmost importance for applications that require parameter exploration such as design optimization and Bayesian calibration. In presence of a large number of parameters, surrogate models that capture correct dependencies may be difficult to construct with traditional techniques. The issue is addressed here with the formulation of the surrogate model constructed via Physics-Informed Neural Networks (PINN) which capture the dependence with respect to the parameters to estimate, while using a limited amount of data. Since forward evaluations of the surrogate model are cheap, parameter exploration is made inexpensive, even when considering a large number of parameters.

Hassanaly, Malik↗

Order conditions for nonlinearly partitioned Runge-Kutta methods

Recently, a new class of nonlinearly partitioned Runge–Kutta (NPRK) methods was proposed for nonlinearly partitioned systems of autonomous ordinary differential equations y' = F(y, y). The target class of problems are those in which different scales, stiffnesses, or physics are coupled in a nonlinear way, wherein the desired partition cannot be written in a classical additive or component-wise fashion. Here we use a rooted-tree analysis to derive full-order conditions for NPRKM methods, where M denotes the number of nonlinear partitions. Due to the nonlinear coupling and thereby the mixed product differentials, it turns out that the standard node-colored rooted tree analysis used in analyzing ODE integrators does not naturally apply. Instead we develop a new edge-colored rooted-tree framework to address the nonlinear coupling. The resulting order conditions are enumerated, are provided directly for up to fourth order with M = 2 and third order with M = 3, and are related to existing order conditions of additive and partitioned RK methods. We conclude with an example that shows how the nonlinear order conditions can be used to obtain an embedded estimate of the state-dependent nonlinear coupling strength in a dynamical system.

97 MATHEMATICS AND COMPUTING↗

Understanding plasma turbulence through exact coherent structures

Plasma turbulence is a key challenge in understanding transport phenomena in magnetically confined plasmas. This work presents a generalized framework to analyze plasma turbulence that utilizes periodic orbit theory. In periodic orbit theory, doubly periodic solutions (coherent structures) of the governing equation(s) serve as building blocks of the considered turbulent dynamics. To illustrate the concept and method, the particularly simple Kuramoto–Sivashinsky (referred to here as LMRT for the original authors: LaQuey, Mahajan, Rutherford, and Tang) trapped-ion mode toy model is used. By applying numerical optimization techniques to the LMRT equation, we extract coherent spacetime patterns that represent the library of allowable fundamental structures of the equation. These structures provide a framework to systematically describe turbulence as a composition of recurrent solutions, revealing an underlying order within chaotic plasma motion. Although illustrated here using the simplified LMRT model for clarity, this framework provides a general strategy that can be extended to more complex and realistic models of plasma turbulence, including gyrokinetic systems. This offers a new method for predicting and potentially controlling transport processes in fusion plasmas by providing a bridge between nonlinear dynamical systems theory and plasma physics in the form of a generalized framework with which to analyze and understand spatially extended nonlinear partial differential equations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

TTDAQ: A Continuous Flow, Timing and Trigger DAQ System

Final Scientific/Technical Report for DOE Award DE-SC0019581, “TTDAQ: A Continuous Flow, Timing and Trigger DAQ System.” The report summarizes Telluric Labs’ Phase II STTR work developing silicon-photonic building blocks for a software-defined, continuous-flow, trigger-less data acquisition system for next-generation high-energy and nuclear-physics detectors. The project focused on radiation-hard photonic integrated circuits, remote optical illumination, dense wavelength-division multiplexing, and a differential microring-resonator transceiver architecture designed to improve high-speed optical link stability and bandwidth. The report describes project objectives, technical accomplishments, AIM Photonics tape-outs, bench characterization, radiation-hardness testing, deferred integration work, and potential applications beyond physics readout.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗