Search NASA⌕ Search

SEARCH · Search NASA

Results for “differentiable programming”

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 217 records · Page 12

Parallel Variable Population Multi-Objective Optimizer (pvpmoo) v1.0

This is a parallel variable population multi-objective optimizer with an adaptive unified differential evolution algorithm or a genetic algorithm. It can also be used for single objective optimization. Some features of this code include: 1) The population size varies from generation to generation to save the total # of objective function evaluations. 2) The population is uniformly distributed to a number of parallel processors for simultaneous objective function evaluation. 3) The objective function evaluation can be attained from an external simulation program with control variables in its input file and objectives calculated from its output files. 4) The optimizer includes an adaptive unified differential evolution algorithm and a real value genetic algorithm. The parameters in the unified differential evolution algorithm can be chosen to attain any mutation schemes in the published literature.

Qiang, Ji↗

Differentiable Multiphysics Codes: A Breakthrough Technology for Simulation and Computing

This document summarizes the findings of a strategic planning exercise commissioned by the Weapons Simulation and Computing, Computational Physics (WSC/CP) program at the Lawrence Livermore National Laboratory (LLNL) in FY24. During the year, the committee met with multiple stakeholder communities to gather input, opinions, suggestions and concerns which have been incorporated throughout this document. The key findings from this exercise are summarized: • The development of multiphysics modelling and simulation (mod/sim) codes and software technologies, their deployment on exascale compute platforms, and their broad adoption across the NNSA is a major success of the Advanced Simulation and Computing (ASC) program and the Exascale Computing Project (ECP). Sustained investment in these core technologies is essential. • Today’s state of the art involves running ensembles of O(100K) simulations to perform uncertainty quantification (UQ) and design studies using multiple statistical methods such as Bayesian optimization to understand sensitivities of our models and explore parameterized design spaces. Even with exascale computing, we are practically limited to O(10) parameters in these studies since the number of simulations required to sample the space scales exponentially with the number of design parameters. • The data from these simulation ensembles is increasingly being used to train machine learned (ML) surrogates (or reduced order models, ROMs) which can then be used for optimization or real time design exploration. However, the trained surrogates are still limited in the number of parameters they can represent due to the sampling limitations previously noted. • Augmenting our suite of integrated multiphysics simulation codes, both current and emerging, with the ability to compute gradients (solution derivatives) of arbitrary simulation outputs with respect to (some or all) simulation inputs would be a breakthrough technology, opening the door to a new era of efficient and automated inverse design based on verified and validated mod/sim capabilities. • This capability, which we refer to as differentiable multiphysics codes (DMCs), would revolutionize both UQ and optimization studies by breaking the curse of dimensionality that presently limits our “gradient-free” ensemble based computing approach. A similar breakthrough occurred in the AI/ML community once the ability to compute gradients of arbitrary loss functions using back-propagation became commonplace. Gradient information from the multiphysics codes can also be used to dramatically improve the efficiency and scale of training of ML/ROM surrogates for rapid assessments. • Achieving this in our suite of codes will be a grand challenge, similar to the amount of effort that was required to transition from CPU to GPU computing. It will require buy-in from the entire WSC/CP program and beyond, including all integrated codes, physics and engineering models, third-party library dependencies and performance portability abstractions. It will also require investment in research and development of numerical methods for computing adjoints of coupled physics across multiple adaptively refined moving meshes and of stochastic (Monte Carlo) and mesh free (SPH) methods. • New software and numerical techniques, largely pioneered by the AI/ML community, make this feasible. Chief among these is automatic differentiation (AD), the ability to employ AD at point-wise locations in a physics calculation (instead of traditional black-box approaches) and the ability to perform “back-propagation in time” (or reverse mode AD) for non-linear partial differential equations (PDEs). Fundamentally, the conclusion of this strategic planning exercise is that the time is right to undertake a large scale effort in WSC, centered on the existing integrated codes, to continue the natural evolution of mod/sim in the age of AI/ML. Instead of attempting to replace mod/sim with purely data driven AI/ML models, we believe the key to success is to integrate AI/ML by building on top of the decades of hard-won knowledge and the verified/validated multiphysics modelling capability that is the hallmark of the ASC program.

97 MATHEMATICS AND COMPUTING↗

The Problem of Package Guidance with Incomplete Information for a Linear Control System with a Delay

The control-package method is applied to solve the problem of constructing guaranteed positional guidance for a linear control dynamical system with a delay under incomplete information about the initial state. The problem is solved for the case when the initial state is contained in a compact functional set. A solvability criterion and a method for the construction of a guidance program package are proposed. The proposed procedure is applied to a particular linear control system of delayed differential equations.

97 MATHEMATICS AND COMPUTING↗

A SEARCH FOR OB ASSOCIATIONS NEAR SOUTHERN LONG-PERIOD CEPHEIDS. V. AQ PUPPIS AND V620 PUPPIS

A photometric UBV survey is presented for 610 stars in a region surrounding the Cepheid AQ Puppis and centered southwest of the variable, based upon photoelectric measures for 14 stars and calibrated iris photometry of photographic plates of the field for 596 stars. An analysis of reddening and distance for program stars indicates that the major dust complex in this direction is {approx}1.8 kpc distant, producing differential extinction described by a ratio of total-to-selective extinction of R = A{sub V} /E{sub B-V} = 3.10 {+-} 0.20. Zero-age main-sequence fitting for the main group of B-type stars along the line of sight yields a distance of 3.21 {+-} 0.19 kpc (V{sub 0} - M{sub V} = 12.53 {+-} 0.13 s.e.). The 29fd97 Cepheid AQ Pup, of field reddening E{sub B-V} = 0.47 {+-} 0.07 (E{sub B-V}(B0) = 0.51 {+-} 0.07), appears to be associated with B-type stars lying within 5' of it as well as with a sparse group of stars, designated Turner 14, centered south of it at J2000.0 = 07:58:37, -29:25:00, with a mean reddening of E{sub B-V} = 0.81 {+-} 0.01. AQ Pup has an inferred luminosity as a cluster member of (M{sub V} ) = -5.40 {+-} 0.25 and an evolutionary age of 3 Multiplication-Sign 10{sup 7} yr. Its observed rate of period increase of +300.1 {+-} 1.2 s yr{sup -1} is an order of magnitude larger than what is observed for Cepheids of comparable period in the third crossing of the instability strip, and may be indicative of a high rate of mass loss or a putative fifth crossing. Another sparse cluster, designated Turner 13, surrounds the newly recognized 2fd59 Cepheid V620 Pup, of space reddening E{sub B-V} = 0.64 {+-} 0.02 (E{sub B-V}(B0) = 0.68 {+-} 0.02), distance 2.88 {+-} 0.11 kpc (V{sub 0} - M{sub V} 12.30 {+-} 0.08 s.e.), evolutionary age 10{sup 8} yr, and an inferred luminosity as a likely cluster member of (M{sub V}) = -2.74 {+-} 0.11. V620 Pup is tentatively identified as a first crosser, pending additional observations.

79 ASTRONOMY AND ASTROPHYSICS↗

Resistive Network Optimal Power Flow: Uniqueness and Algorithms

The optimal power flow (OPF) problem minimizes the power loss in an electrical network by optimizing the voltage and power delivered at the network buses, and is a nonconvex problem that is generally hard to solve. By leveraging a recent development on the zero duality gap of OPF, we propose a second-order cone programming convex relaxation of the resistive network OPF, and study the uniqueness of the optimal solution using differential topology, especially the Poincare-Hopf Index Theorem. We characterize the global uniqueness for different network topologies, e.g., line, radial, and mesh networks. This serves as a starting point to design distributed local algorithms with global behaviors that have low complexity, are computationally fast, and can run under synchronous and asynchronous settings in practical power grids.

Tan, CW↗

First Temperature Profile of a Stellar Flare Using Differential Chromatic Refraction

We present the first derivation of a stellar flare temperature profile from single-band photometry. Stellar flare DWF 030225.574−545707.45129 was detected in 2015 by the Dark Energy Camera as part of the Deeper, Wider, Faster program. The brightness (Δm g = −6.12) of this flare, combined with the high air mass (1.45 ≲ X ≲ 1.75) and blue filter (DES g, 398–548 nm) in which it was observed, provided ideal conditions to measure the zenithward apparent motion of the source due to differential chromatic refraction (DCR) and, from that, infer the effective temperature of the event. We model the flare’s spectral energy distribution as a blackbody to produce the constraints on flare temperature and geometric properties derived from single-band photometry. We additionally demonstrate how simplistic assumptions on the flaring spectrum, as well as on the evolution of flare geometry, can result in solutions that overestimate the effective temperature. Exploiting DCR enables studying chromatic phenomena with ground-based astrophysical surveys, and stellar flares on M dwarfs are a particularly enticing target for such studies due to their ubiquity across the sky and the heightened color contrast between their red quiescent photospheres and the blue flare emission. Our novel method will enable similar temperature constraints for a large sample of objects in upcoming photometric surveys like the Vera C. Rubin Legacy Survey of Space and Time.

Clarke, Riley W. [University of Delaware, Newark, ↗

DOE/NNSA Nuclear Criticality Safety Program – Brief Overview [Slides]

This presentation discusses nuclear data measurements and evaluation work for the NCSP. The objective of the NCSP is to provide measured and evaluated thermal, resonance, unresolved resonance, and fast region cross section data to address the priority NCSP nuclear data needs. It also discusses its vision which is to address multiple Nuclear Data 5- and 10-year goals and attributes identified in the NCSP. The intended final product is for rigorous ENDF/B evaluations to be produced from cross section measurements and analyses. It also discusses NCSP integral experiments and differential experiments in addition to NCSP support of SCALE in the future.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Deep Learning Explicit Differentiable Predictive Control Laws for Buildings

We present a differentiable predictive control (DPC) methodology for learning constrained control laws for unknown nonlinear systems. DPC poses an approximate solution to multiparametric programming problems emerging from explicit nonlinear model predictive control (MPC). Contrary to approximate MPC, DPC does not require supervision by an expert controller. Instead, a system dynamics model is learned from a small dataset of recorded observations of the perturbed system's dynamics and the control law is optimized offline by interaction with the learned system model. The DPC method is based on two sequential steps, i) system identification using a constrained neural state-space model, and ii) optimization of an explicit control law parametrized by another neural network in closed-loop simulation with the identified neural state-space model. The combination of a differentiable closed-loop system and penalty methods for constraint handling of system outputs and inputs allows us to optimize the control law's parameters directly by backpropagating economic MPC loss through the learned system model. By incorporating domain knowledge and leveraging established techniques from optimal control, our method leverages deep neural networks as nonlinear function approximators for system identification and control while avoiding concomitant costs of intractably large datasets, and computationally expensive over-parametrized models. The scalability, data efficiency, and constrained optimal control capability of the proposed DPC method are demonstrated in simulation using a multi-zone building emulator.

Drgona, Jan↗

Cyclosporin A inhibits adipogenic differentiation and regulates immunomodulatory functions of murine mesenchymal stem cells

Highlights: • CsA decreased BM-MSCs proliferation and increased cells apoptosis. • CsA inhibited BM-MSCs adipogenic differentiation. • CsA had no effect on the phenotype and osteogenic differentiation of BM-MSCs. • CsA inhibited BM-MSCs secreting IL-6 and promoted the expression of PD-L2. Aplastic anemia (AA) is generally considered as an immune-mediated bone marrow failure syndrome. Several studies show that bone marrow mesenchymal stem cells (BM-MSCs), as key cellular components of the bone marrow microenvironment, are also involved in the pathogenic mechanism of AA. Cyclosporin A (CsA) is a classic immunosuppressive drug for AA, and it specifically inhibits mammalian T cells by preventing activation of transcription factors involved in cytokine gene expression. However, little is known about the effect of CsA on the BM-MSCs. In this study, murine BM-MSCs were stimulated in the presence of CsA. Further, we found that CsA could inhibit murine BM-MSC proliferation and promote BM-MSC apoptosis, what's more CsA could inhibit adipogenic differentiation. Our study also showed that CsA could inhibit interleukin-6 expression in BM-MSCs, while promoting programmed death-ligand 2 expression. In conclusion, our results proposed that CsA may exert an effect on regulating the bone marrow environment by influencing BM-MSCs, which have a beneficial effect on treating AA.

60 APPLIED LIFE SCIENCES↗

Boundary-Layer-Coupled and Decoupled Clouds in Global Storm-Resolving Models: Comparisons With the ARM Observations

The accurate representation of interactions between clouds and planetary boundary layer (PBL) is a persistent challenge in climate models, critical for simulating surface energy budget. The emergence of kilometer-grid-scale global storm resolving models (GSRMs) offers the potential for enhanced details of PBL processes in these complex interactions. This study evaluates the representation of PBL-coupled and decoupled clouds in nine GSRM simulations against extensive ground-based observations by the Department of Energy Atmospheric Radiation Measurement (ARM) program, across six sites encompassing diverse regimes such as marine and continental environments in tropics and midlatitude. By differentiating coupling based on the relative positions between cloud bases and PBL tops, our analysis focuses on the simulation of PBL height, cloud frequency, position and vertical extent. The GSRMs generally exhibit commendable agreements with observed cloud structures and PBL diurnal cycles across different ARM sites. In contrast to the relatively consistent representation of decoupled clouds, discrepancies exist between the simulated and the observed coupled clouds, particularly in areas of intense convection, for example, over tropical rainforests and mountainous regions. These biases are probably associated with the models' tendency to underestimate the boundary layer humidity and the frequency of coupled clouds within different ranges of PBL heights. This study underscores the importance for continuous improvements in the representation of boundary layer and convection within these global kilometer-grid-scale models.

54 ENVIRONMENTAL SCIENCES↗

Generic Discretization Library

The GenDiL library is a collection of C++ software abstractions designed to discretize and solve partial differential equations (PDEs) for high-performance computing (HPC) applications. Its primary focus is on modern C++ generic programming, which helps ensure portability across various hardware architectures. The central idea behind the library is to provide building blocks for numerical algorithms-such as discretization methods and iteration patterns-so that domain experts can focus on the math, rather than the low-level details of hardware or implementation. By defining abstractions for data types, iteration over computational grids, and scheduling of operations, the library isolates the high-level PDE algorithms from the platform-specific optimizations needed to achieve efficient performance.

Dudouit, Yohann [Lawrence Livermore National Labor↗

A MUC 5 B Gene Polymorphism, rs35705950-T, Confers Protective Effects Against COVID-19 Hospitalization but Not Severe Disease or Mortality

A common MUC5B gene polymorphism, rs35705950-T, is associated with idiopathic pulmonary fibrosis (IPF), but its role in severe acute respiratory syndrome coronavirus 2 infection and disease severity is unclear. To assess whether rs35705950-T confers differential risk for clinical outcomes associated with coronavirus disease (COVID-19) infection among participants in the Million Veteran Program (MVP). The MUC5B rs35705950-T allele was directly genotyped among MVP participants; clinical events and comorbidities were extracted from the electronic health records. Associations between the incidence or severity of COVID-19 and rs35705950-T were analyzed within each ancestry group in the MVP followed by transancestry meta-analysis. Replication and joint meta-analysis were conducted using summary statistics from the COVID-19 Host Genetics Initiative (HGI). Sensitivity analyses with adjustment for additional covariates (body mass index, Charlson comorbidity index, smoking, asbestosis, rheumatoid arthritis with interstitial lung disease, and IPF) and associations with post–COVID-19 pneumonia were performed in MVP subjects. The MUC5B variant rs35705950-T may confer protection in COVID-19 hospitalizations.

59 BASIC BIOLOGICAL SCIENCES↗

A transcriptome software comparison for the analyses of treatments expected to give subtle gene expression responses

Background: In this comparative study we evaluate the performance of four software tools: DNAstar-D (DESeq2), DNAstar-E (edgeR), CLC Genomics and Partek Flow for identification of differentially expressed genes (DEGs) using a transcriptome of E. coli. The RNA-seq data are from the effect of below-background radiation 5.5 nGy total dose (0.2nGy/hr) on E. coli grown shielded from natural radiation 655 m below ground in a pre-World War II steel vault. The gene expression response to three supplemented sources of radiation designed to mimic natural background, 1952 – 5720 nGy in total dose (71–208 nGy/hr), are compared to this “radiation-deprived” treatment. In addition, RNA-seq data of Caenorhabditis elegans nematode from similar radiation treatments was analyzed by three of the software packages. Results: In E. coli, the four software programs identified one of the supplementary sources of radiation (KCl) to evoke about 5 times more transcribed genes than the minus-radiation treatment (69–114 differentially expressed genes, DEGs), and so the rest of the analyses used this KCl vs “Minus” comparison. After imposing a 30-read minimum cutoff, one of the DNAStar options shared two of the three steps (mapping, normalization, and statistic) with Partek Flow (they both used median of ratios to normalize and the DESeq2 statistical package), and these two programs identified the highest number of DEGs in common with each other (53). In contrast, when the programs used different approaches in each of the three steps, between 31 and 40 DEGs were found in common. Regarding the extent of expression differences, three of the four programs gave high fold-change results (15–178 fold), but one (DNAstar’s DESeq2) resulted in more conservative fold-changes (1.5–3.5). In a parallel study comparing three qPCR commercial validation software programs, these programs also gave variable results as to which genes were significantly regulated. Similarly, the C. elegans analysis showed exaggerated fold-changes in CLC and DNAstar’s edgeR while DNAstar-D was more conservative. Conclusions: Regarding the extent of expression (fold-change), and considering the subtlety of the very low level radiation treatments, in E. coli three of the four programs gave what we consider exaggerated fold-change results (15 – 178 fold), but one (DNAstar’s DESeq2) gave more realistic fold-changes (1.5–3.5). When RT-qPCR validation comparisons to transcriptome results were carried out, they supported the more conservative DNAstar-D’s expression results. When another model organism’s (nematode) response to these radiation differences was similarly analyzed, DNAstar-D also resulted in the most conservative expression patterns. Therefore, we would propose DESeq2 (“DNAstar-D”) as an appropriate software tool for differential gene expression studies for treatments expected to give subtle transcriptome responses.

59 BASIC BIOLOGICAL SCIENCES↗

Clean Energy Employment Impacts and Occupational Analyses: Building Envelope & Electrification Upgrades

Under the Communities LEAP (Local Energy Action Program) Pilot, the U.S. Department of Energy (DOE) is working with a group of LEAP communities to bolster workforce development initiatives and support the growth of local employment in the clean energy economy. Through this resource, DOE seeks to provide participating communities with an overview of the types of occupations that could be affected by investments in building envelope and building electrification upgrades in existing building stock - differentiating between the residential and commercial sectors. This occupational analysis is intended to be used by communities to inform workforce initiatives (e.g., local training and education programs) as they implement building energy efficiency and electrification projects at scale.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Insights of Molecular Mechanism of Xylem Development in Five Black Poplar Cultivars

Black poplar (Populus deltoides, P. nigra, and their hybrids) is the main poplar cultivars in China. It offers interesting options of large-scale biomass production for bioenergy due to its rapid growth and high yield. Poplar wood properties were associated with chemical components and physical structures during wood formation. In this study, five poplar cultivars, P. euramericana ‘Zhonglin46’ (Pe1), P. euramericana ‘Guariento’ (Pe2), P. nigra ‘N179’ (Pn1), P. deltoides ‘Danhong’ (Pd1), and P. deltoides ‘Nanyang’ (Pd2), were used to explore the molecular mechanism of xylem development. We analyzed the structural differences of developing xylem in the five cultivars and profiled the transcriptome-wide gene expression patterns through RNA sequencing. The cross sections of the developing xylem showed that the cell wall thickness of developed fiber in Pd1 was thickest and the number of xylem vessels of Pn1 was the least. A total of 10,331 differentially expressed genes were identified among 10 pairwise comparisons of the five cultivars, most of them were related to programmed cell death and secondary cell wall thickening. K-means cluster analysis and Gene Ontology enrichment analysis showed that the genes highly expressed in Pd1 were related to nucleotide decomposition, metabolic process, transferase, and microtubule cytoskeleton; whereas the genes highly expressed in Pn1 were involved in cell wall macromolecule decomposition and polysaccharide binding processes. Based on a weighted gene co-expression network analysis, a large number of candidate regulators for xylem development were identified. And their potential regulatory roles to cell wall biosynthesis genes were validated by a transient overexpression system. This study provides a set of promising candidate regulators for genetic engineering to improve feedstock and enhance biofuel conversion in the bioenergy crop Populus.

59 BASIC BIOLOGICAL SCIENCES↗

Excited Baryon Structure Using Exclusive Reactions with CLAS12

Studying excited nucleon structure through exclusive electroproduction reactions is an important avenue for exploring the nature of the non-perturbative strong interaction. Electrocouplings for $N^*$ states in the mass range below 1.8~GeV have been determined from analyses of CLAS $\pi N$, $\eta N$, and $\pi \pi N$ data. This work made it clear that consistency of independent analyses of exclusive channels with different couplings and non-resonant backgrounds but the same $N^*$ electro-excitation amplitudes, is essential to have confidence in the extracted results. In terms of hadronic coupling, many high-lying $N^*$ states preferentially decay through the $\pi \pi N$ channel instead of $\pi N$. Data from the $KY$ channels will therefore be critical to provide an independent analysis to compare the extracted electrocouplings for the high-lying $N^*$ states against those determined from the $\pi N$ and $\pi \pi N$ channels. A program to study excited $N^*$ decays to non-strange and strange exclusive final states using CLAS12 will measure differential cross sections to be used as input to extract the $\gamma_vNN^*$ transition form factors for the most prominent $N^*$ states in the range of invariant energy $W$ up 3~GeV in the virtually unexplored domain of momentum transfers $Q^2$ up to 12~GeV$^2$.

Carman, Daniel↗

Singular linear-quadratic control problem for systems with linear delay

A singular linear-quadratic optimization problem on the trajectories of non-autonomous linear differential equations with linear delay is considered. The peculiarity of this problem is the fact that this problem has no solution in the class of integrable controls. To ensure the existence of solutions is required to expand the class of controls including controls with impulse components. Dynamical systems with linear delay are used to describe the motion of pantograph from the current collector with electric traction, biology, etc. It should be noted that for practical problems fact singularity criterion of quality is quite commonly occurring, and therefore the study of these problems is surely important. For the problem under discussion optimal programming control contained impulse components at the initial and final moments of time is constructed under certain assumptions on the functional and the right side of the control system.

97 MATHEMATICS AND COMPUTING↗

Validation of transfer map calculation for electrostatic deflectors in the code COSY INFINITY

The code COSY INFINITY uses a beamline coordinate system with a Frenet–Serret frame relative to the reference particle, and calculates differential algebra-valued transfer maps by integrating the ODEs of motion in the respective vector space over a differential algebra (DA). In our work, we described and performed computation of the DA transfer map of an electrostatic spherical deflector in a laboratory coordinate system using two conventional methods: (1) by integrating the ODEs of motion using a numerical integrator and (2) by computing analytically and in closed form the properties of the respective elliptical orbits from Kepler theory. We compared the resulting transfer maps with (3) the DA transfer map of COSY INFINITY’s built-in electrostatic spherical deflector element ESP and (4) the transfer map of the electrostatic spherical deflector computed using the program GIOS, which uses analytic formulas from a paper1 by Hermann Wollnik regarding second-order aberrations. In addition to the electrostatic spherical deflector, we studied an electrostatic cylindrical deflector, where the Kepler theory is not applicable. We computed the DA transfer map by the ODE integration method (1), and we compared it with the transfer maps by (3) COSY INFINITY’s built-in electrostatic cylindrical deflector element ECL and (4) GIOS. The transfer maps of electrostatic spherical and cylindrical deflectors obtained using the direct calculation methods (1) and (2) are in excellent agreement with those computed using (3) COSY INFINITY. On the other hand, we found a significant discrepancy with (4) the program GIOS.

43 PARTICLE ACCELERATORS↗