Source author record

Eric Wolf

Eric Wolf appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.EP math.NA astro-ph.IM astro-ph.SR

Catalog footprint

What is connected

4works

4topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

ExoCAM: A 3D Climate Model for Exoplanet Atmospheres

The TRAPPIST-1 Habitable Atmosphere Intercomparison (THAI) project was initiated to compare 3D climate models that are commonly used for predicting theoretical climates of habitable zone extrasolar planets. One of the core models studied as part of THAI is ExoCAM, an independently curated exoplanet branch of the National Center for Atmospheric Research (NCAR) Community Earth System Model (CESM) version 1.2.1. ExoCAM has been used for studying atmospheres of terrestrial extrasolar planets around a variety of stars. To accompany the THAI project and provide a primary reference, here we describe ExoCAM and what makes it unique from standard configurations of CESM. Furthermore, we also conduct a series of intramodel sensitivity tests of relevant moist physical tuning parameters while using the THAI protocol as our starting point. A common criticism of 3D climate models used for exoplanet modeling is that cloud and convection routines often contain free parameters that are tuned to the modern Earth, and thus may be a source of uncertainty in evaluating exoplanet climates. Here, we explore sensitivities to numerous configuration and parameter selections, including a recently updated radiation scheme, a different cloud and convection physics package, different cloud and precipitation tuning parameters, and a different sea ice albedo. Improvements to our radiation scheme and the modification of cloud particle sizes have the largest effect on global mean temperatures, with variations up to ~10 K, highlighting the requirement for accurate radiative transfer and the importance of cloud microphysics for simulating exoplanetary climates. However for the vast majority of sensitivity tests, climate differences are small. For all cases studied, intramodel differences do not bias general conclusions regarding climate states and habitability.

preprint2019arXiv

Habitable Zone Boundaries for Circumbinary Planets

We use a one-dimensional (1-D) cloud-free climate model to estimate habitable zone (HZ) boundaries for terrestrial planets of masses 0.1 M$_{E}$ and 5 M$_{E}$ around circumbinary stars of various spectral type combinations. Specifically, we consider binary systems with host spectral types F-F, F-G, F-K, F-M, G-G, G-K, G-M, K-K, K-M and M-M. Scaling the background N2 atmospheric pressure with the radius of the planet, we find that the inner edge of the HZ moves inwards towards the star for 5ME compared to 0.1ME planets for all spectral types. This is because the water-vapor column depth is smaller for larger planets and higher temperatures are needed before water vapor completely dominates the outgoing longwave radiation. The outer edge of the HZ changes little due to competing effects of the albedo and greenhouse effect. While these results are broadly consistent with the trend of single star HZ results for different mass planets, there are significant differences between single star and binary star systems for the inner edge of the HZ. Interesting combinations of stellar pairs from our 1-D model results can be used to explore for in-depth climate studies with 3-D climate models. We identify a common HZ stellar flux domain for all circumbinary spectral types

preprint2015arXiv

Method of lines transpose: an efficient A-stable solver for wave propagation

Building upon recent results obtained in [7,8,9], we describe an efficient second order, A-stable scheme for solving the wave equation, based on the method of lines transpose (MOL$^T$), and the resulting semi-discrete (i.e. continuous in space) boundary value problem. In [7], A-stable schemes of high order were derived, and in [9] a high order, fast $\mathcal{O}(N)$ spatial solver was derived, which is matrix-free and is based on dimensional-splitting. In this work, are interested in building a wave solver, and our main concern is the development of boundary conditions. We demonstrate all desired boundary conditions for a wave solver, including outflow boundary conditions, in 1D and 2D. The scheme works in a logically Cartesian fashion, and the boundary points are embedded into the regular mesh, without incurring stability restrictions, so that boundary conditions are imposed without any reduction in the order of accuracy. We demonstrate how the embedded boundary approach works in the cases of Dirichlet and Neumann boundary conditions. Further, we develop outflow and periodic boundary conditions for the MOL$^T$ formulation. Our solver is designed to couple with particle codes, and so special attention is also paid to the implementation of point sources, and soft sources which can be used to launch waves into waveguides.

preprint2013arXiv

Method of Lines Transpose: A Fast Implicit Wave Propagator

As a follow up to \cite{Causley2013}, we provide a detailed description of the numerical implementation of an O(N), A-stable, second order accurate solution of the wave equation, constructed from semi-discrete boundary value problems. We improve on the previous algorithm by replacing the Lax-type correction used in \cite{Causley2013}, which was necessary for convergence when $Δt < Δx/c$, with a more accurate spatial quadrature, which we prove is convergent. We also demonstrate that the resulting solver remains fast even in the case of unstructured meshes, can incorporate domain decomposition, and allows for the implementation of Dirichlet, Neumann, periodic and outflow boundary conditions. Building upon results for the 1d formulation, we utilize alternate direction implicit (ADI) splitting to achieve a fast O(N) solver in higher spatial dimensions. Our solver is built upon line objects and, combined with the flexibility of the integral solver, allows us to solve problems on arbitrary spatial domains, by embedding the boundary in a regular Cartesian mesh. Our solver is designed to couple with particle codes, where scale separation is an issue. We therefore demonstrate the ability of our solver to take time steps well beyond that of the Courant-Friedrichs-Lewy (CFL) stability limit of explicit codes.