Uncertainty Quantification of geochemical and mechanical compaction in layered sedimentary basins

In this work we propose an Uncertainty Quantification methodology for sedimentary basins evolution under mechanical and geochemical compaction processes, which we model as a coupled, time-dependent, non-linear, monodimensional (depth-only) system of PDEs with uncertain parameters. While in previous works (Formaggia et al. 2013, Porta et al., 2014) we assumed a simplified depositional history with only one material, in this work we consider multi-layered basins, in which each layer is characterized by a different material, and hence by different properties. This setting requires several improvements with respect to our earlier works, both concerning the deterministic solver and the stochastic discretization. On the deterministic side, we replace the previous fixed-point iterative solver with a more efficient Newton solver at each step of the time-discretization. On the stochastic side, the multi-layered structure gives rise to discontinuities in the dependence of the state variables on the uncertain parameters, that need an appropriate treatment for surrogate modeling techniques, such as sparse grids, to be effective. We propose an innovative methodology to this end which relies on a change of coordinate system to align the discontinuities of the target function within the random parameter space. The reference coordinate system is built upon exploiting physical features of the problem at hand. We employ the locations of material interfaces, which display a smooth dependence on the random parameters and are therefore amenable to sparse grid polynomial approximations. We showcase the capabilities of our numerical methodologies through two synthetic test cases. In particular, we show that our methodology reproduces with high accuracy multi-modal probability density functions displayed by target state variables (e.g., porosity).Comment: 25 pages, 30 figure

Distinctive rings in the 21 cm signal of the epoch of reionization

It is predicted that sources emitting UV radiation in the Lyman band during the epoch of reionization (EoR) showed a series of discontinuities in their Ly-alpha flux radial profile as a consequence of the thickness of the Lyman line series in the primeval intergalactic medium. Through unsaturated Wouthuysen-Field coupling, these spherical discontinuities are also present in the 21 cm emission of the neutral IGM. In this article, we study the effects these discontinuities have on the differential brightness temperature of the 21 cm signal of neutral hydrogen in a realistic setting including all other sources of fluctuations. We focus on the early phases of the EoR, and we address the question of the detectability by the planned Square Kilometre Array. Such a detection would be of great interest, because these structures could provide an unambiguous diagnostic for the cosmological origin of the signal remaining after the foreground cleaning procedure. Also, they could be used as a new type of standard rulers. We determine the differential brightness temperature of the 21 cm signal in the presence of inhomogeneous Wouthuysen-Field effect using simulations which include (hydro)dynamics and both ionizing and Lyman lines 3D radiative transfer with the code LICORICE. We find that the Lyman horizons are clearly visible on the maps and radial profiles around the first sources of our simulations, but for a limited time interval, typically \Delta z \approx 2 at z \sim 13. Stacking the profiles of the different sources of the simulation at a given redshift results in extending this interval to \Delta z \approx 4. When we take into account the implementation and design planned for the SKA (collecting area, sensitivity, resolution), we find that detection will be challenging. It may be possible with a 10 km diameter for the core, but will be difficult with the currently favored design of a 5 km core.Comment: 10 pages, 10 figures; v2: Section 5.5 rewritten; some new references added; accepted for publication in Astronomy and Astrophysic

Effective Scalar Field Theory for the Electroweak Phase Transition

We investigate an effective model for the finite temperature restoration phase transition of the electroweak theory. It is obtained by dimensional reduction of the $3+1$ dimensional full theory and by subsequent integration over all static gauge degrees of freedom. The resulting theory corresponds to a $3$-dimensional $O(4)$ ferromagnet containing cubic and quartic terms of the field in its potential function. Possible nonperturbative effects of a magnetic screening mass are parametrically included in the potential. We analyse the theory using mean field and numerical Monte Carlo (MC) simulation methods. At the value of the physical Higgs mass, $m_H=37~{\rm GeV}$, considered in the present investigation, we find a discontinuous symmetry restoring phase transition. We determine the critical temperature, order parameter jump, interface tension and latent heat characteristics of the transition. The Monte Carlo results indicate a somewhat weaker first order phase transition as compared to the mean field treatment, demonstrating that non-perturbative fluctuations of the Higgs field are relevant. This effect is especially important for the interface tension. Any observation of hard first order transition could result only from non-perturbative effects related to the gauge degrees of freedom.Comment: 28 pages and 18 figure

Review of high-contrast imaging systems for current and future ground- and space-based telescopes I. Coronagraph design methods and optical performance metrics

The Optimal Optical Coronagraph (OOC) Workshop at the Lorentz Center in September 2017 in Leiden, the Netherlands gathered a diverse group of 25 researchers working on exoplanet instrumentation to stimulate the emergence and sharing of new ideas. In this first installment of a series of three papers summarizing the outcomes of the OOC workshop, we present an overview of design methods and optical performance metrics developed for coronagraph instruments. The design and optimization of coronagraphs for future telescopes has progressed rapidly over the past several years in the context of space mission studies for Exo-C, WFIRST, HabEx, and LUVOIR as well as ground-based telescopes. Design tools have been developed at several institutions to optimize a variety of coronagraph mask types. We aim to give a broad overview of the approaches used, examples of their utility, and provide the optimization tools to the community. Though it is clear that the basic function of coronagraphs is to suppress starlight while maintaining light from off-axis sources, our community lacks a general set of standard performance metrics that apply to both detecting and characterizing exoplanets. The attendees of the OOC workshop agreed that it would benefit our community to clearly define quantities for comparing the performance of coronagraph designs and systems. Therefore, we also present a set of metrics that may be applied to theoretical designs, testbeds, and deployed instruments. We show how these quantities may be used to easily relate the basic properties of the optical instrument to the detection significance of the given point source in the presence of realistic noise.Comment: To appear in Proceedings of the SPIE, vol. 1069

Mechanical testing of metallic foams for 3d model and simulation of cell distribution effects

Cellular materials have a bulk matrix with a larger number of voids named also cells. Metallic foams made by powder technology represent stochastic closed cells. The related inhomogeneity leads to a scattering of results both in terms of stress–strain curves and maximum strength. Scattering is attributed to relative density variations and local cell discontinuities and it is confirmed also in case of dynamic loading. Finite element simulations through geometrical models that are able to capture the void morphology (named “mesoscale models”), confirm these results and some efforts have been already done to quantify the relationship between shape irregularities and mechanical behavior. The aim of this paper is to present the dynamic characterization of an AA7075 closed cell material and to calibrate its mesoscale finite element model according to the related cell shape distribution. Specimens have been derived from a small ingot (45x45x100 mm) divided along sections so that morphological analysis and experimental tests have been carried out. Specimens extracted from a half of the ingot have been used for dynamic compression tests by means of a split Hopkinson bar, meanwhile specimens extracted from the other half of the ingot have been dissected for porosity distribution analyses carried out by means of image analysis. Stress-strain curves obtained from the mechanical tests have been discussed in terms of strain rate and statistical descriptors of the porosity. Successively a 3D-model of the specimen has been generated starting from the Voronoi algorithm, assigning as input the above-mentioned statistical distribution of the porosity. Due to the peculiarity of the cell morphology (e.g. single larger cells), stress-strain localization has been demonstrated as one of the reasons of the scattering found during the experiments. A material model, to reproduce the investigated foam mechanical behavior, has been calibrated. Despite the difference among experiments the material model is able to reproduce all of them. Difference between the model coefficients quantifies roughly the difference due to the local geometry of the cells

First Observational Tests of Eternal Inflation: Analysis Methods and WMAP 7-Year Results

In the picture of eternal inflation, our observable universe resides inside a single bubble nucleated from an inflating false vacuum. Many of the theories giving rise to eternal inflation predict that we have causal access to collisions with other bubble universes, providing an opportunity to confront these theories with observation. We present the results from the first observational search for the effects of bubble collisions, using cosmic microwave background data from the WMAP satellite. Our search targets a generic set of properties associated with a bubble collision spacetime, which we describe in detail. We use a modular algorithm that is designed to avoid a posteriori selection effects, automatically picking out the most promising signals, performing a search for causal boundaries, and conducting a full Bayesian parameter estimation and model selection analysis. We outline each component of this algorithm, describing its response to simulated CMB skies with and without bubble collisions. Comparing the results for simulated bubble collisions to the results from an analysis of the WMAP 7-year data, we rule out bubble collisions over a range of parameter space. Our model selection results based on WMAP 7-year data do not warrant augmenting LCDM with bubble collisions. Data from the Planck satellite can be used to more definitively test the bubble collision hypothesis.Comment: Companion to arXiv:1012.1995. 41 pages, 23 figures. v2: replaced with version accepted by PRD. Significant extensions to the Bayesian pipeline to do the full-sky non-Gaussian source detection problem (previously restricted to patches). Note that this has changed the normalization of evidence values reported previously, as full-sky priors are now employed, but the conclusions remain unchange

Quasi-static imaged-based immersed boundary-finite element model of human left ventricle in diastole

SUMMARY: Finite stress and strain analyses of the heart provide insight into the biomechanics of myocardial function and dysfunction. Herein, we describe progress toward dynamic patient-specific models of the left ventricle using an immersed boundary (IB) method with a finite element (FE) structural mechanics model. We use a structure-based hyperelastic strain-energy function to describe the passive mechanics of the ventricular myocardium, a realistic anatomical geometry reconstructed from clinical magnetic resonance images of a healthy human heart, and a rule-based fiber architecture. Numerical predictions of this IB/FE model are compared with results obtained by a commercial FE solver. We demonstrate that the IB/FE model yields results that are in good agreement with those of the conventional FE model under diastolic loading conditions, and the predictions of the LV model using either numerical method are shown to be consistent with previous computational and experimental data. These results are among the first to analyze the stress and strain predictions of IB models of ventricular mechanics, and they serve both to verify the IB/FE simulation framework and to validate the IB/FE model. Moreover, this work represents an important step toward using such models for fully dynamic fluid–structure interaction simulations of the heart

Space Station Freedom Beta Gimbal Control via Sensitivity Models

Tracking control of the Space Station Freedom solar array beta gimbals is investigated. Of particular interest is the issue of control in the presence of uncertainty in gimbal friction parameters. Sensitivity functions are incorporated into the feedback loop to desensitize the gimbal control law to parameter variations. Simulation results indicated that one such sensitivity function improves the closed-loop performance of the gimbals in the presence of unexpected friction parameter dispersions