Melting Model Research Articles

Geochemistry and U–Pb zircon geochronology of S-type granites in the Karagwe Ankolean Belt, northwestern Tanzania: Implications for rare-metals mineralization

The Mesoproterozoic Karagwe Ankolean Belt in the Central Eastern Africa region hosts suites of granites that intruded supracrustal rocks and are associated with Sn–W mineralization. In northwestern Tanzania, the granites are S-type, peraluminous and calc-alkalic to alkalic-calcic in composition. Geochemical analyses of the granites reveal the following concentrations: SiO2 (71–77 wt%), Al2O3 (11–15 wt%), K2O (3–7 wt%), MgO (0.06–0.86 wt%), TiO2 (0.07–0.36 wt%) and Fe2O3T (0.71–2.40 wt%). The SiO2 contents correlate negatively with MgO, CaO and K2O and positively correlate with Nb, Ta, Zr, Sr, Ba, Hf, Th and the REE (Eu, Ce and Sm). Primitive mantle normalized spider plots show elevated LILE against depletions in Ba, Nb, Ti and La, suggesting inherited arc-magmatic signatures from the protolith. Chondrite-normalized REE spider plots show elevated LREE and relatively flat MREE and HREE as indicated by (La/Yb)cn = 2–9 and negative Eu* anomalies (0.05–0.24). Modelling of nearby Uha group sedimentary rocks partial melts, supports fractional crystallization and upper crustal contamination for generation of the KAB granites. Trace elemental ratios of Y/Nb, Nb/Ta, Zr/Hf and P205 suggest the granites are weakly mineralized to barren in Sn–W.The U–Pb zircon ages point to two successive and overlapping periods: from 1418 ± 15 Ma to 1395 ± 6 Ma (mean = 1406 Ma) and from 1385 ± 8 Ma to 1376 ± 11 Ma (mean = 1381 Ma). Inheritance U–Pb dates of xenocryst zircons, comprise early to late Paleoproterozoic ages ranging from 2376 ± 30 Ma, 2000 ± 27 Ma to 1901 ± 39 Ma, and from 1795 ± 7 Ma to 1738 ± 7 Ma, implying the participation of the Paleoproterozoic sedimentary rocks in the petrogenesis of the granites. A zircon rim overgrowth age of 1015 ± 26 Ma is interpreted to date a metamorphic/anatexis event that coincided with the c.1000 Ma amalgamation of the supercontinent Rodinia and the emplacement of Sn–W enriched granites. The ages of S-type granite magmatism (1418–1376 Ma) obtained from this study are restricted to the Mesoproterozoic within the widespread bimodal magmatism event c.1375 Ma.

Different patterns of changes in glacier mass balance and glacier runoff over the Tarim Basin, Central Asia

AbstractThe climate change impacted glacio‐hydrological regime and thus the availability of water resources in arid region of Central Asia. The effects of climate change in the magnitude or seasonality of regional glacier runoff were still poorly understood in Central Asia. In this study, the glacier runoff, components of glacier runoff, equilibrium line altitude (ELA) and Glacier Mass Balance (GMB) during 1959–2017 are evaluated by elevation‐dependent Degree‐Day glacier melt model in eight sub‐basins of Tarim Basin over Central Asia. The impacts of climate change on glacier and hydrology are assessed. The results suggested that climatic regime shifted to more warm‐wet pattern on glacier zone after 1990 in study area. The ablation and accumulation of glaciers showed different patterns in eight sub‐basins. All sub‐basins showed a glacier mass deficit and GMB displayed a marked decreasing trend, but also exhibiting discrepancy. The mean ELA and rising rate of ELA were higher in the southern region. The glacier runoff increased significantly after 1990 in Tarim Basin, with obviously temporal and spatial variations in sub‐basins. The mean annual volume of glacier runoff was 175.8 × 108 m3. The ice melt was a larger component of glacial runoff in Tarim Basin. The influence of rainfall runoff on glacier runoff was more obviously than snow melt runoff as more precipitation fell as rain in northern region. The larger proportions of snow melt runoff imply more precipitation fell as snowfall in southern region. The elevation‐dependent contributions in glacier runoff showed differences in individual basins. Temperature and precipitation played different role for the glacier runoff increases among the sub‐basins. Differences in sensitivity of GMB and glacier runoff were distinct and vary considerably. A thorough assessment of the spatially and temporally varying melt water originated by glaciers is crucial for the success of water scarcity adaptation under climate change.

Simulation Study on the Energy Utilization Efficiency of a Turbine Impeller Based on a Selective Laser Melting Process

In this paper, a simulation model for Selective Laser Melting (SLM) technology is established to simulate the additive manufacturing process of a turbine impeller for an aerospace engine. By utilizing the simulation model, variations in laser power and scanning speed are employed to obtain simulated results of thermal deformation for the turbine impeller under different laser power and scanning speed conditions. The results indicate that the thermal deformation of the component increases with the augmentation of laser power, decreases with the escalation of scanning speed, and eventually stabilizes. Based on the relationship between thermal deformation and energy, the energy utilization efficiency of the SLM process under different conditions is calculated. The findings demonstrate that, within a certain range of power, the synergistic effect of laser power and scanning speed allows for an increase in energy utilization efficiency and a reduction in processing time while ensuring the mechanical performance of the formed parts. Consequently, this approach proves effective in lowering production costs for complex components based on SLM technology.

The origin of Ni and Mn variations in Hawaiian and MORB olivines and associated basalts

The origins of differences in the Mn and Ni contents of olivine phenocrysts from ocean island basalts (OIB) and mid-ocean ridge basalts (MORB) are disputed. One prevalent hypothesis is mixing between peridotite melts and melts of an olivine-free pyroxenite that forms by reaction between large extents of melting of recycled ocean crust and peridotite. A second is peridotite melting at various pressures without contributions from recycled materials, and a third is mixing of recycled crust into peridotite sources. All models purport to account for olivine data, but an additional requirement is that the models account for corresponding liquid compositions. Here we demonstrate that none of the proposed models alone are able to account for both olivine and liquid data. The pyroxenite model predicts variations in liquid Mn to account for high Fe/Mn ratios, in contrast to the data where the high ratios result primarily from Fe variations. Using an updated model of Ni partitioning, combined with quantitative modeling of polybaric melting, confirms the important role of pressure on olivine compositions, and in addition shows that extent of melting and mantle fertility are also important. Most MORB are consistent with ∼10% -15% melting of peridotite at lower pressures, but low Ni olivines found on ridges adjacent to hot spots require a depleted mantle source. High Ni olivines can be produced by low degree melting of peridotite at high pressure, but such models do not produce the high Ni liquids that are also observed, nor the trace element abundances such as high TiO2. Mixing basalt and peridotite sources is not successful because it fails to account for the low Na/Ti ratios of OIB. A successful model for OIB involves adding ∼3–5% of a low degree eclogite melt with high trace element abundances and low Na/Ti to mantle peridotite, and then melting this source to small extents at high pressure. The need for a low degree melt of eclogite also rules out models that call upon high degrees of melting of eclogite that react with peridotite to produce OIB sources, and shows that this melting does not occur during upwelling beneath hot spots, where high degree melts of eclogite would occur. A firm conclusion of this work is that olivine compositions cannot be used to calculate the amount of recycled crust in mantle sources. Further work to explore variations within OIB and MORB will benefit from the realization that olivine compositions alone cannot be interpreted reliably without quantitative polybaric melting models and attention to co-existing liquids and trace element abundances.

The origin of Na-alkaline lavas revisited: new constraints from experimental melting of amphibole-rich metasomes+lherzolite at uppermost mantle pressure

We present a new experimental dataset for reaction experiments between natural amphibole-clinopyroxene metasomes (hornblendite) and synthetic lherzolite that produced Na-rich alkaline melts. Experiments were conducted at 1, 3 and 4 GPa and 1000–1300 °C. The generated melts range from foidite over basanite to phonotephrite. At 1 GPa between 1000 and 1100 °C amphibole decompression-breakdown products generate a phonotephritic melt. Among the breakdown components rhönite was found to be stable up to 1100 °C and 1 GPa. At 3 and 4 GPa the melt compositions are affected by phlogopite melting and shift to more foiditic compositions. We find that the melting of hornblendites and the reaction of the melt with the lherzolite produce wehrlitic residues with different olivine/clinopyroxene ratios. Wehrlite formation does not always require separate metasomatic processes but can be a direct by-product of alkaline volcanism. We applied a metasome melting model to the magmas of the Kula volcanic province, Turkey, and show that at 1 GPa basanite melts and phonotephrite melts cover the whole range of known Kula lava compositions. The Kula lava compositional trend can be therefore generated by basanite-phonotephrite melt mixing. A comparison of high-pressure (3–4 GPa) melts with natural nephelinite data shows overlap with many major, minor, and trace elements but differences in SiO2, FeO, and TiO2 argue that the natural nephelinite data do not represent primary metasome melts.

Review of the modeling approaches of phase change processes

In recent years, phase change materials have played an important role in the field of energy storage because of their flexibility and high efficiency in energy storage and release. However, most phase change processes are unsteady and highly nonlinear. The ways to obtain exact solutions are urgently needed. This study first summarizes the principles and characteristics of different methods to solve the phase transition process, including analytic methods, numerical methods, and CFD numerical methods. After that, the advantages of employing CFD software to simulate the phase change process are discussed, and the solidification and melting models in FLUENT are introduced. In particular, the influences and weights of natural convection and nanoparticles on the solidification/melting processes are revealed and enumerated in detail. Finally, the challenges and future developments in the solution methods, theoretical models, and numerical simulation applications of phase change materials are prospected. This review shows that most of the phase change heat transfer problems are solved by numerical methods. FLUENT is high-precision software commonly used in CFD software, its solidification and melting model is convenient to simulate the phase change process. Both natural convection and nanoparticles can shorten the melting time of phase change materials.

Melting of glacier ice enhanced by bursting air bubbles

Feedbacks between ice melt, glacier flow and ocean circulation can rapidly accelerate ice loss at tidewater glaciers and alter projections of sea-level rise. At the core of these projections is a model for ice melt that neglects the fact that glacier ice contains pressurized bubbles of air due to its formation from compressed snow. Current model estimates can underpredict glacier melt at termini outside the region influenced by the subglacial discharge plume by a factor of 10–100 compared with observations. Here we use laboratory-scale experiments and theoretical arguments to show that the bursting of pressurized bubbles from glacier ice could be a source of this discrepancy. These bubbles eject air into the seawater, delivering additional buoyancy and impulses of turbulent kinetic energy to the boundary layer, accelerating ice melt. We show that real glacier ice melts 2.25 times faster than clear bubble-free ice when driven by natural convection in a laboratory setting. We extend these results to the geophysical scale to show how bubble dynamics contribute to ice melt from tidewater glaciers. Consequently, these results could increase the accuracy of modelled predictions of ice loss to better constrain sea-level rise projections globally.

Characteristics of Glass in Products of Experimental Modeling of Impact Melts

Characteristics of Glass in Products of Experimental Modeling of Impact Melts

Numerical investigation for concentrated photovoltaic solar system in existence of paraffin equipped with MWCNT nanoparticles

To generate greater electricity and heat saving, a concentrated PVT (Photovoltaic thermal) system was applied in current study. The thermal section of the unit includes tubes containing nanofluid and paraffin containers in existence of MWCNT nanoparticles. Numerical simulation has been incorporated to investigate the impact of installing fins within the paraffin zone. The four mirrors with equal width have been utilized to concentrate the solar irradiation which leads to use of the whole solar spectrum. Laminar 3D flow within the cooling unit has been considered. The verification tests have been done in both views of melting modeling and PVT simulation. The grid independence procedure has been done and also, the best value of time step has been selected. Three configurations of fins (I-shape, short base arrow and long base arrow) have been implemented with four levels of Re (25, 750, 1250 and 1750). The maximum useful heat has been achieved for configuration of short base arrows. Dispersing MWCNT nanoparticles makes useful heat to increase about 9.58% for greatest inlet flow. The minimum temperature of the paraffin zone has been reported for short base arrows and augment of Re leads to reduction of temperature about 1.05%. Among scrutinized configurations, the maximum value of ηele belongs to the configuration of short base arrows and this function enhances with augment of Re about 3.72%.

A Densification prediction Model of Selective Laser Melting Based on BP Neural Network

A Densification prediction Model of Selective Laser Melting Based on BP Neural Network

Numerical Modeling and Simulation of Melting Phenomena for Freeze Valve Analysis in Molten Salt Reactors

In recent years, molten salt reactors (MSRs) have gained new momentum thanks to their potential for innovation in the nuclear industry, and several studies on their compliance with all the expected safety features are currently underway. In terms of passive safety, a strategy currently envisaged in accidental scenarios is to drain by gravity the molten salt, which acts both as fuel and coolant, in an emergency draining tank, ensuring both a subcritical geometry and proper cooling. To activate the draining system, a freeze plug, made of the same salt used in the core, is expected to open when the temperature in the core reaches high values. Up to this point, the freeze valve is still a key concept in the molten salt fast reactor (MSFR), and special attention must be paid to its analysis, given the requirement for passive safety, especially focusing on melting and solidification phenomena related to the molten salt mixture. This work aims to contribute to the macroscale modeling of melting and solidification phenomena relevant to the analysis of the freeze valve behavior. In particular, the focus is on the identification of the numerical models that can be adopted to achieve the quantitative insights needed for the design of the freeze valve. Among the ones available in the literature, the most appropriate models were selected based on a compromise between accuracy and computational efficiency. A critical look at the models allows for a synthetic and consistent formulation of the numerical models and their implementation in the open-source software OpenFOAM. The code was subsequently verified using analytical and numerical solutions already well established in the literature. A good agreement between the results produced by the developed solver and the reference solutions was obtained. In the end, the code was applied to simple case studies related to the freeze valve system, focusing on recognizing whether the developed code can model physical phenomena that can occur in a freeze valve. The results of the simulations are encouraging and show that the code can be used to model single-region melting or solidification problems. As such, this work constitutes a starting point for further development of the code, intending to achieve better quantitative predictions for the design of a freeze valve.

Continuous crustal thinning of the North China Craton in the Early Cretaceous: Evidence from A-type granites

Widespread crust-derived granitoids in North China provide evidence of Mesozoic lithospheric thinning and destruction of the North China Craton. However, the link between lower crust and lithospheric mantle reactivation remains poorly understood. We present whole-rock geochemistry and zircon O-Hf isotopes for Early Cretaceous peralkaline A-type granites (i.e., Xiwanzi, Xiangshan, and Yansehu) from the Yanshan Mountains of the eastern block of the North China Craton. The Xiwanzi granite has a U-Pb age of 134.3 ± 0.4 Ma, with homogeneous whole-rock Nd isotopes [εNd(t) = −7.7 to −7.5] but highly variable zircon O-Hf isotopes [δ18O = 5.1‰ to 7.2‰ and εHf(t) = −7.7 to +9.6]. Such isotopic variations likely resulted from a hybrid juvenile source consisting of partial melting of newly underplated basaltic rocks derived from both the metasomatized ancient and juvenile sub-continental lithospheric mantle. In contrast, the Xiangshan (117.4 ± 1.0 Ma) and Yansehu (116.3 ± 1.4 Ma) granites have relatively uniform isotopic compositions [whole-rock εNd(t) = −12.2 to −10.8, zircon δ18O = 5.2‰ to 6.4‰, and εHf(t) = −17.1 to −4.7], close to those of the Proterozoic mafic granulite xenoliths from the ancient lower crust. The isotopic variations of the Early Cretaceous A-type granites suggest a continuous change in the melting model from the underplated juvenile lowermost crust to the overlying ancient lower crust at high temperatures in an extensional environment. Our results indicate that A-type granites can be used to probe the nature and reactivation of the lower crust of an Archean craton. Significantly, large-scale continuous thinning of lower continental crust provides key insights into the generation of the widespread Mesozoic crust-derived granitoids in North China.

Study of contact melting of plate bundles by molten material in severe reactor accidents

In a severe reactor accident, a crust will form on the surface of the molten material during the core melting process. The crust will have a contact melting with the internal components of the reactor. In this paper, the contact melting process of the molten material on the austenitic stainless steel plate bundles is studied. The contact melting model of parabolic molten material on the plate bundles is proposed, and the rule and main effect factors of the contact melting are analyzed. The results show that the melting velocity is proportional to the slope of the paraboloid, the heat flux and the distance between two plates D. The influence of melt gravity and the plate width on melting velocity is negligible. The thickness of the molten liquid film is proportional to the heat flux and plate width, and it is inversely proportional to the gravity. With the increase of D, the liquid film thickness decreases at first and then increases gradually. The liquid film thickness has a minimum against D. When the width of the plate is small, the width of the plate is the main factor affecting the thickness of the liquid film. The parameters are coupled with each other. In a severe reactor accident, the wider internal components of reactor, which can increase the thickness of the melting liquid film and reduce the net input heat flux from the molten material to the components, are the effective measures to delay the melting process.

Physicochemical Modeling of Melting and Crystallization Processes of Latite and Trachyte

Physicochemical Modeling of Melting and Crystallization Processes of Latite and Trachyte

Rapid Recycling of Subducted Sediments in the Subcontinental Lithospheric Mantle

ABSTRACT Subduction recycling of sediments plays a key role in the geochemical evolution of Earth. The presence of recycled terrigenous sediments in upwelling plumes has been cited to explain the EM2 signature in ocean island volcanics, characterized by particularly high 87Sr/86Sr (>0.706). However, the origin of such isotopic anomalies in continental regions and the role of subducted sediments in the subcontinental lithospheric mantle (SCLM) remain unclear. The Himalaya–Tibet orogen is one of the world’s best places for deciphering continental subduction and the fate of subducted crustal materials in the mantle. Here we present a systematic study of the geochronology, mineral chemistry (especially clinopyroxene), whole-rock chemistry and Sr–Nd–Pb–Hf–O isotopic compositions of Cenozoic potassic–ultrapotassic lavas from the western Kunlun area of northwestern Tibet. New secondary ion mass spectrometry (SIMS) zircon U–Pb dating, coupled with published age results, constrain the timing of volcanism from ~8.3 Ma to the present. These lavas show geochemical characteristics that closely resemble the EM2 mantle end-member represented by the Samoan hotspot. Both whole rocks and individual magmatic clinopyroxenes display arc-like trace-element patterns and remarkably enriched Sr–Nd–Pb–Hf isotope compositions (87Sr/86Sr ≥ 0.7080; εNd ≤ −4.8; 206Pb/204Pb ≥ 18.704; εHf ≤ −2.6). Together with high zircon δ18O values (6.3–10.4‰), the data point to a mantle source enriched by recycled sedimentary materials. Geochemical modeling and geophysical evidence further indicate that the sediments were directly derived from the subducted Indian continental lithosphere during India–Eurasia collision. Partial melting models assuming a hybridized mantle source that contains ~5% Indian continental crust suggest that the primary melts of the potassic–ultrapotassic lavas could be formed by melting of a phlogopite-bearing garnet lherzolite at low melting degrees (1–5%). The magma geochemistry is consistent with the model of mélange melting, implying that the subducted sediments may detach from the downgoing Indian slab and rise up diapirically into the overlying mantle lithosphere. Unlike traditional models of subducted sediments entering the deep mantle, the western Kunlun EM2-like lavas reveal that subducted sediments can be rapidly recycled into the SCLM during continental subduction (probably <50 Myr). We suggest that the SCLM could be an important reservoir for subducted sediments. The findings are important to our understanding of mantle circulation rates and chemical heterogeneities.

Close-contact melting of shear-thinning fluids

This study aims at establishing a model for close-contact melting (CCM) of shear-thinning fluids. We presented a theoretical framework for predicting the variation of liquid melt film thickness and motion of unmelted solid for both Carreau and power-law fluids. We identified the appropriate energy equation considering the convective effect and derived an analytical temperature profile across the liquid film. Using the lubrication approximation, force equilibrium relationships and the corresponding numerical approaches were built. By using laser interferometry and photographic recording methods, we found excellent agreement between numerical solutions and experimental results for Carreau liquids, revealing that the convective effect weakens heat transfer and melting rate. We identified the critical liquid film thickness that determines three situations of CCM in the theoretical model for Carreau fluids. Numerical prediction demonstrated that the CCM of Carreau fluids can be almost equivalent to that of power-law fluids if the initial film thickness is greater than the critical value. Finally, approximate analytical models were developed for both Carreau and power-law models. For the applicability of the approximate analytical solutions, we derived two- and three-dimensional dimensionless phase diagrams of validity range and identified a key dimensionless group $(\varLambda Re)^{4/3}{Re}\left [3\ln (Ste+1)\right ]^{1/3}{Pe}^{-1/3}$ , where $\varLambda$ is dimensionless characteristic time, Re is Reynolds number, Ste is Stefan number and Pe is Peclect number. The reliability of the approximate solutions was verified by comparing with the numerical results. These approximate solutions enable convenient and low-cost computational prediction of the dynamic CCM process of shear-thinning fluids.

Research on Optimization Model of Heat Exchange Fin Structure in Energy Storage System

Aiming at the thermodynamic problem of heat transfer fin in heat storage system, the solidification and melting model and thermodynamic conventional model are established by using simulation software, thermodynamic image processing and other methods. MATLABAB, COMSOL and other software are used to program the optimal rectangular fin spacing Angle and the optimal heat exchange fin structure. Firstly, based on the reason that the simulation software can simulate the heat transfer process of phase change material (PCM) in the actual process to the greatest extent, the plane modeling and thermodynamic equation method are used to analyze the distance Angle of the rectangular fin, and the melting solidification model is established, and the optimal solution is obtained. Finally, combined with the analytic hierarchy model and TOPSIS method, the triangular fin structure and rectangular fin structure were compared from three levels, namely, the heat sink area, temperature rise time and PCM maintenance time of different shapes. And create a heat transfer rate of return to more favorably analyze the data and visually present the comparative results.

The unique solvability of stationary and non-stationary incompressible melt models in the case of their linearization

The unique solvability of stationary and non-stationary incompressible melt models in the case of their linearization

Redefined interface error, 2D verification and validation for pure solid-gallium phase change modeling by enthalpy-porosity methodology

For simulating phase change process of pure solid gallium, previously published literature contains some key problems, including the inappropriate definition & computation for interface errors, and the unsatisfactory numerical verifications & validations. Herein, the 2D verifications and 2D validations for the enthalpy-porosity modeling of pure gallium melting starting from a vertical end wall have been reported comprehensively. The well-known and classic experimental data presented in literature are used to verify and validate the model. The groundbreaking contributions in this work are that: i) The definition and computational method of the interface-position errors between the numerical and experimental results are revisited and proven carefully; ii) The numerical schemes as well as major parameters of 2D modeling gallium melting are selected with overall interface error <9%. However, these schemes and parameters are not appropriate for the solidification modeling of gallium due to the irregular and non-reproducible interface shapes caused by highly anisotropic properties of solid crystal gallium as experimentally reported in literature; iii) The 2D validations of enthalpy-porosity modeling of gallium melting are investigated in detail, and the maximum interface error is within the limits of ±12%. Moreover, the simulated liquid volume fractions and dimensionless heat transfer coefficients are well correlated in terms of related criteria parameters, and they are compared with correlations available in literature. It is concluded that to solidly validate the modeling of pure metal melting, both the global parameters (e.g., volume-averaged liquid fraction) and the local parameters (e.g., interface position) need to be considered together.

An empirical model to calculate snow depth from daily snow water equivalent: SWE2HS 1.0

Abstract. Many methods exist to model snow densification in order to calculate the depth of a single snow layer or the depth of the total snow cover from its mass. Most of these densification models need to be tightly integrated with an accumulation and melt model and need many forcing variables at high temporal resolution. However, when trying to model snow depth (HS) on climatological timescales, which is often needed for winter tourism-related applications, these preconditions can cause barriers. Often, for these types of applications, empirical snow models are used. These can estimate snow accumulation and snowmelt based on daily precipitation and temperature data only. To convert the resultant snow water equivalent (SWE) time series into snow depth, we developed the empirical model SWE2HS. SWE2HS is a multilayer densification model which uses daily snow water equivalent as sole input. A constant new snow density is assumed and densification is calculated via exponential settling functions. The maximum snow density of a single layer changes over time due to overburden and SWE losses. SWE2HS has been calibrated on a data set derived from a network of manual snow stations in Switzerland. It has been validated against independent data derived from automatic weather stations (AWSs) in the European Alps (Austria, France, Germany, Switzerland) and against withheld data from the Swiss manual observer station data set which was not used for calibration. The model fits the calibration data with root mean squared error (RMSE) of 8.4 cm, coefficient of determination (R2) of 0.97, and bias of −0.3 cm; it is able to achieve RMSE of 20.5 cm, R2 of 0.92, and bias of 2.5 cm on the validation data set from automatic weather stations and RMSE of 7.9 cm, R2 of 0.97, and bias of −0.3 cm on the validation data set from manual stations. The temporal evolution of the bulk density can be reproduced reasonably well on all three data sets. Due to its simplicity, the model can be used as post-processing tool for output of any other snow model that provides daily snow water equivalent output. Owing to its empirical nature, SWE2HS should only be used in regions with a similar snow climatology as the European Alps or has to be recalibrated for other snow climatological conditions. The SWE2HS model is available as a Python package which can be easily installed and used.