Search

This paper presents a method of performing topology optimisation of transient heat conduction problems using the parallel-in-time method Parareal. To accommodate the adjoint analysis, the Parareal method was modified to store intermediate time steps. Preliminary tests revealed that Parareal requires many iterations to achieve accurate results and, thus, achieves no appreciable speedup. To mitigate this, a one-shot approach was used, where the time history is iteratively refined over the optimisation process. The method estimates objectives and sensitivities by introducing cumulative objectives and sensitivities and solving for these using a single iteration of Parareal, after which it updates the design using the Method of Moving Asymptotes. The resulting method was applied to a test problem where a power mean of the temperature was minimised. It achieved a peak speedup relative to a sequential reference method of $5\times$ using 16 threads. The resulting designs were similar to the one found by the reference method, both in terms of objective values and qualitative appearance. The one-shot Parareal method was compared to the Parallel Local-in-Time method of topology optimisation. This revealed that the Parallel Local-in-Time method was unstable for the considered test problem, but it achieved a peak speedup of $12\times$ using 32 threads. It was determined that the dominant bottleneck in the one-shot Parareal method was the time spent on computing coarse propagators., Comment: 26 pages, 14 figures. Submitted to SIAM Journal of Scientific Computing

We have conducted extremely ultra-deep pencil beam observations for new satellites around both Uranus and Neptune. Tens of images on several different nights in 2021, 2022 and 2023 were obtained and shifted and added together to reach as faint as 26.9 and 27.2 magnitudes in the r-band around Uranus and Neptune, respectively. One new moon of Uranus, S/2023 U1, and two new moons of Neptune, S/2021 N1 and S/2002 N5, were found. S/2023 U1 was 26.6 mags, is about 7 km in diameter and has a distant, eccentric and inclined retrograde orbit similar to Caliban and Stephano, implying these satellites are fragments from a once larger parent satellite. S/2023 U1 almost completely overlaps Stephano in orbital phase space. S/2021 N1 was 26.9 mags, about 14 km in size and has a retrograde orbit similar to Neso and Psamathe, indicating they are a dynamical family. We find S/2021 N1 is in a Kozai-Lidov orbital resonance. S/2002 N5 was 25.9 mags, is about 23 km in size and it makes a family of distant prograde satellites with Sao and Laomedeia. All three new moons show for the first time dynamical groups of moons exist around both Uranus and Neptune. The creation of these groups likely produced dust that could be the source of red material seen on the leading hemispheres of some larger inner satellites like Titania, Oberon and Umbriel. We also detected all known outer moons of Uranus and Neptune on multiple nights. This survey mostly completes the outer satellites of Uranus to about 8 km and Neptune to about 14 km in diameter. The size distributions of satellite dynamical families around the giant planets shows a strong steepening in the power law size distribution smaller than 5 km in diameter. The satellites of a family become much more common smaller than 5 km and their size distribution is consistent with a collisional break-up of a once larger parent satellite., Comment: Accepted Astronomical Journal

Topology optimization is used to systematically design contact-aided thermo-mechanical regulators, i.e. components whose effective thermal conductivity is tunable by mechanical deformation and contact. The thermo-mechanical interactions are modeled using a fully coupled non-linear thermo-mechanical finite element framework. To obtain the intricate heat transfer response, the components leverage self-contact, which is modeled using a third medium contact method. The effective heat transfer properties of the regulators are tuned by solving a topology optimization problem using a traditional gradient based algorithm. Several designs of thermo-mechanical regulators in the form of switches, diodes and triodes are presented.

Prion-like proteins play crucial parts in biological processes in organisms ranging from yeast to humans. For instance, many neurodegenerative diseases are believed to be caused by the production of prion-like proteins in neural tissue. As such, understanding the dynamics of prion-like protein production is a vital step toward treating neurodegenerative disease. Mathematical models of prion-like protein dynamics show great promise as a tool for predicting disease trajectories and devising better treatment strategies for prion-related diseases. Herein, we investigate a generic model for prion-like dynamics consisting of a class of non-linear ordinary differential equations (ODEs), establishing constraints through a linear stability analysis that enforce the expected properties of mammalian prion-like toxicity. Furthermore, we identify that prion toxicity evolves through three distinct phases for which we provide analytical descriptions using perturbation analyses. Specifically, prion-toxicity is initially characterised by the healthy phase, where the dynamics are dominated by the healthy form of prions, thereafter the system enters the mixed phase, where both healthy and toxic prions interact, and lastly, the system enters the toxic phase, where toxic prions dominate, and we refer to these phases as HeMiTo-dynamics. These findings hold the potential to aid researchers in developing precise mathematical models for prion-like dynamics, enabling them to better understand underlying mechanisms and devise effective treatments for prion-related diseases., Comment: 25 pages and 4 figures

Using the absolute detection calibration and abundant detections of the OSSOS (Outer Solar System Origin Survey) project, we provide population measurements for the main Kuiper Belt. For absolute magnitude $H_r<8.3$, there are 30,000 non-resonant main-belt objects, with twice as many hot-component objects than cold, and with total mass of 0.014 $M_\Earth$, only 1/7 of which is in the cold belt (assuming a cold-object albedo about half that of hot component objects). We show that transneptunian objects with $5.5 < H_r < 8.3$ (rough diameters 400--100~km) have indistinguishable absolute magnitude (size) distributions, regardless of being in the cold classical Kuiper belt (thought to be primordial) or the `hot' population (believed to be implanted after having been formed elsewhere). We discuss how this result was not apparent in previous examinations of the size distribution due to the complications of fitting assumed power-law functional forms to the detections at differing depths. This shared size distribution is surprising in light of the common paradigm that the hot population planetesimals formed in a higher density environment much closer to the Sun, in an environment that also (probably later) formed larger (dwarf planet and bigger) objects. If this paradigm is correct, % << BG added clause our result implies that planetesimal formation was relatively insensitive to the local disk conditions and that the subsequent planet-building process in the hot population did not modify the shape of the planetesimal size distribution in this 50--300~km range., Comment: 27 pages, 3 figures

This article presents a detailed introduction to density-based topology optimisation of fluid flow problems. The goal is to allow new students and researchers to quickly get started in the research area and to skip many of the initial steps, often consuming unnecessarily long time from the scientific advancement of the field. This is achieved by providing a step-by-step guide to the components necessary to understand and implement the theory, as well as extending the supplied MATLAB code. The continuous design representation used and how it is connected to the Brinkman penalty approach, for simulating an immersed solid in a fluid domain, is illustrated. The different interpretations of the Brinkman penalty term and how to chose the penalty parameters are explained. The accuracy of the Brinkman penalty approach is analysed through parametric simulations of a reference geometry. The chosen finite element formulation and the solution method is explained. The minimum dissipated energy optimisation problem is defined and how to solve it using an optimality criteria solver and a continuation scheme is discussed. The included MATLAB implementation is documented, with details on the mesh, pre-processing, optimisation and post-processing. The code has two benchmark examples implemented and the application of the code to these is reviewed. Subsequently, several modifications to the code for more complicated examples are presented through provided code modifications and explanations. Lastly, the computational performance of the code is examined through studies of the computational time and memory usage, along with recommendations to decrease computational time through approximations., Comment: 39 pages, 28 figures

There have been 77 TNOs discovered to be librating in the distant transneptunian resonances (beyond the 2:1 resonance, at semimajor axes greater than 47.7~AU) in four well-characterized surveys: the Outer Solar System Origins Survey (OSSOS) and three similar prior surveys. Here we use the OSSOS Survey Simulator to measure their intrinsic orbital distributions using an empirical parameterized model. Because many of the resonances had only one or very few detections, $j$:$k$ resonant objects were grouped by $k$ in order to have a better basis for comparison between models and reality. We also use the Survey Simulator to constrain their absolute populations, finding that they are much larger than predicted by any published Neptune migration model to date; we also find population ratios that are inconsistent with published models, presenting a challenge for future Kuiper Belt emplacement models. The estimated population ratios between these resonances are largely consistent with scattering-sticking predictions, though further discoveries of resonant TNOs with high-precision orbits will be needed to determine whether scattering-sticking can explain the entire distant resonant population or not., Comment: Accepted for publication in Planetary Sciences Journal (PSJ)

Dynamically excited objects within the Kuiper belt show a bimodal distribution in their surface colors, and these differing surface colors may be a tracer of where these objects formed. In this work we explore radial color distributions in the primordial planetesimal disk and implications for the positions of ice line/color transitions within the Kuiper belt's progenitor populations. We combine a full dynamical model of the Kuiper belt's evolution due to Neptune's migration with precise surface colors measured by the Colours of the Outer Solar System Origins Survey in order to examine the true color ratios within the Kuiper belt and the ice lines within the primordial disk. We investigate the position of a dominant, surface color changing ice-line, with two possible surface color layouts within the initial disk; (1) inner neutral surfaces and outer red, and (2) inner red surfaces and outer neutral. We performed simulations with a primordial disk that truncates at 30 au. By radially stepping the color transition out through 0.5 au intervals we show that both disk configurations are consistent with the observed color fraction. For an inner neutral, outer red primordial disk we find that the color transition can be at $28^{+2}_{-3}$ au at a 95% confidence level. For an inner red, outer neutral primordial disk the color transition can be at $27^{+3}_{-3}$ au at a 95% confidence level., Comment: 18 pages, 8 figures, accepted for publication in PSJ

We present the discovery of 2013 VZ70, the first known horseshoe coorbital companion of Saturn. Observed by the Outer Solar System Origins Survey (OSSOS) for 4.5 years, the orbit of 2013 VZ70 is determined to high precision, revealing that it currently is in `horseshoe' libration with the planet. This coorbital motion will last at least thousands of years but ends ~10 kyr from now; 2013 VZ70 is thus another example of the already-known `transient coorbital' populations of the giant planets, with this being the first known prograde example for Saturn (temporary retrograde coorbitals are known for Jupiter and Saturn). We present a theoretical steady state model of the scattering population of trans-Neptunian origin in the giant planet region (2--34 au), including the temporary coorbital populations of the four giant planets. We expose this model to observational biases using survey simulations in order to compare the model to the real detections made by a set of well-characterized outer Solar System surveys. While the observed number of coorbitals relative to the scattering population is higher than predicted, we show that the number of observed transient coorbitals of each giant planet relative to each other is consistent with a transneptunian source., Comment: 27 pages (double-line spaced manuscript format), 4 figures, 3 tables. Published in The Planetary Science Journal, 2:212 (11pp), 2021 October. Open Access journal article available at https://iopscience.iop.org/article/10.3847/PSJ/ac1c6b

This paper revisits the origin of topology optimisation for fluid flow problems, namely the Poiseuille-based frictional resistance term used to parametrise regions of solid and fluid. The traditional model only works for true topology optimisation, where it is used to approximate solid regions as areas with very small channel height and, thus, very high frictional resistance. It will be shown that if the channel height is allowed to vary continuously and/or the minimum channel height is relatively large and/or meaning is attributed to intermediate design field values, then the predictions of the traditional model are wrong. To remedy this problem, this work introduces an augmentation of the mass conservation equation to allow for continuously varying channel heights. The proposed planar model accurately describes fully-developed flow between two plates of varying channel height. It allows for a significant reduction in the number of degrees-of-freedom, while generally ensuring a high accuracy for low-to-moderate Reynolds numbers in the laminar regime. The accuracy and limitations of both the traditional and proposed models are explored using in-depth parametric studies. The proposed model is used to optimise the height of the fluid channel between two parallel plates and, thus, the topography of the plates for a flow distribution problem. Lastly, it is observed that the proposed model actually produces better topological designs than the traditional model when applied to the topology optimisation of a flow manifold.

Rotation periods of 53 small (diameters $2 < D < 40$ km) Jupiter Trojans (JTs) were derived using the high-cadence light curves obtained by the FOSSIL phase I survey, a Subaru/Hyper Suprime-Cam intensive program. These are the first reported periods measured for JTs with $D < 10$ km. We found a lower limit of the rotation period near 4 hr, instead of the previously published result of 5 hr (Ryan et al. 2017; Szabo et al. 2017, 2020) found for larger JTs. Assuming a rubble-pile structure for JTs, a bulk density of 0.9 gcm$^{-3}$ is required to withstand this spin rate limit, consistent with the value $0.8-1.0$ gcm$^{-3}$ (Marchis et al. 2006; Mueller et al. 2010; Buie et al. 2015; Berthier et al. 2020) derived from the binary JT system, (617) Patroclus-Menoetius system., Comment: Accepted by PSJ on July 9th, 2021. 15 pages, 7 figures, and 3 tables

The cold main classical Kuiper Belt consists of those non-resonant small solar system bodies with low orbital inclinations and orbital semi-major axes between 42.4 and 47.7 au. These objects likely formed \textit{in situ} and the population has experienced minimal collisional modification since formation. Using the Outer Solar System Origins Survey (OSSOS) ensemble sample and characterization, combined with constraints from deeper surveys and supported by evidence from the Minor Planet Center catalog and the Deep Ecliptic Survey, we determine the absolute magnitude $H_r$ distribution of the cold classical belt from $H_r\simeq5$ to 12 (roughly diameters of 400 km to 20 km). We conclude that the cold population's $H_r$ distribution exhibits an exponential cutoff at large sizes. Exponential cutoffs at large sizes are not a natural outcome of pair-wise particle accretion but exponentially tapered power-law size distributions are a feature of numerical simulations of planetesimal formation via a streaming instability. Our observation of an exponential cutoff agrees with previous observational inferences that no large objects ($D \gtrsim 400$ km) exist in the cold population. We note that the asymptotic slope of the $H_r$ distribution is consistent with $\alpha \sim 0.4$ and this asymptotic slope is also found in streaming instability modelling of planetesimal formation and is thus not necessarily associated with achieving collisional equilibrium. Studies of the transneptunian region are providing the parameters that will enable future streaming-instability studies to determine the initial conditions of planetesimal formation in the $\approx$45 au region of the Sun's protoplanetary disk., Comment: 20 pages, 2 figures, accepted in AAS Journals

Here, we present results on the intrinsic collision probabilities, $ P_I$, and range of collision speeds, $V_I$, as a function of the heliocentric distance, $r$, in the trans-Neptunian region. The collision speed is one of the parameters, that serves as a proxy to a collisional outcome e.g., complete disruption and scattering of fragments, or formation of crater, where both processes are directly related to the impact energy. We utilize an improved and de-biased model of the trans-Neptunian object (TNO) region from the "Outer Solar System Origins Survey" (OSSOS). It provides a well-defined orbital distribution model of TNOs, based on multiple opposition observations of more than 1000 bodies. In this work we compute collisional probabilities for the OSSOS models of the main classical, resonant, detached+outer and scattering TNO populations. The intrinsic collision probabilities and collision speeds are computed using the \"{O}pik's approach, as revised and modified by Wetherill for non-circular and inclined orbits. The calculations are carried out for each of the dynamical TNO groups, allowing for inter-population collisions as well as collisions within each TNO population, resulting in 28 combinations in total. Our results indicate that collisions in the trans-Neptunian region are possible over a wide range in ($r, V_I$) phase space. Although collisions are calculated to happen within $r\sim 20 - 200$~AU and $V_I \sim 0.1$~km/s to as high as $V_I\sim9$~km/s, most of the collisions are likely to happen at low relative velocities $V_I<1$~km/s and are dominated by the main classical belt., Comment: 13 pages, 6 figures

The cold classical Kuiper Belt Objects (KBOs) possess a high, $\gtrsim30\%$ binary fraction. Widely separated and dynamically fragile, these binary systems have been useful in tracing the origins of KBOs. A new class of binaries was recently identified by their colours. The so-called blue binaries are unanimously members of the less red compositional class, and exhibit a 100% binary fraction. They appear to be push-out survivors, emplaced in the classical region during Neptune's phases of outward migration. The presence of these binary systems implies that the majority of objects that formed near the cold classical region formed as binaries. Here we present new optical colour measurements of cold classical KBOs from the Colours of the Outer Solar System Origins Survey, including colours of a blue binary discovered by the Solar System Origins Legacy Survey -- 2015 RJ277. The increased size of the colours sample has resulted in order-of-magnitude decrease in the probability that the binaries and singles sample share the same colour distribution. From the Anderson-Darling statistic, this probability is only a 0.3%, while it is only 0.002% when utilizing the difference of means statistic. We find a hint that the blue binaries have inflated free inclinations compared to their red counterparts, consistent with the push-out origin for these bodies., Comment: Accepted for publication in the PSJ. 4 Figures

This paper revisits the optimal shape problem of a single cooling fin using a one-dimensional heat conduction equation with convection boundary conditions. Firstly, in contrast to previous works, we apply an approach using optimality conditions based on requiring stationarity of the Lagrangian functional of the optimisation problem. This yields an optimality condition basis for the commonly touted constant temperature gradient condition. Secondly, we seek to minimise the root temperature for a prescribed thermal power, rather than maximising the heat transfer rate for a constant root temperature as previous works. The optimal solution is shown to be fully equivalent for the two, which may seem obvious but to our knowledge has not been shown directly before. Lastly, it is shown that optimal cooling fins have a Biot number of 1, exhibiting perfect balance between conductive and convective resistances.

A method for density-based topology optimization of heat exchangers with two fluids is proposed. The goal of the optimization process is to maximize the heat transfer from one fluid to the other, under maximum pressure drop constraints for each of the fluid flows. A single design variable is used to describe the physical fields. The solid interface and the fluid domains are generated using an erosion-dilation based identification technique, which guarantees well-separated fluids, as well as a minimum wall thickness between them. Under the assumption of laminar steady flow, the two fluids are modelled separately, but in the entire computational domain using the Brinkman penalization technique for ensuring negligible velocities outside of the respective fluid subdomains. The heat transfer is modelled using the convection-diffusion equation, where the convection is driven by both fluid flows. A stabilized finite element discretization is used to solve the governing equations. Results are presented for two different problems: a two-dimensional example illustrating and verifying the methodology; and a three-dimensional example inspired by shell-and-tube heat exchangers. The optimized designs for both cases show an improved heat transfer compared to the baseline designs. For the shell-and-tube case, the full freedom topology optimization approach is shown to yield performance improvements of up to 113% under the same pressure drop.

The minor planets on orbits that are dynamically stable in Neptune's 1:1 resonance on Gyr timescales were likely emplaced by Neptune's outward migration. We explore the intrinsic libration amplitude, eccentricity, and inclination distribution of Neptune's stable Trojans, using the detections and survey efficiency of the Outer Solar System Origins Survey (OSSOS) and Pan-STARRS1. We find that the libration amplitude of the stable Neptunian Trojan population can be well modeled as a Rayleigh distribution with a libration amplitude width $\sigma_{A_phi}$ of 15$^\circ$. When taken as a whole, the Neptune Trojan population can be acceptably modeled with a Rayleigh eccentricity distribution of width $\sigma_e$ of 0.045 and a typical sin(i) x Gaussian inclination distribution with a width $\sigma_i$ of 14 +/- 2 degrees. However, these distributions are only marginally acceptable. This is likely because, even after accounting for survey detection biases, the known large Hr < 8 and small Hr >= 8 Neptune Trojans appear to have markedly different eccentricities and inclinations. We propose that like the classical Kuiper belt, the stable intrinsic Neptunian Trojan population have dynamically `hot' and dynamically `cold' components to its eccentricity/inclination distribution, with $\sigma_{e-cold}$ ~ 0.02 / $\sigma_{i-cold}$ ~ 6$^\circ$ and $\sigma_{e-hot}$~ 0.05 / $\sigma_{i-hot}$ ~ 18$^\circ$. In this scenario, the `cold' L4 Neptunian Trojan population lacks the Hr >= 8 members and has 13 +11/-6 `cold' Trojans with Hr < 8. On the other hand, the `hot' L4 Neptunian Trojan population has 136 +57/-48 Trojans with Hr < 10 -- a population 2.4 times greater than that of the L4 Jovian Trojans in the same luminosity range., Comment: 15 pages, 7 figures, v3, accepted for publication in Icarus

Observations show that 100-km-class Kuiper belt objects (KBOs) can be divided in (at least) two color groups, hereafter red (R, g-i<1.2) and very red (VR, g-i>1.2), reflecting a difference in their surface composition. This is thought to imply that KBOs formed over a relatively wide range of radial distance, r. The cold classicals at 42r*, with 30

The surface characterization of Trans-Neptunian Binaries (TNBs) is key to understanding the properties of the disk of planetesimals from which these objects formed. In the optical wavelengths, it has been demonstrated that most equal-sized component systems share similar colors, suggesting they have a similar composition. The color homogeneity of binary pairs contrasts with the overall diversity of colors in the Kuiper belt, which was interpreted as evidence that Trans-Neptunian Objects (TNOs) formed from a locally homogeneous and globally heterogeneous protoplanetary disk. In this paradigm, binary pairs must have formed early, before the dynamically hot TNOs were scattered out from their formation location. The latter inferences, however, relied on the assumption that the matching colors of the binary components imply matching composition. Here, we test this assumption by examining the component-resolved photometry of three TNBs found in the Outer Solar System Origins Survey: 505447 (2013 SQ99), 511551 (2014 UD225) and 506121 (2016 BP81), across the visible and J-band near-infrared wavelength range. We report similar colors within 2 sigma for the binary pairs suggestive of similar reflectance spectra and hence surface composition. This advocates for gravitational collapse of pebble clouds as a possible TNO formation route. We however stress that several similarly small TNOs, including at least one binary, have been shown to exhibit substantial spectral variability in the near-infrared, implying color equality of binary pairs is likely to be violated in some cases., Comment: 10 pages, 4 figures, accepted to The Planetary Science Journal

The Latitude Density Search utilized Hyper Suprime-Cam on Subaru Telescope to discover 60 moving objects in the outer Solar System, 54 of which have semi-major axes beyond 30 AU. The images were acquired in exceptional seeing (0.4") and reached a detection limit of m_r~=25.2. The two night arcs were used to calculate orbits which are poorly constrained, however, the distance and inclination are the parameters best constrained by short arcs and a reasonable determination can be made of which objects are cold classical TNOs and which are dynamically excited. We identify 10 objects as likely cold classical objects. We searched all of the detections for binary sources using a trailed Point Spread Function subtraction method, and identified one binary object with a separation of 0.34" and a secondary with 17% the brightness of the primary (2.0 magnitudes fainter). This is the brightest TNO in the sample, the previously known object 471165 (2010 HE79), which has a dynamically excited orbit. Because of the excellent seeing, this search was sensitive to binaries with 0.34" separation and a brightness of >=50% the primary brightness for 7 objects, including one cold classical. This gives an intrinsic binary fraction of ~17% (1 of 6) for the dynamically excited objects within our detection limits. The trailed point spread function subtraction method to identify binaries, fit the two components, and determine the sensitivity limits, used in the Latitude Density Search is a useful tool that could be more broadly applied to identify binary TNOs and track known binary TNO orbits., Comment: 14 pages, 6 figures

With the increasing power density of electronics components, the heat dissipation capacity of heat sinks gradually becomes a bottleneck. Many structural optimization methods, including topology optimization, have been widely used for heat sinks. Due to its high design freedom, topology optimization is suggested for the design of heat sinks using a transient pseudo-3D thermofluid model to acquire better instantaneous thermal performance. The pseudo-3D model is designed to reduce the computational cost and maintain an acceptable accuracy. The model relies on an artificial heat convection coefficient to couple two layers and establish the approximate relationship with the corresponding 3D model. In the model, a constant pressure drop and heat generation rate are treated. The material distribution is optimized to reduce the average temperature of the base plate at the prescribed terminal time. Furthermore, to reduce the intermediate density regions during the density-based topology optimization procedure, a detailed analysis of interpolation functions is made and the penalty factors are chosen on this basis. Finally, considering the engineering application of the model, a practical model with more powerful cooling medium and higher inlet pressure is built. The optimized design shows a better instantaneous thermal performance and provides 66.7% of the pumping power reduction compared with reference design., Comment: 55 pages, 36 figures

Resonant dynamics plays a significant role in the past evolution and current state of our outer Solar System. The population ratios and spatial distribution of Neptune's resonant populations are direct clues to understanding the history of our planetary system. The orbital structure of the objects in Neptune's 2:1 mean-motion resonance (\emph{twotinos}) has the potential to be a tracer of planetary migration processes. Different migration processes produce distinct architectures, recognizable by well-characterized surveys. However, previous characterized surveys only discovered a few twotinos, making it impossible to model the intrinsic twotino population. With a well-designed cadence and nearly 100\% tracking success, the Outer Solar System Origins Survey (OSSOS) discovered 838 trans-Neptunian objects, of which 34 are securely twotinos with well-constrained libration angles and amplitudes. We use the OSSOS twotinos and the survey characterization parameters via the OSSOS Survey Simulator to inspect the intrinsic population and orbital distributions of twotino. The estimated twotino population, 4400$^{+1500}_{-1100}$ with $H_r<8.66$ (diameter$\sim$100km) at 95\% confidence, is consistent with the previous low-precision estimate. We also constrain the width of the inclination distribution to a relatively narrow value of $\sigma_i$=6$^\circ$$^{+1}_{-1}$, and find the eccentricity distribution is consistent with a Gaussian centered on $e_\mathrm{c}=0.275$ with a width $e_\mathrm{w}=0.06$. We find a single-slope exponential luminosity function with $\alpha=0.6$ for the twotinos. Finally, we for the first time meaningfully constrain the fraction of symmetric twotinos, and the ratio of the leading asymmetric islands; both fractions are in a range of 0.2--0.6. These measurements rule out certain theoretical models of Neptune's migration history., Comment: 23 pages, 12 figures, 2 tables, accepted for publication in the Astronomical Journal

While collisional families are common in the asteroid belt, only one is known in the Kuiper belt, linked to the dwarf planet Haumea. The characterization of Haumea's family helps to constrain its origin and, more generally, the collisional history of the Kuiper belt. However, the size distribution of the Haumea family is difficult to constrain from the known sample, which is affected by discovery biases. Here, we use the Outer Solar System Origins Survey (OSSOS) Ensemble to look for Haumea family members. In this OSSOS XVI study we report the detection of three candidates with small ejection velocities relative to the family formation centre. The largest discovery, 2013 UQ15, is conclusively a Haumea family member, with a low ejection velocity and neutral surface colours. Although the OSSOS Ensemble is sensitive to Haumea family members to a limiting absolute magnitude (Hr) of 9.5 (inferred diameter of ~90 km), the smallest candidate is significantly larger, Hr = 7.9. The Haumea family members larger than approximately 20 km in diameter must be characterized by a shallow H-distribution slope in order to produce only these three large detections. This shallow size distribution suggests that the family formed in a graze-and-merge scenario, not a catastrophic collision., Comment: 11 pages, 3 figures, 3 tables

We use published models of the early Solar System evolution with a slow, long-range and grainy migration of Neptune to predict the orbital element distributions and the number of modern-day Centaurs. The model distributions are biased by the OSSOS survey simulator and compared with the OSSOS Centaur detections. We find an excellent match to the observed orbital distribution, including the wide range of orbital inclinations which was the most troublesome characteristic to fit in previous models. A dynamical model, in which the original population of outer disk planetesimals was calibrated from Jupiter Trojans, is used to predict that OSSOS should detect 11+/-4 Centaurs with semimajor axis a<30 au, perihelion distance q>7.5 au and diameter D>10 km (absolute magnitude H_r<13.7 for a 6% albedo). This is consistent with 15 actual OSSOS Centaur detections with H_r<13.7. The population of Centaurs is estimated to be 21,000+/-8,000 for D>10 km. The inner scattered disk at 5010 km bodies and the Oort cloud should contain (5.0+/-1.9)x10^8 D>10 km comets. Population estimates for different diameter cutoffs can be obtained from the size distribution of Jupiter Trojans (N(>D) proportional to D^(-2.1) for 5

Most known trans-Neptunian objects (TNOs) gravitationally scattering off the giant planets have orbital inclinations consistent with an origin from the classical Kuiper belt, but a small fraction of these "scattering TNOs" have inclinations that are far too large (i > 45 deg) for this origin. These scattering outliers have previously been proposed to be interlopers from the Oort cloud or evidence of an undiscovered planet. Here we test these hypotheses using N-body simulations and the 69 centaurs and scattering TNOs detected in the Outer Solar Systems Origins Survey and its predecessors. We confirm that observed scattering objects cannot solely originate from the classical Kuiper belt, and we show that both the Oort cloud and a distant planet generate observable highly inclined scatterers. Although the number of highly inclined scatterers from the Oort Cloud is ~3 times less than observed, Oort cloud enrichment from the Sun's galactic migration or birth cluster could resolve this. Meanwhile, a distant, low-eccentricity 5 Earth-mass planet replicates the observed fraction of highly inclined scatterers, but the overall inclination distribution is more excited than observed. Furthermore, the distant planet generates a longitudinal asymmetry among detached TNOs that is less extreme than often presumed, and its direction reverses across the perihelion range spanned by known TNOs. More complete models that explore the dynamical origins of the planet are necessary to further study these features. With observational biases well-characterized, our work shows that the orbital distribution of detected scattering bodies is a powerful constraint on the unobserved distant solar system., Comment: 17 pages, 11 figures, accepted to AJ

This paper treats topology optimization of natural convection problems. A simplified model is suggested to describe the flow of an incompressible fluid in steady state conditions, similar to Darcy's law for fluid flow in porous media. The equations for the fluid flow are coupled to the thermal convection-diffusion equation through the Boussinesq approximation. The coupled non-linear system of equations is discretized with stabilized finite elements and solved in a parallel framework that allows for the optimization of high resolution three-dimensional problems. A density-based topology optimization approach is used, where a two-material interpolation scheme is applied to both the permeability and conductivity of the distributed material. Due to the simplified model, the proposed methodology allows for a significant reduction of the computational effort required in the optimization. At the same time, it is significantly more accurate than even simpler models that rely on convection boundary conditions based on Newton's law of cooling. The methodology discussed herein is applied to the optimization-based design of three-dimensional heat sinks. The final designs are formally compared with results of previous work obtained from solving the full set of Navier-Stokes equations. The results are compared in terms of performance of the optimized designs and computational cost. The computational time is shown to be decreased to around 5-20% in terms of core-hours, allowing for the possibility of generating an optimized design during the workday on a small computational cluster and overnight on a high-end desktop.

We present variability measurements and partial light curves of Trans-Neptunian Objects (TNOs) from a two-night pilot study using Hyper Suprime-Cam (HSC) on the Subaru Telescope (Maunakea, Hawai'i, USA). Subaru's large aperture (8-m) and HSC's large field of view (1.77 square degrees) allow us to obtain measurements of multiple objects with a range of magnitudes in each telescope pointing. We observed 65 objects with m_r = 22.6--25.5 mag in just six pointings, allowing 20--24 visits of each pointing over the two nights. Our sample, all discovered in the recent Outer Solar System Origins Survey (OSSOS), span absolute magnitudes H_r = 6.2--10.8 mag and thus investigates smaller objects than previous light curve projects have typically studied. Our data supports the existence of a correlation between light curve amplitude and absolute magnitude seen in other works, but does not support a correlation between amplitude and orbital inclination. Our sample includes a number of objects from different dynamical populations within the trans-Neptunian region, but we do not find any relationship between variability and dynamical class. We were only able to estimate periods for 12 objects in the sample and found that a longer baseline of observations is required for reliable period analysis. We find that 31 objects (just under half of our sample) have variability greater than 0.4 magnitudes during all of the observations; in smaller 1.25 hr, 1.85 hr and 2.45 hr windows, the median variability is 0.13, 0.16 and 0.19 mags, respectively. The fact that variability on this scale is common for small TNOs has important implications for discovery surveys (such as OSSOS or the Large Synoptic Survey Telescope) and color measurements., Comment: 49 pages, 14 figures, 7 tables

Both physical and dynamical properties must be considered to constrain the origins of the dynamically excited distant Solar System populations. We present high-precision (g-r) colors for 25 small (Hr>5) dynamically excited Trans-Neptunian Objects (TNOs) and centaurs acquired as part of the Colours of the Outer Solar System Origins Survey (Col-OSSOS). We combine our dataset with previously published measurements and consider a set of 229 colors of outer Solar System objects on dynamically excited orbits. The overall color distribution is bimodal and can be decomposed into two distinct classes, termed `gray' and `red', that each has a normal color distribution. The two color classes have different inclination distributions: red objects have lower inclinations than the gray ones. This trend holds for all dynamically excited TNO populations. Even in the worst-case scenario, biases in the discovery surveys cannot account for this trend: it is intrinsic to the TNO population. Considering that TNOs are the precursors of centaurs, and that their inclinations are roughly preserved as they become centaurs, our finding solves the conundrum of centaurs being the only outer Solar System population identified so far to exhibit this property (Tegler et al. 2016). The different orbital distributions of the gray and red dynamically excited TNOs provide strong evidence that their colors are due to different formation locations in a disk of planetesimals with a compositional gradient., Comment: 24 pages, 5 figures, accepted for publication in AJ

Context. Centaurs are icy objects in transition between the transneptunian region and the inner solar system, orbiting the Sun in the giant planet region. Some Centaurs display cometary activity, which cannot be sustained by the sublimation of water ice in this part of the solar system, and has been hypothesized to be due to the crystallization of amorphous water ice. Aims. In this work, we look at Centaurs discovered by the Outer Solar System Origins Survey (OSSOS) and search for cometary activity. Tentative detections would improve understanding of the origins of activity among these objects. Methods. We search for comae and structures by fitting and subtracting both Point Spread Functions (PSF) and Trailed point-Spread Functions (TSF) from the OSSOS images of each Centaur. When available, Col-OSSOS images were used to search for comae too. Results. No cometary activity is detected in the OSSOS sample. We track the recent orbital evolution of each new Centaur to confirm that none would actually be predicted to be active, and we provide size estimates for the objects. Conclusions. The addition of 20 OSSOS objects to the population of 250 known Centaurs is consistent with the currently understood scenario, in which drastic drops in perihelion distance induce changes in the thermal balance prone to trigger cometary activity in the giant planet region., Comment: 7 pages, 6 figures, Accepted for publication in Astronomy and Astrophysics

The Colours of the Outer Solar System Origins Survey (Col-OSSOS) is acquiring near-simultaneous $g$, $r$, and $J$ photometry of unprecedented precision with the Gemini North Telescope, targeting nearly a hundred trans-Neptunian objects (TNOs) brighter than $m_r=23.6$ mag discovered in the Outer Solar System Origins Survey. Combining the optical and near-infrared photometry with the well-characterized detection efficiency of the Col-OSSOS target sample will provide the first flux-limited compositional dynamical map of the outer Solar System. In this paper, we describe our observing strategy and detail the data reduction processes we employ, including techniques to mitigate the impact of rotational variability. We present optical and near-infrared colors for 35 TNOs. We find two taxonomic groups for the dynamically excited TNOs, the neutral and red classes, which divide at $g-r \simeq 0.75$. Based on simple albedo and orbital distribution assumptions, we find that the neutral class outnumbers the red class, with a ratio of 4:1 and potentially as high as 11:1. Including in our analysis constraints from the cold classical objects, which are known to exhibit unique albedos and $r-z$ colors, we find that within our measurement uncertainty, our observations are consistent with the primordial Solar System protoplanetesimal disk being neutral-class-dominated, with two major compositional divisions in $grJ$ color space., Comment: Accepted to ApJS; on-line supplemental files will be available with the AJS published version of the paper

Topology optimization of natural convection problems is computationally expensive, due to the large number of degrees of freedom (DOFs) in the model and its two-way coupled nature. Herein, a method is presented to reduce the computational effort by use of a reduced-order model governed by simplified physics. The proposed method models the fluid flow using a potential flow model, which introduces an additional fluid property. This material property currently requires tuning of the model by comparison to numerical Navier-Stokes based solutions. Topology optimization based on the reduced-order model is shown to provide qualitatively similar designs, as those obtained using a full Navier-Stokes based model. The number of DOFs is reduced by 50% in two dimensions and the computational complexity is evaluated to be approximately 12.5% of the full model. We further compare to optimized designs obtained utilizing Newton's convection law., Comment: Preprint version. Please refer to final version in Structural Multidisciplinary Optimization https://doi.org/10.1007/s00158-019-02215-9

The migration of Neptune's resonances through the proto-Kuiper belt has been imprinted in the distribution of small bodies in the outer Solar System. Here we analyze five published Neptune migration models in detail, focusing on the high pericenter distance (high-q) trans-Neptunian Objects (TNOs) near Neptune's 5:2 and 3:1 mean-motion resonances, because they have large resonant populations, are outside the main classical belt, and are relatively isolated from other strong resonances. We compare the observationally biased output from these dynamical models with the detected TNOs from the Outer Solar System Origins Survey, via its Survey Simulator. All of the four new OSSOS detections of high-q non-resonant TNOs are on the Sunward side of the 5:2 and 3:1 resonances. We show that even after accounting for observation biases, this asymmetric distribution cannot be drawn from a uniform distribution of TNOs at 2sigma confidence. As shown by previous work, our analysis here tentatively confirms that the dynamical model that uses grainy slow Neptune migration provides the best match to the real high-q TNO orbital data. However, due to extreme observational biases, we have very few high-q TNO discoveries with which to statistically constrain the models. Thus, this analysis provides a framework for future comparison between the output from detailed, dynamically classified Neptune migration simulations and the TNO discoveries from future well-characterized surveys. We show that a deeper survey (to a limiting r-magnitude of 26.0) with a similar survey area to OSSOS could statistically distinguish between these five Neptune migration models., Comment: Accepted for publication in the Astronomical Journal

The Outer Solar System Origins Survey (OSSOS), a wide-field imaging program in 2013-2017 with the Canada-France-Hawaii Telescope, surveyed 155 deg$^{2}$ of sky to depths of $m_r = 24.1$-25.2. We present 838 outer Solar System discoveries that are entirely free of ephemeris bias. This increases the inventory of trans-Neptunian objects (TNOs) with accurately known orbits by nearly 50%. Each minor planet has 20-60 Gaia/Pan-STARRS-calibrated astrometric measurements made over 2-5 oppositions, which allows accurate classification of their orbits within the trans-Neptunian dynamical populations. The populations orbiting in mean-motion resonance with Neptune are key to understanding Neptune's early migration. Our 313 resonant TNOs, including 132 plutinos, triple the available characterized sample and include new occupancy of distant resonances out to semi-major axis $a \sim 130$ au. OSSOS doubles the known population of the non-resonant Kuiper belt, providing 436 TNOs in this region, all with exceptionally high-quality orbits of $a$ uncertainty $\sigma_{a} \leq 0.1\%$; they show the belt exists from $a \gtrsim 37$ au, with a lower perihelion bound of $35$ au. We confirm the presence of a concentrated low-inclination $a\simeq 44$ au "kernel" population and a dynamically cold population extending beyond the 2:1 resonance. We finely quantify the survey's observational biases. Our survey simulator provides a straightforward way to impose these biases on models of the trans-Neptunian orbit distributions, allowing statistical comparison to the discoveries. The OSSOS TNOs, unprecedented in their orbital precision for the size of the sample, are ideal for testing concepts of the history of giant planet migration in the Solar System., Comment: Invited paper, special issue Data: Insights and Challenges in a Time of Abundance. Data tables and example survey simulator are in the supplementary materials (see arXiv source under Downloads > Other formats)

The scattering trans-Neptunian Objects (TNOs) can be measured to smaller sizes than any other distant small-body population. We use the largest sample yet obtained, 68 discoveries, primarily by the Outer Solar System Origins Survey (OSSOS), to constrain the slope of its luminosity distribution, with sensitivity to much fainter absolute $H$ magnitudes than previous work. Using the analysis technique in Shankman et al. (2016), we confirm that a single slope for the $H$-distribution is not an accurate representation of the scattering TNOs and Centaurs, and that a break in the distribution is required, in support of previous conclusions. A bright-end slope of $\alpha_b=0.9$ transitioning to a faint-end slope $\alpha_f$ of 0.4-0.5 with a differential number contrast $c$ from 1 (a knee) to 10 (a divot) provides an acceptable match to our data. We find that break magnitudes $H_b$ of 7.7 and 8.3, values both previously suggested for dynamically hot Kuiper belt populations, are equally non-rejectable for a range of $\alpha_f$ and $c$ in our statistical analysis. Our preferred divot $H$-distribution transitions to $\alpha_f=0.5$ with a divot of contrast $c=3$ at $H_b=8.3$, while our preferred knee $H$-distribution transitions to $\alpha_f=0.4$ at $H_b=7.7$. The intrinsic population of scattering TNOs required to match the OSSOS detections is $3\times10^6$ for $H_r<12$, and $9\times10^4$ for $H_r<8.66$ ($D\gtrsim100$~km), with Centaurs having an intrinsic population two orders of magnitude smaller., Comment: accepted for publication in AJ

We discuss the detection in the Outer Solar System Origins Survey (OSSOS) of two objects in Neptune's distant 9:1 mean motion resonance at semimajor axis $a\approx~130$~au. Both objects are securely resonant on 10~Myr timescales, with one securely in the 9:1 resonance's leading asymmetric libration island and the other in either the symmetric or trailing asymmetric island. These objects are the largest semimajor axis objects with secure resonant classifications, and their detection in a carefully characterized survey allows for the first robust resonance population estimate beyond 100~au. The detection of these objects implies a 9:1 resonance population of $1.1\times10^4$ objects with $H_r<8.66$ ($D~\gtrsim~100$~km) on similar orbits (95\% confidence range of $\sim0.4-3\times10^4$). Integrations over 4~Gyr of an ensemble of clones spanning these objects' orbit fit uncertainties reveal that they both have median resonance occupation timescales of $\sim1$~Gyr. These timescales are consistent with the hypothesis that these objects originate in the scattering population but became transiently stuck to Neptune's 9:1 resonance within the last $\sim1$~Gyr of solar system evolution. Based on simulations of a model of the current scattering population, we estimate the expected resonance sticking population in the 9:1 resonance to be 1000-4500 objects with $H_r<8.66$; this is marginally consistent with the OSSOS 9:1 population estimate. We conclude that resonance sticking is a plausible explanation for the observed 9:1 population, but we also discuss the possibility of a primordial 9:1 population, which would have interesting implications for the Kuiper belt's dynamical history., Comment: accepted for publication in AJ

All surveys include observational biases, which makes it impossible to directly compare properties of discovered trans-Neptunian Objects (TNOs) with dynamical models. However, by carefully keeping track of survey pointings on the sky, detection limits, tracking fractions, and rate cuts, the biases from a survey can be modelled in Survey Simulator software. A Survey Simulator takes an intrinsic orbital model (from, for example, the output of a dynamical Kuiper belt emplacement simulation) and applies the survey biases, so that the biased simulated objects can be directly compared with real discoveries. This methodology has been used with great success in the Outer Solar System Origins Survey (OSSOS) and its predecessor surveys. In this chapter, we give four examples of ways to use the OSSOS Survey Simulator to gain knowledge about the true structure of the Kuiper Belt. We demonstrate how to statistically compare different dynamical model outputs with real TNO discoveries, how to quantify detection biases within a TNO population, how to measure intrinsic population sizes, and how to use upper limits from non-detections. We hope this will provide a framework for dynamical modellers to statistically test the validity of their models., Comment: accepted as a chapter in the "Frontiers in Astronomy and Space Sciences" research topic "From Comets to Pluto and Beyond"

Several different classes of trans-Neptunian objects (TNOs) have been identified based on their optical and near-infrared colors. As part of the Colours of the Outer Solar System Origins Survey, we have obtained $g$, $r$, and $z$ band photometry of 26 TNOs using Subaru and Gemini Observatories. Previous color surveys have not utilized $z$ band reflectance, and the inclusion of this band reveals significant surface reflectance variations between sub-populations. The colors of TNOs in $g-r$ and $r-z$ show obvious structure, and appear consistent with the previously measured bi-modality in $g-r$. The distribution of colors of the two dynamically excited surface types can be modeled using the two-component mixing models from Fraser \& Brown (2012). With the combination of $g-r$ and $r-z$, the dynamically excited classes can be separated cleanly into red and neutral surface classes. In $g - r$ and $r - z$, the two dynamically excited surface groups are also clearly distinct from the cold classical TNO surfaces, which are red, with $g-r\gtrsim$0.85 and $r-z\lesssim$0.6, while all dynamically excited objects with similar $g-r$ colors exhibit redder $r-z$ colors. The $z$ band photometry makes it possible for the first time to differentiate the red excited TNO surfaces from the red cold classical TNO surfaces. The discovery of different $r-z$ colors for these cold classical TNOs makes it possible to search for cold classical surfaces in other regions of the Kuiper belt and to completely separate cold classical TNOs from the dynamically excited population, which overlaps in orbital parameter space., Comment: 11 pages, 2 figures, Accepted to AJ

The accumulating, but small, set of large semi-major axis trans-Neptunian objects (TNOs) shows an apparent clustering in the orientations of their orbits. This clustering must either be representative of the intrinsic distribution of these TNOs, or else arise as a result of observation biases and/or statistically expected variations for such a small set of detected objects. The clustered TNOs were detected across different and independent surveys, which has led to claims that the detections are therefore free of observational bias. This apparent clustering has led to the so-called "Planet 9" hypothesis that a super-Earth currently resides in the distant solar system and causes this clustering. The Outer Solar System Origins Survey (OSSOS) is a large program that ran on the Canada-France-Hawaii Telescope from 2013--2017, discovering more than 800 new TNOs. One of the primary design goals of OSSOS was the careful determination of observational biases that would manifest within the detected sample. We demonstrate the striking and non-intuitive biases that exist for the detection of TNOs with large semi-major axes. The eight large semi-major axis OSSOS detections are an independent dataset, of comparable size to the conglomerate samples used in previous studies. We conclude that the orbital distribution of the OSSOS sample is consistent with being detected from a uniform underlying angular distribution., Comment: Accepted for publication

The Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) is currently the deepest wide- field survey in progress. The 8.2 m aperture of Subaru telescope is very powerful in detect- ing faint/small moving objects, including near-Earth objects, asteroids, centaurs and Tran- Neptunian objects (TNOs). However, the cadence and dithering pattern of the HSC-SSP are not designed for detecting moving objects, making it difficult to do so systematically. In this paper, we introduce a new pipeline for detecting moving objects (specifically TNOs) in a non-dedicated survey. The HSC-SSP catalogs are re-arranged into the HEALPix architecture. Then, the stationary detections and false positive are removed with a machine learning al- gorithm to produce a list of moving object candidates. An orbit linking algorithm and visual inspections are executed to generate the final list of detected TNOs. The preliminary results of a search for TNOs using this new pipeline on data from the first HSC-SSP data release (Mar 2014 to Nov 2015) are also presented., Comment: 32 pages, 10 figures, 2 tables, submitted to HSC special issue in PASJ

The cold classical Kuiper belt objects have low inclinations and eccentricities and are the only Kuiper belt population suspected to have formed in situ. Compared with the dynamically excited populations, which exhibit a broad range of colours and a low binary fraction of ~10% cold classical Kuiper belt objects typically have red optical colours with ~30% of the population found in binary pairs; the origin of these differences remains unclear. We report the detection of a population of blue-coloured, tenuously bound binaries residing among the cold classical Kuiper belt objects. Here we show that widely separated binaries could have survived push-out into the cold classical region during the early phases of Neptune's migration. The blue binaries may be contaminants, originating at ~38 au, and could provide a unique probe of the formative conditions in a region now nearly devoid of objects. The idea that the blue objects, which are predominantly binary, are the products of push-out requires that the planetesimals formed entirely as multiples. Plausible formation routes include planetesimal formation via pebble accretion and subsequent binary production through dynamic friction and binary formation during the collapse of a cloud of solids., Comment: 7 Figures, 3 tables, accepted to Nature Astronomy. Main manuscript and supplement available at http://www.nature.com/articles/s41550-017-0088

We report the discovery of the minor planet 2013 SY$_{99}$, on an exceptionally distant, highly eccentric orbit. With a perihelion of 50.0 au, 2013 SY$_{99}$'s orbit has a semi-major axis of $730 \pm 40$ au, the largest known for a high-perihelion trans-Neptunian object (TNO), well beyond those of (90377) Sedna and 2012 VP$_{113}$. Yet, with an aphelion of $1420 \pm 90$ au, 2013 SY$_{99}$'s orbit is interior to the region influenced by Galactic tides. Such TNOs are not thought to be produced in the current known planetary architecture of the Solar System, and they have informed the recent debate on the existence of a distant giant planet. Photometry from the Canada-France-Hawaii Telescope, Gemini North and Subaru indicate 2013 SY$_{99}$ is $\sim 250$ km in diameter and moderately red in colour, similar to other dynamically excited TNOs. Our dynamical simulations show that Neptune's weak influence during 2013 SY$_{99}$'s perihelia encounters drives diffusion in its semi-major axis of hundreds of astronomical units over 4 Gyr. The overall symmetry of random walks in semi-major axis allow diffusion to populate 2013 SY$_{99}$'s orbital parameter space from the 1000-2000 au inner fringe of the Oort cloud. Diffusion affects other known TNOs on orbits with perihelia of 45 to 49 au and semi-major axes beyond 250 au, providing a formation mechanism that implies an extended population, gently cycling into and returning from the inner fringe of the Oort cloud., Comment: First reviewer report comments incorporated. Comments welcome