Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “smoothing”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Post-Shot Report for OMEGA Double Cylinders (CylDRT 22B)

The direct-drive double cylinder experimental platform is a high-energy-density (HED) science platform designed to image an imploding cylindrical target. The target consists of a directly-driven outer cylinder and a shock-and-collision driven inner cylinder. The purpose of this platform is to study hydrodynamic instability growth on the inner cylinder, the outer surface of which is classically Rayleigh-Taylor unstable during the acceleration phase. We present results from recent experiments at the OMEGA laser facility. These experiments were designed as a proof-of-principle for the platform, using the same cylinder exterior dimensions and direct-drive beam configuration as previous single cylinder experiments. In these experiments, three sets of targets were fielded: no machined perturbations (smooth), a sinusoidal mode-10 perturbation on the outer surface of the inner cylinder, and a mode-20 perturbation on the outer surface of the inner cylinder. The primary diagnostic was a gated x-ray framing camera which imaged the backlit inner cylinder on-axis for sixteen frames over a time window of 1 ns for each shot. A second side-lit radiograph captured one image per shot, diagnosing axial uniformity. In this report we include an overview of the results from both the backlighter and the sidelighter diagnostics. We discuss at length the experimental analysis process. Finally, we present the results of the smooth target implosion trajectory and compare them to post-shot simulations. We see favorable agreement between simulation and experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Exploiting Task Tolerances in Mimicry-based Telemanipulation

We explore task tolerances, i.e., allowable position or rotation inaccuracy, as an important resource to facilitate smooth and effective telemanipulation. Task tolerances provide a robot flexibility to generate smooth and feasible motions; however, in teleoperation, this flexibility may make the user’s control less direct. In this work, we implemented a telema nipulation system that allows a robot to autonomously adjust its configuration within task tolerances. We conducted a user study comparing a telemanipulation paradigm that exploits task tolerances (functional mimicry) to a paradigm that requires the robot to exactly mimic its human operator (exact mimicry), and assess how the choice in paradigm shapes user experience and task performance. Our results show that autonomous adjustments within task tolerances can lead to performance improvements without sacrificing perceived control of the robot. Additionally, we find that users perceive the robot to be more under control, predictable, fluent, and trustworthy in functional mimicry than in exact mimicry.

42 ENGINEERING↗
DOE OSTI2023

Improving Nb3Sn Cavity Performance Using Centrifugal Barrel Polishing

Nb3Sn superconducting radiofrequency (SRF) cavities have been an ongoing research topic for many years motivated by the potential for higher accelerating gradients and quality factors compared to niobium SRF cavities. The highest performing Nb3Sn cavities are manufactured using tin vapor-diffusion coating, which creates a Nb3Sn film with a surface roughness of around 100-200 nm. This is thought to be one of the limiting factors for the accelerating gradient of Nb3Sn cavities due to enhancement of magnetic field near sharp surface features. To smooth Nb3Sn SRF cavities, we have developed a mechanical polishing procedure which uses centrifugal barrel polishing to smooth the surface followed by a secondary tin coating step to repair the surface. We show that the accelerating field of a Nb3Sn SRF cavity is improved by applying this procedure. We also investigate the quench mechanism of the polished cavity by utilizing temperature mapping to measure which regions of the cavity e xperience heating during RF operation. We then cut samples from these regions and analyze the film microstructure and chemical composition in 3D using EDS and EBSD measurements together with a focused ion-beam (FIB) tomography technique.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Improvement of Numerical Methods in Pronghorn

During the fiscal year of 2021, the finite-volume method (FVM) was deployed in Pronghorn to provide improved efficiency, stability, and accuracy for coarse-mesh, thermal-hydraulics problems. While the main goals of fiscal year 2021 were met, several issues emerged from the early deployment of the finite-volume method in Pronghorn. These issues were: The compressible and incompressible formulations that were implemented are inadequate for many nuclear reactor flow problems; Omission of terms accounting for the porosity and Darcy-Forchheimer body force discontinuities in the Rhie-Chow interpolation lead to oscillations in pressure and velocity at these discontinuities; The FVM lacks a correction for non-orthogonal grids for computing accurate pressure gradient leading to loss of accuracy in regions with skewed elements; The FVM currently uses a monolithic solver. Monolithic solvers have issues dealing with the saddle-point nature of the discretized fluid equations. The result are bad convergence if direct factorization is not used and large memory consumption when direct factorization is used. These four issues are addressed in this report. In particular, we report the completion of the following task: Implementation of a weakly compressible formulation in the MOOSE Navier-Stokes module; Implementation of Moukalled’s method for including body forces in the RCI. Additionally, we identified the need to smooth the porosity using Moukalled’s face-cell smoothing operator; Implementation of a non-orthogonal correction for the Green-Gauss gradient computation; Preliminary implementation of a SIMPLE segregated solver.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Advancement of hybrid fluid-kinetic modeling for HEDP and ICF science

We report on the development progress of a hybrid fluid-kinetic code for simulating fluids and plasmas in a wide range of environments, such as laser–matter interactions, inertial confinement fusion, magnetic confinement fusion, and pulsed power. The suite of numerical tools under development utilizes heterogeneous computer architectures and leverages the benefits of particle–based simulation techniques. By working to combine the kinetic particle-in-cell (PIC) model with a particle-based fluid simulation technique, such as smoothed particle hydrodynamics, we are developing a flexible framework capable of accurately modeling complex flows within and between kinetic and fluid regimes. The TriForce code is under development as a C++ framework for parallel, 3D, particle-based, hybrid fluid-kinetic plasma simulations. The fluid half of TriForce will be based upon the meshless smoothed-particle-hydrodynamics (SPH) approach, well-suited for shear, mixing, and turbulence, whereas the kinetic half resembles a traditional particle-in-cell (PIC) code; other particle-based approaches to fluid modeling that do use a mesh are also possible to use and are under investigation. Maxwell’s electromagnetic field equations are solved either via explicit or implicit algorithms or approximated via resistive magnetohydrodynamics (MHD) using an Ohm’s law and resulting induction equation (extended MHD is under development). A primary goal of enabling direct comparisons, from the same code, between results from the variants of MHD and implicit electromagnetic solutions is to improve our fundamental understanding of systems with magnetic fields. The code is under development to recover results from both radiation-MHD and fully kinetic codes in those limits, and is continuing to be developed from other follow-on grants to operate in between where both descriptions may co-exist and interact. For certain applications, it is desired for a simulation to contain fluid ions and electrons as well as kinetic ions and electrons. Typically, it is too computationally intensive to model a full-scale ICF or HEDP experiment fully kinetically since many cycles are expended with very small time steps on modeling the fluid part of a material that is well treated by the fluid approximation. In this case, many traditional PIC particles can be replaced with a single fluid particle representing the thermal part of the distribution function, and there are fewer needed kinetic particles, which describe the non-thermal part and can be sub-cycled relative to the fluid particle advance. Furthermore, a pure fluid code may, depending on the problem, simply lack many physically important details that are beyond the scope of its reduced approximations and assumptions. In this report, we summarize the objectives achieved in the development of the collisional and kinetic half of the code, and the physics problems to which the code has been applied in the areas of advanced and innovative fusion concepts, pulsed power, and magneto-inertial fusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Initial Testing of an In Situ Load Retention Aging Vessel

A thermal aging vessel instrumented with load cells was fabricated. The primary function of the vessel is to continuously monitor the in situ load retention of up to three compressed polymer coupons undergoing thermally accelerated aging under nitrogen. A secondary function is to enable gas sampling of the vessel headspace during thermal aging. Heating of the vessel is achieved using a custom heater jacket. To improve upon our conventional aging study methods which require periodic interruption of aging to perform load testing in an Instron machine at room temperature, this technology aims to automate/facilitate data acquisition/analysis, improve data quality, and enable uninterrupted compression of the polymer which represents the service condition. As an example case to assess functionality of the in situ vessel, the load retention of a siloxane elastomer material additively manufactured by direct-ink-writing (DIW) was measured at three different isothermal aging temperatures for ~1 month. Initial compression of the coupons while near the aging temperature was achieved by temporarily opening the heated vessel to access the interior chamber and manually tightening four nuts to drive the heated compression plate down onto the heated coupons. Initial testing demonstrated achievement of the primary load retention monitoring function. Unfortunately, the vessel leaked which prevented gas sampling; an active purge was used to maintain a nitrogen atmosphere. Welded or otherwise sealed joints, which could be implemented in a future design, would likely eliminate leak paths. To apply time-temperature superposition (TTS), a technique used to provide long-term prediction of the load retention from short-term isothermal data, the load retention needed to be calculated relative to the load at an estimated “equilibrium” time, after most of the transient viscoelastic physical relaxation occurred. The peak load immediately after compression could not be used as the load retention basis for two reasons: (1) age-related changes must be isolated from non-age-related physical relaxation before applying TTS and (2) the manual mechanism used to compress the specimens at the aging temperature was neither smooth nor repeatable which affected the peak load value. To better understand the effect of the mode of initial compression on the measured load, and possibly better estimate “equilibrium” physical relaxation times, systematic stress relaxation experiments were performed using an Instron machine with a thermal chamber. At a given temperature, the DIW polymer was compressed to a fixed strain in either a stepped or continuous manner at two different rates, then held at that strain for 24 hrs. The results indicated that, at a given temperature, the different stress relaxation curves appeared to converge to the same curve at some “equilibrium” time when the non-age-related physical relaxation was mostly complete. Though this observation suggests that the discontinuous manual compression employed by the vessel is feasible, a compression mechanism that is rapid, smooth, and repeatable would enhance its use.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Turbulent Processes that Influence Boundary-Layer Cloud Structure (Final Technical Report)

Many models of the atmosphere simulate vertical profiles of observable fields, such as temperature, moisture, and vertical velocity. To assess the realism of such model output, it is important to evaluate the modeled profiles against observations. However, many extant evaluations consist of comparing time averaged profiles from the model and the observations. Such time averages can smooth out unrealistic features that appear in instantaneous profiles. In this way, the time averages may obscure certain model errors from detection. For instance, many models misrepresent turbulent processes, leading to profiles that are too smoothed out or too noisy. E.g., if a model has coarse vertical resolution and overly strong diffusion, then the model may be unable to produce observed sharp peaks. On the other hand, a model may have artificial approximations that produce spurious jumps in the vertical. Much can be learned by browsing through instantaneous snapshots of model output. But the volume of output produced by a global model is far too great to rely purely on visual inspection. In addition, we need a way to statistically characterize the profiles, so that anomalies can be detected more conveniently.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Precision Polishing of Spheres Via In-Situ Process Monitoring and Machine-Learning-Based Optimization

Precision polishing for high surface finish is a vital step in the making of many high-tech components for optical and semiconductor modules as well as for wide ranging experimental samples for HED science that is at the core of LLNL’s mission. In particular, the fabrication of extremely smooth and round spheres comprised of low atomic number materials (low-Z such as high-density carbon or HDC, beryllium, hydrocarbon polymers, etc.) for Inertial Confinement Fusion (ICF) is entirely reliant on specialized polishing processes to achieve the surface smoothness.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Hybridization Assessment of Trybrid Pumped Storage Hydropower System—Part 1: A Case Study of Corral Summit

This report is a part of the deliverables for technical assistance provided to Cat Creek Energy for the Coral Summit Trybrid (Triple Hybrid-Pumped Storage Hydropower, Battery Energy Storage System, and photovoltaic solar energy) energy project. This report explores the operational benefits and challenges of hybridizing an open loop PSH (200MW) located at Mackay, Custer County, Idaho with solar PV (Ground mount 300MW and floating 40MW) and battery (720MWhr). This document reports two activities performed as a part of the hybridization assessment task 1) optimal resource allocation and energy management strategy, and 2) power quality and reliability assessment. From optimal resource allocation and energy management strategy (activity 1), the following key findings can be observed: • Conventional PSH (CPSH) with two reversible pump turbines and separate penstocks can provide required flexibility equivalent to that from two ternary PSH with separate penstock. With single unit CPSH, upper reservoir head cannot be maintained accurately, the variation of water level is rapid and pump mode flexibility is not available. These disadvantages can be overcome by single unit TPSH. However, using two CPSH units with separate penstocks also overcome these disadvantages with the formation of the hydraulic short circuit between two conventional units. • Flooding of the lower reservoir is a severe concern when considering continuous operation for black start. This limits the duration of continuous operation from PSH alone to around 50 hours. Due to the complementary PV and battery action, the duration of continuous operation and smooth power output can be extended. • An optimization problem is framed that maximizes the power output on an hourly basis while minimizing constraint violations and respecting seasonal variations of solar PV and load profiles . Two value streams, arbitrage and baseload generation are served by this profile. It was uncovered that for smooth power output during regular operation, PV curtailment will be required, or the battery capacity needs to be increased above 90MW to accommodate additional PV. From power quality and reliability assessment (activity 2) the following takeaway points can be observed: • The Trybrid, when integrated at the Lost River bus, and limited to 250MW in generation mode and -150MW in the pump mode, causes no violation of voltage or flow.

13 - HYDRO ENERGY↗
DOE OSTI2026

Electropolishing-induced topographic defects in niobium: Insights and Implications for SRF

Electropolishing (EP) is the premier surface preparation method for high-Q, high-gradient superconducting RF cavities made of Nb. This leaves behind an apparently smooth surface, yet the achievable peak magnetic fields fall well below the superheating field of Nb, in most cases. In this work, the ultimate surface finish of EP was investigated by studying its effect on highly polished Nb samples. EP introduces high slope angle sloped-steps at grain boundaries. The magnetic field enhancement and superheating field suppression factors associated with such a geometry are calculated in the London theory. Despite the by-eye smoothness of electropolished Nb, such defects compromise the stability of the low-loss Meissner state, likely limiting the achievable peak accelerating fields in superconducting RF cavities. Finally, the impact of surface roughness on impurity diffusion is investigated which can link surface roughness to the effectiveness of heat treatments like low-temperature baking or nitrogen infusion in the vortex nucleation or hydride hypotheses. Surface roughness tends to decrease the effective dose of impurities as a result of the expansion of impurities into regions with greater internal angle. The effective dose of impurities can be protected by minimizing slope angles and step heights, ensuring uniformity.

Lechner, Eric [Thomas Jefferson National Accelerat↗
DOE OSTI2023

SPRUCE Vegetation Phenology in Experimental Plots from Phenocam Imagery, 2015-2022

This data set consists of PhenoCam data from the SPRUCE experiment from the beginning of whole ecosystem warming (Hanson et al. 2017) in August 2015 through March 31 of 2023, with start- and end-of-season phenological transition dates derived through the end of autumn 2022. Digital cameras, or phenocams, installed in each SPRUCE enclosure track seasonal variation in vegetation “greenness”, a proxy for vegetation phenology and associated physiological activity. Three separate regions of interest (ROIs) were defined for each camera field of view, corresponding to different vegetation types and demarcating (1) Picea trees (vegetation type EN, for evergreen needleleaf); (2) Larix trees (vegetation type DN, for deciduous needleleaf); and (3) the mixed shrub layer (vegetation type SH). This data set consists of three sets of data files: 3-day summary product files: One file for each camera and each ROI (i.e. vegetation type), characterizing vegetation color at a 3-day time step Contains 35 files in *.csv format inside a compressed (*.zip) file. Transition date file: Estimates “greenness rising” (spring) and “greenness falling” (autumn) transition dates derived from the smoothed daily green chromatic coordinate (GCC) values, for each camera and each ROI (i.e. vegetation type) Contains one file in *.csv format Snow flag files: Indicate days with snow on trees or snow on ground for each experimental enclosure Contains two files in *.csv format, one for snow on trees and one for snow on ground This data set consists of two sets of companion files: Accompanying HTML files show the 90th quantiles of the mean GCC plotted together with transition dates for each vegetation type and plot. Contains three files in HTML format, one for each vegetation type One additional file in HTML format with the transition dates plotted for each vegetation type, by year R files for processing Phenocam files and flags. Contains five files in R file (*.R) format in one compressed (*.zip) file User Note: All imagery is posted in near-real time to the PhenoCam Project web page (http://phenocam.sr.unh.edu/), where it is publicly available. Scroll to “spruce” in the Gallery or link directly to the 29 SPRUCE cameras at https://tinyurl.com/y7z5mau7. The data reported here are based on the complete camera record from SPRUCE and supersedes the previously released phenocam datasets (see Related Data Sets). The estimated transition dates for previously released datasets may differ slightly (in most cases, by ±3 days or less), because following standard PhenoCam processing protocols (Richardson et al. 2018, Scientific Data), smoothing and interpolation, outlier removal, and transition date estimation are always conducted using the full data record.

SPRUCE Experiment, Marcell Experimental Forest, Sp↗
DOE OSTI2024

SPRUCE Vegetation Phenology in Experimental Plots from Phenocam Imagery, 2015-2023

This data set consists of PhenoCam data from the SPRUCE experiment from the beginning of whole ecosystem warming (Hanson et al. 2017) in August 2015 through March 31 of 2024, with start- and end-of-season phenological transition dates derived through the end of autumn 2023. Digital cameras, or phenocams, installed in each SPRUCE enclosure track seasonal variation in vegetation “greenness”, a proxy for vegetation phenology and associated physiological activity. Three separate regions of interest (ROIs) were defined for each camera field of view, corresponding to different vegetation types and demarcating (1) Picea trees (vegetation type EN, for evergreen needleleaf); (2) Larix trees (vegetation type DN, for deciduous needleleaf); and (3) the mixed shrub layer (vegetation type SH). This data set consists of three sets of data files: (1) 3-day summary product files: One file for each camera and each ROI (i.e. vegetation type), characterizing vegetation color at a 3-day time step • Contains 36 files in *.csv format inside a compressed (*.zip) file. (2) Transition date file: Estimates “greenness rising” (spring) and “greenness falling” (autumn) transition dates derived from the smoothed daily green chromatic coordinate (GCC) values, for each camera and each ROI (i.e. vegetation type) • Contains one file in *.csv format (3) Snow flag files: Indicate days with snow on trees or snow on ground for each experimental enclosure • Contains two files in *.csv format, one for snow on trees and one for snow on ground This data set consists of two sets of companion files: (1) Accompanying HTML files show the 90th quantiles of the mean GCC plotted together with transition dates for each vegetation type and plot. • Contains three files in HTML format, one for each vegetation type • One additional file in HTML format with the transition dates plotted for each vegetation type, by year (2) R files for processing Phenocam files and flags. • Contains five files in R file (*.R) format in one compressed (*.zip) file User Note: All imagery is posted in near-real time to the PhenoCam Project web page (http://phenocam.sr.unh.edu/), where it is publicly available. Scroll to “spruce” in the Gallery or link directly to the 29 SPRUCE cameras at https://tinyurl.com/sprucecams. The data reported here are based on the complete camera record from SPRUCE and supersedes the previously released data inclusive of the 2015-2022 data (see Related Data Sets). The estimated transition dates for previously released datasets may differ slightly (in most cases, by ±3 days or less), because following standard PhenoCam processing protocols (Richardson et al. 2018, Scientific Data), smoothing and interpolation, outlier removal, and transition date estimation are always conducted using the full data record.

Spruce and Peatland Responses Under Changing Envir↗
DOE OSTI2025

SPRUCE Vegetation Phenology in Experimental Plots from PhenoCam Imagery, 2015-2024

This data set consists of PhenoCam data from the SPRUCE experiment from the beginning of whole ecosystem warming (Hanson et al. 2017) in August 2015 through March 31 of 2025 (2015-08-24 to 2025-03-31), with start- and end-of-season phenological transition dates derived through the end of autumn 2024. Digital cameras, or phenocams, installed in each SPRUCE enclosure track seasonal variation in vegetation “greenness”, a proxy for vegetation phenology and associated physiological activity. Three separate regions of interest (ROIs) were defined for each camera field of view, corresponding to different vegetation types and demarcating (1) Picea trees (vegetation type EN, for evergreen needleleaf); (2) Larix trees (vegetation type DN, for deciduous needleleaf); and (3) the mixed shrub layer (vegetation type SH). This data set consists of three sets of data files: (1) 3-day summary product files: One file for each camera and each ROI (i.e. vegetation type), characterizing vegetation color at a 3-day time step. • Contains 36 files in *.csv format inside a compressed (*.zip) file. (2) Transition date file: Estimates “greenness rising” (spring) and “greenness falling” (autumn) transition dates derived from the smoothed daily green chromatic coordinate (GCC) values, for each camera and each ROI (i.e., vegetation type). • Contains one file in *.csv format. (3) Snow flag files: Indicate days with snow on trees or snow on ground for each experimental enclosure. • Contains two files in *.csv format, one for snow on trees and one for snow on ground. This data set consists of two sets of companion files: (1) Accompanying HTML files show the 90th quantiles of the mean GCC plotted together with transition dates for each vegetation type and plot. • Contains three files in HTML format, one for each vegetation type. • One additional file in HTML format with the transition dates plotted for each vegetation type, by year. (2) R files for processing PhenoCam files and flags. • Contains five files in R file(*.R) format and the components of the phenocamr package (Version 1.1.4) used for calculating transition dates for 2015-2024. These are contained in a compressed (*.zip) file. User Note: All imagery is posted in near-real time to the PhenoCam Project web page (https://phenocam.nau.edu), where it is publicly available. Scroll to “spruce” in the Gallery or link directly to the 29 SPRUCE cameras at https://tinyurl.com/sprucecams. This data set is based on the complete camera record from SPRUCE and supersedes all previously released PhenoCam datasets (see Related Data Sets). The estimated transition dates for previously released datasets may differ slightly (in most cases, by ±3 days or less), because following standard PhenoCam processing protocols (Richardson et al. 2018, Scientific Data), smoothing and interpolation, outlier removal, and transition date estimation are always conducted using the full data record.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Application of the Approximate Bayesian Computation Algorithm to Gamma-Ray Spectroscopy

Radioisotope identification (RIID) algorithms for gamma-ray spectroscopy aim to infer what isotopes are present and in what amounts in test items. RIID algorithms either use all energy channels in the analysis region or only energy channels in and near identified peaks. Because many RIID algorithms rely on locating peaks and estimating each peak’s net area, peak location and peak area estimation algorithms continue to be developed for gamma-ray spectroscopy. This paper shows that approximate Bayesian computation (ABC) can be effective for peak location and area estimation. Algorithms to locate peaks can be applied to raw or smoothed data, and among several smoothing options, the iterative bias reduction algorithm (IBR) is recommended; the use of IBR with ABC is shown to potentially reduce uncertainty in peak location estimation. Extracted peak locations and areas can then be used as summary statistics in a new ABC-based RIID. ABC allows for easy experimentation with candidate summary statistics such as goodness-of-fit scores and peak areas that are extracted from relatively high dimensional gamma spectra with photopeaks (1024 or more energy channels) consisting of count rates versus energy for a large number of gamma energies.

Burr, Tom↗
DOE OSTI2022

High-Quality Dry Etching of LiNbO3 Assisted by Proton Substitution through H2-Plasma Surface Treatment

The exceptional material properties of Lithium Niobate (LiNbO3) make it an excellent material platform for a wide range of RF, MEMS, phononic and photonic applications; however, nano-micro scale device concepts require high fidelity processing of LN films. Here, we reported a highly optimized processing methodology that achieves a deep etch with nearly vertical and smooth sidewalls. We demonstrated that Ti/Al/Cr stack works perfectly as a hard mask material during long plasma dry etching, where periodically pausing the etching and chemical cleaning between cycles were leveraged to avoid thermal effects and byproduct redeposition. To improve mask quality on X- and Y-cut substrates, a H2-plasma treatment was implemented to relieve surface tension by modifying the top surface atoms. Structures with etch depths as deep as 3.4 µm were obtained in our process across a range of crystallographic orientations with a smooth sidewall and perfect verticality on several crystallographic facets.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Evolution of an Asteroid Family under YORP, Yarkovsky, and Collisions

Any population of asteroids, like asteroid families, will disperse in semimajor axis due to the Yarkovsky effect. The amount of drift is modulated by the asteroid spin state evolution, which determines the balance between the diurnal and seasonal Yarkovsky forces. The asteroid’s spin state is, in turn, controlled in part by the Yarkovsky–O’Keefe–Radzievskii–Paddack (YORP) effect. The otherwise smooth evolution of an asteroid can be abruptly altered by collisions, which can cause impulsive changes in the spin state and can move the asteroid onto a different YORP track. In addition, collisions may also alter the YORP parameters by changing the superficial features and overall shape of the asteroid. Thus, the coupling between YORP and Yarkovsky is also strongly affected by the impact history of each body. To investigate this coupling, we developed a statistical code modeling the time evolution of semimajor axis under YORP–Yarkovsky coupling. It includes the contributions of NYORP (normal YORP), TYORP (tangential YORP), and collisions whose effects are deterministically calculated and not added in a statistical way. We find that both collisions and TYORP increase the dispersion of a family in semimajor axis by making the spin axis evolution less smooth and regular. We show that the evolution of a family’s structure with time is complex and collisions randomize the YORP evolution. In our test families, we do not observe the formation of a “YORP-eye” in the semimajor axis versus diameter distribution, even after a long period of time. If present, the “YORP-eye” might be a relic of an initial ejection velocity pattern of the collisional fragments.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

A Detailed Characterization of HR 8799's Debris Disk with ALMA in Band 7

The exoplanetary system of HR 8799 is one of the rare systems in which multiple planets have been directly imaged. Its architecture is strikingly similar to that of the solar system, with the four imaged giant planets surrounding a warm dust belt analogous to the Asteroid Belt, and themselves being surrounded by a cold dust belt analog to the Kuiper Belt. Previous observations of this cold belt with ALMA in Band 6 (1.3 mm) revealed its inner edge, but analyses of the data differ on its precise location. It was therefore unclear whether or not the outermost planet HR 8799 b was dynamically sculpting it. We present here new ALMA observations of this debris disk in Band 7 (340 GHz, 880 μm). These are the most detailed observations of this disk obtained so far, with a resolution of 1″ (40 au) and sensitivity of 9.8 μJy beam{sup −1}, which allowed us to recover the disk structure with high confidence. In order to constrain the disk morphology, we fit its emission using radiative transfer models combined with a Markov Chain Monte Carlo procedure. We find that this disk cannot be adequately represented by a single power law with sharp edges. It exhibits a smoothly rising inner edge and smoothly falling outer edge, with a peak in between, as expected from a disk that contains a high-eccentricity component, hence confirming previous findings. Whether this excited population and inner edge shape stem from the presence of an additional planet remains, however, an open question.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Minkowski Functionals of SDSS-III BOSS: Hints of Possible Anisotropy in the Density Field?

We present measurements of the Minkowski functionals extracted from the SDSS-III BOSS catalogs. After defining the Minkowski functionals, we describe how an unbiased reconstruction of these statistics can be obtained from a field with masked regions and survey boundaries, validating our methodology with Gaussian random fields and mock galaxy snapshot data. From the BOSS galaxy data, we generate a set of four density fields in three dimensions corresponding to the northern and southern skies of LOWZ and CMASS catalogs, smoothing over large scales (Gaussian smoothing scale of 35 Mpc) such that the field is perturbatively non-Gaussian. We extract the Minkowski functionals from each data set separately, and measure their shapes and amplitudes by fitting a Hermite polynomial expansion. For the shape parameter of the Minkowski functional curves a0, that is related to the bispectrum of the field, we find that the LOWZ-South data presents a systematically lower value of a 0 = –0.080 ± 0.040 than its northern sky counterpart a 0 = 0.032 ± 0.024. Although the significance of this discrepancy is not very high, it potentially indicates some systematics in the data or that the matter density field exhibits anisotropy at low redshift. By assuming a standard isotropic flat ΛCDM cosmology, the amplitudes of the Minkowski functionals from the combination of northern and southern sky data give the constraints Ω c h 2 n s = 0.110 ± 0.006 and 0.111 ± 0.008 for CMASS and LOWZ, respectively, which is in agreement with the Planck ΛCDM best-fit Ω c h 2 n s = 0.116 ± 0.001.

79 ASTRONOMY AND ASTROPHYSICS↗