-
Differentiable Optics with dLux I: Deep calibration of Flat Field and Phase Retrieval with Automatic Differentiation
Authors:
Louis Desdoigts,
Benjamin Pope,
Jordan Dennis,
Peter Tuthill
Abstract:
The sensitivity limits of space telescopes are imposed by uncalibrated errors in the point spread function, photon-noise, background light, and detector sensitivity. These are typically calibrated with specialized wavefront sensor hardware and with flat fields obtained on the ground or with calibration sources, but these leave vulnerabilities to residual time-varying or non-common path aberrations…
▽ More
The sensitivity limits of space telescopes are imposed by uncalibrated errors in the point spread function, photon-noise, background light, and detector sensitivity. These are typically calibrated with specialized wavefront sensor hardware and with flat fields obtained on the ground or with calibration sources, but these leave vulnerabilities to residual time-varying or non-common path aberrations and variations in the detector conditions. It is therefore desirable to infer these from science data alone, facing the prohibitively high dimensional problems of phase retrieval and pixel-level calibration. We introduce a new Python package for physical optics simulation, dLux, which uses the machine learning framework JAX to achieve GPU acceleration and automatic differentiation (autodiff), and apply this to simulating astronomical imaging. In this first of a series of papers, we show that gradient descent enabled by autodiff can be used to simultaneously perform phase retrieval and calibration of detector sensitivity, scaling efficiently to inferring millions of parameters. This new framework enables high dimensional optimization and inference in data analysis and hardware design in astronomy and beyond, which we explore in subsequent papers in this series.
△ Less
Submitted 16 June, 2024; v1 submitted 12 June, 2024;
originally announced June 2024.
-
Fairness and Bias in Algorithmic Hiring: a Multidisciplinary Survey
Authors:
Alessandro Fabris,
Nina Baranowska,
Matthew J. Dennis,
David Graus,
Philipp Hacker,
Jorge Saldivar,
Frederik Zuiderveen Borgesius,
Asia J. Biega
Abstract:
Employers are adopting algorithmic hiring technology throughout the recruitment pipeline. Algorithmic fairness is especially applicable in this domain due to its high stakes and structural inequalities. Unfortunately, most work in this space provides partial treatment, often constrained by two competing narratives, optimistically focused on replacing biased recruiter decisions or pessimistically p…
▽ More
Employers are adopting algorithmic hiring technology throughout the recruitment pipeline. Algorithmic fairness is especially applicable in this domain due to its high stakes and structural inequalities. Unfortunately, most work in this space provides partial treatment, often constrained by two competing narratives, optimistically focused on replacing biased recruiter decisions or pessimistically pointing to the automation of discrimination. Whether, and more importantly what types of, algorithmic hiring can be less biased and more beneficial to society than low-tech alternatives currently remains unanswered, to the detriment of trustworthiness. This multidisciplinary survey caters to practitioners and researchers with a balanced and integrated coverage of systems, biases, measures, mitigation strategies, datasets, and legal aspects of algorithmic hiring and fairness. Our work supports a contextualized understanding and governance of this technology by highlighting current opportunities and limitations, providing recommendations for future work to ensure shared benefits for all stakeholders.
△ Less
Submitted 8 April, 2024; v1 submitted 25 September, 2023;
originally announced September 2023.
-
Just what the doctor ordered: An evaluation of provider preference-based Instrumental Variable methods in observational studies, with application for comparative effectiveness of type 2 diabetes therapy
Authors:
Laura M. Güdemann,
Beverley M. Shields,
John M. Dennis,
Jack Bowden
Abstract:
Instrumental Variables provide a way of addressing bias due to unmeasured confounding when estimating treatment effects using observational data. As instrument prescription preference of individual healthcare providers has been proposed. Because prescription preference is hard to measure and often unobserved, a surrogate measure constructed from available data is often required for the analysis. D…
▽ More
Instrumental Variables provide a way of addressing bias due to unmeasured confounding when estimating treatment effects using observational data. As instrument prescription preference of individual healthcare providers has been proposed. Because prescription preference is hard to measure and often unobserved, a surrogate measure constructed from available data is often required for the analysis. Different construction methods for this surrogate measure are possible, such as simple rule-based methods which make use of the observed treatment patterns, or more complex model-based methods that employ formal statistical models to explain the treatment behaviour whilst considering measured confounders. The choice of construction method relies on aspects like data availability within provider, missing data in measured confounders, and possible changes in prescription preference over time. In this paper we conduct a comprehensive simulation study to evaluate different construction methods for surrogates of prescription preference under different data conditions, including: different provider sizes, missing covariate data, and change in preference. We also propose a novel model-based construction method to address between provider differences and change in prescription preference. All presented construction methods are exemplified in a case study of the relative glucose lowering effect of two type 2 diabetes treatments in observational data. Our study shows that preference-based Instrumental Variable methods can be a useful tool for causal inference from observational health data. The choice of construction method should be driven by the data condition at hand. Our proposed method is capable of estimating the causal treatment effect without bias in case of sufficient prescription data per provider, changing prescription preference over time and non-ignorable missingness in measured confounders.
△ Less
Submitted 15 August, 2023;
originally announced August 2023.
-
Coherent Organizational States in Turbulent Pipe Flow at moderate Reynolds numbers
Authors:
R. Jäckel,
B. Magacho,
B. E. Owolabi,
L. Moriconi,
D. J. C. Dennis,
J. B. R. Loureiro
Abstract:
Turbulent pipe flow is still an essentially open area of research, boosted in the last two decades by considerable progress achieved both on the experimental and numerical frontiers, mainly related to the identification and characterization of coherent structures as basic building blocks of turbulence. It has been a challenging task, however, to detect and visualize these coherent states. We addre…
▽ More
Turbulent pipe flow is still an essentially open area of research, boosted in the last two decades by considerable progress achieved both on the experimental and numerical frontiers, mainly related to the identification and characterization of coherent structures as basic building blocks of turbulence. It has been a challenging task, however, to detect and visualize these coherent states. We address, by means of stereoscopic particle image velocimetry, that issue with the help of a large diameter (6 inches) pipe loop, which allowed us to probe for coherent states at various moderate Reynolds numbers (5300 < Re < 29000)). Although these states have been observed at flow regimes around laminar-turbulent transition (Re $\approx$ 2300) and also at high Reynolds number pipe flow (Re $\approx$ 35000), at moderate Reynolds numbers their existence had not been observed yet by experiment. By conditionally averaging the flow fields with respect to their dominant azimuthal wavenumber of streamwise velocity streaks, we have been able to uncover the existence of ten well-defined coherent flow patterns. It turns out, as a remarkable phenomenon, that their occurrence probabilities and the total number of dominant modes do not essentially change as the Reynolds number is varied. Their occurrence probabilities are noted to be reasonably well described by a Poisson distribution, which suggests that low-speed streaks are created as a Poisson process on the pipe circular geometry.
△ Less
Submitted 26 January, 2023;
originally announced January 2023.
-
Stochastic Model of Organizational State Transitions in a Turbulent Pipe Flow
Authors:
Robert Jäckel,
Bruno Magacho,
Bayode Owolabi,
Luca Moriconi,
David J. C. Dennis,
Juliana B. R. Loureiro
Abstract:
Turbulent pipe flows exhibit organizational states (OSs) that are labelled by discrete azimuthal wavenumber modes and are reminiscent of the traveling wave solutions of low Reynolds number regimes. The discretized time evolution of the OSs, obtained through stereoscopic particle image velocimetry, is shown to be non-Markovian for data acquisition carried out at a structure-resolved sampling rate.…
▽ More
Turbulent pipe flows exhibit organizational states (OSs) that are labelled by discrete azimuthal wavenumber modes and are reminiscent of the traveling wave solutions of low Reynolds number regimes. The discretized time evolution of the OSs, obtained through stereoscopic particle image velocimetry, is shown to be non-Markovian for data acquisition carried out at a structure-resolved sampling rate. In particular, properly defined time-correlation functions for the OS transitions are observed to decay as intriguing power laws, up to a large-eddy time horizon, beyond which they decorrelate at much faster rates. We are able to establish, relying upon a probabilistic description of the creation and annihilation of streamwise streaks, a lower-level {\it{Markovian}} model for the OS transitions, which reproduces their time-correlated behavior with meaningful accuracy. These findings indicate that the OSs are distributed along the pipe as statistically correlated packets of quasi-streamwise vortical structures.
△ Less
Submitted 12 January, 2023;
originally announced January 2023.
-
Modelling Cosmic Radiation Events in the Tree-ring Radiocarbon Record
Authors:
Qingyuan Zhang,
Utkarsh Sharma,
Jordan A. Dennis,
Andrea Scifo,
Margot Kuitems,
Ulf Buentgen,
Mathew J. Owens,
Michael W. Dee,
Benjamin J. S. Pope
Abstract:
Annually-resolved measurements of the radiocarbon content in tree-rings have revealed rare sharp rises in carbon-14 production. These 'Miyake events' are likely produced by rare increases in cosmic radiation from the Sun or other energetic astrophysical sources. The radiocarbon produced is not only circulated through the Earth's atmosphere and oceans, but also absorbed by the biosphere and locked…
▽ More
Annually-resolved measurements of the radiocarbon content in tree-rings have revealed rare sharp rises in carbon-14 production. These 'Miyake events' are likely produced by rare increases in cosmic radiation from the Sun or other energetic astrophysical sources. The radiocarbon produced is not only circulated through the Earth's atmosphere and oceans, but also absorbed by the biosphere and locked in the annual growth rings of trees. To interpret high-resolution tree-ring radiocarbon measurements therefore necessitates modelling the entire global carbon cycle. Here, we introduce 'ticktack', the first open-source Python package that connects box models of the carbon cycle with modern Bayesian inference tools. We use this to analyse all public annual 14C tree data, and infer posterior parameters for all six known Miyake events. They do not show a consistent relationship to the solar cycle, and several display extended durations that challenge either astrophysical or geophysical models.
△ Less
Submitted 25 October, 2022;
originally announced October 2022.
-
Designing Stress-Adaptive Dense Suspensions using Dynamic Covalent Chemistry
Authors:
Grayson L. Jackson,
Joseph M. Dennis,
Neil D. Dolinski,
Michael van der Naald,
Hojin Kim,
Christopher Eom,
Stuart J. Rowan,
Heinrich M. Jaeger
Abstract:
The non-Newtonian behaviors of dense suspensions are central to their use in technological and industrial applications and arise from a network of particle-particle contacts that dynamically adapts to imposed shear. Reported herein are studies aimed at explor-ing how dynamic covalent chemistry between particles and the polymeric solvent can be used to tailor such stress-adaptive contact networks l…
▽ More
The non-Newtonian behaviors of dense suspensions are central to their use in technological and industrial applications and arise from a network of particle-particle contacts that dynamically adapts to imposed shear. Reported herein are studies aimed at explor-ing how dynamic covalent chemistry between particles and the polymeric solvent can be used to tailor such stress-adaptive contact networks leading to their unusual rheological behaviors. Specifically, a room temperature dynamic thia-Michael bond is employed to rationally tune the equilibrium constant (Keq) of the polymeric solvent to the particle interface. It is demonstrated that low Keq leads to shear thinning while high Keq produces antithixotropy, a rare phenomenon where the viscosity increases with shearing time. It is proposed that an increase in Keq increases the polymer graft density at the particle surface and that antithixotropy primar-ily arises from partial debonding of the polymeric graft/solvent from the particle surface and the formation of polymer bridges between particles. Thus, the implementation of dynamic covalent chemistry provides a new molecular handle with which to tailor the macroscopic rheology of suspensions by introducing programable time dependence. These studies open the door to energy absorbing materials that not only sense mechanical inputs and adjust their dissipation as a function of time or shear rate but can switch between these two modalities on demand.
△ Less
Submitted 23 March, 2022;
originally announced March 2022.
-
Triangulating Instrumental Variable, confounder adjustment and Difference-in-Difference methods for comparative effectiveness research in observational data
Authors:
Laura Güdemann,
John M. Dennis,
Andrew P. McGovern,
Lauren R. Rodgers,
Beverley M. Shields,
William Henley,
Jack Bowden
Abstract:
Observational studies can play a useful role in assessing the comparative effectiveness of competing treatments. In a clinical trial the randomization of participants to treatment and control groups generally results in well-balanced groups with respect to possible confounders, which makes the analysis straightforward. However, when analysing observational data, the potential for unmeasured confou…
▽ More
Observational studies can play a useful role in assessing the comparative effectiveness of competing treatments. In a clinical trial the randomization of participants to treatment and control groups generally results in well-balanced groups with respect to possible confounders, which makes the analysis straightforward. However, when analysing observational data, the potential for unmeasured confounding makes comparing treatment effects much more challenging. Causal inference methods such as Instrumental Variable and Prior Even Rate Ratio approaches make it possible to circumvent the need to adjust for confounding factors that have not been measured in the data or measured with error. Direct confounder adjustment via multivariable regression and Propensity score matching also have considerable utility. Each method relies on a different set of assumptions and leverages different aspects of the data. In this paper, we describe the assumptions of each method and assess the impact of violating these assumptions in a simulation study. We propose the prior outcome augmented Instrumental Variable method that leverages data from before and after treatment initiation, and is robust to the violation of key assumptions. Finally, we propose the use of a heterogeneity statistic to decide if two or more estimates are statistically similar, taking into account their correlation. We illustrate our causal framework to assess the risk of genital infection in patients prescribed Sodium-glucose co-transporter-2 inhibitors versus Dipeptidyl peptidase-4 inhibitors as second-line treatment for Type 2 Diabets using observational data from the Clinical Practice Research Datalink.
△ Less
Submitted 27 July, 2022; v1 submitted 18 February, 2022;
originally announced February 2022.
-
C. H. Robinson Uses Heuristics to Solve Rich Vehicle Routing Problems
Authors:
Ehsan Khodabandeh,
Lawrence V. Snyder,
John Dennis,
Joshua Hammond,
Cody Wanless
Abstract:
We consider a wide family of vehicle routing problem variants with many complex and practical constraints, known as rich vehicle routing problems, which are faced on a daily basis by C.H. Robinson (CHR). Since CHR has many customers, each with distinct requirements, various routing problems with different objectives and constraints should be solved. We propose a set partitioning framework with a n…
▽ More
We consider a wide family of vehicle routing problem variants with many complex and practical constraints, known as rich vehicle routing problems, which are faced on a daily basis by C.H. Robinson (CHR). Since CHR has many customers, each with distinct requirements, various routing problems with different objectives and constraints should be solved. We propose a set partitioning framework with a number of route generation algorithms, which have shown to be effective in solving a variety of different problems. The proposed algorithms have outperformed the existing technologies at CHR on 10 benchmark instances and since, have been embedded into the company's transportation planning and execution technology platform.
△ Less
Submitted 30 December, 2019;
originally announced December 2019.
-
Implementation of PFC and RCM for RoCEv2 Simulation in OMNeT++
Authors:
Qian Liu,
Robert D. Russell,
Fabrice Mizero,
Malathi Veeraraghavan,
John Dennis,
Benjamin Jamroz
Abstract:
As traffic patterns and network topologies become more and more complicated in current enterprise data centers and TOP500 supercomputers, the probability of network congestion increases. If no countermeasures are taken, network congestion causes long communication delays and degrades network performance. A congestion control mechanism is often provided to reduce the consequences of congestion. How…
▽ More
As traffic patterns and network topologies become more and more complicated in current enterprise data centers and TOP500 supercomputers, the probability of network congestion increases. If no countermeasures are taken, network congestion causes long communication delays and degrades network performance. A congestion control mechanism is often provided to reduce the consequences of congestion. However, it is usually difficult to configure and activate a congestion control mechanism in production clusters and supercomputers due to concerns that it may negatively impact jobs if the mechanism is not appropriately configured. Therefore, simulations for these situations are necessary to identify congestion points and sources, and more importantly, to determine optimal settings that can be utilized to reduce congestion in those complicated networks. In this paper, we use OMNeT++ to implement the IEEE 802.1Qbb Priority-based Flow Control (PFC) and RoCEv2 Congestion Management (RCM) in order to simulate clusters with RoCEv2 interconnects.
△ Less
Submitted 11 September, 2015;
originally announced September 2015.
-
iSeg: an algorithm for segmentation of genomic data
Authors:
S. B. Girimurugan,
Jonathan Dennis,
Jinfeng Zhang
Abstract:
Identification of functional elements of a genome often requires dividing a sequence of measurements along a genome into segments differing from adjacent segments. In many applications, the mean of the measured values at multiple genomic locations in a segment is used to make inference of the property of interest. The segments with non-zero means often correspond to genomic regions with certain bi…
▽ More
Identification of functional elements of a genome often requires dividing a sequence of measurements along a genome into segments differing from adjacent segments. In many applications, the mean of the measured values at multiple genomic locations in a segment is used to make inference of the property of interest. The segments with non-zero means often correspond to genomic regions with certain biological events, such as changes between two conditions. This problem is often called the segmentation problem in the field of genomics, and the change-point problem in other scientific disciplines. We designed an efficient algorithm, called iSeg, for segmentation of high-throughput genomic profiles. iSeg first utilizes dynamic programming to compute the significance for a large number of candidate segments. It then uses tree-based data structures to detect overlapping significant regions and update them simultaneously. Refinement and merging of significant segments are performed at the end to generate the final segmentation. We evaluate iSeg using both simulated and experimental datasets and show that it performs quite well when compared with existing methods.
△ Less
Submitted 27 June, 2015;
originally announced June 2015.
-
Incorporating Kinetic Physics into a Two-Fluid Solar-Wind Model with Temperature Anisotropy and Low-Frequency Alfven-Wave Turbulence
Authors:
Benjamin D. G. Chandran,
Timothy J. Dennis,
Eliot Quataert,
Stuart D. Bale
Abstract:
We develop a 1D solar-wind model that includes separate energy equations for the electrons and protons, proton temperature anisotropy, collisional and collisionless heat flux, and an analytical treatment of low-frequency, reflection-driven, Alfven-wave turbulence. To partition the turbulent heating between electron heating, parallel proton heating, and perpendicular proton heating, we employ resul…
▽ More
We develop a 1D solar-wind model that includes separate energy equations for the electrons and protons, proton temperature anisotropy, collisional and collisionless heat flux, and an analytical treatment of low-frequency, reflection-driven, Alfven-wave turbulence. To partition the turbulent heating between electron heating, parallel proton heating, and perpendicular proton heating, we employ results from the theories of linear wave damping and nonlinear stochastic heating. We account for mirror and oblique firehose instabilities by increasing the proton pitch-angle scattering rate when the proton temperature anisotropy exceeds the threshold for either instability. We numerically integrate the equations of the model forward in time until a steady state is reached, focusing on two fast-solar-wind-like solutions. These solutions are consistent with a number of observations, supporting the idea that Alfven-wave turbulence plays an important role in the origin of the solar wind.
△ Less
Submitted 29 November, 2011; v1 submitted 13 October, 2011;
originally announced October 2011.
-
Magnetic Nested-wind Scenarios for Bipolar Outflows: Pre-planetary and YSO nebular shaping
Authors:
Timothy J. Dennis,
Adam Frank,
Eric G. Blackman,
Orsola DeMarco,
Bruce Balick,
Sorin Mitran
Abstract:
We present results of a series of magnetohydrodynamic (MHD) and hydro- dynamic (HD) 2.5D simulations of the morphology of outflows driven by nested wide-angle winds - i.e. winds which eminate from a central star as well as from an orbiting accretion disk. While our results are broadly relevent to nested wind systems we have tuned the parameters of the simulations to touch on issues in both Young…
▽ More
We present results of a series of magnetohydrodynamic (MHD) and hydro- dynamic (HD) 2.5D simulations of the morphology of outflows driven by nested wide-angle winds - i.e. winds which eminate from a central star as well as from an orbiting accretion disk. While our results are broadly relevent to nested wind systems we have tuned the parameters of the simulations to touch on issues in both Young Stellar Objects and Planetary Nebula studies. In particular our studies connect to open issues in the early evolution of Planetary Nebulae. We find that nested MHD winds exhibit marked morphological differences from the single MHD wind case along both dimensions of the flow. Nested HD winds on the other hand give rise mainly to geometric distortions of an outflow that is topologically similar to the flow arising from a single stellar HD wind. Our MHD results are insensitive to changes in ambient temperature between ionized and un-ionized circumstellar environments. The results are sensitive to the relative mass-loss rates, and to the relative speeds of the stellar and disk winds. We also present synthetic emission maps of both nested MHD and HD simulations. We find that nested MHD winds show knots of emission appearing on-axis that do not appear in the HD case.
△ Less
Submitted 6 September, 2009;
originally announced September 2009.
-
Hyperfine structure of Sc@C82 from ESR and DFT
Authors:
G. W. Morley,
B. J. Herbert,
S. M. Lee,
K. Porfyrakis,
T. J. S. Dennis,
D. Nguyen-Manh,
R. Scipioni,
J. van Tol,
A. P. Horsfield,
A. Ardavan,
D. G. Pettifor,
J. C. Green,
G. A. D. Briggs
Abstract:
The electron spin g- and hyperfine tensors of the endohedral metallofullerene Sc@C82 are anisotropic. Using electron spin resonance (ESR) and density functional theory (DFT), we can relate their principal axes to the coordinate frame of the molecule, finding that the g-tensor is not axially symmetric. The Sc bond with the cage is partly covalent and partly ionic. Most of the electron spin densit…
▽ More
The electron spin g- and hyperfine tensors of the endohedral metallofullerene Sc@C82 are anisotropic. Using electron spin resonance (ESR) and density functional theory (DFT), we can relate their principal axes to the coordinate frame of the molecule, finding that the g-tensor is not axially symmetric. The Sc bond with the cage is partly covalent and partly ionic. Most of the electron spin density is distributed around the carbon cage, but 5% is associated with the scandium d_yz orbital, and this drives the observed anisotropy.
△ Less
Submitted 11 September, 2008;
originally announced September 2008.
-
Parker/buoyancy instabilities with anisotropic thermal conduction, cosmic rays, and arbitrary magnetic field strength
Authors:
Timothy J. Dennis,
Benjamin D. G. Chandran
Abstract:
We report the results of a local stability analysis for a magnetized, gravitationally stratified plasma containing cosmic rays. We account for cosmic-ray diffusion and thermal conduction parallel to the magnetic field and allow beta to take any value, where p is the plasma pressure and B is the magnetic field strength. We take the gravitational acceleration to be in the -z-direction and the equi…
▽ More
We report the results of a local stability analysis for a magnetized, gravitationally stratified plasma containing cosmic rays. We account for cosmic-ray diffusion and thermal conduction parallel to the magnetic field and allow beta to take any value, where p is the plasma pressure and B is the magnetic field strength. We take the gravitational acceleration to be in the -z-direction and the equilibrium magnetic field to be in the y-direction, and we derive the dispersion relation for small-amplitude instabilities and waves in the large-|k_x| limit. We use the Routh-Hurwitz criterion to show analytically that the necessary and sufficient criterion for stability in this limit is n k_B dT/dz + dp_cr/dz + (1/8pi)dB^2/dz > 0, where T is the temperature, n is the number density of thermal particles, and p_cr is the cosmic-ray pressure. We present approximate analytical solutions for the normal modes in the low- and high-diffusivity limits, show that they are consistent with the derived stability criterion, and compare them to numerical results obtained from the full, unapproximated, dispersion relation. Our results extend earlier analyses of buoyancy instabilities in galaxy-cluster plasmas to the beta <= 1 regime. Our results also extend earlier analyses of the Parker instability to account for anisotropic thermal conduction, and show that the interstellar medium is more unstable to the Parker instability than was predicted by previous studies in which the thermal plasma was treated as adiabatic.
△ Less
Submitted 3 September, 2008;
originally announced September 2008.
-
PPN as Explosions: Bullets vs Jets and Nebular Shaping
Authors:
Timothy J. Dennis,
Andrew J. Cunningham,
Adam Frank,
Bruce Balick,
Eric G. Blackman,
Sorin Mitran
Abstract:
Many proto-planetary nebulae (PPN) appear as narrow collimated structures sometimes showing multiple, roughly aligned lobes. In addition, many PPN flows have been shown to have short acceleration times. In this paper we explore whether jet or ``bullet'' (a massive clump) models fit the observations of individual collimated lobes adequately by comparing simulations of both radiatively cooled (sta…
▽ More
Many proto-planetary nebulae (PPN) appear as narrow collimated structures sometimes showing multiple, roughly aligned lobes. In addition, many PPN flows have been shown to have short acceleration times. In this paper we explore whether jet or ``bullet'' (a massive clump) models fit the observations of individual collimated lobes adequately by comparing simulations of both radiatively cooled (stable) jets and bullets. We find that the clump model is somewhat favored over jets because (1) it leads to greater collimation of outflows (2) it accounts better and more naturally for ring-like structures observed in the PPN CRL 618, and (3) it is more successful in reproducing the Hubble-flow character of observed kinematics in some PPN. In addition, bullets naturally account for observed multipolar flows, since the likely MHD launch mechanisms required to drive outflows make multiple non-aligned jets unlikely. Thus we argue that PPN outflows may be driven by explosive MHD launch mechanisms such as those discussed in the context of supernovae (SNe) and gamma-ray bursts(GRB).
△ Less
Submitted 13 July, 2007; v1 submitted 11 July, 2007;
originally announced July 2007.
-
A Comparison of Two Shallow Water Models with Non-Conforming Adaptive Grids: classical tests
Authors:
Amik St-Cyr,
Christiane Jablonowski,
John M. Dennis,
Henry M. Tufo,
Stephen J. Thomas
Abstract:
In an effort to study the applicability of adaptive mesh refinement (AMR) techniques to atmospheric models an interpolation-based spectral element shallow water model on a cubed-sphere grid is compared to a block-structured finite volume method in latitude-longitude geometry. Both models utilize a non-conforming adaptation approach which doubles the resolution at fine-coarse mesh interfaces. The…
▽ More
In an effort to study the applicability of adaptive mesh refinement (AMR) techniques to atmospheric models an interpolation-based spectral element shallow water model on a cubed-sphere grid is compared to a block-structured finite volume method in latitude-longitude geometry. Both models utilize a non-conforming adaptation approach which doubles the resolution at fine-coarse mesh interfaces. The underlying AMR libraries are quad-tree based and ensure that neighboring regions can only differ by one refinement level.
The models are compared via selected test cases from a standard test suite for the shallow water equations. They include the advection of a cosine bell, a steady-state geostrophic flow, a flow over an idealized mountain and a Rossby-Haurwitz wave. Both static and dynamics adaptations are evaluated which reveal the strengths and weaknesses of the AMR techniques. Overall, the AMR simulations show that both models successfully place static and dynamic adaptations in local regions without requiring a fine grid in the global domain. The adaptive grids reliably track features of interests without visible distortions or noise at mesh interfaces. Simple threshold adaptation criteria for the geopotential height and the relative vorticity are assessed.
△ Less
Submitted 15 February, 2007;
originally announced February 2007.