scholarly journals A two-dimensional glacier–fjord coupled model applied to estimate submarine melt rates and front position changes of Hansbreen, Svalbard

2018 ◽  
Vol 64 (247) ◽  
pp. 745-758 ◽  
Author(s):  
E. DE ANDRÉS ◽  
J. OTERO ◽  
F. NAVARRO ◽  
A. PROMIŃSKA ◽  
J. LAPAZARAN ◽  
...  
Keyword(s):  
Coupled Model ◽  
Small Scale ◽  
Two Dimensional ◽  
Time Step ◽  
Software Packages ◽  
Front Position ◽  
Circulation Patterns ◽  
Glacier Front ◽  
Order Of Magnitude ◽  
Frontal Ablation

ABSTRACTWe have developed a two-dimensional coupled glacier–fjord model, which runs automatically using Elmer/Ice and MITgcm software packages, to investigate the magnitude of submarine melting along a vertical glacier front and its potential influence on glacier calving and front position changes. We apply this model to simulate the Hansbreen glacier–Hansbukta proglacial–fjord system, Southwestern Svalbard, during the summer of 2010. The limited size of this system allows us to resolve some of the small-scale processes occurring at the ice–ocean interface in the fjord model, using a 0.5 s time step and a 1 m grid resolution near the glacier front. We use a rich set of field data spanning the period April–August 2010 to constrain, calibrate and validate the model. We adjust circulation patterns in the fjord by tuning subglacial discharge inputs that best match observed temperature while maintaining a compromise with observed salinity, suggesting a convectively driven circulation in Hansbukta. The results of our model simulations suggest that both submarine melting and crevasse hydrofracturing exert important controls on seasonal frontal ablation, with submarine melting alone not being sufficient for reproducing the observed patterns of seasonal retreat. Both submarine melt and calving rates accumulated along the entire simulation period are of the same order of magnitude, ~100 m. The model results also indicate that changes in submarine melting lag meltwater production by 4–5 weeks, which suggests that it may take up to a month for meltwater to traverse the englacial and subglacial drainage network.

scholarly journals Calving controlled by melt-under-cutting: detailed calving styles revealed through time-lapse observations

2019 ◽  
Vol 60 (78) ◽  
pp. 20-31 ◽  
Author(s):  
Penelope How ◽  
Kristin M. Schild ◽  
Douglas I. Benn ◽  
Riko Noormets ◽  
Nina Kirchner ◽  
...  
Keyword(s):  
Tidal Cycle ◽  
Time Lapse ◽  
Prognostic Models ◽  
Small Scale ◽  
Cycle Plays ◽  
Tidewater Glacier ◽  
Glacier Front ◽  
Detailed Record ◽  
Frontal Ablation

AbstractWe present a highly detailed study of calving dynamics at Tunabreen, a tidewater glacier in Svalbard. A time-lapse camera was trained on the terminus and programmed to capture images every 3 seconds over a 28-hour period in August 2015, producing a highly detailed record of 34 117 images from which 358 individual calving events were distinguished. Calving activity is characterised by frequent events (12.8 events h−1) that are small relative to the spectrum of calving events observed, demonstrating the prevalence of small-scale calving mechanisms. Five calving styles were observed, with a high proportion of calving events (82%) originating at, or above, the waterline. The tidal cycle plays a key role in the timing of calving events, with 68% occurring on the falling limb of the tide. Calving activity is concentrated where meltwater plumes surface at the glacier front, and a ~ 5 m undercut at the base of the glacier suggests that meltwater plumes encourage melt-under-cutting. We conclude that frontal ablation at Tunabreen may be paced by submarine melt rates, as suggested from similar observations at glaciers in Svalbard and Alaska. Using submarine melt rate to calculate frontal ablation would greatly simplify estimations of tidewater glacier losses in prognostic models.

scholarly journals Glacier–plume or glacier–fjord circulation models? A 2-D comparison for Hansbreen–Hansbukta system, Svalbard

2021 ◽  
pp. 1-14
Author(s):  
Eva De Andrés ◽  
Jaime Otero ◽  
Francisco J. Navarro ◽  
Waldemar Walczowski
Keyword(s):  
Water Depth ◽  
Seasonal Variations ◽  
Stress Fields ◽  
Two Dimensional ◽  
Coupled Models ◽  
Fjord Circulation ◽  
Glacier Front ◽  
Iceberg Calving ◽  
Best Fit ◽  
Frontal Ablation

Abstract Up to 30% of the current tidewater mass loss in Svalbard corresponds to frontal ablation through submarine melting and calving. We developed two-dimensional (2-D) glacier–line–plume and glacier–fjord circulation coupled models, both including subglacial discharge, submarine melting and iceberg calving, to simulate Hansbreen–Hansbukta system, SW Svalbard. We ran both models for 20 weeks, throughout April–August 2010, using different scenarios of subglacial discharge and crevasse water depth. Both models showed large seasonal variations of submarine melting in response to transient fjord temperatures and subglacial discharges. Subglacial discharge intensity and crevasse water depth influenced calving rates. Using the best-fit configuration for both parameters our two coupled models predicted observed front positions reasonably well (±10 m). Although the two models showed different melt-undercutting front shapes, which affected the net-stress fields near the glacier front, no significant effects on the simulated glacier front positions were found. Cumulative calving (91 and 94 m) and submarine melting (108 and 118 m) along the simulated period showed in both models (glacier–plume and glacier–fjord) a 1:1.2 ratio of linear frontal ablation between the two mechanisms. Overall, both models performed well on predicting observed front positions when best-fit subglacial discharges were imposed, the glacier–plume model being 50 times computationally faster.

scholarly journals The Influence of an Integration Time Step on Dynamic Calculation of a Vehicle-Track-Bridge under High-Speed Railway

Mathematics ◽  
2021 ◽  
Vol 9 (4) ◽  
pp. 431
Author(s):  
Junjie Ye ◽  
Hao Sun
Keyword(s):  
High Speed ◽  
Coupled Model ◽  
Short Wave ◽  
Integration Time ◽  
Vertical Acceleration ◽  
Dynamic Calculation ◽  
Time Step ◽  
High Speed Railway ◽  
Track Irregularity ◽  
Track Bridge

In order to study the influence of an integration time step on dynamic calculation of a vehicle-track-bridge under high-speed railway, a vehicle-track-bridge (VTB) coupled model is established. The influence of the integration time step on calculation accuracy and calculation stability under different speeds or different track regularity states is studied. The influence of the track irregularity on the integration time step is further analyzed by using the spectral characteristic of sensitive wavelength. According to the results, the disparity among the effect of the integration time step on the calculation accuracy of the VTB coupled model at different speeds is very small. Higher speed requires a smaller integration time step to keep the calculation results stable. The effect of the integration time step on the calculation stability of the maximum vertical acceleration of each component at different speeds is somewhat different, and the mechanism of the effect of the integration time step on the calculation stability of the vehicle-track-bridge coupled system is that corresponding displacement at the integration time step is different. The calculation deviation of the maximum vertical acceleration of the car body, wheel-sets and bridge under the track short wave irregularity state are greatly increased compared with that without track irregularity. The maximum vertical acceleration of wheel-sets, rails, track slabs and the bridge under the track short wave irregularity state all show a significant declining trend. The larger the vibration frequency is, the smaller the range of integration time step is for dynamic calculation.

scholarly journals Statistical Uncertainty of DNS in Geometries without Homogeneous Directions

2021 ◽  
Vol 11 (4) ◽  
pp. 1399
Author(s):  
Jure Oder ◽  
Cédric Flageul ◽  
Iztok Tiselj
Keyword(s):  
Channel Flow ◽  
A Priori ◽  
Time Integration ◽  
Navier Stokes ◽  
Statistical Uncertainty ◽  
Time Step ◽  
Order Of Magnitude ◽  
Averaging Time ◽  
The Individual ◽  
Time Required

In this paper, we present uncertainties of statistical quantities of direct numerical simulations (DNS) with small numerical errors. The uncertainties are analysed for channel flow and a flow separation case in a confined backward facing step (BFS) geometry. The infinite channel flow case has two homogeneous directions and this is usually exploited to speed-up the convergence of the results. As we show, such a procedure reduces statistical uncertainties of the results by up to an order of magnitude. This effect is strongest in the near wall regions. In the case of flow over a confined BFS, there are no such directions and thus very long integration times are required. The individual statistical quantities converge with the square root of time integration so, in order to improve the uncertainty by a factor of two, the simulation has to be prolonged by a factor of four. We provide an estimator that can be used to evaluate a priori the DNS relative statistical uncertainties from results obtained with a Reynolds Averaged Navier Stokes simulation. In the DNS, the estimator can be used to predict the averaging time and with it the simulation time required to achieve a certain relative statistical uncertainty of results. For accurate evaluation of averages and their uncertainties, it is not required to use every time step of the DNS. We observe that statistical uncertainty of the results is uninfluenced by reducing the number of samples to the point where the period between two consecutive samples measured in Courant–Friedrichss–Levy (CFL) condition units is below one. Nevertheless, crossing this limit, the estimates of uncertainties start to exhibit significant growth.

MAGNETIC FIELD AMPLIFICATION BY RELATIVISTIC SHOCKS IN AN INHOMOGENEOUS MEDIUM

2012 ◽  
Vol 08 ◽  
pp. 364-367
Author(s):  
YOSUKE MIZUNO ◽  
MARTIN POHL ◽  
JACEK NIEMIEC ◽  
BING ZHANG ◽  
KEN-ICHI NISHIKAWA ◽  
...  
Keyword(s):  
Magnetic Field ◽  
Inhomogeneous Medium ◽  
Magnetic Energy ◽  
Small Scale ◽  
Two Dimensional ◽  
Filamentary Structure ◽  
Field Amplification ◽  
Relativistic Shocks ◽  
Field Lines ◽  
Magnetohydrodynamic Simulations

We perform two-dimensional relativistic magnetohydrodynamic simulations of a mildly relativistic shock propagating through an inhomogeneous medium. We show that the postshock region becomes turbulent owing to preshock density inhomogeneity, and the magnetic field is strongly amplified due to the stretching and folding of field lines in the turbulent velocity field. The amplified magnetic field evolves into a filamentary structure in two-dimensional simulations. The magnetic energy spectrum is flatter than the Kolmogorov spectrum and indicates that the so-called small-scale dynamo is occurring in the postshock region. We also find that the amplitude of magnetic-field amplification depends on the direction of the mean preshock magnetic field.

A Numerical Study of Unsteady Natural Convection in a Rectangular Enclosure: The Effect of Compressibility

2004 ◽  
Author(s):  
K. M. Akyuzlu ◽  
Y. Pavri ◽  
A. Antoniou
Keyword(s):  
Mathematical Model ◽  
Stokes Equations ◽  
Numerical Study ◽  
Vertical Wall ◽  
Ideal Gas ◽  
Working Fluid ◽  
Two Dimensional ◽  
Algebraic Equations ◽  
Rectangular Enclosure ◽  
Circulation Patterns

A two-dimensional, mathematical model is adopted to investigate the development of buoyancy driven circulation patterns and temperature contours inside a rectangular enclosure filled with a compressible fluid (Pr=1.0). One of the vertical walls of the enclosure is kept at a higher temperature then the opposing vertical wall. The top and the bottom of the enclosure are assumed insulated. The physics based mathematical model for this problem consists of conservation of mass, momentum (two-dimensional Navier-Stokes equations) and energy equations for the enclosed fluid subjected to appropriate boundary conditions. The working fluid is assumed to be compressible through a simple ideal gas relation. The governing equations are discretized using second order accurate central differencing for spatial derivatives and first order forward finite differencing for time derivatives where the computation domain is represented by a uniform orthogonal mesh. The resulting nonlinear equations are then linearized using Newton’s linearization method. The set of algebraic equations that result from this process are then put into a matrix form and solved using a Coupled Modified Strongly Implicit Procedure (CMSIP) for the unknowns (primitive variables) of the problem. A numerical experiment is carried out for a benchmark case (driven cavity flow) to verify the accuracy of the proposed solution procedure. Numerical experiments are then carried out using the proposed compressible flow model to simulate the development of the buoyancy driven circulation patterns for Rayleigh numbers between 103 and 105. Finally, an attempt is made to determine the effect of compressibility of the working fluid by comparing the results of the proposed model to that of models that use incompressible flow assumptions together with Boussinesq approximation.

scholarly journals Small-scale spatio-temporal characteristics of accumulation rates in western Dronning Maud Land, Antarctica

2008 ◽  
Vol 54 (185) ◽  
pp. 315-323 ◽  
Author(s):  
Helgard Anschütz ◽  
Daniel Steinhage ◽  
Olaf Eisen ◽  
Hans Oerter ◽  
Martin Horwath ◽  
...  
Keyword(s):  
Accumulation Rate ◽  
Temporal Variations ◽  
Spatial Variations ◽  
Small Scale ◽  
Accumulation Rates ◽  
Order Of Magnitude ◽  
Dronning Maud Land ◽  
Spatio Temporal ◽  
Ground Penetrating ◽  
Firn Core

AbstractSpatio-temporal variations of the recently determined accumulation rate are investigated using ground-penetrating radar (GPR) measurements and firn-core studies. The study area is located on Ritscherflya in western Dronning Maud Land, Antarctica, at an elevation range 1400–1560 m. Accumulation rates are derived from internal reflection horizons (IRHs), tracked with GPR, which are connected to a dated firn core. GPR-derived internal layer depths show small relief along a 22 km profile on an ice flowline. Average accumulation rates are about 190 kg m−2 a−1 (1980–2005) with spatial variability (1σ) of 5% along the GPR profile. The interannual variability obtained from four dated firn cores is one order of magnitude higher, showing 1σ standard deviations around 30%. Mean temporal variations of GPRderived accumulation rates are of the same magnitude or even higher than spatial variations. Temporal differences between 1980–90 and 1990–2005, obtained from two dated IRHs along the GPR profile, indicate temporally non-stationary processes, linked to spatial variations. Comparison with similarly obtained accumulation data from another coastal area in central Dronning Maud Land confirms this observation. Our results contribute to understanding spatio-temporal variations of the accumulation processes, necessary for the validation of satellite data (e.g. altimetry studies and gravity missions such as Gravity Recovery and Climate Experiment (GRACE)).

Half division two-dimensional method for simple sinusoidal quantities parameters accurate measurement

2021 ◽  
pp. 16-21
Author(s):  
Kirill Yu. Solomentsev ◽  
Vyacheslav I. Lachin ◽  
Aleksandr E. Pasenchuk
Keyword(s):  
Measuring Instruments ◽  
Measured Signal ◽  
Two Dimensional ◽  
Measuring Device ◽  
New Approach ◽  
Precision Measuring ◽  
Useful Signal ◽  
Midpoint Method ◽  
Order Of Magnitude ◽  
Current Measuring

Several variants of half division two-dimensional method are proposed, which is the basis of a fundamentally new approach for constructing measuring instruments for sinusoidal or periodic electrical quantities. These measuring instruments are used in the diagnosis of electric power facilities. The most general variant, called midpoint method, is considered. The proposed midpoint method allows you to measure much smaller than using widespread methods, alternating currents or voltages, especially when changing the amplitude of the measured signal in very wide ranges, by 1–2 orders of magnitude. It is shown that using the midpoint method it is possible to suppress sinusoidal or periodic interference in the measuring path, in particular, to measure small alternating current when sinusoidal or periodic interference is 1–2 orders of magnitude higher than the useful signal. Based on the results of comparative tests, it was found that the current measuring device implementing the midpoint method is an order of magnitude more sensitive than the currently used high-precision measuring instruments.

Experiments on wind-driven heat exchange processes over melting snow

2021 ◽  
Author(s):  
Michael Haugeneder ◽  
Tobias Jonas ◽  
Dylan Reynolds ◽  
Michael Lehning ◽  
Rebecca Mott
Keyword(s):  
Snow Cover ◽  
Surface Wind ◽  
Infrared Camera ◽  
Snow Accumulation ◽  
Internal Boundary ◽  
Small Scale ◽  
Snow Melt ◽  
Two Dimensional ◽  
Atmospheric Conditions ◽  
Near Surface

<p>Snowmelt runoff predictions in alpine catchments are challenging because of the high spatial variability of t<span>he snow cover driven by </span>various snow accumulation and ablation processes. In spring, the coexistence of bare and snow-covered ground engages a number of processes such as the enhanced lateral advection of heat over partial snow cover, the development of internal boundary layers, and atmospheric decoupling effects due to increasing stability at the snow cover. The interdependency of atmospheric conditions, topographic settings and snow coverage remains a challenge to accurately account for these processes in snow melt models.<br>In this experimental study, we used an Infrared Camera (VarioCam) pointing at thin synthetic projection screens with negligible heat capacity. Using the surface temperature of the screen as a proxy for the air temperature, we obtained a two-dimensional instantaneous measurement. Screens were installed across the transition between snow-free and snow-covered areas. With IR-measurements taken at 10Hz, we capture<span> the dynamics of turbulent temperature fluctuations</span><span> </span>over the patchy snow cover at high spatial and temporal resolution. From this data we were able to obtain high-frequency, two-dimensional windfield estimations adjacent to the surface.</p><p>Preliminary results show the formation of a stable internal boundary layer (SIBL), which was temporally highly variable. Our data suggest that the SIBL height is very shallow and strongly sensitive to the mean near-surface wind speed. Only strong gusts were capable of penetrating through this SIBL leading to an enhanced energy input to the snow surface.</p><p>With these type of results from our experiments and further measurements this spring we aim to better understand small scale energy transfer processes over patch snow cover and it’s dependency on the atmospheric conditions, enabling to improve parameterizations of these processes in coarser-resolution snow melt models.</p>

Microdeformations in Sands by Digital Image Processing and Analysis

1996 ◽  
Vol 1548 (1) ◽  
pp. 31-37 ◽  
Author(s):  
Scott A. Raschke ◽  
Roman D. Hryciw ◽  
Gregory W. Donohoe
Keyword(s):  
Image Processing ◽  
Digital Image Processing ◽  
Digital Image ◽  
Small Scale ◽  
Two Dimensional ◽  
Localized Deformation ◽  
Displacement Fields ◽  
Particulate Media ◽  
Image Processing And Analysis ◽  
Video Images

Laboratory experiments are typically performed on particulate media to study stress-deformation behavior and to verify or calibrate computer models from controlled or measured boundary stresses and displacements. However, such data do not permit the formation of shear bands, displacement fields within flowing granular media, and other small-scale localized deformation phenomena to be identified. Described are two semiautomated computer vision techniques for accurately determining the two-dimensional displacement field in granular soils from video images obtained through a transparent planar viewing window. The techniques described are applicable for studying the behavior of particulate media under plane strain and certain axisymmetric test conditions. Digital image processing and analysis routines are used in two different computer programs, Tracker and Tracer, Tracker uses a graphical user interface that allows individual particles to be selected and tracked through a sequence of digital video images. A contrast edge detection algorithm delineates the two-dimensional projected boundaries of particles. The location of the centroid of each particle selected for tracking is determined from the boundary to quantify the trajectory of each particle. Tracer maps the trace or trajectory of specially dyed fluorescent particles in a sequence of video frames. A thresholding technique segments individual particle trajectories. Together, Tracker and Tracer provide a set of tools for identifying small-scale displacement fields in particulate assemblies deforming under either quasi-static or rapid loading (such as gravity flow).

Sign in / Sign up

Export Citation Format

Share Document