Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “coordinate computing”

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 1,207 records · Page 67

WM26 Paper Multi-Robot Collaboration for Hazardous Environments

Hazardous nuclear and industrial facilities are rarely designed for robots. Work in these domains demand precise manipulation and robust mobility in cluttered, constrained spaces where off-the-shelf platforms struggle and “one-size-fits-all” machines become costly and complex. Idaho National Laboratory (INL) is developing an autonomous, multi-robot inspection system that coordinates task-specific platforms rather than relying on a single omni-tool robot. An electric truck serves as a power and compute hub for a custom manipulator co-developed with Florida International University (FIU), a commercial mini crawler, a pan–tilt–zoom camera, and a Nexxis Argus LiDAR mapping system. Working in concert, these robots generate spatial, radiation, and temperature maps of the pit environments at the Hanford Waste Tank Farms. These systems will capture visual records and environmental telemetry to allow for analysis post inspection. The system architecture uses Robot Operating System 2 (ROS 2) for publish/subscribe integration, NVIDIA Isaac Sim and Unity for simulation and visualization, and algorithms such as NVBlox to fuse data into unified 3D overlays. This robot-agnostic approach reduces operator burden by enabling autonomy across heterogeneous platforms and lets each robot be used where it is strongest. Having autonomous functions means operators don’t have to fully control multiple different components. The ease of use could allow for more widespread adoption of advanced robotics at waste management sites that see continued use. By coordinating simpler, purpose-built mechanisms, the approach lowers design and manufacturing complexity, reduces capital risk in contaminated settings, and improves controllability for complex inspection and manipulation tasks. We present the architecture, early results, and lessons learned from building and deploying this coordinated multi-robot system, with the goal of accelerating safe, cost-effective adoption of advanced robotics at waste-management sites.

42 - ENGINEERING↗

Empowering Scientific Discovery Through Computing at the Advanced Photon Source

This paper explores the challenges and solutions for managing and processing the vast amount of data generated by the Advanced Photon Source (APS), a synchrotron light source facility producing ultra-bright x-rays for diverse scientific domains. With 68 experimental beamlines covering materials research, biology, and more, the APS serves a wide user base across academia, government, and industry. The ongoing upgrade of the APS storage ring and installation of new instruments will amplify data generation and processing demands. This paper discusses the approach to address these demands through automated data processing using standardized workflows that produce faster scientific insights. The APS Data Management System coordinates various data related tasks to manage storage, data transfer, metadata cataloging, data processing, and interfaces with tools provided by Globus. Through integration with the Argonne Leadership Computing Facility (ALCF), APS users can efficiently access high-performance computing resources. Standardized workflows have led to reduced computational burdens on scientists and greater accessibility of high performance computing resources. We demonstrate how standardization and collaboration enable scientists to rapidly convert raw data into meaningful scientific results, establishing a streamlined path from data collection to analysis and ultimately to publication.

Parraga, Hannah↗

Coordinated observations of X-ray bright BL Lacertae objects

No new International Ultraviolet Explorer (IUE) observations were scheduled. The work on the spectral extraction routine was continued. The computer code was adapted from the RDAF PDP 11/44 computer to a VAX 11/750, in order that it run more quickly. A chapter on ultraviolet observations of blazars for a book celebrating the scientific accomplishment of IUE was written in colloboration. A copy of the final draft is attached.

Urry, C. M.↗

Algorithm for Wavefront Sensing Using an Extended Scene

A recently conceived algorithm for processing image data acquired by a Shack-Hartmann (SH) wavefront sensor is not subject to the restriction, previously applicable in SH wavefront sensing, that the image be formed from a distant star or other equivalent of a point light source. That is to say, the image could be of an extended scene. (One still has the option of using a point source.) The algorithm can be implemented in commercially available software on ordinary computers. The steps of the algorithm are the following: 1. Suppose that the image comprises M sub-images. Determine the x,y Cartesian coordinates of the centers of these sub-images and store them in a 2xM matrix. 2. Within each sub-image, choose an NxN-pixel cell centered at the coordinates determined in step 1. For the ith sub-image, let this cell be denoted as si(x,y). Let the cell of another subimage (preferably near the center of the whole extended-scene image) be designated a reference cell, denoted r(x,y). 3. Calculate the fast Fourier transforms of the sub-sub-images in the central NxN portions (where N < N and both are preferably powers of 2) of r(x,y) and si(x,y). 4. Multiply the two transforms to obtain a cross-correlation function Ci(u,v), in the Fourier domain. Then let the phase of Ci(u, v) constitute a phase function, phi(u,v). 5. Fit u and v slopes to phi (u,v) over a small u,v subdomain. 6. Compute the fast Fourier transform, Si(u,v) of the full NxN cell si(x,y). Multiply this transform by the u and phase slopes obtained in step 4. Then compute the inverse fast Fourier transform of the product. 7. Repeat steps 4 through 6 in an iteration loop, cumulating the u and slopes, until a maximum iteration number is reached or the change in image shift becomes smaller than a predetermined tolerance. 8. Repeat steps 4 through 7 for the cells of all other sub-images.

Sidick, Erkin↗

Latitude dependence of eddy variances

The eddy variance of a meteorological field must tend to zero at high latitudes due solely to the nature of spherical polar coordinates. The zonal averaging operator defines a length scale: the circumference of the latitude circle. When the circumference of the latitude circle is greater than the correlation length of the field, the eddy variance from transient eddies is the result of differences between statistically independent regions. When the circumference is less than the correlation length, the eddy variance is computed from points that are well correlated with each other, and so is reduced. The expansion of a field into zonal Fourier components is also influenced by the use of spherical coordinates. As is well known, a phenomenon of fixed wavelength will have different zonal wavenumbers at different latitudes. Simple analytical examples of these effects are presented along with an observational example from satellite ozone data. It is found that geometrical effects can be important even in middle latitudes.

Bowman, Kenneth P.↗

New Rover Conops with High-Performance Onboard Computing: Give Up Raw Data to Reduce Ops Cost and Do More Science

A major portion of time during the tactical operation of Mars rovers is spent for selecting, prioritizing, and coordinating sciences and engineering activities such that they fit within resource constraints, including the downlink data volume, energy, and time. In particular, the downlink data volume constraint is getting particularly tighter in recent missions because modern instruments produce increasingly high data volume while the communication bandwidth is essentially bounded by the law of physics. Tactical operation would be substantially simplified, hence the operation cost could be reduced, if the data volume constraint is relaxed or even removed. In this abstract, we propose a new operation paradigm for achieving this goal. The key observation is that, both in science and engineering applications, the bit size of raw data is typically much greater than the volume of processed information that is needed for scientific or engineering analysis. For example, a full-resolution image from Mastcam-Z, the main science camera on Perseverance, is about 700 kB in volume and we downlinked 29,685 images up to Sol 243, totaling ~20 GB of data. But of course, scientists do not use every pixel of these images; what they really look for in the images are geological features, typically represented by specific geometric configurations or textures. An end product after processing hundreds of Mascam-Z images could be a single geological map summarizing the spatial distribution of the features. For another example, a 100-meter drive of Perseverance produces 7-12 MB of drive telemetry, which records every detail of the rover's motion at 8 Hz, including position, attitude, steering angles, encoder readings, motor currents and many other information. But what the ground engineers eventually pay attention to is the signs of anomaly, such as excessive motor currents or high slip; if a drive is nominal, the vast majority of this data is unused. What if, then, we process the raw data onboard and only downlink the processed data that is relevant to scientific or engineering analyses, such as a list of detected science features (with cropped images) or a list of potential signs of anomaly while driving? A major roadblock for such onboard, high-level information processing has been the onboard computational resource. RAD750, the main onboard computer of Perseverance, is obviously not sufficient for performing complex image or signal processing such as object detection, semantic segmentation, or anomaly detection. Interestingly, RAD750 is not the best processor that Perseverance has; Qualcomm's Snapdragon 801, a modern mobile processor, is on her Heli Base Station, a device for communicating with Mars Helicopter Ingenuity; also, Intel's Atom E3845 processors are on engineering cameras. In the reminder of this paper, we will introduce two particular uses cases of these high-performance co-processors (meaning auxiliary CPU, GPU, or other types of processors that are separate from the main processor that runs the main flight software) for lowering operation cost and accommodating more science activities for a given communication constraint.

Didier, A.↗

Three-dimensional unsteady viscous flow analysis over airfoil sections

A three-dimensional solution procedure for the approximate form of the Navier-Stokes equation was exercised in the two- and three-dimensional modes to compute the unsteady turbulent boundary layer on a flat plate corresponding to the data of Karlsson. The procedure is based on the use of a consistently split Linearized Block Implicit technique in conjunction with a QR operator scheme. New time-dependent upstream boundary conditions were developed that yielded realistic solutions for the interior in the vicinity of the upstream boundary. Comparisons of the computation employing these boundary conditions with the data indicate that both qualitative and quantitative agreement was obtained for the mean velocity and the in phase and out of phase components of the first harmonic of the velocity. In addition, the calculation gave results for the skin friction phase angle that had expected physical behavior for large distances downstream of the inflow boundary. For the three-dimensional case, the two-dimensional data of Karlsson was considered, but in a coordinate system skewed at 45 deg to the free stream direction. The results of the calculations were in excellent agreement with the data and the two-dimensional computations.

Weinberg, B. C.↗

Computational study of duct and pipe flows using the method of pseudocompressibility

A viscous, three-dimensional, incompressible, Navier-Stokes Computational Fluid Dynamics code employing pseudocompressibility is used for the prediction of laminar primary and secondary flows in two 90-degree bends of constant cross section. Under study are a square cross section duct bend with 2.3 radius ratio and a round cross section pipe bend with 2.8 radius ratio. Sensitivity of predicted primary and secondary flow to inlet boundary conditions, grid resolution, and code convergence is investigated. Contour and velocity versus spanwise coordinate plots comparing prediction to experimental data flow components are shown at several streamwise stations before, within, and after the duct and pipe bends. Discussion includes secondary flow physics, computational method, computational requirements, grid dependence, and convergence rates.

Williams, Robert W.↗

Static and dynamic properties of atomic nuclei with high-resolution potentials

Here, we compute ground-state and dynamical properties of 4 He and 16 O nuclei using as input high-resolution, phenomenological nucleon-nucleon and three-nucleon forces that are local in coordinate space. The nuclear Schrodinger equation for both nuclei is accurately solved employing the auxiliary-field diffusion Monte Carlo approach. For the 4 He nucleus, detailed benchmarks are carried out with the hyperspherical harmonics method. In addition to presenting results for the binding energies and radii, we also analyze the momentum distributions of these nuclei and their Euclidean response function corresponding to the isoscalar density transition. The latter quantity is particularly relevant for lepton-nucleus scattering experiments, as it paves the way to quantum Monte Carlo calculations of electroweak response functions of 16 O.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Development of a three-dimensional supersonic inlet flow analysis

A method for computing three dimensional flow in supersonic inlets is described. An approximate set of governing equations is given for viscous flows which have a primary flow direction. The governing equations are written in general orthogonal coordinates. These equations are modified in the subsonic region of the flow to prevent the phenomenon of branching. Results are presented for the two sample cases: a Mach number equals 2.5 flow in a square duct, and a Mach number equals 3.0 flow in a research jet engine inlet. In the latter case the computed results are compared with the experimental data. A users' manual is included.

Buggeln, R. C.↗

Quantum Distributed Unit Commitment: An Application in Microgrids

The dawn of quantum computing brings on a revolution in the way combinatorially complex power system problems such as Unit Commitment are solved. The Unit Commitment problem complexity is expected to increase in the future because of the trend toward the increase of penetration of intermittent renewables. Even though quantum computing has proven effective for solving a host of problems, its applications for power systems’ problems have been rather limited. Here, in this paper, a quantum unit commitment is innovatively formulated and the quantum version of the decomposition and coordination alternate direction method of multipliers (ADMM) is established. The above is achieved by devising quantum algorithms and by exploiting the superposition and entanglement of quantum bits (qubits) for solving subproblems, which are then coordinated through ADMM to obtain feasible solutions. The main contributions of this paper include: 1) the innovative development of a quantum model for Unit Commitment; 2) development of decomposition and coordination-supported framework which paves the way for the utilization of limited quantum resources to potentially solve the large-scale discrete optimization problems; 3) devising the novel quantum distributed unit commitment (QDUC) to solve the problem in a larger scale than currently available quantum computers are capable of solving. The QDUC results are compared with those from its classical counterpart, which validate the efficacy of quantum computing.

97 MATHEMATICS AND COMPUTING↗

Analysis and computation of three-dimensional flow in strongly curved ducts

Numerical solutions are presented for three-dimensional laminar and turbulent flow in curved ducts of rectangular cross section and significant curvature. The analysis is based on a primary-secondary velocity decomposition in a given coordinate system, and leads to approximate governing equations which correct an a priori inviscid solution for viscous effects, secondary flows, total pressure distortion, heat transfer, and internal flow blockage and losses. Solution of the correction equations is accomplished as an initial-value problem in space using an implicit forward-marching technique. The overall solution procedure requires significantly less computational effort than Navier-Stokes algorithms. The present solution procedure is effective even with the extreme local mesh resolution which is necessary to resolve near-wall sublayer regions in turbulent flow calculations. Computed solutions for both laminar and turbulent flow compare very favorably with available analytical and experimental results.

Kreskovsky, J. P.↗

Dimetalloylene (M‐E‐M) Complexes of Heavier Main Group Elements Ge, Sn, Pb, Bi via Cleavage of E‐X Bonds (X=N(SiMe 3 ) 2 , O t Bu) with an Iridium Hydride

Abstract Reactions of the Ir V hydride [ Me BDI Dipp ]IrH 4 {BDI=(Dipp)NC(Me)CH(Me)CN(Dipp); Dipp=2,6‐ i Pr 2 C 6 H 3 } with E[N(SiMe 3 ) 2 ] 2 (E=Sn, Pb) afforded the unusual dimeric dimetallotetrylenes ([ Me BDI Dipp ]IrH) 2 ( μ 2 ‐E) 2 in good yields. Moreover, ([ Me BDI Dipp ]IrH) 2 ( μ 2 ‐Ge) 2 was formed in situ from thermal decomposition of [ Me BDI Dipp ]Ir(H) 2 Ge[N(SiMe 3 ) 2 ] 2 . These reactions are accompanied by liberation of HN(SiMe 3 ) 2 and H 2 through the apparent cleavage of an E−N(SiMe 3 ) 2 bond by Ir−H. In a reversal of this process, ([ Me BDI Dipp ]IrH) 2 ( μ 2 ‐E) 2 reacted with excess H 2 to regenerate [ Me BDI Dipp ]IrH 4 . Varying the concentrations of reactants led to formation of the trimeric ([ Me BDI Dipp ]IrH 2 ) 3 ( μ 2 ‐E) 3 . The further scope of this synthetic route was investigated with group 15 amides, and ([ Me BDI Dipp ]IrH) 2 ( μ 2 ‐Bi) 2 was prepared by the reaction of [ Me BDI Dipp ]IrH 4 with Bi(NMe 2 ) 3 or Bi(O t Bu) 3 to afford the first example of a “naked” two‐coordinate Bi atom bound exclusively to transition metals. A viable mechanism that accounts for the formation of these products is proposed. Computational investigations of the Ir 2 E 2 (E=Sn, Pb) compounds characterized them as open‐shell singlets with confined nonbonding lone pairs at the E centers. In contrast, Ir 2 Bi 2 is characterized as having a closed‐shell singlet ground state.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Optimizing power system restoration with damaged communications

Utility procedures for power system blackstart and restoration typically assume that energization decisions can be reliably communicated across the grid. In reality, the communications and control network would likely also be affected in power outages, such as those caused by extreme weather events or cyber-attacks. This paper studies the effect of damage to the power system communications and control infrastructure on restoration operations following a blackout. We model the communications infrastructure as a graph, overlaying the power grid, and imposing the requirement that every energized element in the power grid be observable from a control center. We expand on a specialized branch-and-bound algorithm from the literature to optimize the restoration process and devise an initialization heuristic and a rounding heuristic to improve solution speed. We perform numerical experiments on synthetic systems for Illinois and Texas with outages based on a solar flare or hurricane. We compare the results of our specialized branch-and-bound algorithm to the results from (i) the initialization heuristic alone, (ii) a variation of this heuristic that we use as a baseline, and (iii) the restoration optimization for the power system without communications constraints. Here, we find that damage to the communications infrastructure significantly increases the time required to re-energize the grid. Moreover, by simultaneously optimizing communications repairs and grid energization decisions, we are able to re-energize the grid significantly faster than if communications repairs and energization decisions were made independently or with partial coordination, motivating improvements to current industry practice.

97 MATHEMATICS AND COMPUTING↗

Cyclic Peptides for Lanthanide Binding

Lanthanide ions are difficult to separate from one another due to their similar chemical properties. The discovery of lanthanide-binding peptides and proteins in nature has led to an increased interest in the possibility of utilizing the strong binding of peptides to lanthanide ions for their separations; as such, there has been an effort to identify or design peptides with improved lanthanide binding and selectivity toward particular lanthanide ions. Here, in this study, we designed and characterized lanthanide-binding cyclic peptides (LBCPs) with molecular dynamics simulations, electronic structure calculations, and emission spectroscopy. Luminescent decay measurements were done to determine the number of water molecules coordinated to the Eu 3+ ion in Eu-LBCP complexes and compare to the predicted number of water molecules by computation to assess the lanthanide-binding affinity of LBCPs. Measured stability constants show binding of the LBCPs to the Eu 3+ ion with stronger than micromolar affinity. We were able to identify multiple peptides that selectively bind to middle lanthanides. We describe the structural basis of the lanthanide-binding selectivity trend with strongest binding to the middle lanthanides, followed by the heavier lanthanides, and finally to the lighter ions.

ions↗

Distributed Rules-Based Deconfliction of ADMS Applications: Part 1: Requirements & Decomposition

This work introduces a rules-based deconfliction methodology for resolving conflicting device control commands issued by advanced power applications considering a range of technical, economic, environmental, and social objectives. The methodology is designed to serve as one of multiple alternative implementations (along with application cooperation and global optimization) for the numerical component of the Deconfliction Pipeline. Development of the methodology is divided into two parts. This first document introduces the requirements, context, and methods for decomposing the deconfliction problem using the Laminar Coordination Framework and Variable Grid Structures. The deconfliction problem is decomposed into a distributed optimization problem based on the concept of quasi-static grid segments, which form independent distributed areas for control and coordination. Selection of the optimal number of decompositions of the deconfliction problem should be made based on a tradeoff analysis between computational speed and global optimality. This second document will define an initial set of technical, economic, and environmental criteria, as well as thirty specific qualitative rules that are used as part of the deconfliction methodology to eliminate non-viable setpoint alternatives. The deconfliction optimization problem is converted into a ranking of individual discrete setpoints, which are scored by the extent to which they satisfy specific decision criteria. The ranking is determined through the concepts of deconfliction exclusivity, priority, and preference. Several multi-criteria decision-making frameworks are examined with the simple multi-attribute rating technique exploiting ranks (SMARTER) recommended as a simple implementation alternative that aligns with the steps of the rules-based deconfliction methodology.

97 MATHEMATICS AND COMPUTING↗

New Methods for Predicting Non-Born-Oppenheimer Chemistry

Current methods for modeling non-adiabatic molecular dynamics face fundamental limitations when treating geometric phase effects: quantum mechanical phenomena where nuclear wavepackets acquire phase shifts when encircling conical intersections. Existing approaches either neglect these effects entirely or rely on potential energy surfaces arising from the Born-Oppenheimer approximation, which introduce artificial singularities and can overestimate geometric phase contributions. This project developed a new theoretical framework based on exact factorization (XF) methods to overcome these limitations. We derived mathematical formulations for hybrid quantum-classical XF dynamics that selectively treat critical nuclear degrees of freedom quantum mechanically while propagating others classically. This approach addresses the computational intractability that has previously limited exact methods to toy systems. Key innovations include a new approach to systematically identifying nuclear coordinates requiring quantum treatment, as well as novel implementation strategies that interface with existing quantum chemistry codes. The project also developed a proof-of-concept code for treating Jahn-Teller systems and creation of educational materials on non-adiabatic dynamics geared at the graduate level. The theoretical framework developed will enable future systematically improvable calculations of nuclear quantum effects in realistic molecular systems, filling a critical gap in non-adiabatic dynamics methods. This foundation supports future development of predictive tools for designing energy-relevant photochemical processes where quantum coherence effects may be exploited to control reaction outcomes.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Asymmetric hypersonic flow

A general method for the analysis of the inviscid asymmetric hypersonic flow fields enveloping smooth bodies of general shape is given. The method is based on the assumption of a thin shock layer which yields an expression for pressure in generalized Mises coordinates. Numerical results for elliptic cones at angle of attack are shown to compare well with experiments and other theories. The computing logic for a blunt body is described, and a limiting solution at the stagnation point is presented.

Maslen, S. H.↗