SEARCH · Search NASA

Results for “Angular Reference”

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 559 records · Page 31

All sources NASA NTRS DOE OSTI

DOE OSTI2017

Evaluation of different inertial control methods for variable-speed wind turbines simulated by fatigue, aerodynamic, structures and turbulence (FAST)

To mitigate the degraded power system inertia and undesirable primary frequency response caused by large-scale wind power integration, the frequency support capabilities of variable-speed wind turbines is studied in this work. This is made possible by controlled inertial response, which is demonstrated on a research turbine - controls advanced research turbine, 3-bladed (CART3). Two distinct inertial control (IC) methods are analysed in terms of their impacts on the grids and the response of the turbine itself. The released kinetic energy in the IC methods are determined by the frequency measurement or shaped active power reference in the turbine speed-power plane. The wind turbine model is based on the high-fidelity turbine simulator fatigue, aerodynamic, structures and turbulence, which constitutes the aggregated wind power plant model with the simplified power converter model. The IC methods are implemented over the baseline CART3 controller, evaluated in the modified 9-bus and 14-bus testing power grids considering different wind speeds and different wind power penetration levels. The simulation results provide various insights on designing such kinds of ICs. The authors calculate the short-term dynamic equivalent loads and give a discussion about the turbine structural loadings related to the inertial response.

17 WIND ENERGY↗

NASA NTRS2020

Performance of NOAA-20 VIIRS Solar Diffuser and Solar Diffuser Stability Monitor

Visible Infrared Imaging Radiometer Suite (VIIRS) radiometrically calibrates its reflective solar bands (RSBs) primarily through a sunlit onboard solar diffuser (SD). The sunlit SD provides a known radiance under the condition that the absolute product of the SD screen transmittance and the bidirectional reflectance distribution function (BRDF) along the SD-to-telescope direction is accurately known. The BRDF changes due to solar exposure. The change, referred to as the H-factor, is monitored by the onboard SD stability monitor (SDSM). The accuracy of the retrieved H-factor propagates to the retrieved F-factor which corrects the scene spectral radiance. High accuracy of the retrieved H-factor relies on high accuracies in the SDSM screen relative effective transmittance and the relative product of the SD screen effective transmittance and the BRDF at the mission start, and a high SDSM detector signal-to-noise ratio (SNR). This article briefly reviews the algorithms used for the NOAA-20 (N20) VIIRS RSB on-orbit radiometric calibration. Additionally, we show the performance of the N20 VIIRS SDSM, giving the SDSM detector SNRs and the SDSM detector gain temporal changes. We develop a model for the SNRs. The model shows that the decreased SNRs in time are due to the detector gain decreases. We also show the N20 VIIRS SD on-orbit performance, measured by the retrieved H-factor and the estimated standard deviation of its error. The H-factor for the telescope view is obtained from the H-factor for the SDSM view, multiplied by an H-factor angular dependence term. We use an innovative method to determine the angular dependence, using the dependence obtained for the Suomi National Polar-orbiting Partnership (SNPP) VIIRS.

VIIRS↗

NASA NTRSDate not supplied

New Features of the NEQAIR Radiation Code

The longest-lived code for predicting shock layer radiation, NEQAIR, is now in its 5th decade of service. Substantial changes to the code have been made over the previous decade, the most recent report of which was at the 5th Workshop on Radiation in High Temperature Gases in 2014, for the version referred to as NEQAIR14. This paper will review some of the improvements made to the NEQAIR code since then, which is now at v15.2. Some of these features are discussed briefly below. NEQAIR15 and subsequent versions have enabled parallel evaluation of multiple lines of sight. This is accomplished by utilizing the HDF5 file format and placing multiple lines into a single file, LOS.h5, which is used for both input and output. This approach enables straightforward parallel execution both over the number of lines of sight and the number of points per line. For large problems, runtime reduces linearly with the number of nodes deployed since each line is processed independently by a subset of MPI ranks. Three applications of the multi-line solver are discussed. The first has to do with performing loosely coupled radiation-flowfield solutions. In this case the computed absorption and emission coefficients are used to evaluate the total energy absorbed or emitted at each point, allowing evaluation of the volumetric source term in the flowfield. The second computation is for obtaining heat flux from nonuniform flows, which require integration over spherical co-ordinates. These are of particular interest for evaluating radiation on the vehicle backshell. This 3D option improves the angular integration scheme and allows adaptive line selection that together reduce the number of lines required by about an order of magnitude. The final application is for remote observation, which is essentially the 3D integration problem over a small solid angle. For all three of these computations, data can be stored in the HDF5 file which allows a NEQAIR run to be restarted when it times out, or to add atmospheric absorption or instrument scan functions. An additional level of parallelism is enabled in NEQAIR15.2 using GPU routines. The GPU parallelism has realized up to 8x speed-up when running on a single core but diminishes as CPU parallelism is increased. For running multi-line simulations, it may be easier to reserve a large number of CPU nodes than to obtain the number of GPU nodes required for similar performance. A GUI, known as NEQTPY, allows for reading and creating input files, running NEQAIR, and displaying results. A significant feature of NEQTPY is the ability to perform spectral fits to data. The fits can operate on a single line spectrum (radiance vs. wavelength) or a 3D input file with multiple columns of data. Other new features include improved constants, additional species, more detailed non-Boltzmann modelling, advanced user controls, the ability to read and calculate spectra from HITRAN datafiles, photodissociation and photoionization cross-sections. A “fast” automatic grid option may reduce the size and time of spectral calculations while still maintaining good accuracy for total heat flux.

Brett A Cruden↗

NASA NTRSDate not supplied

An Updated Secondary Lunar Meteoroid Ejecta Model for Engineering Design

Introduction: The surface of the Moon is constantly being bombarded by a flux of meteoroids of various sizes. Impacts due to these meteoroids produce secondary ejecta material at much lower speeds but with a total mass larger than the original impactor. Details about the secondary ejecta are important for planning missions on the lunar surface. In this work, an updated ejecta model is presented called the Meteoroid Model of Secondary Ejecta (MeMoSeE), to re-place the Apollo-era ejecta model, NASA SP-8013 [1], in the SLS-SPEC 159 Design Specification for Natural Environments (DSNE) [2]. The model pro-duces secondary ejecta flux environments for a user-specified location on the lunar surface, and sorts the incoming secondary flux by angular direction and speed. Methods: MeMoSeE is separated into three parts: the inputs, the conversion step, and the integration of fluxes. Inputs to the model include the primary meteoroid fluxes and the primary near-Earth object (NEO) fluxes. Meteoroid fluxes, both asteroidal and cometary, are calculated using the Meteoroid Engineering Model (MEM3) [3] for different locations on the Moon. For each surface location, an ephemeris is generated using the JPL HORIZONS System [4] that feeds into MEM3. The NEO fluxes are approximated by the high-density population of MEM3 (i.e., only the directionality), where the speed distribution of the NEO fluxes is renormalized to match observations [5]. The regolith properties are used as defined in the DSNE [2]. The conversion step utilizes scaling laws given by Housen & Holsapple 2011 [6] to convert the primary impactor flux to the total mass of secondary ejecta. The ejecta distribution, at the point-of-impact (POI), is separated into a zenith angle and azimuthal angle distribution. The zenith angle distribution follows a beta distribution where the peak depends on the impact altitude angle and the impact azimuth [7]. We employ an ejecta azimuth distribution that is based on Rival & Mandeville 1999 [8] which focuses ejecta in the downstream direction for more oblique impacts. Finally, during the integration step, we sum secondary ejecta number fluxes at a particular region-of-interest (ROI) that originated from many POI locations over the entire surface of the Moon. We keep track of both altitude and azimuth angle bins as well as a range of speed bins, following the igloo gridding as done in MEM3 [9]. The ejecta particle size distribution and density is assumed to be the same as the lunar regolith [10, 2]. Results: The primary fluxes are computed for one Metonic cycle (19 years) for various locations over the lunar surface with a fixed orientation. Both the angular and speed distributions of the primary fluxes are dependent on the latitude and longitude. The over-all primary fluxes show a roughly 13% increase from the eastern limb to the western limb. In general, the speed distribution of the secondary fluxes span from a user-defined minimum speed to the escape speed of the Moon (2.38 km/s), roughly following a power-law relation [6]. Different parts of the speed distribution come from different primary impact locations on the Moon. The secondary ejecta is dominated by the slowest speeds, where these particles originate nearby the ROI. For speeds around 71% of the escape speed, the secondary ejecta originates from locations near the antipodal point. On the other hand, for secondary ejecta speeds that exceed roughly 90% the escape speed, the ejecta particles come from all over the lunar surface to the ROI. Comparing the secondary ejecta fluxes from MeMoSeE with NASA SP-8013 [1], there is a reduction by about 2-3 orders of magnitude for secondary ejecta particles greater than 1 µg. These estimates agree with recent findings from Bjorkman & Christiansen 2019 [11]. The secondary ejecta fluxes are also compared with the primary fluxes, where the ejecta fluxes are roughly an order of magnitude greater than the primary fluxes. References: [1] Cour-Palais, B. G., (1969) NASA SP-8013. [2] NASA SLS-SPEC-159 Rev. H (2020). [3] Moorhead, A. V., et al. (2019) JS&R, 1-17. [4] Giorgini, J. D., (2015) IAUGA, 29, 2256293. [5] Moorhead, A. V., (2020) Memo OSMA/MEO/Lunar-001. [6] Housen, K. R., and Holsapple, K. A., (2011) Icarus, 211(1), 856-875. [7] Gault, D. E., and Wedekind, J. A., (1978) L&PSCP, 9, 3843-3875. [8] Rival, M., and Mandeville, J., (1999) Space Debris, 1(1), 45-57. [9] Moorhead, A. V., (2019) MEM3 User Guide. [10] Carrier, W. D., (1973) The Moon, 6(3-4), 250-263. [11] Bjorkman, M. D. & Christiansen, E. L., (2019) ODC

impacts↗

DOE OSTI2017

Galaxy And Mass Assembly (GAMA): Gas Fueling of Spiral Galaxies in the Local Universe. I. The Effect of the Group Environment on Star Formation in Spiral Galaxies

We quantify the effect of the galaxy group environment (for group masses of 10{sup 12.5}–10{sup 14.0} M {sub ⊙}) on the current star formation rate (SFR) of a pure, morphologically selected sample of disk-dominated (i.e., late-type spiral) galaxies with redshift ≤0.13. The sample embraces a full representation of quiescent and star-forming disks with stellar mass M {sub *} ≥ 10{sup 9.5} M {sub ⊙}. We focus on the effects on SFR of interactions between grouped galaxies and the putative intrahalo medium (IHM) of their host group dark matter halos, isolating these effects from those induced through galaxy–galaxy interactions, and utilizing a radiation transfer analysis to remove the inclination dependence of derived SFRs. The dependence of SFR on M {sub *} is controlled for by measuring offsets Δlog(ψ {sub *}) of grouped galaxies about a single power-law relation in specific SFR, ψ{sub ∗}∝M{sub ∗}{sup −0.45±0.01}, exhibited by non-grouped “field” galaxies in the sample. While a small minority of the group satellites are strongly quenched, the group centrals and a large majority of satellites exhibit levels of ψ {sub *} statistically indistinguishable from their field counterparts, for all M {sub *}, albeit with a higher scatter of 0.44 dex about the field reference relation (versus 0.27 dex for the field). Modeling the distributions in Δlog(ψ {sub *}), we find that (i) after infall into groups, disk-dominated galaxies continue to be characterized by a similar rapid cycling of gas into and out of their interstellar medium shown prior to infall, with inflows and outflows of ∼1.5–5 x SFR and ∼1–4 x SFR, respectively; and (ii) the independence of the continuity of these gas flow cycles on M {sub *} appears inconsistent with the required fueling being sourced from gas in the circumgalactic medium on scales of ∼100 kpc. Instead, our data favor ongoing fueling of satellites from the IHM of the host group halo on ∼Mpc scales, i.e., from gas not initially associated with the galaxies upon infall. Consequently, the color–density relation of the galaxy population as a whole would appear to be primarily due to a change in the mix of disk- and spheroid-dominated morphologies in the denser group environment compared to the field, rather than to a reduced propensity of the IHM in higher-mass structures to cool and accrete onto galaxies. We also suggest that the required substantial accretion of IHM gas by satellite disk-dominated galaxies will lead to a progressive reduction in the specific angular momentum of these systems, thereby representing an efficient secular mechanism to transform morphology from star-forming disk-dominated types to more passive spheroid-dominated types.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

AND YET IT MOVES: THE DANGERS OF ARTIFICIALLY FIXING THE MILKY WAY CENTER OF MASS IN THE PRESENCE OF A MASSIVE LARGE MAGELLANIC CLOUD

Motivated by recent studies suggesting that the Large Magellanic Cloud (LMC) could be significantly more massive than previously thought, we explore whether the approximation of an inertial Galactocentric reference frame is still valid in the presence of such a massive LMC. We find that previous estimates of the LMC’s orbital period and apocentric distance derived assuming a fixed Milky Way (MW) are significantly shortened for models where the MW is allowed to move freely in response to the gravitational pull of the LMC. Holding other parameters fixed, the fraction of models favoring first infall is reduced. Due to this interaction, the MW center of mass within the inner 50 kpc can be significantly displaced in phase-space in a very short period of time that ranges from 0.3 to 0.5 Gyr by as much as 30 kpc and 75 km s{sup −1}. Furthermore, we show that the gravitational pull of the LMC and response of the MW are likely to significantly affect the orbit and phase space distribution of tidal debris from the Sagittarius dwarf galaxy (Sgr). Such effects are larger than previous estimates based on the torque of the LMC alone. As a result, Sgr deposits debris in regions of the sky that are not aligned with the present-day Sgr orbital plane. In addition, we find that properly accounting for the movement of the MW around its common center of mass with the LMC significantly modifies the angular distance between apocenters and tilts its orbital pole, alleviating tensions between previous models and observations. While these models are preliminary in nature, they highlight the central importance of accounting for the mutual gravitational interaction between the MW and LMC when modeling the kinematics of objects in the MW and Local Group.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2001

High Order Difference Method for Low Mach Number Aeroacoustics

A high order finite difference method with improved accuracy and stability properties for computational aeroacoustics (CAA) at low Mach numbers is proposed. The Euler equations are split into a conservative and a symmetric non- conservative portion to allow the derivation of a generalized energy estimate. Since the symmetrization is based on entropy variables, that splitting of the flux derivatives is referred to as entropy splitting. Its discretization by high order central differences was found to need less numerical dissipation than conventional conservative schemes. Owing to the large disparity of acoustic and stagnation quantities in low Mach number aeroacoustics, the split Euler equations are formulated in perturbation form. The unknowns are the small changes of the conservative variables with respect to their large stagnation values. All nonlinearities and the conservation form of the conservative portion of the split flux derivatives can be retained, while cancellation errors are avoided with its discretization opposed to the conventional conservative form. The finite difference method is third-order accurate at the boundary and the conventional central sixth-order accurate stencil in the interior. The difference operator satisfies the summation by parts property analogous to the integration by parts in the continuous energy estimate. Thus, strict stability of the difference method follows automatically. Spurious high frequency oscillations are suppressed by a characteristic-based filter similar to but without limiter. The time derivative is approximated by a 4-stage low-storage second-order explicit Runge-Kutta method. The method has been applied to simulate vortex sound at low Mach numbers. We consider the Kirchhoff vortex, which is an elliptical patch of constant vorticity rotating with constant angular frequency in irrotational flow. The acoustic pressure generated by the Kirchhoff vortex is governed by the 2D Helmholtz equation, which can be solved analytically using separation of variables.

Mueller, B.↗

DOE OSTI2014

Planck intermediate results: XVII. Emission of dust in the diffuse interstellar medium from the far-infrared to microwave frequencies

The dust-Hi correlation is used to characterize the emission properties of dust in the diffuse interstellar medium (ISM) from far infrared wavelengths to microwave frequencies. The field of this encompasses the part of the southern sky best suited to study the cosmic infrared and microwave backgrounds. We cross-correlate sky maps from Planck , the Wilkinson Microwave Anisotropy Probe ( WMAP ), and the diffuse infrared background experiment ( DIRBE ), at 17 frequencies from 23 to 3000 GHz, with the Parkes survey of the 21 cm line emission of neutral atomic hydrogen, over a contiguous area of 7500 deg 2 centred on the southern Galactic pole. We present a general methodology to study the dust-Hi correlation over the sky, including simulations to quantify uncertainties. Our analysis yields four specific results. (1) We map the temperature, submillimetre emissivity, and opacity of the dust per H-atom. The dust temperature is observed to be anti-correlated with the dust emissivity and opacity. We interpret this result as evidence of dust evolution within the diffuse ISM. The mean dust opacity is measured to be (7.1 ± 0.6) × 10 -27 cm 2 H -1 × (v/353 GHz) 1.53 ± 0.03 for 100 ≤ v ≤ 353 GHz. This is a reference value to estimate hydrogen column densities from dust emission at submillimetre and millimetre wavelengths. (2) We map the spectral index β mm of dust emission at millimetre wavelengths (defined here as v ≤ 353 GHz), and find it to be remarkably constant at β mm = 1.51 ± 0.13. We compare it with the far infrared spectral index β FIR derived from greybody fits at higher frequencies, and find a systematic difference, β mm -β FIR = -0.15, which suggests that the dust spectral energy distribution (SED) flattens at v ≤ 353 GHz. (3) We present spectral fits of the microwave emission correlated with Hi from 23 to 353 GHz, which separate dust and anomalous microwave emission (AME). We show that the flattening of the dust SED can be accounted for with an additional component with a blackbody spectrum. This additional component, which accounts for (26 ± 6)% of the dust emission at 100GHz, could represent magnetic dipole emission. Alternatively, it could account for an increasing contribution of carbon dust, or a flattening of the emissivity of amorphous silicates, at millimetre wavelengths. These interpretations make different predictions for the dust polarization SED. (4) We analyse the residuals of the dust-Hi correlation. We identify a Galactic contribution to these residuals, which we model with variations of the dust emissivity on angular scales smaller than that of our correlation analysis. This model of the residuals is used to quantify uncertainties of the CIB power spectrum in a companion Planck paper.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2018

Iterative map-making with two-level preconditioning for polarized cosmic microwave background data sets: A worked example for ground-based experiments

Context. An estimation of the sky signal from streams of time ordered data (TOD) acquired by the cosmic microwave background (CMB) experiments is one of the most important steps in the context of CMB data analysis referred to as the map-making problem. The continuously growing CMB data sets render the CMB map-making problem progressively more challenging in terms of computational cost and memory in particular in the context of ground-based experiments with their operational limitations as well as the presence of contaminants. Aims. We study a recently proposed, novel class of the Preconditioned Conjugate Gradient (PCG) solvers which invoke two-level preconditioners in the context of the ground-based CMB experiments. We compare them against the PCG solvers commonly used in the map-making context considering their precision and time-to-solution. Methods. We compare these new methods on realistic, simulated data sets reflecting the characteristics of current and forthcoming CMB ground-based experiments. We develop a divide-and-conquer implementation of the approach where each processor performs a sequential map-making for a subset of the TOD. Results. We find that considering the map level residuals, the new class of solvers permits us to achieve a tolerance that is better than the standard approach by up to three orders of magnitude, where the residual level often saturates before convergence is reached. This often corresponds to an important improvement in the precision of the recovered power spectra in particular on the largest angular scales. The new method also typically requires fewer iterations to reach a required precision and therefore shorter run times are required for a single map-making solution. However, the construction of an appropriate two-level preconditioner can be as costly as a single standard map-making run. Nevertheless, if the same problem needs to be solved multiple times, for example, as in Monte Carlo simulations, this cost is incurred only once, and the method should be competitive, not only as far as its precision is concerned but also its performance.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

Theoretical calculations of rotationally inelastic collisions of He with NaK(A {sup 1}Σ{sup +}): Transfer of population, orientation, and alignment

We have performed extensive calculations to investigate thermal energy, rotationally inelastic collisions of NaK (A{sup 1}Σ{sup +}) with He. We determined a potential energy surface using a multi-reference configuration interaction wave function as implemented by the GAMESS electronic structure code, and we have performed coupled channel scattering calculations using the Arthurs and Dalgarno formalism. We also calculate the Grawert coefficients B{sub λ}(j, j′) for each j → j′ transition. These coefficients are used to determine the probability that orientation and alignment are preserved in collisions taking place in a cell environment. The calculations include all rotational levels with j or j′ between 0 and 50, and total (translational and rotational) energies in the range 0.0002–0.0025 a.u. (∼44–550 cm{sup −1}). The calculated cross sections for transitions with even values of Δj tend to be larger than those for transitions with odd Δj, in agreement with the recent experiments of Wolfe et al. (J. Chem. Phys. 134, 174301 (2011)). The calculations of the energy dependence of the cross sections and the calculations of the fraction of orientation and alignment preserved in collisions also exhibit distinctly different behaviors for odd and even values of Δj. The calculations also indicate that the average fraction of orientation or alignment preserved in a transition becomes larger as j increases. We interpret this behavior using the semiclassical model of Derouard, which also leads to a simple way of visualizing the distribution of the angles between the initial and final angular momentum vectors j and j′. Finally, we compare the exact quantum results for j → j′ transitions with results based on the simpler, energy sudden approximation. That approximation is shown to be quite accurate.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2015

ASTEROID LIGHT CURVES FROM THE PALOMAR TRANSIENT FACTORY SURVEY: ROTATION PERIODS AND PHASE FUNCTIONS FROM SPARSE PHOTOMETRY

We fit 54,296 sparsely sampled asteroid light curves in the Palomar Transient Factory survey to a combined rotation plus phase-function model. Each light curve consists of 20 or more observations acquired in a single opposition. Using 805 asteroids in our sample that have reference periods in the literature, we find that the reliability of our fitted periods is a complicated function of the period, amplitude, apparent magnitude, and other light-curve attributes. Using the 805-asteroid ground-truth sample, we train an automated classifier to estimate (along with manual inspection) the validity of the remaining ∼53,000 fitted periods. By this method we find that 9033 of our light curves (of ∼8300 unique asteroids) have “reliable” periods. Subsequent consideration of asteroids with multiple light-curve fits indicates a 4% contamination in these “reliable” periods. For 3902 light curves with sufficient phase-angle coverage and either a reliable fit period or low amplitude, we examine the distribution of several phase-function parameters, none of which are bimodal though all correlate with the bond albedo and with visible-band colors. Comparing the theoretical maximal spin rate of a fluid body with our amplitude versus spin-rate distribution suggests that, if held together only by self-gravity, most asteroids are in general less dense than ∼2 g cm{sup −3}, while C types have a lower limit of between 1 and 2 g cm{sup −3}. These results are in agreement with previous density estimates. For 5–20 km diameters, S types rotate faster and have lower amplitudes than C types. If both populations share the same angular momentum, this may indicate the two types’ differing ability to deform under rotational stress. Lastly, we compare our absolute magnitudes (and apparent-magnitude residuals) to those of the Minor Planet Center’s nominal (G = 0.15, rotation-neglecting) model; our phase-function plus Fourier-series fitting reduces asteroid photometric rms scatter by a factor of ∼3.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

Effects on image quality of a 2D antiscatter grid in x-ray digital breast tomosynthesis: Initial experience using the dual modality (x-ray and molecular) breast tomosynthesis scanner

Purpose: Radiation scattered from the breast in digital breast tomosynthesis (DBT) causes image degradation, including loss of contrast between cancerous and background tissue. Unlike in 2-dimensional (2D) mammography, an antiscatter grid cannot readily be used in DBT because changing alignment between the tube and detector during the scan would result in unacceptable loss of primary radiation. However, in the dual modality breast tomosynthesis (DMT) scanner, which combines DBT and molecular breast tomosynthesis, the tube and detector rotate around a common axis, thereby maintaining a fixed tube-detector alignment. This C-arm geometry raises the possibility of using a 2D (cellular) focused antiscatter grid. The purpose of this study is to assess change in image quality when using an antiscatter grid in the DBT portion of a DMT scan under conditions of fixed radiation dose. Methods: Two 2D focused prototype grids with 80 cm focal length were tested, one stack-laminated from copper (Cu) and one cast from a tungsten-polymer (W-poly). They were reciprocated using a motion scheme designed to maximize transmission of primary x-ray photons. Grid-in and grid-out scatter-to-primary ratios (SPRs) were measured for rectangular blocks of material simulating 30%, 50%, and 70% glandular tissue compositions. For assessment of changes in image quality through the addition of a grid, the Computerized Imaging Reference Systems, Inc., phantom Model 011A containing a set of 1 cm thick blocks simulating a range of glandular/adipose ratios from 0/100 to 100/0 was used. To simulate 6.5 and 8.5 cm thick compressed breasts, 1 cm thick slices of PMMA were added to the Model 011A phantom. DBT images were obtained with and without the grid, with exposure parameters fixed for a given compressed thickness. Signal-difference-to-noise ratios (SDNRs), contrast, and voxel value-based attenuation coefficients (μ) were measured for all blocks from reconstructed phantom images. Results: For 4, 6, and 8 cm tissue-equivalent block phantom thicknesses, the inclusion of the W-poly grid reduced the SPR by factors of 5, 6, and 5.8, respectively. For the same thicknesses, the copper grid reduced the SPR by factors of 3.9, 4.5, and 4.9. For the 011A phantom, the W-poly grid raised the SDNR of the 70/30 block from 0.8, −0.32, and −0.72 to 0.9, 0.76, and 0.062 for the 4.5, 6.5, and 8.5 cm phantoms, respectively. It raised the SDNR of the 100/0 block from 3.78, 1.95, and 1.0 to 3.79, 3.67, and 3.25 for the 4.5, 6.5, and 8.5 cm phantoms, respectively. Inclusion of the W-poly grid improved the accuracy of image-based μ values for all block compositions. However, smearing of attenuation across slices due to limited angular sampling decreases the sensitivity of voxel values to changing composition compared to theoretical μ values. Conclusions: Under conditions of fixed radiation dose to the breast, use of a 2D focused grid increased contrast, SDNR, and accuracy of estimated attenuation for mass-simulating block compositions in all phantom thicknesses tested, with the degree of improvement depending upon material composition. A 2D antiscatter grid can be usefully incorporated in DBT systems that employ fully isocentric tube-detector rotation.

60 APPLIED LIFE SCIENCES↗

NASA NTRSDate not supplied

End-to-end Modeling of Lyot Coronagraph for LUVOIR: Sensitivity to Segment-Level Aberrations

Direct imaging and spectroscopic characterization of atmospheres of Earth-like exoplanets require starlight suppression to detect a faint planetary signal at small angular separations. Future space telescopes will require segmentation to survey a statistically significant population of nearby stars systems by increasing aperture size beyond the limits set by traditional monolithic apertures; however, segmentation also represents an additional challenge due to the introduction of temporal pupil discontinuities in terms of both phase and amplitude. For this study, a baseline Lyot coronagraph is used with the Large UV / Optical /Infrared (LUVOIR) mission’s 15-m on-axis, segmented optical telescope assembly (OTA). We present a high-fidelity, end-to-end modeling approach developed as part of the Segmented Coronagraph Design and Analysis (SCDA) group that includes diffractive optical propagations at the 1e-10 contrast level across the entire optical train including the OTA and intermediate re-imaging optics inside the coronagraph instrument. We perform a tolerancing analysis to assess the contrast sensitivity to wavefront perturbations for a set of key telescope terms including segment-level low-order and mid-spatial frequency aberrations. As an input to the end-to-end telescope model, we use two industry-provided segment-level wavefront aberration datasets pertaining: (1) thermal gradients and accelerations (courtesy of Ball-Aerospace and L3-Harris), and (2) telescope dynamics including line-of-sight and segment jitter (courtesy of Lockheed-Martin). Contrast performance outputs of the coronagraph are fed into a mission design reference simulator to assess the scientific yield sensitivity. End-to-end modeling will incorporate performance improvements due to active wavefront sensing and control and/or segment-level metrology to improve system performance.

Optics↗

NASA NTRS2001

Improved Spacecraft Materials for Radiation Protection

Methods by which radiation shielding is optimized need to be developed and materials of improved shielding characteristics identified and validated. The galactic cosmic rays (GCR) are very penetrating and the energy absorbed by the astronaut behind the shield is nearly independent of shield composition and even the shield thickness. However, the mix of particles in the transmitted beam changes rapidly with shield material composition and thickness. This results in part from the breakup of the high-energy heavy ions of the GCR which make contributions to biological effects out of proportion to their deposited energy. So the mixture of particles in the radiation field changes with shielding and the control of risk contributions from dominant particle types is critical to reducing the hazard to the astronaut. The risk of biological injury for a given particle type depends on the type of biological effect and is specific to cell or tissue type. Thus, one is faced with choosing materials which may protect a given tissue against a given effect but leave unchanged or even increase the risk of other effects in the same tissue or increase the risks to other adjacent tissues of a different type in the same individual. The optimization of shield composition will then be tied to a specific tissue and risk to that tissue. Such peculiarities arise from the complicated mixture of particles, the nature of their biological response, and the details of their interaction with material constituents. Aside from the understanding of the biological response to specific components, one also needs an accurate understanding of the radiation emerging from the shield material. This latter subject has been a principal element of this project. In the past ten years our understanding of space radiation interactions with materials has changed radically, with a large impact on shield design. For example, the NCRP estimated that only 2 g/sq cm. of aluminum would be required to meet the annual 500 mSv limit for the exposure of the blood forming organs (this limit is strictly for LEO but can be used as a guideline for the Mars mission analysis). The current estimates require aluminum shield thicknesses above 50 g/sq cm., which is impractical. In such a heavily shielded vehicle, the neutrons produced throughout the vehicle also contribute significantly to the exposure and this demands greater care in describing the angular dependence of secondary particle production processes. As such the continued testing of databases and transport procedures in laboratory and spaceflight experiments has continued. This has been the focus of much of the last year's activity and has resulted in improved neutron prediction capability. These new methods have also improved our understanding of the surface environment of Mars. The Mars 2003 NRA HEDS related surface science requirements were driven by the need to validate predictions on the upward flux of neutrons produced in the Martian regolith and bedrock made by the codes developed under this project. The codes used in the surface environment definition are also being used to look at in situ resources for the development of construction material for Martian surface facilities. For example, synthesis of polyimides and polyethylene as binders of regolith for developing basic structural elements has been studied and targets built for accelerator beam testing of radiation shielding properties. Preliminary mechanical tests have also been promising. Improved spacecraft materials have been identified (using the criteria reported by this project at the last conference) as potentially important for future shielding materials. These are liquid hydrogen, hydrogenated nanofibers, liquid methane, LiH, Polyethylene, Polysulfone, and Polyetherimide (in order of decreasing shield performance). Some of the materials are multifunctional and are required for other onboard systems. We are currently preparing software for trade studies with these materials relative to the Mars Reference Mission as required in the project's final year.

Wilson, John W.↗

DOE OSTI2016

Advanced Energy Harvesting Control Schemes for Marine Renewable Energy Devices

The Advanced Energy Harvesting Control Schemes for Marine Renewable Energy Devices (Project) investigated, analyzed and modeled advanced turbine control schemes with the objective of increasing the energy harvested by hydrokinetic turbines in turbulent flow. Ocean Renewable Power Company (ORPC) implemented and validated a feedforward controller to increase power capture; and applied and tested the controls on ORPC’s RivGen® Power Systems in Igiugig, Alaska. Assessments of performance improvements were made for the RivGen® in the Igiugig environment and for ORPC’s TidGen® Power System in a reference tidal environment. Annualized Energy Production (AEP) and Levelized Cost of Energy (LCOE) improvements associated with implementation of the recommended control methodology were made for the TidGen® Power System in the DOE reference tidal environment. System Performance Advancement (SPA) goals were selected for the project. SPA targets were to improve Power to Weight Ratio (PWR) and system Availability, with the intention of reducing Levelized Cost of Electricity (LCOE). This project focused primarily reducing in PWR. Reductions in PWR of 25.5% were achieved. Reductions of 20.3% in LCOE were achieved. This project evaluated four types of controllers which were tested in simulation, emulation, a laboratory flume, and the field. The adaptive Kω2 controller performs similarly to the non-adaptive version of the same controller and may be useful in tidal channels where the mean velocity is continually evolving. Trends in simulation were largely verified through experiments, which also provided the opportunity to test assumptions about turbine responsiveness and control resilience to varying scales of turbulence. Laboratory experiments provided an essential stepping stone between simulation and implementation on a field-scale turbine. Experiments also demonstrated that using “energy loss” as a metric to differentiate between well-designed controllers operating at an optimal tip-speed ratio set-point is difficult, which anticipated the outcome from field experiments. The clear message is that the feedforward Kω2 controller out-performs the feedback controllers in almost all aspects and modes of evaluation. The controllers proved a substantial improvement over the baseline performance of the TidGen® turbine, in terms of energy capture. The effects of noise-contaminated angular velocity signals were investigated and validated by simulation as an explanation for the performance limitations observed for TidGen® turbine operations in Eastport, Maine. Measurements of loads performed as part of the laboratory testing indicate that there are limited differences in average axial thrust force between control architectures. This suggests that none of the control strategies are likely to substantially affect loads on the turbine support structure. Velocity measurements during the ORPC RivGen® turbine deployment at Igiugig, Alaska, in 2014 were used to assess the variability of the river flow. Results suggest that the river flow is approximately steady, in the mean sense, at any particular location in the river, with random turbulent fluctuations that are around 10% of the mean flow. The mean flow in the center channel of the river is 2.5 m/s, with reductions near the riverbanks and in the shallows. Spectral analysis and lagged correlation results indicate that temporal fluctuations at a given point are dominated by large scale fluctuations, such that measurements at the turbine location are just as useful for inflow control implementation as upstream measurements. At this site, and likely at many other river sites, flow is generally steady at a given location, but flow varies dramatically between locations, particularly laterally across the river. The primary result is that a lateral change in position of a few meters results in changes to flow speed that far exceed the turbulence fluctuations at any given location. The turbulence is dominated by long time scales. Following final system tests, the RivGen® device was submerged and each evaluated controller was tested across a range of gain/set point values and filter configurations for a minimum of 5 minutes, with longer runs attempted for well-performing cases. In addition to testing controllers during the 2015 deployment season, LGL Alaska Research Associates Inc. (LGL) performed a fish monitoring study in compliance with Alaska Department of Fish and Game fisheries habitat permit for the project. During the season, LGL reviewed 10% of 179 one-hour blocks of footage (6 minutes on the hour) and documented a total of over 1200 fish in the vicinity of the RivGen® device, including over 800 salmon smolt and over 350 adult salmon. No evidence of adverse effects including passage delay by upstream migrating salmon was noted. This is an important result for future deployments and has a direct impact on commercial system designs.

16 TIDAL AND WAVE POWER↗

NASA NTRS1995

Investigation of computational and spectral analysis methods for aeroacoustic wave propagation

Most computational fluid dynamics (CFD) schemes are not adequately accurate for solving aeroacoustics problems, which have wave amplitudes several orders of magnitude smaller yet with frequencies larger than the flow field variations generating the sound. Hence, a computational aeroacoustics (CAA) algorithm should have minimal dispersion and dissipation features. A dispersion relation preserving (DRP) scheme is, therefore, applied to solve the linearized Euler equations in order to simulate the propagation of three types of waves, namely: acoustic, vorticity, and entropy waves. The scheme is derived using an optimization procedure to ensure that the numerical derivatives preserve the wave number and angular frequency of the partial differential equations being discretized. Consequently, simulated waves propagate with the correct wave speeds and exhibit their appropriate properties. A set of radiation and outflow boundary conditions, compatible with the DRP scheme and derived from the asymptotic solutions of the governing equations, are also implemented. Numerical simulations are performed to test the effectiveness of the DRP scheme and its boundary conditions. The computed solutions are shown to agree favorably with the exact solutions. The major restriction appears to be that the dispersion relations can be preserved only for waves with wave lengths longer than four or five spacings. The boundary conditions are found to be transparent to the outgoing disturbances. However, when the disturbance source is placed closer to a boundary, small acoustic reflections start appearing. CAA generates enormous amounts of temporal data which needs to be reduced to understand the physical problem being simulated. Spectral analysis is one approach that helps us in extracting information which often can not be easily interpreted in the time domain. Thus, three different methods for the spectral analysis of numerically generated aeroacoustic data are studied. First, the capabilities of two traditional methods for spectral analysis, namely, the Blackman-Tukey method and periodogram method, are compared in estimating the spectra of a simple-periodic process. The periodogram is then applied to analyze transitory-deterministic processes. Finally, these two methods are compared with a more recent method, referred as the Weighted-Overlapped-Segment-Averaging (WOSA) method, in estimating the spectra of a chaotic (random-like) process. From the demonstrative case for the spectral analyses of data generated by simple-periodic process, the periodogram method is found to give a better estimate of the steep-sloped spectra than the Blackman-Tukey method. Also, for this problem, the Hanning window is found to perform better with the periodogram method than with the Blackman-Tukey method. Finally, for the spectral analysis of data generated by the chaotic process, the periodogram method does not perform well, whereas, the WOSA and Blackman-Tukey methods give equivalently good results.

Vanel, Florence O.↗

DOE OSTI2015

THE ACCRETION OF DARK MATTER SUBHALOS WITHIN THE COSMIC WEB: PRIMORDIAL ANISOTROPIC DISTRIBUTION AND ITS UNIVERSALITY

The distribution of galaxies displays anisotropy on different scales and it is often referred to as galaxy alignment. To understand the origin of galaxy alignments on small scales, one must investigate how galaxies were accreted in the early universe and quantify their primordial anisotropy at the time of accretion. In this paper we use N-body simulations to investigate the accretion of subhalos, focusing on their alignment with halo shape and the orientation of mass distribution on the large scale, defined using the Hessian matrix of the density field. The large/small (e1/e3) eigenvalues of the Hessian matrix define the fast/slow collapse direction of matter on the large scale. We find that: (1) the halo major axis is well aligned with the e3 (slow collapse) direction, and it is stronger for massive halos; (2) subhalos are predominantly accreted along the major axis of the host halo, and the alignment increases with the host halo mass. Most importantly, this alignment is universal; (3) accretion of subhalos with respect to the e3 direction is not universal. In massive halos, subhalos are accreted along the e3 (even more strongly than the alignment with the halo major axis), but in low-mass halos subhalos are accreted perpendicular to e3. The transitional mass is lower at high redshift. The last result well explains the puzzling correlation (both in recent observations and simulations) that massive galaxies/halos have their spin perpendicular to the filament, and the spin of low-mass galaxies/halos is slightly aligned with the filament, under the assumption that the orbital angular momentum of subhalos is converted to halo spin.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Predicting magnetic properties of single-molecule magnets from self-interaction-free density-functional theory (Final Report)

In this project we investigated electronic structure of an intermediate-sized copper-based molecule and magnetic and hyperfine properties of several small non-magnetic and magnetic molecules including transition-metal elements by applying self-interaction corrections to density functional theory (DFT). In a sufficient number of cases, DFT-calculated electronic structure and magnetic properties of single-molecule magnets qualitatively differ from corresponding experimental data. This is partly due to self-interacting electrons within the DFT formalism. Recently, an efficient method to correct the self-interactions was proposed, i.e., Fermi-Lowdin prbital (FLO) based self-interaction corrected (SIC) methodology, within DFT. Henceforth, this method is referred to as FLO-SIC method which exists in FLOSIC code. We used this FLO-SIC method for our studies of electronic structure and magnetic and hyperfine properties of small magnetic molecules and non-magnetic molecules. Our study will provide insight into predictions of magnetic properties of single-molecule magnets where self-interaction corrections play a critical role. There are two components of this project. In the first work, we studied the electronic structure of a planar mononuclear Cu-based molecule in two oxidation states, using DFT with the FLO-SIC method. We chose this system because it is small enough and it includes a transition metal element. We found that the standard FLO-SIC method takes too much compute time even for the small transition-metal molecule and so we slightly modified the method in order to expedite the process. In the dianionic state, we found that the FLO-SIC spin density agrees quantitatively with accurate quantum chemistry methods, while DFT spin density without self-interaction corrections are severely deviated from the quantum chemistry methods. We also showed that the energy gap between the highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) of the dianionic state is larger than that of the monoanionic state. This result is consistent with experimental data. In the second work, we investigated how the interaction between the electron spin and nuclear spin is affected by electron self-interactions within small non-magnetic and magnetic molecules. Such an interaction is called hyperfine interaction. For molecules without significant orbital angular momentum, the hyperfine interaction consists of Fermi contact and dipolar interaction terms. Since the Fermi contact term depends on electron spin density at the nuclear site, it would be highly affected by self-interaction corrections. Therefore, we calculated the hyperfine interaction for the small molecules using DFT with the slightly modified expedited FLO-SIC method which was obtained in the first work, and compared the results to experimental data and DFT calculations without self-interaction corrections. We found significant improvement of the Fermi contact term computed using the FLO-SIC method for small magnetic molecules. Overall, the first and second work provided positive outlook of application of the FLO-SIC method to magnetic molecules and systems including transition-metal elements.

36 MATERIALS SCIENCE↗