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 379 records · Page 21

Real-Time Adaptive Machine Learning Platform

As part of the Cyclotron Road program, Dauntless.io, Inc. extended its existing ultrafast (0.00001 – 0.001 second) real-time, adaptive machine learning platform for modeling & control of complex, high-value physical systems. Objectives included ability to scale ubiquitously & remain differentiated on a 5+ year horizon in the highly competitive, AI / machine learning market — a rapidly growing market that (overall) is estimated to be worth $191 B by 2024 at a CAGR of 37%.

97 MATHEMATICS AND COMPUTING↗

Agricultural Impacts on Nitrogen Cycling: Climate and Air Pollution (Final Technical Report)

We have successfully added a parameterization of ammonia emissions from synthetic fertilizer and manure from agriculture to Earth System Models. We have evaluated how emissions change in the future and the processes driving these changes. We have examined the impact of agricultural ammonia emissions on present-day air quality and have evaluated the impact of resolving the emissions with high temporal resolution. We have also examined the impact of how Earth System Models treat the cycling of emissions in soils. Our code, the Flows of Agricultural Nitrogen – version 2 – model (FANv2) has been run in both the CLM framework and the ELM framework. To my knowledge the code has not yet been fully implemented into ELM. Volatilization of ammonia (NH3) from fertilizers and livestock wastes forms a significant pathway of nitrogen losses in agricultural ecosystems and constitutes the largest source of atmospheric emissions of NH3. We have provided a major update to the process model FAN (Flow of Agricultural Nitrogen), which evaluates NH3 emissions interactively within an Earth system model; in this work, the Community Earth System Model (CESM) is used. The updated version (FANv2) includes a more detailed treatment of both physical and agricultural processes, which allows the model to differentiate between the volatilization losses from animal housings, manure storage, grazed pastures, and the application of manure and different types of mineral fertilizers. The modeled ammonia emissions are first evaluated at a local scale against experimental data for various types of fertilizers and manure, and they are subsequently run globally to evaluate NH3 emissions for 2010–2015 based on gridded datasets of fertilizer use and livestock populations. Comparison of regional emissions shows that FANv2 agrees with previous inventories for North America and Europe and is within the range of previous inventories for China. However, due to higher NH3 emissions in Africa, India, and Latin America, the global emissions simulated by FANv2(48 Tg N) are 30 %–40 % higher than in the existing inventories.

54 ENVIRONMENTAL SCIENCES↗

Report for the ASCR Workshop on Basic Research Needs in Quantum Computing and Networking - 2023

Employing quantum mechanical resources in computing and networking opens the door to new computation and communication models and potential disruptive advantages over classical counterparts. However, quantifying and realizing such advantages face extensive scientific and engineering challenges. Investments by the Department of Energy (DOE) have driven progress toward addressing such challenges. Quantum algorithms have been recently developed, in some cases offering asymptotic exponential advantages in speed or accuracy, for fundamental scientific problems such as simulating physical systems, solving systems of linear equations, or solving differential equations. Empirical demonstrations on nascent quantum hardware suggest better performance than classical analogs on specialized computational tasks favorable to the quantum computing systems. However, demonstration of an end-to-end, substantial and rigorously quantifiable quantum performance advantage over classical analogs remains a grand challenge, especially for problems of practical value. The definition of requirements for quantum technologies to exhibit scalable, rigorous, and transformative performance advantages for practical applications also remains an outstanding open question, namely, what will be required to ultimately demonstrate practical quantum advantage?

97 MATHEMATICS AND COMPUTING↗

Report for the ASCR Workshop on Basic Research Needs in Quantum Computing and Networking - 2023

Employing quantum mechanical resources in computing and networking opens the door to new computation and communication models and potential disruptive advantages over classical counterparts. However, quantifying and realizing such advantages face extensive scientific and engineering challenges. Investments by the Department of Energy (DOE) have driven progress toward addressing such challenges. Quantum algorithms have been recently developed, in some cases offering asymptotic exponential advantages in speed or accuracy, for fundamental scientific problems such as simulating physical systems, solving systems of linear equations, or solving differential equations. Empirical demonstrations on nascent quantum hardware suggest better performance than classical analogs on specialized computational tasks favorable to the quantum computing systems. However, demonstration of an end-to-end, substantial and rigorously quantifiable quantum performance advantage over classical analogs remains a grand challenge, especially for problems of practical value. The definition of requirements for quantum technologies to exhibit scalable, rigorous, and transformative performance advantages for practical applications also remains an outstanding open question, namely, what will be required to ultimately demonstrate practical quantum advantage?

97 MATHEMATICS AND COMPUTING↗

An improved mechanistic model for ammonia volatilization in Earth system models: Flow of Agricultural Nitrogen version 2 (FANv2)

Abstract. Volatilization of ammonia (NH 3 ) from fertilizers and livestock wastes forms a significant pathway of nitrogen losses in agricultural ecosystems and constitutes the largest source of atmospheric emissions of NH 3 . This paper describes a major update to the process model FAN (Flow of Agricultural Nitrogen), which evaluates NH 3 emissions interactively within an Earth system model; in this work, the Community Earth System Model (CESM) is used. The updated version (FANv2) includes a more detailed treatment of both physical and agricultural processes, which allows the model to differentiate between the volatilization losses from animal housings, manure storage, grazed pastures, and the application of manure and different types of mineral fertilizers. The modeled ammonia emissions are first evaluated at a local scale against experimental data for various types of fertilizers and manure, and they are subsequently run globally to evaluate NH3 emissions for 2010–2015 based on gridded datasets of fertilizer use and livestock populations. Comparison of regional emissions shows that FANv2 agrees with previous inventories for North America and Europe and is within the range of previous inventories for China. However, due to higher NH 3 emissions in Africa, India, and Latin America, the global emissions simulated by FANv2 (48 Tg N) are 30 %–40 % higher than in the existing inventories.

Vira, Julius↗

Automated controller for liquid-cooled garments

Automated controller is governed by inlet temperature of coolant and temperature differential across garment. Controller eliminates restrictive body attachments of physical sensors, while regulating temperature more efficiently than manual adjustment.

Kuznetz, L. H.↗

Mass and energy balance in the 1973 August 9 flare

The mass and energy balance of the thermal plasma during the decay phase of the solar flare of August 9, 1973, are studied. The analysis is based on observationally determined values for the differential emission measure, density, turbulent and bulk velocities, and physical dimensions. The total particle content and total thermal energy content of the flare plasmas with temperatures above 100,000 K and their variation with time are calculated. The particle loss and the energy losses through radiation, conduction, and convection are evaluated. The decrease in total particle content can be accounted for by the convective losses through the loop footprints at 100,000 K. Radiation is the dominant energy loss mechanism although convective losses at 100,000 K can be important. Conductive losses at 100,000 K into cooler chromospheric material appear to be negligible. The decrease in the total energy content during the decay phase is equal to the sum of the energy losses over the period of observation. No requirement is found for continued heating during the decay phase.

Dere, K. P.↗

Chaos and nonlinear dynamics of single-particle orbits in a magnetotaillike magnetic field

The properties of charged-particle motion in Hamiltonian dynamics are studied in a magnetotaillike magnetic field configuration. It is shown by numerical integration of the equation of motion that the system is generally nonintegrable and that the particle motion can be classified into three distinct types of orbits: bounded integrable orbits, unbounded stochastic orbits, and unbounded transient orbits. It is also shown that different regions of the phase space exhibit qualitatively different responses to external influences. The concept of 'differential memory' in single-particle distributions is proposed. Physical implications for the dynamical properties of the magnetotail plasmas and the possible generation of non-Maxwellian features in the distribution functions are discussed.

Chen, J.↗

Optimization methods and silicon solar cell numerical models

The goal of this project is the development of an optimization algorithm for use with a solar cell model. It is possible to simultaneously vary design variables such as impurity concentrations, front junction depth, back junctions depth, and cell thickness to maximize the predicted cell efficiency. An optimization algorithm has been developed and interfaced with the Solar Cell Analysis Program in 1 Dimension (SCAPID). SCAPID uses finite difference methods to solve the differential equations which, along with several relations from the physics of semiconductors, describe mathematically the operation of a solar cell. A major obstacle is that the numerical methods used in SCAPID require a significant amount of computer time, and during an optimization the model is called iteratively until the design variables converge to the value associated with the maximum efficiency. This problem has been alleviated by designing an optimization code specifically for use with numerically intensive simulations, to reduce the number of times the efficiency has to be calculated to achieve convergence to the optimal solution. Adapting SCAPID so that it could be called iteratively by the optimization code provided another means of reducing the cpu time required to complete an optimization. Instead of calculating the entire I-V curve, as is usually done in SCAPID, only the efficiency is calculated (maximum power voltage and current) and the solution from previous calculations is used to initiate the next solution.

Girardini, K.↗

Optimization methods and silicon solar cell numerical models

An optimization algorithm for use with numerical silicon solar cell models was developed. By coupling an optimization algorithm with a solar cell model, it is possible to simultaneously vary design variables such as impurity concentrations, front junction depth, back junction depth, and cell thickness to maximize the predicted cell efficiency. An optimization algorithm was developed and interfaced with the Solar Cell Analysis Program in 1 Dimension (SCAP1D). SCAP1D uses finite difference methods to solve the differential equations which, along with several relations from the physics of semiconductors, describe mathematically the performance of a solar cell. A major obstacle is that the numerical methods used in SCAP1D require a significant amount of computer time, and during an optimization the model is called iteratively until the design variables converge to the values associated with the maximum efficiency. This problem was alleviated by designing an optimization code specifically for use with numerically intensive simulations, to reduce the number of times the efficiency has to be calculated to achieve convergence to the optimal solution.

Girardini, K.↗

Dusty cometary atmospheres

Governing equations describing the gas and dust production and the accelerating dusty gas flow in a cometary atmosphere are compiled and compared with the information available on the inner coma of the Halley comet. Cometary dusty gas flow calculations, based on the Probstein (1968) approach, yield a first order differential equation for the gas velocity with only one physical solution possible, which is the one of a transonic accelerating gas 'wind'. According to this model, the comet wind starts subsonically at the nucleus, goes through the singular sonic point, and then accelerates further.

Gombosi, Tamas I.↗

The electron beam instability and turbulence theories

Extensions and practical applications of recent observations of electron beam-plasma interactions are investigated for the range of turbulence theories, extending from quasi-linear to strong turbulence theory, which have been developed on the basis of the Langmuir-wave excitation model. Electron foreshock observations have indicated that linear instability theory must encompass the excitation of waves whose frequencies are substantially different from those of the plasma frequency; the point of departure for such extensions should be a quantitative test of existing theories, and particle simulations conducive to such testing are presented. A step-by-step addition of physical considerations is used in such simulation studies to differentiate among nonlinear turbulence effects.

Dum, C. T.↗

Shear-induced opening of the coronal magnetic field

This work describes the evolution of a model solar corona in response to motions of the footpoints of its magnetic field. The mathematics involved is semianalytic, with the only numerical solution being that of an ordinary differential equation. This approach, while lacking the flexibility and physical details of full MHD simulations, allows for very rapid computation along with complete and rigorous exploration of the model's implications. We find that the model coronal field bulges upward, at first slowly and then more dramatically, in response to footpoint displacements. The energy in the field rises monotonically from that of the initial potential state, and the field configuration and energy appraoch asymptotically that of a fully open field. Concurrently, electric currents develop and concentrate into a current sheet as the limiting case of the open field is approached. Examination of the equations shows rigorously that in the asymptotic limit of the fully open field, the current layer becomes a true ideal MHD singularity.

Wolfson, Richard↗

Rapid Switching and Modulation by use of Coupled VCSELs

Devices incorporating coupled vertical-cavity surface-emitting lasers (VCSELs) have been proposed as means to effect unprecedentedly rapid modulation and/or switching of light beams in advanced optical communication and information-processing systems. A basic device according to the proposal would comprise two or more VCSELs (1) positioned so close to each other that they are coupled and (2) connected to a current source that biases them at a steady current above the threshold current for lasing. The operation of VCSELs under the conditions stated above has been simulated computationally with an algorithm that solves an approximation of the Maxwell-Bloch equations. (These are nonlinear coupled partial differential equations that model the relevant aspects of the physics of semiconductor lasers.) The results of the simulations show that the intensities of the beams oscillate, out of phase relative to each other, at a frequency that is typically of the order of tens of gigahertz. In particular, one simulation was performed for two coupled VCSELs, each having a circular aperture 5.6 microns in diameter, operating at a wavelength of 980 nm. The results of the simulation showed that the intensities of the two beams oscillated in opposite phase (see Figure 1) at a frequency of 42 GHz. In the far field, the beams were found to be separated by an angle of approx.=8deg. The intensities of the beams were also found to oscillate in the near field, 90deg out of phase with the oscillation in the far field. Figure 2 depicts a simple near-field switching device in which an assembly of microlenses would couple the two oscillating beams to two receivers.

Goorjian, Peter M.↗

Simulation Apps Bring Us Closer to Mars

My work with mathematical modeling and computer simulations began in earnest in 1987 when I signed on with the National Aeronautics and Space Administration, better known as NASA. At the time I was developing computer simulations to support the design of carbon dioxide (CO2) removal systems for the International Space Station (ISS) life support system. CO2 is a human metabolic waste, produced at a rate of a kilogram per day per crewmember. CO2 must be removed from the crew cabin, as it will quickly become detrimental to crew health. The systems used to remove CO2 are typically based on sorption processes, which include complex interactions of heat transfer, mass transfer, and gas flow through porous media. At this time, there were no commercial options for solving a set of coupled mathematical models such as partial differential equations (PDEs); you either force-fit the physics into a thermal analysis package or wrote your own code, complete with discretization, meshing, and solution algorithms. Unfortunately, coding CO2 removal processes from scratch did not allow an appropriate focus on the challenging yet crucial task of understanding and capturing the underlying physics via appropriate mathematical models. Configuration control was often simply ignored due to tight schedules, resulting in a code that would be quicker for a second party to rewrite than modify. In the early 2000s, I decided to move CO2 removal simulations to a platform with built-in meshing, solvers, and postprocessing capabilities, and that could solve user-defined multiphysics PDEs. The program chosen to meet these needs was COMSOL Multiphysics®. Along with freeing the engineer to focus on the underlying physics, a degree of configuration control was automatically achieved via a consistent user interface, thus allowing COMSOL users to share computer models. My team has developed simulations that are already providing a valuable debugging capability for the ISS CO2 removal system, and will provide guidance in upgrades to that system.

Knox, Jim↗

Differentiable, Learnable, Regionalized Process-Based Models With Multiphysical Outputs can Approach State-Of-The-Art Hydrologic Prediction Accuracy

Predictions of hydrologic variables across the entire water cycle have significant value for water resources management as well as downstream applications such as ecosystem and water quality modeling. Recently, purely data-driven deep learning models like long short-term memory (LSTM) showed seemingly insurmountable performance in modeling rainfall runoff and other geoscientific variables, yet they cannot predict untrained physical variables and remain challenging to interpret. Here, we show that differentiable, learnable, process-based models (called δ models here) can approach the performance level of LSTM for the intensively observed variable (streamflow) with regionalized parameterization. We use a simple hydrologic model HBV as the backbone and use embedded neural networks, which can only be trained in a differentiable programming framework, to parameterize, enhance, or replace the process-based model's modules. Without using an ensemble or post-processor, δ models can obtain a median Nash-Sutcliffe efficiency of 0.732 for 671 basins across the USA for the Daymet forcing data set, compared to 0.748 from a state-of-the-art LSTM model with the same setup. For another forcing data set, the difference is even smaller: 0.715 versus 0.722. Meanwhile, the resulting learnable process-based models can output a full set of untrained variables, for example, soil and groundwater storage, snowpack, evapotranspiration, and baseflow, and can later be constrained by their observations. Both simulated evapotranspiration and fraction of discharge from baseflow agreed decently with alternative estimates. The general framework can work with models with various process complexity and opens up the path for learning physics from big data.

54 ENVIRONMENTAL SCIENCES↗

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↗