Classification and reduction of pilot error

Human error is a primary or contributing factor in about two-thirds of commercial aviation accidents worldwide. With the ultimate goal of reducing pilot error accidents, this contract effort is aimed at understanding the factors underlying error events and reducing the probability of certain types of errors by modifying underlying factors such as flight deck design and procedures. A review of the literature relevant to error classification was conducted. Classification includes categorizing types of errors, the information processing mechanisms and factors underlying them, and identifying factor-mechanism-error relationships. The classification scheme developed by Jens Rasmussen was adopted because it provided a comprehensive yet basic error classification shell or structure that could easily accommodate addition of details on domain-specific factors. For these purposes, factors specific to the aviation environment were incorporated. Hypotheses concerning the relationship of a small number of underlying factors, information processing mechanisms, and error types types identified in the classification scheme were formulated. ASRS data were reviewed and a simulation experiment was performed to evaluate and quantify the hypotheses

The fragmentation of protostellar discs: the Hill criterion for spiral arms

We present a new framework to explain the link between cooling and fragmentation in gravitationally unstable protostellar discs. This framework consists of a simple model for the formation of spiral arms, as well as a criterion, based on the Hill radius, to determine if a spiral arm will fragment. This detailed model of fragmentation is based on the results of numerical simulations of marginally stable protostellar discs, including those found in the literature, as well as our new suite of 3-D radiation hydrodynamics simulations of an irradiated, optically-thick protostellar disc surrounding an A star. Our set of simulations probes the transition to fragmentation through a scaling of the physical opacity. This model allows us to directly calculate the critical cooling time of Gammie (2001), with results that are consistent with those found from numerical experiment. We demonstrate how this model can be used to predict fragmentation in irradiated protostellar discs. These numerical simulations, as well as the model that they motivate, provide strong support for the hypothesis that gravitational instability is responsible for creating systems with giant planets on wide orbits.Comment: 11 page, 10 figures, submitted to MNRA

Chemistry in a gravitationally unstable protoplanetary disc

Until now, axisymmetric, alpha-disc models have been adopted for calculations of the chemical composition of protoplanetary discs. While this approach is reasonable for many discs, it is not appropriate when self-gravity is important. In this case, spiral waves and shocks cause temperature and density variations that affect the chemistry. We have adopted a dynamical model of a solar-mass star surrounded by a massive (0.39 Msun), self-gravitating disc, similar to those that may be found around Class 0 and early Class I protostars, in a study of disc chemistry. We find that for each of a number of species, e.g. H2O, adsorption and desorption dominate the changes in the gas-phase fractional abundance; because the desorption rates are very sensitive to temperature, maps of the emissions from such species should reveal the locations of shocks of varying strengths. The gas-phase fractional abundances of some other species, e.g. CS, are also affected by gas-phase reactions, particularly in warm shocked regions. We conclude that the dynamics of massive discs have a strong impact on how they appear when imaged in the emission lines of various molecular species.Comment: 10 figures and 3 tables, accepted for publication in MNRA

ALMA and VLA Observations of the HD 141569 System

We present VLA 9 mm (33 GHz) observations of the HD 141569 system from semester 16A. The observations achieve a resolution of 0.25 arcsec ($\sim28$ au) and a sensitivity of $4.7~\mu \rm Jy~beam^{-1}$. We find (1) a $52\pm 5~\mu$Jy point source at the location of HD 141569A that shows potential variability, (2) the detected flux is contained within the SED-inferred central clearing of the disc meaning the spectral index of the dust disc is steeper than previously inferred, and (3) the M dwarf companions are also detected and variable. Previous lower-resolution VLA observations (semester 14A) found a higher flux density, interpreted as solely dust emission. When combined with ALMA observations, the VLA 14A observations suggested the spectral index and grain size distribution of HD 141569's disc was shallow and an outlier among debris systems. Using archival ALMA observations of HD 141569 at 0.87 mm and 2.9 mm we find a dust spectral index of $\alpha_{\rm mm} = 1.81\pm 0.20$. The VLA 16A flux corresponds to a brightness temperature of $\sim5\times10^{6}$ K, suggesting strong non-disc emission is affecting the inferred grain properties. The VLA 16A flux density of the M2V companion HD 141569B is $149\pm9~\mu$Jy, corresponding to a brightness temperature of $\sim2\times10^{8}$ K and suggesting significant stellar variability when compared to the VLA14A observations, which are smaller by a factor of $\sim6$.Comment: Accepted for publication in MNRAS, 8 pages, 6 figure

On the accuracy of solving confluent Prony systems

In this paper we consider several nonlinear systems of algebraic equations which can be called "Prony-type". These systems arise in various reconstruction problems in several branches of theoretical and applied mathematics, such as frequency estimation and nonlinear Fourier inversion. Consequently, the question of stability of solution with respect to errors in the right-hand side becomes critical for the success of any particular application. We investigate the question of "maximal possible accuracy" of solving Prony-type systems, putting stress on the "local" behavior which approximates situations with low absolute measurement error. The accuracy estimates are formulated in very simple geometric terms, shedding some light on the structure of the problem. Numerical tests suggest that "global" solution techniques such as Prony's algorithm and ESPRIT method are suboptimal when compared to this theoretical "best local" behavior

Numerical determination of the material properties of porous dust cakes

The formation of planetesimals requires the growth of dust particles through collisions. Micron-sized particles must grow by many orders of magnitude in mass. In order to understand and model the processes during this growth, the mechanical properties, and the interaction cross sections of aggregates with surrounding gas must be well understood. Recent advances in experimental (laboratory) studies now provide the background for pushing numerical aggregate models onto a new level. We present the calibration of a previously tested model of aggregate dynamics. We use plastic deformation of surface asperities as the physical model to bring critical velocities for sticking into accordance with experimental results. The modified code is then used to compute compression strength and the velocity of sound in the aggregate at different densities. We compare these predictions with experimental results and conclude that the new code is capable of studying the properties of small aggregates.Comment: Accepted for publication in A&

Activity Identification and Local Linear Convergence of Douglas--Rachford/ADMM under Partial Smoothness

Convex optimization has become ubiquitous in most quantitative disciplines of science, including variational image processing. Proximal splitting algorithms are becoming popular to solve such structured convex optimization problems. Within this class of algorithms, Douglas--Rachford (DR) and alternating direction method of multipliers (ADMM) are designed to minimize the sum of two proper lower semi-continuous convex functions whose proximity operators are easy to compute. The goal of this work is to understand the local convergence behaviour of DR (resp. ADMM) when the involved functions (resp. their Legendre-Fenchel conjugates) are moreover partly smooth. More precisely, when both of the two functions (resp. their conjugates) are partly smooth relative to their respective manifolds, we show that DR (resp. ADMM) identifies these manifolds in finite time. Moreover, when these manifolds are affine or linear, we prove that DR/ADMM is locally linearly convergent. When $J$ and $G$ are locally polyhedral, we show that the optimal convergence radius is given in terms of the cosine of the Friedrichs angle between the tangent spaces of the identified manifolds. This is illustrated by several concrete examples and supported by numerical experiments.Comment: 17 pages, 1 figure, published in the proceedings of the Fifth International Conference on Scale Space and Variational Methods in Computer Visio