Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Physical quantities”

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 19 records

Calculation of surface and top of atmosphere radiative fluxes from physical quantities based on ISCCP data sets. 1: Method and sensitivity to input data uncertainties

The largest uncertainty in upwelling shortwave (SW) fluxes (approximately equal 10-15 W/m(exp 2), regional daily mean) is caused by uncertainties in land surface albedo, whereas the largest uncertainty in downwelling SW at the surface (approximately equal 5-10 W/m(exp 2), regional daily mean) is related to cloud detection errors. The uncertainty of upwelling longwave (LW) fluxes (approximately 10-20 W/m(exp 2), regional daily mean) depends on the accuracy of the surface temperature for the surface LW fluxes and the atmospheric temperature for the top of atmosphere LW fluxes. The dominant source of uncertainty is downwelling LW fluxes at the surface (approximately equal 10-15 W/m(exp 2)) is uncertainty in atmospheric temperature and, secondarily, atmospheric humidity; clouds play little role except in the polar regions. The uncertainties of the individual flux components and the total net fluxes are largest over land (15-20 W/m(exp 2)) because of uncertainties in surface albedo (especially its spectral dependence) and surface temperature and emissivity (including its spectral dependence). Clouds are the most important modulator of the SW fluxes, but over land areas, uncertainties in net SW at the surface depend almost as much on uncertainties in surface albedo. Although atmospheric and surface temperature variations cause larger LW flux variations, the most notable feature of the net LW fluxes is the changing relative importance of clouds and water vapor with latitude. Uncertainty in individual flux values is dominated by sampling effects because of large natrual variations, but uncertainty in monthly mean fluxes is dominated by bias errors in the input quantities.

Zhang, Y.-C.↗

Calculation of surface and top of atmosphere radiative fluxes from physical quantities based on ISCCP data sets. 2: Validation and first results

We use global, multiyear observations of the properties of clouds, the atmosphere, and the surface to calculate global shortwave (SW) and longwave (LW) fluxes at the top of the atmosphere and at the surface at a resolution of 280 km and 3 hours for every third month from April 1985 to January 1989. Our validation studies suggest that the specification of cloud effects is no longer the dominant uncertainty in reconstructing the radiative fluxes at the top of atmosphere and at the surface. Rather cloud property uncertainties are now roughly equal contributors to the flux uncertainty, along with surface and atmospheric properties. The resulting SW and LW flux data sets suggest the following conclusions: (1) The net SW heating of Earth appears predominantly at the surface, whereas the net LW cooling arises predominantly from the atmosphere. The net cooling effect of clouds on top of atmospheric radiation appears primarily at the surface rather than in the atmosphere. (2) Clouds have almost no net effect on the global mean radiation balance of the atmosphere, but they enhance the latitudinal gradient in the LW cooling and reinforce the radiative forcing for the mean atmospheric circulation. Clouds act to mute seasonal contrasts however. (3) Clouds enhance the land-ocean contrasts of the atmospheric cooling, reinforcing the growth of standing eddy motions; but reduce land-ocean contrasts of the surface heating.

Rossow, W. B.↗

Quantifying the unknown impact of segmentation uncertainty on image-based simulations

Abstract Image-based simulation, the use of 3D images to calculate physical quantities, relies on image segmentation for geometry creation. However, this process introduces image segmentation uncertainty because different segmentation tools (both manual and machine-learning-based) will each produce a unique and valid segmentation. First, we demonstrate that these variations propagate into the physics simulations, compromising the resulting physics quantities. Second, we propose a general framework for rapidly quantifying segmentation uncertainty. Through the creation and sampling of segmentation uncertainty probability maps, we systematically and objectively create uncertainty distributions of the physics quantities. We show that physics quantity uncertainty distributions can follow a Normal distribution, but, in more complicated physics simulations, the resulting uncertainty distribution can be surprisingly nontrivial. We establish that bounding segmentation uncertainty can fail in these nontrivial situations. While our work does not eliminate segmentation uncertainty, it improves simulation credibility by making visible the previously unrecognized segmentation uncertainty plaguing image-based simulation.

42 ENGINEERING↗

Measuring Economy-Wide Circularity of the United States: An Input-Output Model in Mass Units

The ideal of creating closed material cycles by transforming the way we make, use and repurpose goods has become known as the circular economy. Besides its visionary appeal, material efficiency strategies central to circular economy can allow us to meet decarbonization goals that are otherwise out of reach. Production of basic materials such as cement, iron and steel and petrochemicals is one of the largest drivers of greenhouse gas emissions. Measuring circularity, and understanding its relationship with other sustainability metrics, however, is still difficult due to lack of data in mass units covering the entire economic system. Input-output (I-O) tables were originally created as a means of tracking the monetary flows that represent exchanges of goods and services within an economy, but have found additional uses in life cycle assessment and material flow accounts. In the traditional approach of environmentally extended I-O tables, emissions and energy data augment monetary flows. This approach can lead to price effects that distort physical quantities. Also, circular economy strategies and their goals relate to mass, not monetary flows. Therefore, it is best to simulate them using physical quantities. With these issues in mind, we have developed I-O tables in mass units to measure the flow of goods in the U.S. economy. This allows us to better measure how circular the economy really is, and how policy changes to affect this circularity may also affect decarbonization goals. Improving knowledge of these linkages could allow manufacturers to make changes in their production to achieve their sustainability and decarbonization targets. Our tool also provides a standard approach and public repository for data in physical units, making such data more available, useful, and meaningful. Following an established set of material flow metrics used by the European Union, we have used our tool to calculate material footprints over time differentiated by oil and gas versus other extractive industries. This approach also shows the relative trade balances of the U.S. in materials. We have developed a case study for the iron and steel sector, showing how different decarbonization scenarios affect not only economy-wide greenhouse gas emissions, but also total material use.

circular economy↗

Demonstration of Quantum Energy Teleportation on Superconducting Quantum Hardware

Teleporting physical quantities to remote locations remains a key challenge for quantum information science and technology. Quantum teleportation has enabled the transfer of quantum information, but teleportation of quantum physical quantities has not yet been realized. Here we report the realization and observation of quantum energy teleportation on real superconducting quantum hardware. We achieve this by using several IBM superconducting quantum computers. The results are consistent with the exact solution of the theory and are improved by the mitigation of measurement error. Quantum energy teleportation requires only local operations and classical communication. Furthermore, our results provide a realistic benchmark that is fully achievable with current quantum computing and communication technologies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Correlation function distributions for O ( N ) lattice field theories in the disordered phase

Numerical computations in strongly interacting quantum field theories are often performed using Monte Carlo sampling methods. A key task in these calculations is to estimate the value of a given physical quantity from the distribution of stochastic samples that are generated using the Monte Carlo method. Typically, the sample mean and sample variance are used to define the expectation values and uncertainties of computed quantities. However, the Monte Carlo sample distribution contains more information than these basic properties, and it is useful to investigate it more generally. In this work, the exact form of the probability distributions of two-point correlation functions at zero momentum in O ( N ) lattice field theories in the disordered phase and in infinite volume are determined. These distributions allow for a robust investigation of the efficacy of the Monte Carlo sampling procedure and are shown also to allow for improved estimators of the target physical quantity to be constructed. The theoretical expectations are shown to agree with numerical calculations in the O ( 2 ) model. Published by the American Physical Society 2024

Astronomy & Astrophysics↗

A study of the use of abstract types for the representation of engineering units in integration and test applications

Physical quantities using various units of measurement can be well represented in Ada by the use of abstract types. Computation involving these quantities (electric potential, mass, volume) can also automatically invoke the computation and checking of some of the implicitly associable attributes of measurements. Quantities can be held internally in SI units, transparently to the user, with automatic conversion. Through dimensional analysis, the type of the derived quantity resulting from a computation is known, thereby allowing dynamic checks of the equations used. The impact of the possible implementation of these techniques in integration and test applications is discussed. The overhead of computing and transporting measurement attributes is weighed against the advantages gained by their use. The construction of a run time interpreter using physical quantities in equations can be aided by the dynamic equation checks provided by dimensional analysis. The effects of high levels of abstraction on the generation and maintenance of software used in integration and test applications are also discussed.

Johnson, Charles S.↗

Strictly Enforcing Invertibility and Conservation in CNN-Based Super Resolution for Scientific Datasets

Abstract Recently, deep convolutional neural networks (CNNs) have revolutionized image “super resolution” (SR), dramatically outperforming past methods for enhancing image resolution. They could be a boon for the many scientific fields that involve imaging or any regularly gridded datasets: satellite remote sensing, radar meteorology, medical imaging, numerical modeling, and so on. Unfortunately, while SR-CNNs produce visually compelling results, they do not necessarily conserve physical quantities between their low-resolution inputs and high-resolution outputs when applied to scientific datasets. Here, a method for “downsampling enforcement” in SR-CNNs is proposed. A differentiable operator is derived that, when applied as the final transfer function of a CNN, ensures the high-resolution outputs exactly reproduce the low-resolution inputs under 2D-average downsampling while improving performance of the SR schemes. The method is demonstrated across seven modern CNN-based SR schemes on several benchmark image datasets, and applications to weather radar, satellite imager, and climate model data are shown. The approach improves training time and performance while ensuring physical consistency between the super-resolved and low-resolution data. Significance Statement Recent advancements in using deep learning to increase the resolution of images have substantial potential across the many scientific fields that use images and image-like data. Most image super-resolution research has focused on the visual quality of outputs, however, and is not necessarily well suited for use with scientific data where known physics constraints may need to be enforced. Here, we introduce a method to modify existing deep neural network architectures so that they strictly conserve physical quantities in the input field when “super resolving” scientific data and find that the method can improve performance across a wide range of datasets and neural networks. Integration of known physics and adherence to established physical constraints into deep neural networks will be a critical step before their potential can be fully realized in the physical sciences.

54 ENVIRONMENTAL SCIENCES↗

A unified neural-network framework for nucleon imaging from numerical simulations of QCD

Parton distributions encode the momentum-space structure and, in their generalizations, the spatial tomography of quarks and gluons inside hadrons, the building blocks of visible matter. We present a unified neural-network approach that learns these distributions directly from matrix elements calculated via numerical simulations of quantum chromodynamics (QCD) on the lattice by fitting two complementary inputs simultaneously: data matched to physical quantities via known momentum-space and coordinate-space formalisms. Utilizing data from both methods stabilizes the extraction and mitigates biases that can arise when either is used alone. We validate the method on controlled mock data and apply it to lattice-QCD matrix elements to extract parton distribution functions (PDFs). We show benefits of such an approach for determining the physical quantities. We further extend the framework to zero-skewness generalized parton distributions and demonstrate nucleon tomography within the same neural-network parameterization. Our results provide an adaptable and systematically improvable approach for extracting partonic distributions from Euclidean correlators. It can incorporate polarization, additional channels, and future experimental constraints from current and future facilities, such as the Electron-Ion Collider.

Hadronic Spectroscopy↗

Effects of Different Closure Choices in Core-collapse Supernova Simulations

The two-moment method is widely used to approximate the full neutrino transport equation in core-collapse supernova (CCSN) simulations, and different closures lead to subtle differences in the simulation results. In this paper, we compare the effects of closure choices on various physical quantities in 1D and 2D time-dependent CCSN simulations with our multigroup radiation hydrodynamics code FORNAX. We find that choices of the third-order closure relations influence the time-dependent simulations only slightly. Choices of the second-order closure relation have larger consequences than choices of the third-order closure, but these are still small compared to the remaining variations due to ambiguities in some physical inputs such as the nuclear equation of state. We also find that deviations in Eddington factors are not monotonically related to deviations in physical quantities, which means that simply comparing the Eddington factors does not inform one concerning which closure is better.

79 ASTRONOMY AND ASTROPHYSICS↗

Distorted-Toroidal Flux Rope Model

The 3D characterization of magnetic flux ropes observed in the heliosphere has been a challenging task for decades. This is mainly due to the limitations on inferring the 3D global topology and physical properties from the 1D time series from any spacecraft. To advance our understanding of magnetic flux ropes whose configuration departs from the typical stiff geometries, here we present an analytical solution for a 3D flux rope model with an arbitrary cross section and a toroidal global shape. This constitutes the next level of complexity following the elliptic-cylindrical (EC) geometry. The mathematical framework was established by Nieves-Chinchilla et al. with the EC flux rope model, which describes a magnetic topology with an elliptical cross section as a first approach to changes in the cross section. In the distorted-toroidal flux rope model, the cross section is described by a general function. The model is completely described by a nonorthogonal geometry and the Maxwell equations can be consistently solved to obtain the magnetic field and relevant physical quantities. As a proof of concept, this model is generalized in terms of the radial dependence of current density components. The last part of this paper is dedicated to a specific function, F (φ) = δ (1 - λ cosφ), to illustrate possibilities of the model. This model paves the way toward the investigation of complex distortions of magnetic structures in the solar wind. Future investigations will explore these distortions in depth by analyzing specific events; studying implications for physical quantities, such as magnetic fluxes, helicity, or energy; and evaluating the force balance with the ambient solar wind that allows such distortions.

Solar wind↗

The relation between quasars’ optical spectra and variability

Abstract Brightness variation is an essential feature of quasars, but its mechanism and relationship to other physical quantities are not understood well. We aimed to find the relationship between the optical variability and spectral features to reveal the regularity behind the random variation. It is known that a quasar’s Fe ii/Hβ flux ratio and equivalent width of [O iii]5007 are negatively correlated; this is called Eigenvector 1. In this work, we visualized the relationship between the position on this Eigenvector 1 (EV1) plane and how the brightness of the quasars had changed after ∼10 yr. We conducted three analyses, using a different quasar sample in each. The first analysis showed the relation between the quasars’ distributions on the EV1 plane and how much they had changed brightness, using 13438 Sloan Digital Sky Survey quasars. This result shows how brightness changes later are clearly related to the position on the EV1 plane. In the second analysis, we plotted the sources reported as “changing-look quasars” (or “changing-state quasars”) on the EV1 plane. This result shows that the position on the EV1 plane corresponds to the activity level of each source, and the bright or dim states of them are distributed on the opposite sides divided by the typical quasar distribution. In the third analysis, we examined the transition vectors on the EV1 plane using sources with multiple-epoch spectra. This result shows that the brightening and dimming sources move on a similar path and they reach a position corresponding to the opposite activity level. We also found this trend is opposite to the empirical rule that $R_{\rm {Fe\, \small {II}}}$ positively correlated with the Eddington ratio, which has been proposed based on the trends of a large number of quasars. From all these analyses, it is indicated that quasars tend to oscillate between both sides of the distribution ridge on the EV1 plane; each of them corresponds to a dim state and a bright state. This trend in optical variation suggests that significant brightness changes, such as changing-look quasars, are expected to repeat.

Astronomy & Astrophysics↗

SYMBA: symbolic computation of squared amplitudes in high energy physics with machine learning

The cross section is one of the most important physical quantities in high-energy physics and the most time consuming to compute. While machine learning has proven to be highly successful in numerical calculations in high-energy physics, analytical calculations using machine learning are still in their infancy. In this work, we use a sequence-to-sequence model, specifically, a transformer, to compute a key element of the cross section calculation, namely, the squared amplitude of an interaction. We show that a transformer model is able to predict correctly 97.6% and 99% of squared amplitudes of quantum chromodynamics and quantum electrodynamics processes, respectively, at a speed that is up to orders of magnitude faster than current symbolic computation frameworks. We discuss the performance of the current model, its limitations and possible future directions for this work.

97 MATHEMATICS AND COMPUTING↗

Comparison of Local Scale Measured and Modeled Brightness Temperatures and Snow Parameters from the CLPX 2003 by Means of a Dense Medium Radiative Transfer Theory Model

Microwave remote sensing offers distinct advantages for observing the cryosphere. Solar illumination is not required, and spatial and temporal coverage are excellent from polar-orbiting satellites. Passive microwave measurements are sensitive to the two most useful physical quantities for many hydrological applications: physical temperature and water content/state. Sensitivity to the latter is a direct result of the microwave sensitivity to the dielectric properties of natural media, including snow, ice, soil (frozen or thawed), and vegetation. These considerations are factors motivating the development of future cryospheric satellite remote sensing missions, continuing and improving on a 26-year microwave measurement legacy. Perhaps the biggest issues regarding the use of such satellite measurements involve how to relate parameter values at spatial scales as small as a hectare to observations with sensor footprints that may be up to 25 x 25 km. The NASA Cold-land Processes Field Experiment (CLPX) generated a dataset designed to enhance understanding of such scaling issues. CLPX observations were made in February (dry snow) and March (wet snow), 2003 in Colorado, USA, at scales ranging from plot scale to 25 x 25 km satellite footprints. Of interest here are passive microwave observations from ground-based, airborne, and satellite sensors, as well as meteorological and snowpack measurements that will enable studies of the effects of spatial heterogeneity of surface conditions on the observations. Prior to performing such scaling studies, an evaluation of snowpack forward modelling at the plot scale (least heterogeneous scale) is in order. This is the focus of this paper. Many forward models of snow signatures (brightness temperatures) have been developed over the years. It is now recognized that a dense medium radiative transfer (DMRT) treatment represents a high degree of physical fidelity for snow modeling, yet dense medium models are particularly sensitive to snowpack structural parameters such as grain size, density, and depth---parameters that may vary substantially within a snowpack. Microwave radiometric data and snow pit measurements collected at the Local-Scale Observation Site (LSOS) during the third Intensive Observation Period (IOP3) of the CLPX have been used to test the capabilities of a DMRT model using the Quasi Crystalline Approximation with Coherent Potential (QCA-CP). The radiometric measurements were made by the University of Tokyo s Ground Based Microwave Radiometer-7 (GBMR-7) system. We evaluate the degree to which a DMRT-based model can accurately reproduce the GBMR-7 brightness temperatures at different frequencies and incidence angles.

Tedescol, Marco↗

Physics constrained learning for data-driven inverse modeling from sparse observations

Deep neural networks (DNN) have been used to model nonlinear relations between physical quantities. Those DNNs are embedded in physical systems described by partial differential equations (PDE) and trained by minimizing a loss function that measures the discrepancy between predictions and observations in some chosen norm. This loss function often includes the PDE constraints as a penalty term when only sparse observations are available. As a result, the PDE is only satisfied approximately by the solution. However, the penalty term typically slows down the convergence of the optimizer for stiff problems. We present a new approach that trains the embedded DNNs while numerically satisfying the PDE constraints. We develop an algorithm that enables differentiating both explicit and implicit numerical solvers in reverse-mode automatic differentiation. This allows the gradients of the DNNs and the PDE solvers to be computed in a unified framework. We demonstrate that our approach enjoys faster convergence and better stability in relatively stiff problems compared to the penalty method. Furthermore, our approach allows for the potential to solve and accelerate a wide range of data-driven inverse modeling, where the physical constraints are described by PDEs and need to be satisfied accurately.

97 MATHEMATICS AND COMPUTING↗

Modeling of nonequilibrium space plasma flows

Godunov-type numerical solution of the 20 moment plasma transport equations. One of the centerpieces of our proposal was the development of a higher order Godunov-type numerical scheme to solve the gyration dominated 20 moment transport equations. In the first step we explored some fundamental analytic properties of the 20 moment transport equations for a low b plasma, including the eigenvectors and eigenvalues of propagating disturbances. The eigenvalues correspond to wave speeds, while the eigenvectors characterize the transported physical quantities. In this paper we also explored the physically meaningful parameter range of the normalized heat flow components. In the second step a new Godunov scheme type numerical method was developed to solve the coupled set of 20 moment transport equations for a quasineutral single-ion plasma. The numerical method and the first results were presented at several national and international meetings and a paper describing the method has been published in the Journal of Computational Physics. To our knowledge this is the first numerical method which is capable of producing stable time-dependent solutions to the full 20 (or 16) moment set of transport equations, including the full heat flow equation. Previous attempts resulted in unstable (oscillating) solutions of the heat flow equations. Our group invested over two man-years into the development and implementation of the new method. The present model solves the 20 moment transport equations for an ion species and thermal electrons in 8 domain extending from a collision dominated to a collisionless region (200 km to 12,000 km). This model has been applied to study O+ acceleration due to Joule heating in the lower ionosphere.

Gombosi, Tamas↗

TPSAS-NF1676L-30451-DND

The fundamental laws of the physics of fluid flow in the continuum regime are embedded in a somewhat "simple" set of partial differential equations describing the transport of mass, momentum, energy, and/orentropy. These basic equations can be manipulated to construct other important physical quantities such as vorticity that help further elucidate flow physics. Additionally, these equations can be nondimensionalized to identify important nondimensional groups governing compressible, reacting fluid flows, most prominently Reynolds, Mach, Froude, Prandtl, Lewis, Schmidt, Damkohler, and Zeldovich numbers. In this lecture, we present these governing equations and discuss the significance of the different terms and nondimensional groups. Furthermore, since these equations cannot be solved directly on current supercomputers for turbulent flows at Reynolds numbers of engineering interest, a scale reduction by averaging or filtering is introduced, which leads to a "closure" problem. The closure is then provided by introducing some common models. The current introduction is by no means complete, but it and the references herein may serve as a good starting point for learning about the most common approaches to turbulent reacting flows.

Tomasz (Tom) G Drozda↗