Search NASA⌕ Search

SEARCH · Search NASA

Results for “Physical Algorithm”

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 577 records · Page 32

Resilient information and inference networks under mixed-trust sensing

With ubiquitous digitization, sensing, and computational intelligence deployed in increasingly more and broader domains, including critical infrastructure, potentially misleading and destabilizing effects of multimodal anomalies and adversarial behavior are growing in importance. Here, we develop randomized and reinforcement learning-based strategies for strategically recruiting and utilizing deployed (and, thus, vulnerable and potentially faulty and/or compromised) nodes from information and inference networks, while defending against adversaries that attempt to misguide assessments of inferred variables. Recognizing that, besides communication and other costs, sampling from any observable node can either provide true data or dangerously expose our inference to misinformation (without being easily distinguishable what actually happens), the proposed strategies proceed by progressively recruiting nodes and cautiously scaling their information contribution based on assumed, or, in our reinforcement learning approach, intelligently weighed trustworthiness, with the learning approach also considering network-wide, threat-inclusive risk/value tradeoffs. While avoiding the hardware, communication, analytical and computational burden of explicit redundancy, the proposed defensive schemes enable on-the-fly assessments of underlying processes, and system-wide situational awareness with demonstrable resilience against adversarial activities.

97 - MATHEMATICS AND COMPUTING↗

Static Subspace Approximation for Random Phase Approximation Correlation Energies: Implementation and Performance

Developing theoretical understanding of complex reactions and processes at interfaces requires using methods that go beyond semilocal density functional theory to accurately describe the interactions between solvent, reactants and substrates. Methods based on many-body perturbation theory, such as the random phase approximation (RPA), have previously been limited due to their computational complexity. However, this is now a surmountable barrier due to the advances in computational power available, in particular through modern GPU-based supercomputers. In this work, we describe the implementation of RPA calculations within BerkeleyGW and show its favorable computational performance on large complex systems relevant for catalysis and electrochemistry applications. Our implementation builds off of the static subspace approximation which, by employing a compressed representation of the frequency dependent polarizability, enables the evaluation of the RPA correlation energy with significant acceleration and systematically controllable accuracy. We find that the computational cost of calculating the RPA correlation energy scales only linearly with system size for systems containing up to 50 thousand bands, and is expected to scale quadratically thereafter. We also show excellent strong scaling results across several supercomputers, demonstrating the performance and portability of this implementation.

algorithmic development↗

Shadow molecular dynamics for flexible multipole models

Shadow molecular dynamics provide an efficient and stable atomistic simulation framework for flexible charge models with long-range electrostatic interactions. Shadow molecular dynamics simulations are driven by approximate “shadow” Born–Oppenheimer potentials for which the exact charges and forces are directly accessible without relying on costly (and approximate) iterative solvers. While previous implementations have been limited to atomic monopole charge distributions, we extend this approach to flexible multipole models. We derive detailed expressions for the shadow energy functions, potentials, and force terms, explicitly incorporating monopole–monopole, dipole–monopole, and dipole–dipole interactions. In our formulation, both atomic monopoles and atomic dipoles are treated as extended dynamical variables alongside the propagation of the nuclear degrees of freedom. We demonstrate that introducing the additional dipole degrees of freedom preserves the stability and accuracy previously seen in monopole-only shadow molecular dynamics simulations. In addition, we present a shadow molecular dynamics scheme where the monopole charges are held fixed while the dipoles remain flexible. Our extended shadow dynamics provide a framework for stable, computationally efficient, and versatile molecular dynamics simulations involving long-range interactions between flexible multipoles. This is of particular current interest in combination with machine-learned interatomic potentials, including long-range electrostatic interactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Performance of the ATLAS forward proton Time-of-Flight detector in Run 2

We present performance studies of the Time-of-Flight (ToF) subdetector of the ATLAS Forward Proton (AFP) detector at the LHC. Efficiencies and resolutions are measured using high-statistics data samples collected at low and moderate pile-up in 2017, the first year when the detectors were installed on both sides of the interaction region. While low efficiencies are observed, of the order of a few percent, the resolutions of the two ToF detectors measured individually are 21 ps and 28 ps, yielding an expected resolution of the longitudinal position of the interaction, z vtx , in the central ATLAS detector of 5.3 ± 0.6 mm. This is in agreement with the observed width of the distribution of the difference between z vtx , measured independently by the central ATLAS tracker and by the ToF detector, of 6.0 ± 2.0 mm.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Application of deep learning methods for beam size control during user operation at the Advanced Light Source

Past research at the Advanced Light Source (ALS) provided a proof-of-principle demonstration that deep learning methods could be effectively employed to compensate for the significant perturbations to the transverse electron beam size induced by user-controlled adjustments of the insertion devices. However, incorporating these methods into the ALS’ daily operations has faced notable challenges. The complexity of the system’s operational requirements and the significant upkeep demands has restricted their sustained application during user operation. Here, we introduce the development of a more robust neural network (NN)-based algorithm that utilizes a novel online fine-tuning approach and its systematic integration into the day-to-day machine operations. Our analysis emphasizes the process of NN model selection, demonstrates the superior performance of the NN-based method over traditional feedback methods, and examines the effectiveness and resilience of the new algorithm during user-operation scenarios. Published by the American Physical Society 2024

43 PARTICLE ACCELERATORS↗

Optimization of X-ray event screening using ground and in-orbit data for the Resolve instrument onboard the XRISM satellite

The X-Ray Imaging and Spectroscopy Mission (XRISM) satellite was successfully launched and put into a low-Earth orbit on September 6, 2023 (UT). The Resolve instrument onboard XRISM hosts an X-ray microcalorimeter detector, which was designed to achieve a high-resolution ( ≤ 7 eV FWHM at 6 keV), high-throughput, and non-dispersive spectroscopy over a wide energy range. It also excels in a low background with a requirement of < 2 × 10 -3 s -1 keV -1 (0.3 to 12.0 keV), which is equivalent to only one background event per spectral bin per 100-ks exposure. Event screening to discriminate X-ray events from background is a key to meeting the requirement. We present the result of the Resolve event screening using data sets recorded on the ground and in orbit based on the heritage of the preceding X-ray microcalorimeter missions, in particular, the Soft X-ray Spectrometer onboard ASTRO-H. We optimize and evaluate 19 screening items of three types based on (1) the event pulse shape, (2) relative arrival times among multiple events, and (3) good time intervals. We show that the initial screening, which is applied for science data products in the performance verification phase, reduces the background rate to 1.8 × 10 -3 s -1 keV -1 meeting the requirement. We further evaluate the additional screening utilizing the correlation among some pulse shape properties of X-ray events and show that it further reduces the background rate, particularly in the < 2 keV band. Over 0.3 to 12 keV, the background rate becomes 1.0 × 10 -3 s -1 keV -1 .

47 OTHER INSTRUMENTATION↗

Digital Twin Framework for PIP-II Linac: AI-Driven Multi-Scale Modeling from Ion Source to 800 MeV

The PIP-II linac will enable >1.2 MW beam power for DUNE, requiring unprecedented operational reliability across its warm front-end (RFQ, MEBT) and five distinct SRF sections operating at 162.5/325/650 MHz. We present a comprehensive digital twin framework uniquely combining a fully differentiable fast beam transport code with neural network surrogates trained on high-fidelity PIC simulations, capturing space charge and nonlinear dynamics beyond traditional envelope codes while achieving 10⁴ speedup at <1% accuracy. End-to-end differentiability enables gradient-based optimization across 500+ parameters simultaneously previously impossible with conventional tools while the model incorporates static/dynamic errors and serves as a virtual commissioning platform for diverse hardware integration. The framework facilitates reinforcement learning for pulsed/CW mode transitions, predictive maintenance through anomaly detection, and autonomous tuning algorithm development with real-time execution capability. Validation against physics simulations shows excellent agreement for the front-end, with initial results demonstrating potential for 30% commissioning time reduction and proactive fault mitigation, providing a scalable blueprint for operating next-generation high-intensity accelerators.

Pathak, Abhishek [Fermilab] (ORCID:000000021704208↗

Sensitivity of Simulations of Double-detonation Type Ia Supernovae to Integration Methodology

Abstract We study the coupling of hydrodynamics and reactions in simulations of the double-detonation model for Type Ia supernovae. When assessing the convergence of simulations, the focus is usually on spatial resolution; however, the method of coupling the physics together as well as the tolerances used in integrating a reaction network also play an important role. In this paper, we explore how the choices made in both coupling and integrating the reaction portion of a simulation (operator/Strang splitting versus the simplified spectral deferred corrections method we introduced previously) influences the accuracy, efficiency, and nucleosynthesis of simulations of double detonations. We find no need to limit reaction rates or reduce the simulation time step to the reaction timescale. The entire simulation methodology used here is GPU-accelerated and made freely available as part of the Castro simulation code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Antarctic sea ice concentrations inferred from Nimbus 5 ESMR and Landsat imagery

Landsat Multispectral Scanner near-infrared band images (MSS-7) are used to calculate Antarctic ice concentration values and these results are quantitatively compared to those derived from Nimbus 5 electrically scanning microwave radiometer (ESMR) brightness temperature data. The set of images includes regions near the ice edge and near the continental boundary during late winter or spring. The percentage of ice cover is derived from the microwave data, using an algorithm that incorporates the brightness temperature emissivity and interpolated climatological physical temperatures and that is not dependent on spatial resolution. Cloud-free Landsat images are used to provide an alternative determination of large-scale ice concentrations, but this method is limited in its ability to resolve individual ice flows, especially near the ice edge. A proportional classification procedure is applied to Landsat data based on observed reflection and the results are correlated with similarly-derived ESMR ice concentration values. The corresponding ice concentration values are found to agree to + or - 15%, with the largest uncertainties occurring in regions of low ice concentration.

Comiso, J. C.↗

Aerothermal modeling

The objectives, approach, and status of a program to develop the computational fluid dynamics tools needed to improve combustor design and analysis are outlined. The calculation procedure selected consists of a finite difference solution of the time averaged, steady state, primitive variable, elliptic form of the Reynolds equations. Standard TEACH type numerics are used to solve the resulting equations. These include hybrid differencing, SIMPLE algorithm for the pressure field, line by line iterative solution using the ADI method and the tridiagonal matrix algorithm (TDMA). Convergence is facilitated by using under relaxation. The physical processes are modeled by a two equation eddy viscosity model for turbulence; combustion is represented by a simple, irreversible, one step chemical reaction whose rate is influenced only by the time scale of the turbulence. The model evaluation procedure is also described.

Sturgess, G. J.↗

A general concurrent algorithm for plasma particle-in-cell simulation codes

The general concurrent particle-in-cell (GCPIC) algorithm has been used to implement an electrostatic particle-in-cell code on a 32-node hypercube parallel computer. The GCPIC algorithm decomposes the PIC code by dividing the particle simulation physical domain into subdomains that are equal in number to the number of processors; all subdomains will accordingly possess approximately equal numbers of particles. The portion of the code which updates particle positions and velocities is nearly 100 percent efficient when the number of particles increases linearly with that of hypercube processors.

Liewer, Paulett C.↗

Computational fluid dynamics research and applications at NASA Langley Research Center

Information on computational fluid dynamics (CFD) research and applications carried out at the NASA Langley Research Center is given in viewgraph form. The Langley CFD strategy, the five-year plan in CFD and flow physics, 3-block grid topology, the effect of a patching algorithm, F-18 surface flow, entropy and vorticity effects that improve accuracy of unsteady transonic small disturbance theory, and the effects of reduced frequency on first harmonic components of unsteady pressures due to airfoil pitching are among the topics covered.

South, Jerry C., Jr.↗

Dynamics and control of a tumbling satellite retrieval device

The problem of capturing a satellite in an uncontrolled three-axis tumble is addressed. The proposed solution is a truss which opens and engulfs the satellite, resulting in coincident centers of mass. Torques are then applied to the truss so that it matches the tumble rates of the satellite and remains motionless in the satellite's reference frame. Equations are derived which show how to vary the mass distribution to match any set of inertia moment ratios. The physical system is shown to be controllable. An optimal control algorithm for tumble matching is then derived and a computer simulation is developed.

Neff, J. M.↗

Adaptive Control Strategies for Flexible Robotic Arm

The control problem of a flexible robotic arm has been investigated. The control strategies that have been developed have a wide application in approaching the general control problem of flexible space structures. The following control strategies have been developed and evaluated: neural self-tuning control algorithm, neural-network-based fuzzy logic control algorithm, and adaptive pole assignment algorithm. All of the above algorithms have been tested through computer simulation. In addition, the hardware implementation of a computer control system that controls the tip position of a flexible arm clamped on a rigid hub mounted directly on the vertical shaft of a dc motor, has been developed. An adaptive pole assignment algorithm has been applied to suppress vibrations of the described physical model of flexible robotic arm and has been successfully tested using this testbed.

Bialasiewicz, Jan T.↗

Retrieval of Soil Moisture and Roughness from the Polarimetric Radar Response

The main objective of this investigation was the characterization of soil moisture using imaging radars. In order to accomplish this task, a number of intermediate steps had to be undertaken. In this proposal, the theoretical, numerical, and experimental aspects of electromagnetic scattering from natural surfaces was considered with emphasis on remote sensing of soil moisture. In the general case, the microwave backscatter from natural surfaces is mainly influenced by three major factors: (1) the roughness statistics of the soil surface, (2) soil moisture content, and (3) soil surface cover. First the scattering problem from bare-soil surfaces was considered and a hybrid model that relates the radar backscattering coefficient to soil moisture and surface roughness was developed. This model is based on extensive experimental measurements of the radar polarimetric backscatter response of bare soil surfaces at microwave frequencies over a wide range of moisture conditions and roughness scales in conjunction with existing theoretical surface scattering models in limiting cases (small perturbation, physical optics, and geometrical optics models). Also a simple inversion algorithm capable of providing accurate estimates of soil moisture content and surface rms height from single-frequency multi-polarization radar observations was developed. The accuracy of the model and its inversion algorithm is demonstrated using independent data sets. Next the hybrid model for bare-soil surfaces is made fully polarimetric by incorporating the parameters of the co- and cross-polarized phase difference into the model. Experimental data in conjunction with numerical simulations are used to relate the soil moisture content and surface roughness to the phase difference statistics. For this purpose, a novel numerical scattering simulation for inhomogeneous dielectric random surfaces was developed. Finally the scattering problem of short vegetation cover above a rough soil surface was considered. A general scattering model for grass-blades of arbitrary cross section was developed and incorporated in a first order random media model. The vegetation model and the bare-soil model are combined and the accuracy of the combined model is evaluated against experimental observations from a wheat field over the entire growing season. A complete set of ground-truth data and polarimetric backscatter data were collected. Also an inversion algorithm for estimating soil moisture and surface roughness from multi-polarized multi-frequency observations of vegetation-covered ground is developed.

Sarabandi, Kamal↗

High Altitude Measurements of Radiance at High Spectral and Spatial Resolution for SIMBIOS Sensor Calibration, Validation, and Intercomparisons

The successful combination of data from different ocean color sensors depends on the correct interpretation of signal from each of these sensors. Ideally, the sensor measured signals are calibrated to geophysical units of spectral radiance, and sensor artifacts are removed and corrected. The calibration process resamples the signal into a common radiometric data space so that subsequent ocean color algorithms that are applied to the data are based on physical processes and are inherently sensor independent. The objective of this project is to calibrate and validate the on-orbit radiometric characteristics of Sea-Viewing Wide Field-of-view Sensor (SeaWiFS) with underflights of NASA's calibrated Airborne Visible/Infrared Imaging Spectrometer (AVIRIS). This objective is feasible because AVIRIS measures the same spectral range as SeaWIFS at higher spectral resolution. In addition to satellite sensor underflights, the AVIRIS project has supported comparison and analysis of the radiometric calibration standards used for AVIRIS and SeaWIFS. To date, both the OCTS and SeaWIFS satellite sensors have been underflown by AVIRIS with matching spectral, spatial, geometric, radiometric, and temporal domains. The calibration and validation objective of this project is pursued for the following reasons: (1) Calibration is essential for the quantitative use of SeaWIFS and other SIMBIOS (Sensor Intercomparison and Merger for Biological and Interdisciplinary Oceanic Studies) sensor data; (2) Calibration in the laboratory of spaceborne sensors is challenging; (3) Satellite sensors are subjected aging on the ground and to trauma during launch; (4) The Earth orbit environment is significantly different than the laboratory calibration environment; (5) Through years of effort AVIRIS has been demonstrated to be well calibrated; and (6) AVIRIS can match the spectral and spatial observation characteristics near the top of the atmosphere at the time of SeaWIFS measurements.

Green, Robert O.↗

The Three-Dimensional Combustion of Heterogeneous Propellants

A numerical framework is described which permits the calculation of the 3-D combustion field supported by a heterogeneous propellant, allowing for complete coupling between the condensed phase physics, the gas-phase physics, and the unsteady, uneven, regressing surface. A random packing algorithm is used to construct models of ammonium-perchlorate in hydroxyl-terminated-polybutadiene propellants which mimic experimental propellants designed by R. Miller, and these are numerically burnt. Mean burning rates are compared with experimental data for four packs, over a pressure range of 7-200atm.

Massa, L.↗

Diabatic Initialization of Mesoscale Models in the Southeastern United States: Can 0 to 12h Warm Season QPF be Improved?

It is well known that numerical warm season quantitative precipitation forecasts lack significant skill for numerous reasons. Some are related to the model--it may lack physical processes required to realistically simulate convection or the numerical algorithms and dynamics employed may not be adequate. Others are related to initialization-mesoscale features play an important role in convective initialization and atmospheric observation systems are incapable of properly depicting the three-dimensional stability structure at the mesoscale. The purpose of this study is to determine if a mesoscale model initialized with a diabatic initialization scheme can improve short-term (0 to 12h) warm season quantitative precipitation forecasts in the Southeastern United States. The Local Analysis and Prediction System (LAPS) developed at the Forecast System Laboratory is used to diabatically initialize the Pennsylvania State University/National center for Atmospheric Research (PSUNCAR) Mesoscale Model version 5 (MM5). The SPORT Center runs LAPS operationally on an hourly cycle to produce analyses on a 15 km covering the eastern 2/3 of the United States. The 20 km National Centers for Environmental Prediction (NCEP) Rapid Update Cycle analyses are used for the background fields. Standard observational data are acquired from MADIS with GOES/CRAFT Nexrad data acquired from in-house feeds. The MM5 is configured on a 140 x 140 12 km grid centered on Huntsville Alabama. Preliminary results indicate that MM5 runs initialized with LAPS produce improved 6 and 12h QPF threat scores compared with those initialized with the NCEP RUC.

Lapenta, William M.↗