Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DECOUPLING”

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,711 records · Page 95
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Data-Driven RANS Turbulence Closures for Forced Convection Flow in Reactor Downcomer Geometry

Recent progress in data-driven turbulence modeling has shown its potential to enhance or replace traditional equation-based Reynolds-averaged Navier-Stokes (RANS) turbulence models. Here, this work utilizes invariant neural network (NN) architectures to model Reynolds stresses and turbulent heat fluxes in forced convection flows (when the models can be decoupled). As the considered flow is statistically one dimensional, the invariant NN architecture for the Reynolds stress model reduces to the linear eddy viscosity model. To develop the data-driven models, direct numerical and RANS simulations in vertical planar channel geometry mimicking a part of the reactor downcomer are performed. Different conditions and fluids relevant to advanced reactors (sodium, lead, unitary-Prandtl-number fluid, and molten salt) constitute the training database. The models enabled accurate predictions of velocity and temperature, and compared to the baseline k–τ turbulence model with the simple gradient diffusion hypothesis, do not require tuning of the turbulent Prandtl number. The data-driven framework is implemented in the open-source graphics processing unit–accelerated spectral element solver nekRS and has shown the potential for future developments and consideration of more complex mixed convection flows.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2018

Limiting Surface Density Method Adapted to Large Arrays of Heterogeneous Shipping Packages with Nonlinear Responses

The classic Limiting Surface Density (LSD) method is an empirical calculation technique for analyzing and setting mass limits for fissile items in storage arrays. LSD is a desirable method because it can reduce or eliminate the need for lengthy detailed Monte Carlo models of storage arrays. The original (or classic) method was developed based on idealized arrays of bare spherical metal items in air-spaced cubic units in a water-reflected cubic array. In this case, the geometric and material-based surface densities were acceptably correlated by linear functions. Later updates to the method were made to allow for concrete reflection rather than water, cylindrical masses rather than spheres, different material forms, and noncubic arrays. However, in the intervening four decades since those updates, little work has been done to update the method, especially for use with contemporary highly heterogeneous shipping packages that are noncubic and stored in noncubic arrays. In this work, the LSD method is reevaluated for application to highly heterogeneous shipping packages for fissile material. The package modeled is the 9975 shipping package, currently the primary package used to store fissile material at Savannah River Site’s K-Area Complex. The package is neither cubic nor rectangular but resembles nested cylinders of stainless steel, lead, aluminum, and Celotex. The fissile content is assumed to be a cylinder of plutonium metal. The packages may be arranged in arrays with both an equal number of packages per side (package cubic) and an unequal number of packages per side (noncubic). The cubic arrangements are used to derive the 9975-specific material and geometry constants for the classic linear form LSD method. The linear form of the LSD, with noncubic array adjustment, is applied and evaluated against computational models for these packages to determine the critical unit fissile mass. Sensitivity equations are derived from the classic method, and these are also used to make projections of the critical unit fissile mass. It was discovered that the heterogeneous packages have a nonlinear surface density versus critical mass relationship compared to the acceptably linear response of bare spherical fissile masses. Methodology is developed to address the nonlinear response. In so doing, the solution to the nonlinear LSD method becomes decoupled from the critical mass of a single unit, adding to its flexibility. The ability of the method to predict changes in neutron multiplication due to perturbations in a parameter is examined to provide a basis for analyzing upset conditions. In conclusion, a full rederivation of the classic LSD method from diffusion theory is also included as this was found to be lacking in the available literature.

42 ENGINEERING↗
DOE OSTI2020

Nonlinear Elimination Applied to Radiation Diffusion

We apply a nonlinearly preconditioned, quasi-Newton framework to accelerate the numerical solution of the thermal radiative transfer (TRT) equations. This framework was inspired by the unpublished method that has existed for years in Teton, Lawrence Livermore National Laboratory’s deterministic TRT code. In this paper, we cast this iteration scheme within a formal nonlinear preconditioning framework and compare its performance against other iteration schemes in the framework. With proper choices of iteration controls for the various levels of the solver, we can recover the standard linearized one-step method, a full nonlinear Newton scheme, as well as the method in Teton. In brief, the nonlinear preconditioning TRT scheme formally eliminates the material temperature equation from the nonlinear system in a nonlinear analog of a Schur complement. This nonlinear elimination step involves solving a decoupled nonlinear equation for each spatial degree of freedom and is therefore inexpensive. By applying a quasi-Newton iteration scheme on the new system, we obtain a three-level iteration scheme that is at least as efficient as commonly used TRT schemes. The new method allows full convergence to the nonlinear backward Euler time-discretized system, increasing accuracy and robustness, while using a similar number of linear iterations as the more common linearized one-step methods Eq. (4).

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2014

Evaluating the Accuracy and Precision of Multiple Abundance Estimators Using State-Space Models: A Case Study for a Threatened Population of Chinook Salmon in Johnson Creek, Idaho

Abstract Over the last century, Chinook Salmon Oncorhynchus tshawytscha populations in the Pacific Northwest have experienced dramatic declines, leading to many of them being listed as threatened under the Endangered Species Act. The abundance of these threatened populations relative to the thresholds for delisting remains the primary metric for assessing recovery, yet determining the true population abundances from multiple survey types with unknown levels of accuracy and precision remains difficult. The abundance of the spring–summer Chinook Salmon population in Johnson Creek, Idaho, has been measured using a mark–recapture survey and three different redd count surveys (RCSs) that vary temporally and spatially. Using a state-space model, we determined the accuracy and precision of each survey type by decoupling the observation error of the survey from the process error describing the annual variability in the true population abundance. We then extended the results of the model to determine the risk of managers’ incorrectly delisting the population (a type I error) or incorrectly keeping it listed (a type II error). Finally, we show that salmon managers with data-limited populations (primarily those with only single-pass index RCSs) might use the results of our risk analysis to determine whether expanding survey efforts to minimize management risks is appropriate when they are confronted with dwindling financial resources. For example, we determined that although both the multiple-pass extended RCS (CV = 0.06) and mark–recapture surveys (CV = 0.14) provide unbiased estimates of salmon abundance in Johnson Creek, the mark–recapture study can have annual costs that are 30–100 times greater. Managers may determine that directing research funds toward acquiring information unique to weir-based mark–recapture surveys (i.e., migration timing, good genetics samples, etc.) may not be justified for all populations. Received October 7, 2013; accepted March 5, 2014

Chasco, Brandon E.↗
DOE OSTI2017

Vacuum System and Modeling for the Materials Plasma Exposure Experiment

Understanding the science of plasma-material interactions (PMI) is essential for the future development of fusion facilities. The design of divertors and first walls for the next generation of long-pulse fusion facilities, such as a Fusion Nuclear Science Facility (FNSF) or a DEMO, requires significant PMI research and development. In order to meet this need, a new linear plasma facility, the Materials Plasma Exposure Experiment (MPEX) is proposed, which will produce divertor relevant plasma conditions for these next generation facilities. The device will be capable of handling low activation irradiated samples and be able to remove and replace samples without breaking vacuum. A Target Exchange Chamber (TEC) which can be disconnected from the high field environment in order to perform in-situ diagnostics is planned for the facility as well. The vacuum system for MPEX must be carefully designed in order to meet the requirements of the different heating systems, and to provide conditions at the target similar to those expected in a divertor. An automated coupling-decoupling (“autocoupler”) system is designed to create a high vacuum seal, and will allow the TEC to be disconnected without breaking vacuum in either the TEC or the primary plasma materials interaction chamber. This autocoupler, which can be actuated remotely in the presence of the high magnetic fields, has been designed and prototyped, and shows robustness in a variety of conditions. The vacuum system has been modeled using a simplified finite element analysis, and indicates that the design goals for the pressures in key regions of the facility are achievable.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Final Design of the Material Plasma Exposure eXperiment

The Material Plasma Exposure eXperiment (MPEX) has completed its design phase. MPEX will be a unique facility to investigate plasma material interactions (PMIs) under fusion prototypic divertor conditions in steady state. This includes plasma exposure conditions expected in a fusion reactor divertor. Materials to be investigated will include solids, liquids, and neutron pre-irradiated materials. Electron and ion heating will allow for a large operational domain ranging from erosion-dominated PMI conditions to deposition-dominated PMI conditions. An overview of the final design for all MPEX systems is given. Here, in particular, it is shown how mission-driven project requirements have led to detailed design choices with innovative solutions. Examples are the water-cooled helicon window, the electron cyclotron heating launcher, the target holder and manipulator, and the autonomous decouplers.

Fusion reactor↗
DOE OSTI2013

CAPTURE OF PLANETESIMALS BY GAS DRAG FROM CIRCUMPLANETARY DISKS

Growing giant planets have circumplanetary disks around them in the late stage of their formation if their mass is sufficiently large. We examine capture of relatively large planetesimals that are decoupled from the gas inflow, due to gas drag from a circumplanetary disk of a growing giant planet. Assuming that the structure of the circumplanetary disk is axisymmetric, and solving the three-body problem including gas drag, we perform analytic and numerical calculations for capture of planetesimals. When planetesimal random velocity is small, planetesimals approaching in the retrograde direction are more easily captured, owing to their larger velocity relative to the gas. Planetesimals with large orbital inclinations interact with the disk for a short period of time and show lower capture rates. The effect of ablation on capture rates seems insignificant, although mass loss due to ablation would be significant in the case of high random velocity. We also examine the effect of non-uniform radial distribution of planetesimals in the protoplanetary disk due to gap opening by the planet. When the random velocity of planetesimals is small, the planetesimal capture rate decreases rapidly as the half width of the gap in the planetesimal disk increases from two planetary Hill radii to three planetary Hill radii; planetesimals with low random velocities cannot approach the planet in the case of a sufficiently wide gap. Our results show that the radial distribution and random velocity of planetesimals in the protoplanetary disk are essentially important for the understanding of capture of planetesimals by circumplanetary disks.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

THE HCN-WATER RATIO IN THE PLANET FORMATION REGION OF DISKS

We find a trend between the mid-infrared HCN/H{sub 2}O flux ratio and submillimeter disk mass among T Tauri stars in Taurus. While it may seem puzzling that the molecular emission properties of the inner disk (<few AU) are related to the properties of the outer disk (beyond {approx}20 AU) probed by the submillimeter continuum, an interesting possible interpretation is that the trend is a result of planetesimal and protoplanet formation. Because objects this large are decoupled from the accretion flow, when they form, they can lock up water (and oxygen) beyond the snow line, thereby enhancing the C/O ratio in the inner disk and altering the molecular abundances there. We discuss the assumptions that underlie this interpretation, a possible alternative explanation, and related open questions that motivate future work. Whatever its origin, understanding the meaning of the relation between the HCN/H{sub 2}O ratio and disk mass is of interest as trends like this among T Tauri disk properties are relatively rare.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

UBIQUITOUS NON-THERMALS IN ASTROPHYSICAL PLASMAS: RESTATING THE DIFFICULTY OF MAINTAINING MAXWELLIANS

This paper outlines the rather narrow conditions on a radiatively decoupled plasma where a Maxwell-Boltzmann (MB) distribution can be assumed with confidence. The complementary non-thermal distribution with non-perturbative kurtosis is argued to have a much broader purview than has previously been accepted. These conditions are expressed in terms of the electron Knudsen number, K{sub e} , the ratio of the electron mean free path to the scale length of electron pressure. Rather generally, f(v < v{sub 2}(K{sub e} )) will be Gaussian, so that MB atomic or wave particle effects controlled by speeds v < v{sub 2} {identical_to} w(15/8K{sub e} ){sup 1/4} will remain defensible, where w is the most probable speed. The sufficient condition for Spitzer-Braginskii plasma fluid closure at the energy equation requires globally K{sub e} (s) {<=} 0.01; this global condition pertains to the maximum value of K{sub e} along the arc length s of the magnetic field (to its extremities) provided that contiguous plasma remains uncoupled from the radiation field. The non-thermal regime K{sub e} > 0.01 is common in all main-sequence stellar atmospheres above approximately 0.05 stellar radii from the surface. The entire solar corona and wind are included in this regime where non-thermal distributions with kurtosis are shown to be ubiquitous, heat flux is not well modeled by Spitzer-Braginskii closure, and fluid modeling is qualitative at best.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

SPECTROSCOPIC OBSERVATIONS OF SN 2012fr: A LUMINOUS, NORMAL TYPE Ia SUPERNOVA WITH EARLY HIGH-VELOCITY FEATURES AND A LATE VELOCITY PLATEAU

We present 65 optical spectra of the Type Ia SN 2012fr, 33 of which were obtained before maximum light. At early times, SN 2012fr shows clear evidence of a high-velocity feature (HVF) in the Si II {lambda}6355 line that can be cleanly decoupled from the lower velocity ''photospheric'' component. This Si II {lambda}6355 HVF fades by phase -5; subsequently, the photospheric component exhibits a very narrow velocity width and remains at a nearly constant velocity of {approx}12,000 km s{sup -1} until at least five weeks after maximum brightness. The Ca II infrared triplet exhibits similar evidence for both a photospheric component at v Almost-Equal-To 12,000 km s{sup -1} with narrow line width and long velocity plateau, as well as an HVF beginning at v Almost-Equal-To 31,000 km s{sup -1} two weeks before maximum. SN 2012fr resides on the border between the ''shallow silicon'' and ''core-normal'' subclasses in the Branch et al. classification scheme, and on the border between normal and high-velocity Type Ia supernovae (SNe Ia) in the Wang et al. system. Though it is a clear member of the ''low velocity gradient'' group of SNe Ia and exhibits a very slow light-curve decline, it shows key dissimilarities with the overluminous SN 1991T or SN 1999aa subclasses of SNe Ia. SN 2012fr represents a well-observed SN Ia at the luminous end of the normal SN Ia distribution and a key transitional event between nominal spectroscopic subclasses of SNe Ia.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

AN IN-DEPTH VIEW OF THE MID-INFRARED PROPERTIES OF POINT SOURCES AND THE DIFFUSE ISM IN THE SMC GIANT H II REGION, N66

The focus of this work is to study mid-infrared point sources and the diffuse interstellar medium (ISM) in the low-metallicity ({approx}0.2 Z{sub Sun }) giant H II region N66 in order to determine properties that may shed light on star formation in these conditions. Using the Spitzer Space Telescope's Infrared Spectrograph, we study polycyclic aromatic hydrocarbon (PAH), dust continuum, silicate, and ionic line emission from 14 targeted infrared point sources as well as spectra of the diffuse ISM that is representative of both the photodissociation regions (PDRs) and the H II regions. Among the point source spectra, we spectroscopically confirm that the brightest mid-infrared point source is a massive embedded young stellar object, we detect silicates in emission associated with two young stellar clusters, and we see spectral features of a known B[e] star that are commonly associated with Herbig Be stars. In the diffuse ISM, we provide additional evidence that the very small grain population is being photodestroyed in the hard radiation field. The 11.3 {mu}m PAH complex emission exhibits an unexplained centroid shift in both the point source and ISM spectra that should be investigated at higher signal-to-noise and resolution. Unlike studies of other regions, the 6.2 {mu}m and 7.7 {mu}m band fluxes are decoupled; the data points cover a large range of I{sub 7.7}/I{sub 11.3} PAH ratio values within a narrow band of I{sub 6.2}/I{sub 11.3} ratio values. Furthermore, there is a spread in PAH ionization, being more neutral in the dense PDR where the radiation field is relatively soft, but ionized in the diffuse ISM/PDR. By contrast, the PAH size distribution appears to be independent of local ionization state. Important to unresolved studies of extragalactic low-metallicity star-forming regions, we find that emission from the infrared-bright point sources accounts for only 20%-35% of the PAH emission from the entire region. These results make a comparative data set to other star-forming regions with similarly hard and strong radiation fields.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

HOT ELECTROMAGNETIC OUTFLOWS. II. JET BREAKOUT

We consider the interaction between radiation, matter, and a magnetic field in a compact, relativistic jet. The entrained matter accelerates outward as the jet breaks out of a star or other confining medium. In some circumstances, such as gamma-ray bursts (GRBs), the magnetization of the jet is greatly reduced by an advected radiation field while the jet is optically thick to scattering. Where magnetic flux surfaces diverge rapidly, a strong outward Lorentz force develops and radiation and matter begin to decouple. The increase in magnetization is coupled to a rapid growth in Lorentz factor. We take two approaches to this problem. The first examines the flow outside the fast magnetosonic critical surface, and calculates the flow speed and the angular distribution of the radiation field over a range of scattering depths. The second considers the flow structure on both sides of the critical surface in the optically thin regime, using a relaxation method. In both approaches, we find how the terminal Lorentz factor and radial profile of the outflow depend on the radiation intensity and optical depth at breakout. The effect of bulk Compton scattering on the radiation spectrum is calculated by a Monte Carlo method, while neglecting the effects of internal dissipation. The peak of the scattered spectrum sits near the seed peak if radiation pressure dominates the acceleration, but is pushed to a higher frequency if the Lorentz force dominates. The unscattered seed radiation can form a distinct, low-frequency component of the spectrum, especially if the magnetic Poynting flux dominates.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

PULSAR STATE SWITCHING FROM MARKOV TRANSITIONS AND STOCHASTIC RESONANCE

Markov processes are shown to be consistent with metastable states seen in pulsar phenomena, including intensity nulling, pulse-shape mode changes, subpulse drift rates, spin-down rates, and X-ray emission, based on the typically broad and monotonic distributions of state lifetimes. Markovianity implies a nonlinear magnetospheric system in which state changes occur stochastically, corresponding to transitions between local minima in an effective potential. State durations (though not transition times) are thus largely decoupled from the characteristic timescales of various magnetospheric processes. Dyadic states are common but some objects show at least four states with some transitions forbidden. Another case is the long-term intermittent pulsar B1931+24 that has binary radio-emission and torque states with wide, but non-monotonic duration distributions. It also shows a quasi-period of 38 ± 5 days in a 13 yr time sequence, suggesting stochastic resonance in a Markov system with a forcing function that could be strictly periodic or quasi-periodic. Nonlinear phenomena are associated with time-dependent activity in the acceleration region near each magnetic polar cap. The polar-cap diode is altered by feedback from the outer magnetosphere and by return currents from the equatorial region outside the light cylinder that may also cause the neutron star to episodically charge and discharge. Orbital perturbations of a disk or current sheet provide a natural periodicity for the forcing function in the stochastic-resonance interpretation of B1931+24. Disk dynamics may introduce additional timescales in observed phenomena. Future work can test the Markov interpretation, identify which pulsar types have a propensity for state changes, and clarify the role of selection effects.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Global numerical modeling of energetic proton acceleration in a coronal mass ejection traveling through the solar corona

The acceleration of protons and electrons to high (sometimes GeV/nucleon) energies by solar phenomena is a key component of space weather. These solar energetic particle (SEP) events can damage spacecraft and communications, as well as present radiation hazards to humans. In-depth particle acceleration simulations have been performed for idealized magnetic fields for diffusive acceleration and particle propagation, and at the same time the quality of MHD simulations of coronal mass ejections (CMEs) has improved significantly. However, to date these two pieces of the same puzzle have remained largely decoupled. Such structures may contain not just a shock but also sizable sheath and pileup compression regions behind it, and may vary considerably with longitude and latitude based on the underlying coronal conditions. In this work, we have coupled results from a detailed global three-dimensional MHD time-dependent CME simulation to a global proton acceleration and transport model, in order to study time-dependent effects of SEP acceleration between 1.8 and 8 solar radii in the 2005 May 13 CME. We find that the source population is accelerated to at least 100 MeV, with distributions enhanced up to six orders of magnitude. Acceleration efficiency varies strongly along field lines probing different regions of the dynamically evolving CME, whose dynamics is influenced by the large-scale coronal magnetic field structure. We observe strong acceleration in sheath regions immediately behind the shock.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Transformation of a Virgo Cluster dwarf irregular galaxy by ram pressure stripping: IC3418 and its fireballs

We present optical imaging and spectroscopy and H I imaging of the Virgo Cluster galaxy IC 3418, which is likely a 'smoking gun' example of the transformation of a dwarf irregular into a dwarf elliptical galaxy by ram pressure stripping. IC 3418 has a spectacular 17 kpc length UV-bright tail comprised of knots, head-tail, and linear stellar features. The only Hα emission arises from a few H II regions in the tail, the brightest of which are at the heads of head-tail UV sources whose tails point toward the galaxy ('fireballs'). Several of the elongated tail sources have Hα peaks outwardly offset by ∼80-150 pc from the UV peaks, suggesting that gas clumps continue to accelerate through ram pressure, leaving behind streams of newly formed stars which have decoupled from the gas. Absorption line strengths, measured from Keck DEIMOS spectra, together with UV colors, show star formation stopped 300 ± 100 Myr ago in the main body, and a strong starburst occurred prior to quenching. While neither Hα nor H I emission are detected in the main body of the galaxy, we have detected 4 × 10{sup 7} M {sub ☉} of H I from the tail with the Very Large Array. The velocities of tail H II regions, measured from Keck LRIS spectra, extend only a small fraction of the way to the cluster velocity, suggesting that star formation does not happen in more distant parts of the tail. Stars in the outer tail have velocities exceeding the escape speed, but some in the inner tail should fall back into the galaxy, forming halo streams.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Plasma physical parameters along coronal-mass-ejection-driven shocks. I. Ultraviolet and white-light observations

In this work, UV and white-light (WL) coronagraphic data are combined to derive the full set of plasma physical parameters along the front of a shock driven by a coronal mass ejection. Pre-shock plasma density, shock compression ratio, speed, and inclination angle are estimated from WL data, while pre-shock plasma temperature and outflow velocity are derived from UV data. The Rankine-Hugoniot (RH) equations for the general case of an oblique shock are then applied at three points along the front located between 2.2 and 2.6 R {sub ☉} at the shock nose and at the two flanks. Stronger field deflection (by ∼46°), plasma compression (factor ∼2.7), and heating (factor ∼12) occur at the nose, while heating at the flanks is more moderate (factor 1.5-3.0). Starting from a pre-shock corona where protons and electrons have about the same temperature (T{sub p} ∼ T{sub e} ∼ 1.5 × 10{sup 6} K), temperature increases derived with RH equations could better represent the proton heating (by dissipation across the shock), while the temperature increase implied by adiabatic compression (factor ∼2 at the nose, ∼1.2-1.5 at the flanks) could be more representative of electron heating: the transit of the shock causes a decoupling between electron and proton temperatures. Derived magnetic field vector rotations imply a draping of field lines around the expanding flux rope. The shock turns out to be super-critical (sub-critical) at the nose (at the flanks), where derived post-shock plasma parameters can be very well approximated with those derived by assuming a parallel (perpendicular) shock.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Dynamical model for the toroidal sporadic meteors

More than a decade of radar operations by the Canadian Meteor Orbit Radar have allowed both young and moderately old streams to be distinguished from the dispersed sporadic background component. The latter has been categorized according to broad radiant regions visible to Earth-based observers into three broad classes: the helion and anti-helion source, the north and south apex sources, and the north and south toroidal sources (and a related arc structure). The first two are populated mainly by dust released from Jupiter-family comets and new comets. Proper modeling of the toroidal sources has not to date been accomplished. Here, we develop a steady-state model for the toroidal source of the sporadic meteoroid complex, compare our model with the available radar measurements, and investigate a contribution of dust particles from our model to the whole population of sporadic meteoroids. We find that the long-term stable part of the toroidal particles is mainly fed by dust released by Halley type (long period) comets (HTCs). Our synthetic model reproduces most of the observed features of the toroidal particles, including the most troublesome low-eccentricity component, which is due to a combination of two effects: particles' ability to decouple from Jupiter and circularize by the Poynting-Robertson effect, and large collision probability for orbits similar to that of the Earth. Our calibrated model also allows us to estimate the total mass of the HTC-released dust in space and check the flux necessary to maintain the cloud in a steady state.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Accreting planets as dust dams in 'transition' disks

We investigate under what circumstances an embedded planet in a protoplanetary disk may sculpt the dust distribution such that it observationally presents as a 'transition' disk. We concern ourselves with 'transition' disks that have large holes (≳ 10 AU) and high accretion rates (∼10{sup –9}-10{sup –8} M {sub ☉} yr{sup –1}), particularly, those disks which photoevaporative models struggle to explain. Adopting the observed accretion rates in 'transition' disks, we find that the accretion luminosity from the forming planet is significant, and can dominate over the stellar luminosity at the gap edge. This planetary accretion luminosity can apply a significant radiation pressure to small (s ≲ 1 μm) dust particles provided they are suitably decoupled from the gas. Secular evolution calculations that account for the evolution of the gas and dust components in a disk with an embedded, accreting planet, show that only with the addition of the radiation pressure can we explain the full observed characteristics of a 'transition' disk (NIR dip in the spectral energy distribution (SED), millimeter cavity, and high accretion rate). At suitably high planet masses (≳ 3-4 M{sub J} ), radiation pressure from the accreting planet is able to hold back the small dust particles, producing a heavily dust-depleted inner disk that is optically thin to infrared radiation. The planet-disk system will present as a 'transition' disk with a dip in the SED only when the planet mass and planetary accretion rate are high enough. At other times, it will present as a disk with a primordial SED, but with a cavity in the millimeter, as observed in a handful of protoplanetary disks.

79 ASTRONOMY AND ASTROPHYSICS↗