Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “joint inversion”

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 235 records · Page 13

Scientific and Technical Assistance for the Deployment of a Flexible Airborne Spectrometer System During C-MAPExp and COMEX: COMEX Version - 1.9T

The COMEX (CO2 and MEthane eXperiment) campaign supports the mission definition of CarbonSat and HyspIRI (Hyperspectral Infrared Imager) by providing representative airborne remote sensing data ‐ MAMAP (Methane Airborne MAPper) for CarbonSat; the Airborne Visual InfraRed Imaging Spectrometer (Classic & Next Generation) AVIRIS‐C/AVIRIS‐NG for HyspIRI ‐ as well as ground-based and airborne in‐situ data. The objectives of the COMEX campaign activities are (see Campaign Implementation Plan (RD‐4)): 1. Investigate spatial/spectral resolution trade‐offs for CH4 anomaly detection and flux inversion by comparison of MAMAP‐derived emission estimates with AVIRIS/AVIRIS‐NG derived data. 2. Evaluate sun‐glint observation geometry on CH4 retrievals for marine sources. 3. Characterize the effect of Surface Spectral Reflectance (SSR) heterogeneity on trace gas retrievals of CO2 and CH4 for medium and low‐resolution spectrometry. 4. Identify benefits from joint SWIR/TIR (ShortWave InfraRed/Thermal InfraRed ) data for trace gas detection and retrieval by comparison of MAMAP and AVIRIS/AVIRIS‐NG NIR/SWIR data with MAKO (Aerospace Corp.)TIR data. The ability to derive emission source strength for a range of strong emitting targets by remote sensing will be evaluated from combined AVIRIS‐NG and MAMAP data, adding significant value to the HyspIRI campaign AVIRIS‐NG dataset. The data will be used to quantify anomalies in atmospheric CO2 and CH4 from strong local greenhouse gas sources e.g. localized industrial complexes, landfills, etc. and to derive CO2 and CH4 emissions estimates from atmospheric gradient measurements. The original campaign concept was developed by University of Bremen and BRI. The COMEX campaign is funded bilaterally by NASA and ESA (European Space Agency). Whereas NASA funds the US part of the project via a contract with Dr. Ira Leifer, BRI (Bubbleology Research International), the contribution of MAMAP to the COMEX campaign is funded by ESA within the COMEX‐E project and NASA with respect to a 50 percent contribution to the flight-related costs of flying MAMAP on a US aircraft. The Data Acquisition Report (RD‐9) describes the instrumentation used, the measurements made by the team during the COMEX campaign in May/June 2014 and August/September 2014 in California, and an initial assessment of the data quality.

Point Sources↗

Computational neural learning formalisms for manipulator inverse kinematics

An efficient, adaptive neural learning paradigm for addressing the inverse kinematics of redundant manipulators is presented. The proposed methodology exploits the infinite local stability of terminal attractors - a new class of mathematical constructs which provide unique information processing capabilities to artificial neural systems. For robotic applications, synaptic elements of such networks can rapidly acquire the kinematic invariances embedded within the presented samples. Subsequently, joint-space configurations, required to follow arbitrary end-effector trajectories, can readily be computed. In a significant departure from prior neuromorphic learning algorithms, this methodology provides mechanisms for incorporating an in-training skew to handle kinematics and environmental constraints.

Gulati, Sandeep↗

Estimating Soil Moisture from Satellite Microwave Observations

Cooperative research in microwave remote sensing between the Hydrological Sciences Branch of the NASA Goddard Space Flight Center and the Earth Sciences Faculty of the Vrije Universiteit Amsterdam began with the Botswana Water and Energy Balance Experiment and has continued through a series of highly successful International Research Programs. The collaboration between these two research institutions has resulted in significant scientific achievements, most notably in the area of satellite-based microwave remote sensing of soil moisture. The Botswana Program was the first joint research initiative between these two institutions, and provided a unique data base which included historical data sets of Scanning Multifrequency Microwave Radiometer (SN4NM) data, climate information, and extensive soil moisture measurements over several large experimental sites in southeast Botswana. These data were the basis for the development of new approaches in physically-based inverse modelling of soil moisture from satellite microwave observations. Among the results from this study were quantitative estimates of vegetation transmission properties at microwave frequencies. A single polarization modelling approach which used horizontally polarized microwave observations combined with monthly composites of Normalized Difference Vegetation Index was developed, and yielded good results. After more precise field experimentation with a ground-based radiometer system, a dual-polarization approach was subsequently developed. This new approach realized significant improvements in soil moisture estimation by satellite. Results from the Botswana study were subsequently applied to a desertification monitoring study for the country of Spain within the framework of the European Community science research programs EFEDA and RESMEDES. A dual frequency approach with only microwave data was used for this application. The Microwave Polarization Difference Index (MPDI) was calculated from 37 GHz data and used to derive the one-way canopy transmissivity. Using a simple radiative transfer model, this information was combined with horizontally polarized 6.6 GHz SMMR observations to derive a 9-year time series of soil moisture for all of Spain at a one quarter degree spatial scale. Both day and night SMMR observations were used independently, in order to check the consistency of the results. A first order Fourier Transform was performed on the mean monthly soil moisture values to identify major characteristics of time series such as trend, amplitude, and phase shift.

Owe, M.↗

Combined Pre-supernova Alert System with KamLAND and Super-Kamiokande

Preceding a core-collapse supernova (CCSN), various processes produce an increasing amount of neutrinos of all flavors characterized by mounting energies from the interior of massive stars. Among them, the electron antineutrinos are potentially detectable by terrestrial neutrino experiments such as KamLAND and Super-Kamiokande (SK) via inverse beta decay interactions. Once these pre-supernova (pre-SN) neutrinos are observed, an early warning of the upcoming CCSN can be provided. In light of this, KamLAND and SK, both located in the Kamioka mine in Japan, have been monitoring pre-SN neutrinos since 2015 and 2021, respectively. Recently, we performed a joint study between KamLAND and SK on pre-SN neutrino detection. A pre-SN alert system combining the KamLAND detector and the SK detector was developed and put into operation, which can provide a supernova alert to the astrophysics community. Fully leveraging the complementary properties of these two detectors, the combined alert is expected to resolve a pre-SN neutrino signal from a 15 M ⊙ star within 510 pc of the Earth at a significance level corresponding to a false alarm rate of no more than 1 per century. For a Betelgeuse-like model with optimistic parameters, it can provide early warnings up to 12 hr in advance.

79 ASTRONOMY AND ASTROPHYSICS↗

Optimizing Retrieval Spaces of Bio-Optical Models for Remote Sensing of Ocean Color

We investigated the optimal number of independent parameters required to accurately represent spectral remote sensing reflectances (𝑅 rs ) by performing principal component analysis on quality controlled in situ and synthetic 𝑅 rs data. We found that retrieval algorithms should be able to retrieve no more than four free parameters from 𝑅 rs spectra for most ocean waters. In addition, we evaluated the performance of five different bio-optical models with different numbers of free parameters for the direct inversion of in-water inherent optical properties (IOPs) from in situ and synthetic 𝑅 rs data. The multi-parameter models showed similar performances regardless of the number of parameters. Considering the computational cost associated with larger parameter spaces, we recommend bio-optical models with three free parameters for the use of IOP or joint retrieval algorithms.

Ocean Color↗

Diagnostics: Chapter 8 of the special issue: on the path to tokamak burning plasma operation

This chapter presents the activity conducted by the ITPA topical group (TG) on Diagnostics over about the last 15 years. Following a general introduction of the ITER Diagnostics led by their measurement roles, the document is organized in several subchapters detailing the design support, research and development activity conducted by each of the specialist working groups (WGs) of the TG. Please note that the magnetic diagnostics were supported at the TG without a specific WG. Their status is included in the general introduction. In the following some highlights of the subchapter’s contents are provided. Recent advances in ITER first wall (FW) diagnostics for the measurements of plasma-metallic wall interaction in support of the ITER research plan are reported. An InfraRed imaging Video Bolometer for ITER has been developed and tested on several tokamaks to measure the radiated power loss. A laser-induced breakdown spectroscopy (LIBS) technique which utilizes a pulsed laser beam to ablate locally by forming a crater, will measure local tritium inventory in the FW material. Real-time Residual Gas Analyzers will measure the neutral gas composition in a divertor port and an equatorial port during plasma operation. Due to the full metallic FW environment, the plasma-wall interaction in ITER will face several challenges such as the compromised radiated power and divertor heat flux measurements by reflection. Ray tracing and analysis codes have been developed to eliminate and correct the effects of reflection in the measurements. The characteristics of the reflecting surfaces depending on the roughness and angle of the incidence have been measured by dedicated experiments, and the results were applied to the reflection elimination. For the measurement of the metallic impurity radiation induced by eroded metallic atoms, a vacuum ultraviolet spectrometer has been developed and tested. An extensive thermonuclear diagnostic suite will be required to support the operation of ITER and the planned experimental program for future burning plasma experiments. Due to the harsh environmental conditions, the implementation of diagnostic systems in ITER is a major challenge. These conditions include high levels of neutron and gamma fluxes, neutron heating, particle bombardment. Therefore, the selection and design of diagnostic systems must take into account a number of phenomena previously unseen in diagnostic design. For this reason, the measurement of neutrons and confined or lost fast ions, with particular emphasis on alpha particles, is critical to ITER. The diagnostics associated with these measurements will be important for future plasma-burning experiments at ITER. The high neutron emission and very large plasma size in ITER make neutron diagnostics the main diagnostic method used to measure plasma parameters such as fusion power, fusion power density, ion temperature, energy of fast ions and their spatial distributions in the plasma core. Active spectroscopy techniques are methods where a neutral particle beam is injected into the plasma and information on plasma parameters is extracted from the measurement of line emission resulting from the beam-plasma interaction, either by plasma ions or by beam atoms. Spatial localization is achieved by crossing the beamline and multiple observation lines. The ITER plasma will be a high temperature, moderately dense, fully ionized collisional plasma. The plasma facing surfaces are principally metallic being fashioned from beryllium or tungsten but many other elements, arising from either structural or from operational needs, may enter this plasma. The energy range of the emitted photons range from meV (infra-red) to multi keV (x-rays) and originate from all areas of the plasma volume. The primary role of passive emission diagnostics is to identify what is in the plasma from spectral signatures. Extracting quantitative information from these measurements such as impurity content, ion temperature, rotation, degree of detachment and radiated power depends on calibrated instruments, a physics model of the atomic and molecular processes and plasma transport and an analysis workflow that takes into account environmental effects such as reflections. The particular needs for ITER have prompted a multi-machine, many-year effort to address all these aspects and this chapter reviews the work on diagnostic design, experiments and new analysis techniques. An overview of the laser diagnostics to be implemented on ITER is also provided in this paper. This includes descriptions of the Thomson scattering in the core, edge and divertor regions, polarimetry and interferometry diagnostics used for measuring plasma density and also measurements of helium density in the divertor using Laser Induced Flourescence. Techniques which can allow improvements on current measurements are also addressed in particular expanding poloidal polarimetry measurements to measure field fluctuations and proposed use of dispersion interferometery which has a number of advantages over existing methods. This paper identifies particular areas where further research and testing on existing tokamaks is useful even at this advanced stage to inform the design of diagnostics for ITER. Outstanding areas of concern for the implementation of laser diagnostics, in particular with a view to reliable operation are identified. An overview of the latest developments of microwave diagnostic systems and techniques is given. The primary focus is the contributions for ITER—the next step burning plasma experiment—which is supplemented by describing recent progress of techniques applicable for fusion experiments beyond ITER. The contributions are intentionally kept concise, and are being supplemented by a rich list of references for further studies. Radiation induced effects are receiving continuous and well-deserved attention of the ITER diagnostic community and they are in many cases one of the primary design drivers of the ITER diagnostic systems. The paper summarizes recent progress in this area focusing primarily on the ITER diagnostics but in some cases provides also outlook for the possible solutions for even more demanding radiation environment of fusion reactors beyond ITER. Despite advancements in the area of modeling and simulation of various radiation induced effects, experimental testing in a nuclear environment as close as possible to the target one is still seen as unavoidable for proper qualification of particular diagnostic functional elements. Recent advancement within three diagnostic areas: optical diagnostics, magnetics and bolometers is covered. Encouraging results on qualification of silica glass vacuum window assemblies are presented. In the area of magnetic sensors, progress of irradiation tests performed on ITER in-vessel LTCC inductive sensors is presented with outlook for novel technological approaches to inductive sensors utilizing thick printing and photolithography technologies being highlighted. Summary of advancements in the area of steady state magnetic field sensors based on Hall effect is given. New results of neutron irradiation test of the ITER borosilicate glass inserts for vacuum electrical feedthroughs are summarized finding negligible swelling at target level of neutron fluence. Off-line irradiation tests of fiber optic current sensors for plasma current measurement demonstrated that both for gamma doses up to 5 MGy and a total neutron fluence up to 10 15 cm −2 , radiation induced changes are still compatible with required measurement accuracy on ITER. The ITER bolometers are given as an example how considering radiation effects may influence the diagnostic design. Finally, outlook for future main R&D directions is outlined. All optical and laser-based diagnostics in ITER will be using mirrors to guide plasma radiation toward detectors, cameras and sensors. In the hostile plasma, radiation and particle environment the optical characteristics of diagnostic mirrors will degrade directly affecting the entire performance of involved diagnostic systems. An assessment of factors affecting mirror performance is provided. Among the prime adverse factors are deposition of plasma impurities, sputtering of mirror surface and steam ingress in the vicinity of mirrors. Within the International Tokamak Physics Activity with active support by ITER central team and domestic agencies, the structured research and development (R&D) program on mitigation of risks for diagnostic mirrors is underway. Within this program the mirror material development, the passive mitigation of mirror degradation by using diagnostic ducts and shutters along with an active mirror recovery program comprising the in-situ mirror cleaning and calibration is underway. Recent developments in diagnostic mirror R&D are described in this Chapter along with an example of their implementation of R&D solutions in ITER Infrared Thermography diagnostic. An assessment of still open engineering and physics questions, considerations on mirror risks during an early phase of ITER operation are given along with an overview of diagnostic mirror evolution in the late ITER operation stage toward the demonstration fusion power plant. Several crucial areas of diagnostic R&D outlined in ITER Research Plan are addressed. The basic control groups in a fusion reactor can be broken-down in five categories: (1) plasma position, magnetic configuration, and plasma current control, (2) profile control and confinement optimization, (3) MHD control and suppression, (4) edge dissipation control, radiation and plasma exhaust control and (5) break-down optimization. These categories are coupled via the physics (a control action in one domain will affect the other domains) and via shared actuators (e.g. ECRH for impurity accumulation avoidance, current density distribution control and MHD suppression). Consequently, a supervisory control system should determine the priority of the various control tasks, their couplings, and the interfaces with the safety and interlock system. For the systematic development of the various controllers taking the complexity of the plasma and the control system into account, a model-based approach is required. A short historical overview is given of the developments in systems and control theory and control engineering with special emphasis on those developments that are most relevant for Nuclear Fusion research and operation. An overview is given of the state of the field of fusion plasma control for the control categories. It will be shown how synthetic diagnostics are being developed in ITER and how they are used in diagnostic design and design validation and how they can be in model-based controller synthesis using relatively simple models. In modern control methods, multiple diagnostics are used to constrain relatively simple models. The constrained models provide an estimate for the state. This opens the route to state controllers, such as model predictive control. A major challenge in nuclear fusion research is the coherent combination of data from heterogeneous diagnostics and modeling codes for machine control and safety as well as physics studies. Measured data from different diagnostics often provide information about the same subset of physical parameters. Additionally, information provided by some diagnostics might be needed for the analysis of other diagnostics. A joint analysis of complementary and redundant data allows, e.g. to improve the reliability of parameter estimation, to increase the spatial and temporal resolution of profiles, to obtain synergistic effects, to consider diagnostics interdependencies and to find and resolve data inconsistencies. Physics-based modeling and parameter relationships provide additional information improving the treatment of ill-posed inversion problems. A coherent combination of all kind of available information within a probabilistic framework allows for improved data analysis results. The concept of integrated data analysis (IDA) in the framework of Bayesian probability theory is outlined and contrasted with conventional data analysis. Components of the probabilistic approach are summarized and specific ingredients beneficial for data analysis at fusion devices are discussed.

ITER↗

An adaptive inverse kinematics algorithm for robot manipulators

An adaptive algorithm for solving the inverse kinematics problem for robot manipulators is presented. The algorithm is derived using model reference adaptive control (MRAC) theory and is computationally efficient for online applications. The scheme requires no a priori knowledge of the kinematics of the robot if Cartesian end-effector sensing is available, and it requires knowledge of only the forward kinematics if joint position sensing is used. Computer simulation results are given for the redundant seven-DOF robotics research arm, demonstrating that the proposed algorithm yields accurate joint angle trajectories for a given end-effector position/orientation trajectory.

Colbaugh, R.↗

Simplified Dynamic Control Of Redundant Manipulators

Simplified scheme proposed for dynamic control of robotic manipulator having redundant joints; that is, extra degrees of freedom beyond needed to perform task, to position and orient end effector at specified position and/or move it along specified trajectory. Extra degrees of freedom used to perform simultaneous subtask. Scheme adaptive and based on observed performance of manipulator. Involves neither complicated mathematical model of dynamics of manipulator nor time-consuming inverse kinematic transformation.

Seraji, Homayoun↗

Contact control for advanced applications of light weight arms

Many applications of robotic and teleoperated manipulator arms require operation in contact and non-contact regimes. This paper deals with both regimes and the transition between them with special attention given to problems of flexibility in the links and drives. This is referred to as contact control. Inverse dynamics is used to plan the tip motion of the flexible link so that the free motion can stop very near the contact surface without collision due to overshoot. Contact must occur at a very low speed since the high frequency impact forces are too sudden to be affected by any feedback generated torques applied to a joint at the other end of the link. The effect of approach velocity and surface properties are discussed. Force tracking is implemented by commands to the deflection states of the link and the contact force. This enables a natural transition between tip position and tip force control that is not possible when the arm is treated as rigid. The effect of feedback gain, force trajectory, and desired final force are of particular interest and are studied. Experimental results are presented on a one link arm and the system performance in the overall contact task is analyzed. Extension to multi-link cases with potential applications are discussed.

Book, Wayne J.↗

A Simple Stochastic Model for Generating Broken Cloud Optical Depth and Top Height Fields

A simple and fast algorithm for generating two correlated stochastic twodimensional (2D) cloud fields is described. The algorithm is illustrated with two broken cumulus cloud fields: cloud optical depth and cloud top height retrieved from Moderate Resolution Imaging Spectrometer (MODIS). Only two 2D fields are required as an input. The algorithm output is statistical realizations of these two fields with approximately the same correlation and joint distribution functions as the original ones. The major assumption of the algorithm is statistical isotropy of the fields. In contrast to fractals and the Fourier filtering methods frequently used for stochastic cloud modeling, the proposed method is based on spectral models of homogeneous random fields. For keeping the same probability density function as the (first) original field, the method of inverse distribution function is used. When the spatial distribution of the first field has been generated, a realization of the correlated second field is simulated using a conditional distribution matrix. This paper is served as a theoretical justification to the publicly available software that has been recently released by the authors and can be freely downloaded from http://i3rc.gsfc.nasa.gov/Public codes clouds.htm. Though 2D rather than full 3D, stochastic realizations of two correlated cloud fields that mimic statistics of given fields have proved to be very useful to study 3D radiative transfer features of broken cumulus clouds for better understanding of shortwave radiation and interpretation of the remote sensing retrievals.

Prigarin, Sergei M.↗

Joint state-parameter estimation for the reduced fracture model via the united filter

Here, in this paper, we introduce an effective United Filter method for jointly estimating the solution state and physical parameters in flow and transport problems within fractured porous media. Fluid flow and transport in fractured porous media are critical in subsurface hydrology, geophysics, and reservoir geomechanics. Reduced fracture models, which represent fractures as lower-dimensional interfaces, enable efficient multi-scale simulations. However, reduced fracture models also face accuracy challenges due to modeling errors and uncertainties in physical parameters such as permeability and fracture geometry. To address these challenges, we propose a United Filter method, which integrates the Ensemble Score Filter (EnSF) for state estimation with the Direct Filter for parameter estimation. EnSF, based on a score-based diffusion model framework, produces ensemble representations of the state distribution without deep learning. Meanwhile, the Direct Filter, a recursive Bayesian inference method, estimates parameters directly from state observations. The United Filter combines these methods iteratively: EnSF estimates are used to refine parameter values, which are then fed back to improve state estimation. Numerical experiments demonstrate that the United Filter method surpasses the state-of-the-art Augmented Ensemble Kalman Filter, delivering more accurate state and parameter estimation for reduced fracture models. This framework also provides a robust and efficient solution for PDE-constrained inverse problems with uncertainties and sparse observations.

Bayesian inference↗

Adaptive impedance control of redundant manipulators

A scheme for controlling the mechanical impedance of the end-effector of a kinematically redundant manipulator is presented. The proposed control system consists of two subsystems: an adaptive impedance controller which generates the Cartesian-space control input F (is a member of Rm) required to provide the desired end-effector impedance characteristics, and an algorithm that maps this control input to the joint torque T (is a member of Rn). The F to T map is constructed so that the robot redundancy is utilized to improve either the kinematic or dynamic performance of the robot. The impedance controller does not require knowledge of the complex robot dynamic model or parameter values for the robot, the payload, or the environment, and is implemented without calculation of the robot inverse kinematic transformation. As a result, the scheme is very general and is computationally efficient for on-line implementation.

Colbaugh, R.↗

Completion Effects on Diagnosing Multistage Facture Treatments with Distributed Temperature Sensing

Summary Distributed temperature sensing (DTS) is a valuable tool to diagnose multistage hydraulic fracture treatments. When a stage interval is shut in, the clusters that take more fluid during pumping warm up more slowly. Therefore, the fluid volume injected into each cluster can be quantitatively interpreted by numerical inversion of the warm-back temperature behavior. This general concept assumes that the different warm-back behavior is controlled by only the injected fluid volume; however, recent observations of DTS data indicate that completion configurations significantly influence the warm-back behavior. This paper investigates the completion effects on the DTS interpretation. In ideal conditions, when a stage is fractured, the upstream stage intervals should show an almost uniform temperature that is close to the injected fluid temperature. This is due to the high fluid velocity of injected fluid in the wellbore, and the upstream intervals have not been perforated (noncommunicating intervals), so the only heat transfer is heat conduction between the wellbore fluid and the surrounding reservoir. But the field DTS data show considerably irregular variations in temperature along the upstream stage intervals. These variations are caused by the completion effects. The nonuniform temperature profile is caused by different heat transfer behavior induced by completion hardware along the production casing string, such as joints, clamps, and blast protectors, and by the sensing cable location in the cement, as well as the cement quality. Because the varying heat transfer behavior impacts the warm-back behavior as well as the temperature profile, the completion effects need to be considered in DTS interpretation. A method of DTS interpretation considering the completion effects to diagnose multistage fracture treatments was developed. Because the heat transfer between a wellbore and a reservoir depends on the overall heat transfer coefficient describing heat conduction through the completion in a forward model, this parameter needs to be tuned all along the wellbore. To calibrate the completion effect, the temperature inversion is conducted using the temperature measured at a stage interval that is upstream of a stage interval currently being treated. Because the interpreted stage interval is not perforated at that time, the thermal behavior at the noncommunicating interval is governed by only the heat conduction through the completion environment. Once the effective values of the overall heat transfer coefficient are estimated along the interpreted stage interval, they can be assumed to be constant physical parameters. Then, the fluid volume distribution is interpreted by using the effective overall heat transfer coefficient profile along each interval. This study provides a field application of the developed interpretation method. The new interpretation method provides more accurate diagnosis of fracture treatments by DTS interpretation.

Engineering↗

Planck 2018 results XII. Galactic astrophysics using polarized dust emission

Observations of the submillimetre emission from Galactic dust, in both total intensity I and polarization, have received tremendous interest thanks to the Planck full-sky maps. In this paper we make use of such full-sky maps of dust polarized emission produced from the third public release of Planck data. As the basis for expanding on astrophysical studies of the polarized thermal emission from Galactic dust, we present full-sky maps of the dust polarization fraction p, polarization angle ψ, and dispersion function of polarization angles S . The joint distribution (one-point statistics) of p and NH confirms that the mean and maximum polarization fractions decrease with increasing NH. The uncertainty on the maximum observed polarization fraction, pmax = 22.0 -1.4 +3.5% at 353 GHz and 80' resolution, is dominated by the uncertainty on the Galactic emission zero level in total intensity, in particular towards diffuse lines of sight at high Galactic latitudes. Furthermore, the inverse behaviour between p and S found earlier is seen to be present at high latitudes. This follows the S ∝ p -1 relationship expected from models of the polarized sky (including numerical simulations of magnetohydrodynamical turbulence) that include effects from only the topology of the turbulent magnetic field, but otherwise have uniform alignment and dust properties. Thus, the statistical properties of p, ψ, and S for the most part reflect the structure of the Galactic magnetic field. Nevertheless, we search for potential signatures of varying grain alignment and dust properties. First, we analyse the product map S × p, looking for residual trends. While the polarization fraction p decreases by a factor of 3-4 between NH = 10 20 cm -2 and NH = 2 x 10 22 cm -2 , out of the Galactic plane, this product S × p only decreases by about 25%. Because S is independent of the grain alignment efficiency, this demonstrates that the systematic decrease in p with N H is determined mostly by the magnetic-field structure and not by a drop in grain alignment. This systematic trend is observed both in the diffuse interstellar medium (ISM) and in molecular clouds of the Gould Belt. Second, we look for a dependence of polarization properties on the dust temperature, as we would expect from the radiative alignment torque (RAT) theory. We find no systematic trend of S × p with the dust temperature T d , whether in the diffuse ISM or in the molecular clouds of the Gould Belt. In the diffuse ISM, lines of sight with high polarization fraction p and low polarization angle dispersion S tend, on the contrary, to have colder dust than lines of sight with low p and high S . We also compare the Planck thermal dust polarization with starlight polarization data in the visible at high Galactic latitudes. The agreement in polarization angles is remarkable, and is consistent with what we expect from the noise and the observed dispersion of polarization angles in the visible on the scale of the Planck beam. The two polarization emission-to-extinction ratios, R P/p and R S/V , which primarily characterize dust optical properties, have only a weak dependence on the column density, and converge towards the values previously determined for translucent lines of sight. We also determine an upper limit for the polarization fraction in extinction, p V /E(B - V), of 13% at high Galactic latitude, compatible with the polarization fraction p ≈ 20% observed at 353 GHz. Taken together, these results provide strong constraints for models of Galactic dust in diffuse gas.

79 ASTRONOMY AND ASTROPHYSICS↗

Yield Estimation of the August 2020 Beirut Explosion by Using Physics‐Based Propagation Simulations of Regional Infrasound

Abstract The August 2020 Beirut explosion is the largest single‐fired ammonium nitrate explosion documented in history. The massive explosion excited loud infrasound in the atmosphere, and clear waveforms were recorded by a regional infrasound array at an epicentral distance of 100 km, allowing for accurate measurements of explosion energy. We estimate the explosion size based on the infrasound waveform inversion. Unlike conventional inversions using empirical models, we perform full 3‐D finite‐difference simulations to obtain a physics‐based propagation model for the inversion. Accurate numerical modeling of infrasound is challenging as the propagation is substantially affected by the turbulent atmosphere. Instead of a single deterministic prediction, we provide a range of waveform predictions by running multiple simulations with stochastic weather forecast models, which allows for comprehensive uncertainty analysis of numerical modeling and estimated yields. Finally, we expand the yield estimation technique for seismoacoustic analysis and demonstrate the substantial advantage of the joint approach.

58 GEOSCIENCES↗

Scaled Vecchia Approximation for Fast Computer-Model Emulation

Many scientific phenomena are studied using computer experiments consisting of multiple runs of a computer model while varying the input settings. Gaussian processes (GPs) are a popular tool for the analysis of computer experiments, enabling interpolation between input settings, but direct GP inference is computationally infeasible for large datasets. We adapt and extend a powerful class of GP methods from spatial statistics to enable the scalable analysis and emulation of large computer experiments. Specifically, we apply Vecchia’s ordered conditional approximation in a transformed input space, with each input scaled according to how strongly it relates to the computer-model response. The scaling is learned from the data by estimating parameters in the GP covariance function using Fisher scoring. Our methods are highly scalable, enabling estimation, joint prediction, and simulation in near-linear time in the number of model runs. In several numerical examples, our approach substantially outperformed existing methods.

97 MATHEMATICS AND COMPUTING↗

Efficient Method for Analyzing Hyperspectral Remote Sensing Data

The current and next generation satellite infrared sensors such as AIRS, IASI, CrIS, and CLARREO all have one thing in common: large number of spectral channels. In order to fully utilize the information content of these sensors, a large number of Radiative Transfer (RT) calculations through the inhomogeneous atmosphere are needed. It is also computationally intensive to invert atmospheric and surface properties using all the channels. Usually, only subsets of channels are used to perform physical inversions for atmospheric profiles. We will explore ways to speed up RT calculations and inversions and test the algorithms using data from the Infrared Atmospheric Sounding Interferometer (IASI) satellite instrument. We will describe a retrieval algorithm using a Principal Component-based Radiative Transfer Model (PCRTM) for generating atmospheric temperature/moisture/ozone profiles and surface properties. The retrieval algorithm performs both the radiative transfer calculations and inversions in the Principal Component (PC) domain. The inversion algorithm is based on a non-linear Levenberg- Marquardt method with climatology covariance matrices and a priori information as constraints. One advantage of this approach is that it uses all information content from the ultraspectral data so that the retrieval is less sensitive to instrument noise and eliminates the need for selecting a sub-set of the channels. We will also use data collected during the Joint Airborne IASI Validation Experiment (JAIVEx) field campaign to validate the algorithm and IASI retrievals.

Liu, Xu↗

Tilt/Tip/Piston Manipulator with Base-Mounted Actuators

A proposed three-degree-of-freedom (tilt/tip/piston) manipulator, suitable for aligning an optical or mechanical component, would offer several advantages over prior such manipulators: Unlike in some other manipulators, no actuator would support the weight of another actuator: All of the actuators would be mounted on a base. Hence, there would be less manipulated weight. The basic geometry of the manipulator would afford mechanical advantage: that is, actuator motions would be larger than the motions they produce in the manipulated object. Mechanical advantage inherently increases the accuracy and resolution of manipulation. Unlike in some other manipulators, it would not be necessary to route power and/or data lines through manipulator joints. The proposed manipulator (see figure) would include three prismatic actuators (T1N1, T2N2, and T3N3) mounted on the base and operating in the same plane. Examples of suitable prismatic actuators include lead-screw mechanisms, linear hydraulic motors, piezoelectric linear drives, inchworm-movement linear stepping motors, and linear flexure drives. The actuators would control the lengths of links R1T1, R2T2, and R3T3. Three spherical joints (P1, P2, and P3) would be located at the corners of an equilateral triangle of side length q on the platform holding the object to be manipulated. Three inextensible limbs (R1P1, R2P2, and R3P3) having length r would connect the spherical joints on the platform to revolute joints (R1, R2, and R3) at the ends of the actuator-controlled links R1T1, R2T2, and R3T3. By varying the lengths of these links, one could control the tilt, tip, and piston coordinates of the platform. Closed-form equations for direct or forward kinematics of the manipulator (given the lengths of the variable links, find the tilt, tip, and piston coordinates) have been derived. The equations of inverse kinematics (find the variable link lengths needed to obtain the desired tilt, tip, and piston coordinates) have also been derived.

Tahmasebi, Farhad↗