SEARCH · Search NASA

Results for “Discrete Rotation”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 541 records · Page 30

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Algorithm to extract direction in 2D discrete distributions and a continuous Frobenius norm

In this study, we present a novel algorithm for determining directionality in 2D distributions of discrete data. We compare a reference dataset with a known direction to a measured dataset with an unknown direction by the Frobenius norm of the difference (FND) to find the unknown direction. To generalize this concept, we develop a continuous Frobenius norm of the difference (CFND) as a continuous analog of the FND and derive its analytical expression. By relating fitted and normalized 2D Gaussian distributions, we show that the CFND approximates the FND, and we validate this relationship with computer simulations. We find that a first-order approximation of the CFND between two similar Gaussian distributions takes the form of an absolute sine function, offering a simple analytical form with potential applications in specialized areas such as segmented inverse beta decay neutrino detectors, astronomy, machine learning, and more. Our methodology consists of modeling a 2D Gaussian distribution, binning the data into a histogram, and encoding it as a square matrix. Rotating this matrix around its geometric center and comparing it to a measured dataset using the FND gives us rotational data that we fit with an absolute sine function. The location of the minimum of this fit is the angle closest to the true angle of the direction in the measured dataset. We present the derivation and discuss initial applications of the CFND in our novel algorithm, demonstrating its success in approximating directionality in 2D distributions.

Physics↗

DOE OSTI2023

Solute effects upon dislocation motion and recovery in Mg alloys (Final Report)

The objective the research was to develop a firmer understanding of the interactions between substitutional solute atoms and dislocations within Mg alloys. These interactions govern the absolute and relative mobilities of various types of dislocations in Mg (e.g., basal < a>, non-basal < a>, and pyramidal < c+a>). Furthermore, they have an impact on dislocation recovery processes (e.g., cross-glide, climb, rearrangement, and annihilation). Ultimately, solute-dislocation interactions strongly impact a) strain hardening, b) strain rate sensitivity, c) plastic anisotropy of textured polycrystals, d) texture evolution, e) dislocation substructures evolution, and even f) recrystallization behavior. It is broadly known that Mg alloys exhibit poor low temperature formability. Since formability is largely governed by properties (a) – (c) in this list, it is critical to develop a better understanding of solute-dislocation interactions, if one hopes to improve the situation. While the theory of static solute strengthening is well developed, especially for alloys with face centered cubic (FCC) crystal structures, there are outstanding questions related to applications to hexagonal close packed (HCP) crystal structures and dynamic strain aging (DSA) in materials of various structures. DSA has far-reaching implications for metal formability and, in the case of Mg alloys, appears to correlate with the so-called rare-earth (RE) texture which has been shown to benefit formability. The Portevin-Le Chatelier (PLC) effect, and associated negative strain rate sensitivity, occur at higher temperatures (>100℃) in Mg alloys as compared with similar Al alloys, even though they have similar melting points and solute diffusivities. Our preliminary research has shown that modern, physics-based models of DSA can be tuned to describe the behavior of Mg alloys if a rather higher activation enthalpy is assumed for cross-core diffusivity. While this partially explains the delay in DSA to higher temperatures, it is also hypothesized that this delay is due in part to the intrinsically more thermally activated (rate sensitive) nature of non-basal < a> dislocation motion which is required for macroscopic flow of Mg alloys, whereas octahedral slip in many FCC metals like aluminum is essentially athermal at room temperature. It was originally proposed to employ a combination of in-situ diffraction-based experimental characterization to validate existing theory. It was envisioned to perform in-situ transmission electron microscopy (TEM) to assess individual dislocation behavior and in-situ high-energy X-ray diffraction (HEXRD) techniques which were showing great promise for elucidating collective dislocation behavior, 2 including recovery, especially if the contributions to various forms of diffraction peak broadening (𝜂𝜂,𝜔𝜔,and 2𝜃𝜃) can be effectively integrated. Finally, it was envisioned to perform discrete dislocation dynamics (DDD) modeling approaches to aide in the interpretation of both TEM and HEXRD experiments. In the end, mechanical tests were performed on more complex Mg alloys which exhibited evidence of dynamic strain aging, and this led to the establishment of another project. Mechanical test data obtained at McMaster University served as the basis of an assessment of the applicability of Bazinski’s “stress equivalence” theory of solute strengthening to polycrystalline alloys of Mg. Although we did not succeed in applying the approach to Mg alloys, we did develop expertise with the HEXRD approach using a BCC, β-Ti alloy and demonstrated numerous new capabilities that may be applied to any polycrystalline material in collaboration with researchers at CHESS and around the world: (1) assessment of details of the elastoplastic transition (yielding) using a combination of HEXRD and full-field polycrystal plasticity modeling, (2) the first-ever experimental observation of strong stress rotation within the individual grains of a polycrystalline material, and (3) a comprehensive analysis of the grain-level dislocation density evolution based upon diffraction peak broadening along 𝜂𝜂,𝜔𝜔,and 2𝜃𝜃 directions. This final aspect allowed us to confirm that dislocations were gliding on multiple plane types and not restricted to {110} type planes, and it also provided clues as to why some grains were unloading during straining, with surprising implications for our understanding of the effects of geometrically necessary dislocations (GNDs). Finally, graduate student, Mohammed Shabana, developed a MATLAB code which confirmed the conclusions of Prof. Catalin Picu (Rensallear Polytechnic Institute, RPI) regarding the effect of solute-trapped, forest dislocations on the breaking stress of Lomer lock junctions in FCC metal alloys. He applied the same anisotropic line-tension model to a variety of dislocation junction configurations and found an inconsistency in the widely cited results of Dupuy and Fivel regarding the Hirth Lock, and he outlined an approach to extend these finding to HCP Mg alloys that we are still pursuing with discretionary fundings at UVA.

36 MATERIALS SCIENCE↗

DOE OSTI2016

Hydrogen segregation to inclined Σ3 < 110 > twin grain boundaries in nickel

Low-mobility twin grain boundaries dominate the microstructure of grain boundary-engineered materials and are critical to understanding their plastic deformation behaviour. The presence of solutes, such as hydrogen, has a profound effect on the thermodynamic stability of the grain boundaries. This work examines the case of a Σ3 grain boundary at inclinations from 0° ≤ Φ ≤ 90°. The angle Φ corresponds to the rotation of the Σ3 (1 1 1) < 1 1 0 > (coherent) into the Σ3 (1 1 2) < 1 1 0 > (lateral) twin boundary. To this end, atomistic models of inclined grain boundaries, utilising empirical potentials, are used to elucidate the finite-temperature boundary structure while grand canonical Monte Carlo models are applied to determine the degree of hydrogen segregation. In order to understand the boundary structure and segregation behaviour of hydrogen, the structural unit description of inclined twin grain boundaries is found to provide insight into explaining the observed variation of excess enthalpy and excess hydrogen concentration on inclination angle, but the explanatory power is limited by how the enthalpy of segregation is affected by hydrogen concentration. At higher concentrations, the grain boundaries undergo a defaceting transition. In order to develop a more complete mesoscale model of the interfacial behaviour, an analytical model of boundary energy and hydrogen segregation that relies on modelling the boundary as arrays of discrete 1/3 < 1 1 1 > disconnections is constructed. Lastly, the complex interaction of boundary reconstruction and concentration-dependent segregation behaviour exhibited by inclined twin grain boundaries limits the range of applicability of such an analytical model and illustrates the fundamental limitations for a structural unit model description of segregation in lower stacking fault energy materials.

36 MATERIALS SCIENCE↗

DOE OSTI2015

Electron-impact vibrational excitation of cyclopropane

We report a very detailed test of the ab initio discrete momentum representation (DMR) method of calculating vibrational excitation of polyatomic molecules by electron impact, by comparison of its results with an extensive set of experimental data, covering the entire range of scattering angles from 10{sup ∘} to 180{sup ∘} and electron energies from 0.4 to 20 eV. The DMR calculations were carried out by solving the two-channel Lippmann-Schwinger equation in the momentum space, and the interaction between the scattered electron and the target molecule was described by exact static-exchange potential corrected by a density functional theory (DFT) correlation-polarization interaction that models target’s response to the field of incoming electron. The theory is found to quantitatively reproduce the measured spectra for all normal modes, even at the difficult conditions of extreme angles and at low energies, and thus provides full understanding of the excitation mechanism. It is shown that the overlap of individual vibrational bands caused by limited experimental resolution and rotational excitation must be properly taken into account for correct comparison of experiment and theory. By doing so, an apparent discrepancy between published experimental data could be reconciled. A substantial cross section is found for excitation of the non-symmetric HCH twisting mode ν{sub 4} of A{sub 1}{sup ″} symmetry by the 5.5 eV A{sub 2}{sup ′} resonance, surprisingly because the currently accepted selection rules predict this process to be forbidden. The DMR theory shows that the excitation is caused by an incoming electron in an f-wave of A{sub 2}{sup ′} symmetry which causes excitation of the non-symmetric HCH twisting mode ν{sub 4} of the A{sub 1}{sup ″} symmetry and departs in p- and f-waves of A{sub 2}{sup ″} symmetry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1994

The reflection for dense plant canopies from the one-angle radiative transfer equation

An essential component of remote sensing of vegetation canopies from satellites is fundamental understanding. Since passive remote is driven by photons, the modeling of photon interactions with vegetation is a basic building block in that understanding. Several such photon transport models have been developed during the past two decades and continue to be developed. Different approaches have been followed including monte carlo, radiosity methods, geometric shadowing, and radiative transfer. In general, each approach has application for canopies with specific attributes. This presentation concerns the application of radiative transfer to dense vegetation canopies in which the soil does not participate. The approach taken here is novel in that a consistent theory for photon transport for non-rotationally invariant leaf scattering is developed in a canopy with a general leaf angle distribution (LAD). The theory is limited to the one-angle approximation (azimuthally averaged radiance) and is based on Chandrasekhar's analytical theory. While such a model is admittedly only approximate, it does fulfill a unique function in our search for understanding. In particular, the model is simple in its construct yet contains the essential features of canopy architecture that are mainly responsible for observed responses. Thus, this model will not only be a predictive tool but also an educational one. The mathematical setting is the radiative transfer equation in a dense (semiinfinite) canopy. The leaf scattering phase function is assumed to be Lambertian with different reflectance and transmittance. In addition, abaxial and adaxial differentiation is allowed which effectively destroys optical reciprocity. The analytical solution for the canopy BRDF is obtained by manipulation of the integral transport equation (a la Chandrasekhar) for a general LAD. With discretization of the. leaf angle, the resulting set of integral equations are solved iteratively including an acceleration procedure when the single scatter albedo is near one (in the NIR). Results will be compared to the LARS soybean canopy radiances as well as to broadleaf results from a recent Ames experiment.

Ganapol, B. D.↗

DOE OSTI2017

Anatomy of an Asteroid Breakup: The Case of P/2013 R3

We present an analysis of new and published data on P/2013 R3, the first asteroid detected while disintegrating. Thirteen discrete components are measured in the interval between UT 2013 October 01 and 2014 February 13. We determine a mean, pair-wise velocity dispersion among these components of Δv = 0.33 ± 0.03 m s{sup −1} and find that their separation times are staggered over an interval of ∼5 months. Dust enveloping the system has, in the first observations, a cross-section of ∼30 km{sup 2} but fades monotonically at a rate consistent with the action of radiation pressure sweeping. The individual components exhibit comet-like morphologies and also fade except where secondary fragmentation is accompanied by the release of additional dust. We find only upper limits to the radii of any embedded solid nuclei, typically ∼100–200 m (geometric albedo 0.05 assumed). Combined, the components of P/2013 R3 would form a single spherical body with a radius of ≲400 m, which is our best estimate of the size of the precursor object. The observations are consistent with rotational disruption of a weak (cohesive strength of ∼50 to 100 N m{sup −2}) parent body, ∼400 m in radius. Estimated radiation (YORP) spin-up times of this parent are ≲1 Myr, shorter than the collisional lifetime. If present, water ice sublimating at as little as 10{sup −3} kg s{sup −1} could generate a torque on the parent body rivaling the YORP torque. Under conservative assumptions about the frequency of similar disruptions, the inferred asteroid debris production rate is ≳10{sup 3} kg s{sup −1}, which is at least 4% of the rate needed to maintain the Zodiacal Cloud.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2011

Convergence of quantum random walks with decoherence

In this paper, we study the discrete-time quantum random walks on a line subject to decoherence. The convergence of the rescaled position probability distribution p(x,t) depends mainly on the spectrum of the superoperator L{sub kk}. We show that if 1 is an eigenvalue of the superoperator with multiplicity one and there is no other eigenvalue whose modulus equals 1, then P(({nu}/{radical}(t)),t) converges to a convex combination of normal distributions. In terms of position space, the rescaled probability mass function p{sub t}(x,t){identical_to}p({radical}(t)x,t), x is an element of Z/{radical}(t), converges in distribution to a continuous convex combination of normal distributions. We give a necessary and sufficient condition for a U(2) decoherent quantum walk that satisfies the eigenvalue conditions. We also give a complete description of the behavior of quantum walks whose eigenvalues do not satisfy these assumptions. Specific examples such as the Hadamard walk and walks under real and complex rotations are illustrated. For the O(2) quantum random walks, an explicit formula is provided for the scaling limit of p(x,t) and their moments. We also obtain exact critical exponents for their moments at the critical point and show universality classes with respect to these critical exponents.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1998

Stratified Diffractive Optic Approach for Creating High Efficiency Gratings

Gratings with high efficiency in a single diffracted order can be realized with both volume holographic and diffractive optical elements. However, each method has limitations that restrict the applications in which they can be used. For example, high efficiency volume holographic gratings require an appropriate combination of thickness and permittivity modulation throughout the bulk of the material. Possible combinations of those two characteristics are limited by properties of currently available materials, thus restricting the range of applications for volume holographic gratings. Efficiency of a diffractive optic grating is dependent on its approximation of an ideal analog profile using discrete features. The size of constituent features and, consequently, the number that can be used within a required grating period restricts the applications in which diffractive optic gratings can be used. These limitations imply that there are applications which cannot be addressed by either technology. In this paper we propose to address a number of applications in this category with a new method of creating high efficiency gratings which we call stratified diffractive optic gratings. In this approach diffractive optic techniques are used to create an optical structure that emulates volume grating behavior. To illustrate the stratified diffractive optic grating concept we consider a specific application, a scanner for a space-based coherent wind lidar, with requirements that would be difficult to meet by either volume holographic or diffractive optic methods. The lidar instrument design specifies a transmissive scanner element with the input beam normally incident and the exiting beam deflected at a fixed angle from the optical axis. The element will be rotated about the optical axis to produce a conical scan pattern. The wavelength of the incident beam is 2.06 microns and the required deflection angle is 30 degrees, implying a grating period of approximately 4 microns. Creating a high efficiency volume grating with these parameters would require a grating thickness that cannot be attained with current photosensitive materials. For a diffractive optic grating, the number of binary steps necessary to produce high efficiency combined with the grating period requires feature sizes and alignment tolerances that are also unattainable with current techniques. Rotation of the grating and integration into a space-based lidar system impose the additional requirements that it be insensitive to polarization orientation, that its mass be minimized and that it be able to withstand launch and space environments.

Chambers, Diana M.↗

DOE OSTI2023

Topological Circular Dichroism in Chiral Multifold Semimetals

Uncovering the physical contents of the nontrivial topology of quantum states is a critical problem in condensed matter physics. Here, we study the topological circular dichroism in chiral semimetals using linear response theory and first-principles calculations. We show that, when the low-energy spectrum respects emergent SO(3) rotational symmetry, topological circular dichroism is forbidden for Weyl fermions, and thus is unique to chiral multifold fermions. This is a result of the selection rule that is imposed by the emergent symmetry under the combination of particle-hole conjugation and spatial inversion. Using first-principles calculations, we predict that topological circular dichroism occurs in CoSi for photon energy below about 0.2 eV. In conclusion, our Letter demonstrates the existence of a response property of unconventional fermions that is fundamentally different from the response of Dirac and Weyl fermions, motivating further study to uncover other unique responses.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS1975

Numerical analysis and parametric studies of the buckling of composite orthotropic compression and shear panels

A computer program is presented which was developed for the combined compression and shear of stiffened variable thickness orthotropic composite panels on discrete springs: boundary conditions are general and include elastic boundary restraints. Buckling solutions are obtained by using a newly developed trigonometric finite difference procedure which improves the solution convergence rate over conventional finite difference methods. The classical general shear buckling results which exist only for simply supported panels over a limited range of orthotropic properties, were extended to the complete range of these properties for simply supported panels and, in addition, to the complete range of orthotropic properties for clamped panels. The program was also applied to parametric studies which examine the effect of filament orientation upon the buckling of graphite-epoxy panels. These studies included an examination of the filament orientations which yield maximum shear or compressive buckling strength for panels having all four edges simply supported or clamped over a wide range of aspect ratios. Panels with such orientations had higher buckling loads than comparable, equal weight, thin skinned aluminum panels. Also included among the parameter studies were examinations of combined axial compression and shear buckling and examinations of panels with rotational elastic edge restraints.

Housner, J. M.↗

DOE OSTI2025

Algorithm to extract direction in 2D discrete distributions and a continuous Frobenius norm

In this study, we present a novel algorithm for determining directionality in 2D distributions of discrete data. We compare a reference dataset with a known direction to a measured dataset with an unknown direction by the Frobenius norm of the difference (FND) to find the unknown direction. To generalize this concept, we develop a continuous Frobenius norm of the difference (CFND) as a continuous analog of the FND and derive its analytical expression. By relating fitted and normalized 2D Gaussian distributions, we show that the CFND approximates the FND, and we validate this relationship with computer simulations. We find that a first-order approximation of the CFND between two similar Gaussian distributions takes the form of an absolute sine function, offering a simple analytical form with potential for specialized applications in segmented inverse beta decay (IBD) neutrino detectors, astronomy, machine learning, and more. Although this method may easily extend to 3D scalar fields, our focus here is on 2D real-valued fields as it directly applies to directionality. Our methodology consists of modeling a 2D Gaussian distribution, binning the data into a histogram, and encoding it as a square matrix. Rotating this matrix around its geometric center and comparing it to a measured dataset using the FND gives us rotational data that we fit with an absolute sine function. The location of the minimum of this fit is the angle closest to the true angle of the direction in the measured dataset. We present the derivation and discuss initial applications of the CFND in our novel algorithm, demonstrating its success in approximating directionality in 2D distributions.

Data Analysis, Statistics and Probability (physics↗

DOE OSTI2019

Validation of Numerical Models of the Offshore Wind Turbine From the Alpha Ventus Wind Farm Against Full-Scale Measurements Within OC5 Phase III

The main objective of the Offshore Code Comparison Collaboration Continuation, with Correlation (OC5) project is validation of aero-hydro-servo-elastic simulation tools for offshore wind turbines (OWTs) through comparison of simulated results to the response data of physical systems. Phase III of the OC5 project validates OWT models against the measurements recorded on a Senvion 5M wind turbine supported by the OWEC Quattropod from the alpha ventus offshore wind farm. The following operating conditions of the wind turbine were chosen for the validation: (1) Idling below the cut-in wind speed; (2) Rotor-nacelle assembly (RNA) rotation maneuver below the cut-in wind speed; (3) Power production below and above the rated wind speed; and (4) Shutdown. A number of validation load cases were defined based on these operating conditions. The following measurements were used for validation: (1) Strains and accelerations recorded on the support structure; (2) Pitch, yaw, and azimuth angles, generator speed, and electrical power recorded from the RNA. Strains were not directly available from the majority of the OWT simulation tools. Therefore, strains were calculated based on out-of-plane bending moments, axial forces, and cross-sectional properties of the structural members. Also, a number of issues arose during the validation: (1) The need for a thorough quality check of sensor measurements; (2) The sensitivity of the turbine loads to the controller and airfoil properties, which were only approximated in the modeling approach; (3) The importance of estimating and applying an appropriate damping value for the structure; and (4) The importance of wind characteristics beyond turbulence on the loads. The simulation results and measurements were compared in terms of time series, discrete Fourier transforms, power spectral densities, probability density functions of strains and accelerometers. A good match was achieved between the measurements and models set up by OC5 Phase III participants.

17 WIND ENERGY↗

NASA NTRS1997

Mechanical signals in plant development: a new method for single cell studies

Cell division, which is critical to plant development and morphology, requires the orchestration of hundreds of intracellular processes. In the end, however, cells must make critical decisions, based on a discrete set of mechanical signals such as stress, strain, and shear, to divide in such a way that they will survive the mechanical loads generated by turgor pressure and cell enlargement within the growing tissues. Here we report on a method whereby tobacco protoplasts swirled into a 1.5% agarose entrapment medium will survive and divide. The application of a controlled mechanical load to agarose blocks containing protoplasts orients the primary division plane of the embedded cells. Photoelastic analysis of the agarose entrapment medium can identify the lines of principal stress within the agarose, confirming the hypothesis that cells divide either parallel or perpendicular to the principal stress tensors. The coincidence between the orientation of the new division wall and the orientation of the principal stress tensors suggests that the perception of mechanical stress is a characteristic of individual plant cells. The ability of a cell to determine a shear-free orientation for a new partition wall may be related to the applied load through the deformation of the matrix material. In an isotropic matrix a uniaxial load will produce a rotationally symmetric strain field, which will define a shear-free plane. Where high stress intensities combine with the loading geometry to produce multiaxial loads there will be no axis of rotational symmetry and hence no shear free plane. This suggests that two mechanisms may be orienting the division plane, one a mechanism that works in rotationally symmetrical fields, yielding divisions perpendicular to the compressive tensor, parallel to the long axis of the cell, and one in asymmetric fields, yielding divisions parallel to the short axis of the cell and the compressive tensor.

NASA Discipline Number 40-50↗

DOE OSTI2019

Validation of Numerical Models of the Offshore Wind Turbine from the Alpha Ventus Wind Farm Against Full-Scale Measurements Within OC5 Phase III: Preprint

The main objective of the Offshore Code Comparison Collaboration Continuation, with Correlation (OC5) project, is validation of aero-hydro-servo-elastic simulation tools for offshore wind turbines (OWTs) through comparison of simulated results to the response data of physical systems. Phase III of the OC5 project validates OWT models against the measurements recorded on Senvion 5M wind turbine supported by the OWEC Quattropod from the alpha ventus offshore wind farm. The following operating conditions of the wind turbine were chosen for the validation: (1) Idling below the cut-in wind speed; (2) Rotornacelle assembly (RNA) rotation maneuver below the cut-in wind speed; (3) Power production below and above the rated wind speed; and (4) Shut-down. A number of validation load cases were defined based on these operating conditions. The following measurements were used for validation: (1) Strains and accelerations recorded on the support structure; (2) Pitch, yaw, and azimuth angles, generator speed, and electrical power recorded from the RNA. Strains were not directly available from the majority of the OWT simulation tools. Therefore, strains were calculated based on out-of-plane bending moments, axial forces, and geometrical properties of the structural members. Also, a number of issues arose during the validation: (1) The need for a thorough quality check of sensor measurements; (2) The sensitivity of the turbine loads to the controller and airfoil properties, which were only approximated in the modeling approach; (3) The importance of estimating and applying an appropriate damping value for the structure; and (4) The importance of wind characteristics beyond turbulence on the loads. The simulation results and measurements were compared in terms of time series, discrete Fourier transforms (DFTs), power spectral densities (PSDs), probability density functions (PDFs) of strains and accelerometers. A good match was achieved between the measurements and models set up by OC5 Phase III participants.

17 WIND ENERGY↗

NASA NTRS1990

A computational procedure for the dynamics of flexible beams within multibody systems

The dynamic analysis of three dimensional elastic beams which experience large rotational and large deformational motions are examined. The beam motion is modeled using an inertial reference for the translational displacements and a body-fixed reference for the rotational quantities. Finite strain rod theories are then defined in conjunction with the beam kinematic description which accounts for the effects of stretching, bending, torsion, and transverse shear deformations. A convected coordinate representation of the Cauchy stress tensor and a conjugate strain definition is introduced to model the beam deformation. To treat the beam dynamics, a two-stage modification of the central difference algorithm is presented to integrate the translational coordinates and the angular velocity vector. The angular orientation is then obtained from the application of an implicit integration algorithm to the Euler parameter/angular velocity kinematical relation. The combined developments of the objective internal force computation with the dynamic solution procedures result in the computational preservation of total energy for undamped systems. The present methodology is also extended to model the dynamics of deployment/retrieval of the flexible members. A moving spatial grid corresponding to the configuration of a deployed rigid beam is employed as a reference for the dynamic variables. A transient integration scheme which accurately accounts for the deforming spatial grid is derived from a space-time finite element discretization of a Hamiltonian variational statement. The computational results of this general deforming finite element beam formulation are compared to reported results for a planar inverse-spaghetti problem.

Downer, Janice Diane↗

DOE OSTI2021

Computational optimal transport for molecular spectra: The fully discrete case

In this study, the use of computational optimal transport is investigated as a tool for comparing two molecular spectra. Unlike other techniques for comparing molecular spectra in a pattern-recognition framework, transport distances simultaneously encode information about line positions and intensities. In addition, it is shown that transport distances are a useful alternative to Euclidean distances as Euclidean distances are based on line-by-line comparisons, while transport distances reflect broader features of molecular spectra and adequately compare spectra with different resolutions. This paper includes a tutorial on the use of optimal transport and investigates several well-chosen examples to illustrate the utility of computational optimal transport for comparing molecular spectra.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1998

Investigating the Near-Infrared Properties of Planetary Nebulae II. Medium Resolution Spectra: Medium Resolution Spectra - 2

We present medium-resolution (R approximately 700) near-infrared (lambda = 1 - 2.5 micrometers) spectra of a sample of planetary nebulae (PNe). A narrow slit was used which sampled discrete locations within the nebulae; observations were obtained at one or more positions in the 41 objects included in the survey. The PN spectra fall into one of four general categories: H1 emission line-dominated PNe, H1 and H2 emission line PNe, H2 emission line-dominated PNe, and continuum-dominated PNe. These categories correlate with morphological type, with the elliptical PNe falling into the first group, and the bipolar PNe primarily in the H2 and continuum emission groups. The categories also correlate with C/O ratio, with the O-rich objects falling into the first group and the C-rich objects in the groups. Other spectral features were observed in all catagories, such as continuum emission from the central star, and warm dust continuum emission towards the long wavelength end of the spectra. H2 was detected in four PNe in this survey for the first time. An analysis was performed using the H2 line ratios in all of the PN spectra in the survey where a sufficient number of lines were observed to determine the ortho-to-para ratio and the rotational and vibrational excitation temperatures of the H-2 in those objects. One unexpected result from this analysis is that the H-2 is excited by absorption of ultraviolet photons in most of the PNe, although there are several PNe in which collisional excitation plays an important role. The correlation between bipolar morphology and H2 emission has been strengthened with the new detections of H2 in this survey.

Hora, Joseph L.↗

DOE OSTI2016

Accelerated iterative beam angle selection in IMRT

Purpose: Iterative methods for beam angle selection (BAS) for intensity-modulated radiation therapy (IMRT) planning sequentially construct a beneficial ensemble of beam directions. In a naïve implementation, the nth beam is selected by adding beam orientations one-by-one from a discrete set of candidates to an existing ensemble of (n − 1) beams. The best beam orientation is identified in a time consuming process by solving the fluence map optimization (FMO) problem for every candidate beam and selecting the beam that yields the largest improvement to the objective function value. This paper evaluates two alternative methods to accelerate iterative BAS based on surrogates for the FMO objective function value. Methods: We suggest to select candidate beams not based on the FMO objective function value after convergence but (1) based on the objective function value after five FMO iterations of a gradient based algorithm and (2) based on a projected gradient of the FMO problem in the first iteration. The performance of the objective function surrogates is evaluated based on the resulting objective function values and dose statistics in a treatment planning study comprising three intracranial, three pancreas, and three prostate cases. Furthermore, iterative BAS is evaluated for an application in which a small number of noncoplanar beams complement a set of coplanar beam orientations. This scenario is of practical interest as noncoplanar setups may require additional attention of the treatment personnel for every couch rotation. Results: Iterative BAS relying on objective function surrogates yields similar results compared to naïve BAS with regard to the objective function values and dose statistics. At the same time, early stopping of the FMO and using the projected gradient during the first iteration enable reductions in computation time by approximately one to two orders of magnitude. With regard to the clinical delivery of noncoplanar IMRT treatments, we could show that optimized beam ensembles using only a few noncoplanar beam orientations often approach the plan quality of fully noncoplanar ensembles. Conclusions: We conclude that iterative BAS in combination with objective function surrogates can be a viable option to implement automated BAS at clinically acceptable computation times.

61 RADIATION PROTECTION AND DOSIMETRY↗