scholarly journals Coupling ice flow models of varying orders of complexity with the Tiling method

2012 ◽  
Vol 58 (210) ◽  
pp. 776-786 ◽  
Author(s):  
Helene Seroussi ◽  
Hachmi Ben Dhia ◽  
Mathieu Morlighem ◽  
Eric Larour ◽  
Eric Rignot ◽  
...  
Keyword(s):  
Numerical Models ◽  
Computational Cost ◽  
Higher Order ◽  
Ice Flow ◽  
Flow Models ◽  
Continental Scale ◽  
Critical Areas ◽  
Warming Climate ◽  
Computationally Intensive ◽  
A New Technique

AbstractIce flow numerical models are essential for predicting the evolution of ice sheets in a warming climate. Recent research emphasizes the need for higher-order and even full-Stokes flow models, instead of the traditional shallow-ice approximation, whose assumptions are not valid in certain critical areas. These higher-order models are, however, computationally intensive and difficult to use at the continental scale. Here we present a new technique, the Tiling method, to couple ice flow models of varying orders of complexity. The goal of the method is to limit the spatial extent of where higherorder models are applied to reduce the computational cost, while maintaining the model precision. We apply this method on synthetic geometries to demonstrate its practical use. We first use a geometry for which all models yield the same results to check the consistency of the method. Then we apply our method to a geometry for which a full-Stokes model is required in the vicinity of the ice front. Our results show that the hybrid models present significant improvements over mono-model approaches and reduce computational times.

Investigating basal thaw as a potential driver of ice flow acceleration in Antarctica

2020 ◽  
Author(s):  
Eliza Dawson ◽  
Dustin Schroeder ◽  
Winnie Chu ◽  
Elisa Mantelli ◽  
Helene Seroussi
Keyword(s):  
Numerical Models ◽  
Ice Sheet ◽  
Ice Streams ◽  
Ice Flow ◽  
Flow Acceleration ◽  
Critical Areas ◽  
Potential Impact ◽  
Basal Melting ◽  
The Impact ◽  
Selection Of

<p>Glacial thermal processes exert a fundamental control on ice flow, governing viscosity and frozen-to-thawed transitions at the ice-bed interface. Across Antarctica, frozen bed regions characterized by numerical models and geophysical observations, can also reduce ice flow by increasing basal traction. Some frozen bed regions can separate or confine fast-flowing glaciers and ice streams. Others separate inland catchments with thawed beds from the grounding zone of marine ice-sheet sectors. If regions with frozen bed experienced thawing, such a transition may lead to ice-sheet acceleration, reconfiguration, or retreat. To investigate the potential impact of such a thermal transition, we use the JPL/UCI Ice Sheet System Model (ISSM) to identify vulnerable regions across Antarctica that are close to the basal melting point. We assess the impact of thawing these regions by quantifying resulting volume changes and surface expressions. This allows us to identify the areas of the ice sheet where the thermal regime at the ice-bed interface has the largest potential impact on ice-sheet stability and sea-level contribution. We also examine the potential basal temperature and thaw-propagation thresholds governing this process. We then compare the ISSM results to a selection of ice-penetrating radar sounding observations to refine our constraints of the configuration, distribution, and extent of these thermally critical areas.</p>

scholarly journals Efficiently Simulating an Endograft Deployment: A Methodology for Detailed CFD Analyses

2020 ◽  
Vol 48 (10) ◽  
pp. 2449-2465
Author(s):  
Faidon Kyriakou ◽  
Craig Maclean ◽  
William Dempster ◽  
David Nash
Keyword(s):  
Numerical Models ◽  
Flow Analysis ◽  
Endovascular Aneurysm Repair ◽  
Computational Cost ◽  
Element Analysis ◽  
Fea Model ◽  
Order Of Magnitude ◽  
Modelling Methodology ◽  
Computationally Intensive ◽  
Cfd Analyses

Abstract Numerical models of endografts for the simulation of endovascular aneurysm repair are increasingly important in the improvement of device designs and patient outcomes. Nevertheless, current finite element analysis (FEA) models of complete endograft devices come at a high computational cost, requiring days of runtime, therefore restricting their applicability. In the current study, an efficient FEA model of the Anaconda™ endograft (Terumo Aortic, UK) was developed, able to yield results in just over 4 h, an order of magnitude less than similar models found in the literature. The model was used to replicate a physical device that was deployed in a 3D printed aorta and comparison of the two shapes illustrated a less than 5 mm placement error of the model in the regions of interest, consistent with other more computationally intensive models in the literature. Furthermore, the final goal of the study was to utilize the deployed fabric model in a hemodynamic analysis that would incorporate realistic fabric folds, a feature that is almost always omitted in similar simulations. By successfully exporting the deployed graft geometry into a flow analysis, it was illustrated that the inclusion of fabric wrinkles enabled clinically significant flow patterns such as flow stagnation and recirculation to be detected, paving the way for this modelling methodology to be used in future for stent design optimisation.

scholarly journals Hydrodynamic Simulation of an Irregularly Meandering Gravel-Bed River: Comparison of MIKE 21 FM and Delft3D Flow models

2018 ◽  
Vol 40 ◽  
pp. 02004 ◽  
Author(s):  
Parna Parsapour-moghaddam ◽  
Colin D. Rennie ◽  
Jonathan Slaney
Keyword(s):  
Flow Model ◽  
Numerical Models ◽  
Computational Cost ◽  
Grid Cell ◽  
Flood Management ◽  
Water Levels ◽  
Test Analysis ◽  
Flow Models ◽  
Gravel Bed ◽  
Mike 21

This study aims at hydrodynamic modelling of Bow River, which passes through the City of Calgary, Canada. Bow River has a mobile gravel bed. Erosion and deposition processes were exacerbated by a catastrophic flood in 2013. Channel banks were eroded at various locations, and large gravel bars formed, which could lead to water level changes and accordingly more flooding. This study investigates the performance of Delft3D-Flow and MIKE 21 FM to simulate the hydrodynamics of the river during the 2013 flood. MIKE 21FM employs unstructured triangular mesh while Delft3D-Flow model uses curvilinear structured grids. Performance of each model was evaluated by the available historical water levels. The results of this study demonstrated that, with approximately the same averaged grid resolution, MIKE 21 FM resulted in more accurate results with a higher computational cost compared to the Delft3DFlow model. It was shown that Delft3D-Flow model may require higher grid cell resolution to result in comparably same depth-averaged velocities throughout the study area. However, considering the balance between the computational cost and the accuracy of the results, both models were capable to adequately replicate the hydrodynamics of the river during the 2013 flood. Results of statistical KS and ANOVA test analysis showed that the model predictions were sensitive to the horizontal eddy viscosity and the Manning roughness. This confirms the necessity of an appropriate calibration of the generated numerical models. The findings of this study shed light on the Bow River flood modelling, which can guide flood management.

scholarly journals Sensitivity of Greenland ice sheet projections to spatial resolution in higher-order simulations: the AWI contribution to ISMIP6-Greenland using ISSM

2020 ◽  
Author(s):  
Martin Rückamp ◽  
Heiko Goelzer ◽  
Angelika Humbert
Keyword(s):  
Sea Level Rise ◽  
Sea Level ◽  
Ice Sheet ◽  
Greenland Ice Sheet ◽  
Higher Order ◽  
Initial State ◽  
Ice Flow ◽  
Outlet Glacier ◽  
Flow Models ◽  
High Emission

Abstract. Projections of the contribution of the Greenland ice sheet to future sea-level rise include uncertainties primarily due to the imposed climate forcing and the initial state of the ice sheet model. Several state-of-the-art ice flow models are currently being employed on various grid resolutions to estimate future mass changes in the framework of the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6). Here we investigate the sensitivity to grid resolution on centennial sea-level contributions from the Greenland ice sheet and study the mechanism at play. To this end, we employ the finite-element higher-order ice flow model ISSM and conduct experiments with four different horizontal resolutions, namely 4, 2, 1 and 0.75 km. We run the simulation based on the ISMIP6 core GCM MIROC5 under the high emission scenario RCP8.5 and consider both atmospheric and oceanic forcing in full and separate scenarios. Under the full scenarios, finer simulations unveil up to ~5 % more sea-level rise compared to the coarser resolution. The sensitivity depends on the magnitude of outlet glacier retreat, which is implemented as a series of retreat masks following the ISMIP6 protocol. Without imposed retreat under atmosphere-only forcing, the resolution dependency exhibits an opposite behaviour with about ~ 5 % more sea-level contribution in the coarser resolution. The sea-level contribution indicates a converging behaviour

PoLIM: an open source 2D higher-order thermomechanically coupled mountain glacier flow model

2020 ◽  
Author(s):  
Yuzhe Wang ◽  
Tong Zhang
Keyword(s):  
Flow Model ◽  
Order Approximation ◽  
Water Pressure ◽  
Computational Cost ◽  
Higher Order ◽  
Thermomechanical Coupling ◽  
High Mountain ◽  
Climate Data ◽  
Mountain Glacier ◽  
Ice Flow

<p>The worldwide glacier is retreating and is expected to continue shrinking in a warming climate. Understanding the dynamics of glaciers is essential for the knowledge of sea-level rise, water resources in high mountain and arid regions, and the potential glacier hazards. Over the past decades, various 3D higher-order and full-Stokes ice flow models including thermomechanical coupling have been developed, and some have opened their source codes. However, such 3D modeling requires detailed datasets about surface and bedrock topography, variable climatic conditions, and high computational cost. Due to difficulties in measuring glacier thickness, only a small minority of glaciers around the globe have ice thickness observations. It is also a challenge to downscale the climate data (e.g., air temperature, precipitation) to the glacier surface, particularly, in rugged high-mountain terrains. In contrast to 3D models, flowline models only require inputs along the longitudinal profile and are thus computationally efficient. They continue to be useful tools for simulating the evolution of glaciers and studying the particular phenomena related to glacier dynamics. In this study, we present a two-dimensional thermomechanically coupled ice flow model named PoLIM (Polythermal Land Ice Model). The velocity solver of PoLIM is developed based on the higher-order approximation (Blatter-Pattyn type). It includes three critical features for simulating the dynamics of mountain glaciers: 1) an enthalpy-based thermal model to describe the heat transfer, which is particularly convenient to simulate the polythermal structures; 2) a drainage model to simulate the water transport in the temperate ice layer driven by gravity; 3) a subglacial hydrology model to simulate the subglacial water pressure for the coupling with the basal sliding law. We verify PoLIM with several standard benchmark experiments (e.g., ISMIP-HOM, enthalpy, SHMIP) in the glacier modeling community. PoLIM shows a good performance and agrees well with these benchmark results, indicating its reliable and robust capability of simulating the thermomechanical behaviors of glaciers.</p>

scholarly journals Using the unstable manifold correction in a Picard iteration to solve the velocity field in higher-order ice-flow models

2010 ◽  
Vol 56 (196) ◽  
pp. 257-261 ◽  
Author(s):  
Bert De Smedt ◽  
Frank Pattyn ◽  
Pieter De Groen
Keyword(s):  
Velocity Field ◽  
Unstable Manifold ◽  
Real Data ◽  
Higher Order ◽  
Picard Iteration ◽  
Ice Flow ◽  
Flow Models ◽  
Original Algorithm ◽  
Fast Solution ◽  
Higher Order Models

AbstractWe address the usefulness of the unstable manifold correction (UMC) in a Picard iteration for the solution of the velocity field in higher-order ice-flow models. We explain under- and overshooting and how one can remedy them. We then discuss the rationale behind the UMC, initially developed to remedy overshooting, and how it was previously introduced in a Picard iteration to calculate the velocity field in higher-order models. Using a laminar-flow experiment with two higher-order model implementations, we demonstrate that it is not overshooting, but undershooting, that is the main problem when using a proper implementation for the calculation of the velocity field in higher-order models. We also consider a variant of the original UMC algorithm that often enables a relatively fast solution, but is theoretically less sound. Therefore, neither the variant nor the original algorithm is suited for these problems. We present a more appropriate, stable and simple relaxed Picard algorithm and demonstrate that, compared to the true Picard iteration and the variant of the original UMC algorithm, it results in a faster solution of the velocity field in higher-order models for problems with real data.

scholarly journals Ice flow models and glacial erosion over multiple glacial–interglacial cycles

2014 ◽  
Vol 2 (1) ◽  
pp. 389-428
Author(s):  
R. M. Headley ◽  
T. A. Ehlers
Keyword(s):  
Time Scales ◽  
Landscape Evolution ◽  
Numerical Models ◽  
Higher Order ◽  
Glacial Erosion ◽  
Distinctive Features ◽  
Flow Models ◽  
Erosion Rates ◽  
Long Time ◽  
Order Of Magnitude

Abstract. Mountain topography is constructed through a variety of interacting processes. Over glaciological time scales, even simple representations of glacial-flow physics can reproduce many of the distinctive features formed through glacial erosion. However, detailed comparisons at orogen time and length scales hold potential for quantifying the influence of glacial physics in landscape evolution models. We present a comparison using two different numerical models for glacial flow over single and multiple glaciations, within a modified version of the ICE-Cascade landscape evolution model. This model calculates not only glaciological processes but also hillslope and fluvial erosion and sediment transport, isostasy, and temporally and spatially variable orographic precipitation. We compare the predicted erosion patterns using a modified SIA as well as a nested, 3-D Stokes-flow model calculated using COMSOL Multiphysics. Both glacial-flow models predict different patterns in time-averaged erosion rates. However, these results are sensitive to the climate and the ice temperature. For warmer climates with more sliding, the higher-order model has a larger impact on the erosion rate, with variations of almost an order of magnitude. As the erosion influences the basal topography and the ice deformation affects the ice thickness and extent, the higher-order glacial model can lead to variations in total ice-covered that are greater than 30%, again with larger differences for temperate ice. Over multiple glaciations and long-time scales, these results suggest that consideration of higher-order glacial physics may be necessary, particularly in temperate, mountainous settings.

scholarly journals Ice flow models and glacial erosion over multiple glacial–interglacial cycles

2015 ◽  
Vol 3 (1) ◽  
pp. 153-170 ◽  
Author(s):  
R. M. Headley ◽  
T. A. Ehlers
Keyword(s):  
Landscape Evolution ◽  
Numerical Models ◽  
Higher Order ◽  
Glacial Erosion ◽  
Distinctive Features ◽  
Flow Models ◽  
Erosion Rates ◽  
Covered Area ◽  
Order Of Magnitude ◽  
Long Timescales

Abstract. Mountain topography is constructed through a variety of interacting processes. Over glaciological timescales, even simple representations of glacial-flow physics can reproduce many of the distinctive features formed through glacial erosion. However, detailed comparisons at orogen time and length scales hold potential for quantifying the influence of glacial physics in landscape evolution models. We present a comparison using two different numerical models for glacial flow over single and multiple glaciations, within a modified version of the ICE-Cascade landscape evolution model. This model calculates not only glaciological processes but also hillslope and fluvial erosion and sediment transport, isostasy, and temporally and spatially variable orographic precipitation. We compare the predicted erosion patterns using a modified SIA as well as a nested, 3-D Stokes flow model calculated using COMSOL Multiphysics. Both glacial-flow models predict different patterns in time-averaged erosion rates. However, these results are sensitive to the climate and the ice temperature. For warmer climates with more sliding, the higher-order model yields erosion rates that vary spatially and by almost an order of magnitude from those of the SIA model. As the erosion influences the basal topography and the ice deformation affects the ice thickness and extent, the higher-order glacial model can lead to variations in total ice-covered area that are greater than 30% those of the SIA model, again with larger differences for temperate ice. Over multiple glaciations and long timescales, these results suggest that higher-order glacial physics should be considered, particularly in temperate, mountainous settings.

Two-Dimensional Local Deep-Level Transient Spectroscopy Imaging Using Super-Higher-Order Scanning Nonlinear Dielectric Microscopy

2016 ◽  
Author(s):  
N. Chinone ◽  
Y. Cho ◽  
R. Kosugi ◽  
Y. Tanaka ◽  
S. Harada ◽  
...  
Keyword(s):  
Deep Level Transient Spectroscopy ◽  
Deep Level ◽  
Higher Order ◽  
Trap Density ◽  
New Technique ◽  
Two Dimensional ◽  
Nonlinear Dielectric ◽  
Annealing Conditions ◽  
Transient Spectroscopy ◽  
A New Technique

Abstract A new technique for local deep level transient spectroscopy (DLTS) imaging using super-higher-order scanning nonlinear dielectric microscopy is proposed. Using this technique. SiCVSiC structure samples with different post oxidation annealing conditions were measured. We observed that the local DLTS signal decreases with post oxidation annealing (POA), which agrees with the well-known phenomena that POA reduces trap density. Furthermore, obtained local DLTS images had dark and bright areas, which is considered to show the trap distribution at/near SiCVSiC interface.

scholarly journals A Fast Estimator for Binary Choice Models with Spatial, Temporal, and Spatio-Temporal Interdependence

2021 ◽  
pp. 1-7
Author(s):  
Julian Wucherpfennig ◽  
Aya Kachi ◽  
Nils-Christian Bormann ◽  
Philipp Hunziker
Keyword(s):  
Computational Cost ◽  
Choice Models ◽  
Reduced Form ◽  
Binary Choice ◽  
Monte Carlo Experiments ◽  
Binary Choice Models ◽  
Spatio Temporal ◽  
Computationally Intensive ◽  
Pseudo Maximum Likelihood ◽  
Parameter Values

Abstract Binary outcome models are frequently used in the social sciences and economics. However, such models are difficult to estimate with interdependent data structures, including spatial, temporal, and spatio-temporal autocorrelation because jointly determined error terms in the reduced-form specification are generally analytically intractable. To deal with this problem, simulation-based approaches have been proposed. However, these approaches (i) are computationally intensive and impractical for sizable datasets commonly used in contemporary research, and (ii) rarely address temporal interdependence. As a way forward, we demonstrate how to reduce the computational burden significantly by (i) introducing analytically-tractable pseudo maximum likelihood estimators for latent binary choice models that exhibit interdependence across space and time and by (ii) proposing an implementation strategy that increases computational efficiency considerably. Monte Carlo experiments show that our estimators recover the parameter values as good as commonly used estimation alternatives and require only a fraction of the computational cost.

Sign in / Sign up

Export Citation Format

Share Document