Heat Flow Model Research Articles

Study on hybrid blue-IR laser welding on AZ31B magnesium alloy

How to perform a laser welding is crucial to assure the quality of welds on magnesium alloys, and the challenges are the suppression of bubble generation and the reduction of pore defects. In this paper, hybrid blue-IR laser welding (HBI-LW) is proposed to improve the stability of a molten pool and keyhole; HBI-LW is produced by using dual coaxial lasers that have different wavelengths and spot diameters. A multi-physics heat flow model is developed to analyze the effect of composite heat sources on porosity quantitatively. The combined volumetric heat source model has been verified experimentally. The experiments show that when the laser power is 1800 W and the power of HBI is increased to 600 W, the porosity is formed at a lower level of about 0.5 %. The results are aligned well with that from numerical simulations on the proposed volumetric heat source model. In comparison with LW, HBI-LW is able to generate a larger melt pool and plasma plume and to smoothen the rear wall of a keyhole. This regulates the flow of the molten pool, suppresses the keyhole closure, and reduces the occurrence of keyhole collapse effectively. The dynamics of the keyhole and the forming mechanism of pores are investigated. This research provides theoretical and practical guides to suppress pore defects in laser welding on magnesium alloy.

Analysis of the effect of the heat stabilizer outer surface finning on soil cooling efficiency

Relevance. The construction of buildings and structures in permafrost conditions is associated with the problem of soil thawing. To prevent this phenomenon, heat stabilizer have proven themselves good in practice. Their successful use is difficult without preliminary modeling of heat flows in the soil–heat stabilizer system. The most effective design of a passive heat stabilizer is an installation with two-phase of refrigerant. An additional increase in the efficiency of using a two-phase heat stabilizer is possible due to the finning of the outer surface of the outer pipe of it. Aim. To calculate the efficiency of using a heat stabilizer with and without finning of the outer surface. Objects. Heat stabilizer, refrigerant, heat transfer, frozen soil. Methods. The mathematical model for describing the processes of heat and mass transfer in the heat stabilizer–frozen soil system is considered in an axisymmetric formulation. There are three subtasks: the movement of liquid refrigerant in the inner tube of the heat stabilizer, the two-phase convective flow of refrigerant in the gap between the inner and outer tubes and conductive heat exchange in the system of the heat stabilizer–frozen soil. The influence of the finning of the outer surface of the heat stabilizer is considered within the framework of the modified concept of the skin factor. Results and conclusions. The authors have obtained the pressure and temperature distributions along the length of the heat stabilizer. It was found that the presence of fins on the outer surface of the external pipe of the heat stabilizer increases the heat flow from the soil by 10%. It is determined that the soil is effectively cooled within 1 m from the heat stabilizer; this distance is recommended as optimal for the placement of an adjacent heat stabilizer.

Performance evaluation and finite element modeling of heat, mass, and fluid flow inside a hybrid solar dryer during drying of paddy grains

Performance evaluation and finite element modeling of heat, mass, and fluid flow inside a hybrid solar dryer during drying of paddy grains

A Thermal Model for a Hybrid Gas Bearing System in Oil-Free Turbochargers

Abstract With little drag friction, gas bearings operate at high rotor speed and high temperature and thus enable long operating life along with material damping for mechanical energy dissipation. However, most computational tools for modeling gas bearings are not accessible to bearing manufacturers or potential end-users, and the models are restricted to specific geometries and operating conditions or are missing important features, thus limiting their utility. This paper presents the thermal energy transport in a hybrid gas bearing (HGB) for oil-free turbochargers (TCs) with integrated heat and fluid flow models to produce a comprehensive thermohydrodynamic analysis predictive tool and to design and evaluate air-lubricated gas bearings by predictions and experiments. An efficient algorithm couples the solution of the Reynolds equations and the thermal energy transport equations in the film between the rotor and a pad of the hybrid gas bearing and updates the temperature-dependent viscosity and film thicknesses during the iterative process. The balance of the bulk-flow thermal energy transport is that the summation of drag power losses plus heat convection into the film equals to the thermal energy advected by the film. The predictions show that the film temperature varies along circumferential and axial directions by following the balance of the energy flows and shows high film temperatures corresponding to the small film thicknesses and thus large film pressures. The predicted result also shows that the peak temperature of the film occurs at the exit side of the film near the trailing edge of each pad. Thermocouples installed axially at the leading and trailing edges of each pad measure temperature of the film to validate the thermal module. The test repeats and averaged for each rotor speed from 2 krpm to 22 krpm at intervals of 2 krpm, for supply pressure varying along 25 psi to 100 psi, and predicted film temperature matches well with the measured one. The parameters of interest in this analysis include the gas supply condition, bearing configuration, and the gas properties at high altitudes and high rotation speeds. The parametric study results show that the density and kinematic viscosity are the most dominant properties of the gas lubrication for the film temperature.

Multi-Criteria Calibration of a Thermo-Mechanical Model of Steel Plate Welding in Vacuum

This paper proposes a procedurefor multi-criteria calibration of a thermo-mechanical model for numerical simulation of welding in the space vacuum. A finite-element model of a steel plate is created. Experimental and computational data are obtained. An inverse problem is formulated for the vector identification of five calibration parameters from the heat-flow model. They are evaluated for adequacy with controlled accuracy according to four criteria. An optimization problem is solved using a two-step interactive procedure. The parameter space studying method (PSI) has been applied to the study of multidimensional regions by means of quasi-uniform sounding. A Pareto-optimal set is defined. It is used to determine reduced ranked Pareto subsets by μ-selection. Salukvadze optimum is also determined.

Coupled heat and fluid transport in pulled extrusion of optical fiber preforms

In the fabrication of optical fibers, the viscosity of the glass varies dramatically with temperature so that heat transfer plays an important role in the deformation of the fiber geometry. To date, models of fiber drawing processes have either assumed isothermal conditions or that advection is the dominant heat-transfer mechanism, with conduction therefore neglected. While advective heat transfer and neglect of conduction are valid in the fast process of drawing a preform to a fiber, heat conduction is also important in the slow process of extruding a fiber preform. In this paper, we derive and solve, in the context of pulled extrusion, a coupled heat and fluid flow model that incorporates heat conduction. The dramatic variation of viscosity with temperature means that this problem is extremely challenging to solve via standard numerical techniques and we therefore develop a novel finite-difference numerical solution method that proves to be highly robust. We use this method to show that conduction significantly affects the size of internal holes at the exit of the device and hence has important consequences for the process of manufacturing optical fibers.

A dynamic migration route planning optimization strategy based on real-time energy state observation considering flexibility and energy efficiency of thermal power unit

New energy resources compromising intermittent and fluctuating natures have been integrated into the power grid on a large scale, and thermal power units are obliged to promote the depth of deep peak shaving and flexible response-ability to cope with this new scenario. The safety, flexibility, and high efficiency of coal-fired units are the ultimate goals of its transformation, and there is a mutual influence and mutual restriction relationship among these indicators. To exploit the flexibility of thermal power units, take into account efficiency and safety, and play a supporting role in the new power system, a control strategy that weighs various indicators needs to be developed. Based on heat flow model and dynamic state space model, and updated matrix in real time, the key process parameters affecting the flexibility of heat transfer can be fully characterized. By using the observed state and calculation, the load response index and the energy efficiency index are optimized under the boundary constraint. The results show that the multi-objective optimization methods can avoid over regulation, reduce the flue gas flow by 3.23%–3.59%, and reduce the fluctuation of state quantity under the premise of satisfying the load response rate, about 2.32%–2.52% away from the safety boundary of the design temperature parameter to achieve the flexibility of frequency modulation, high efficiency of energy transmission and operation safety.

Integrated framework of multi−physics mechanism models for gas−powder flow, molten pool evolution and part deformation in high alloy steel additive manufacturing process

Multi−physics numerical models for metal additive manufacturing (MAM) could replace the expensive trial−and−error experiments, reveal the mechanisms of printing process, and provide an effective means of optimizing process parameters and mitigating part defects. However, the existing simulations of MAM process seldom consider the coupling effects between multiple recursive processes. There is an urgent need to overcome the dilemma of limiting MAM process simulation to a localized single sub−process. For directed energy deposition (DED) additive manufacturing process, the integrated simulation framework proposed in this study consists of a gas−powder flow model, a heat and fluid flow model, and a sequentially coupled thermal–mechanical model, which can be used to predict the feedstock feeding process, the evolution of molten pool, and the residual deformations of metal part, respectively. The accurate results of the three recursive sub−processes can be obtained by delivering the validated parameters between models. Corresponding experiments are conducted to verify the accuracy of the integrated framework based on 12CrNi2 alloy powders. The deviation of the predicted focal distance of powder streams is 5.14 %, and the deviations of the predicted molten pool height and width are all less than 8.5 %. The integrated framework could rapidly and accurately predict the residual deformation tendency and the position where the maximum deformation occurs. For single−track and eight−layer line deposition, the deviations of the predicted maximum deformation are all within 0.1 mm. The integrated framework could efficiently provide accurate and comprehensive solutions for optimizing DED process.

JIRAM Observations of Volcanic Flux on Io: Distribution and Comparison to Tidal Heat Flow Models

AbstractJuno has allowed clear, high‐resolution imaging of Io's polar volcanoes using the Jovian Infrared Auroral Mapper (JIRAM) instrument. We have used data from JIRAM's M‐band (4.78 μm) imager from 11 Juno orbits to construct a global map of volcanic flux. This map provides short‐term insight into the spatial distribution of volcanoes and the ways in which high‐ and low‐latitude volcanoes differ. Using spherical harmonic analysis, we quantitatively compare our volcanic flux map to the surface heat flow distribution expected from models of Io's tidal heat deposition (summarized in de Kleer, Park, et al. (2019, https://doi.org/10.26206/d4wc‐6v82). Our observations confirm previously detected systems of bright volcanoes at high latitudes. Our study finds that both poles are comparably active and that the observed flux distribution is inconsistent with an asthenospheric heating model, although the south pole is viewed too infrequently to establish reliable trends.

Influence of Mg Content on Microstructure Coarsening, Molten Pool Size, and Hardness of Laser Remelted Al(-x)-Mg-Sc Alloys.

Investigations leading to high-quality surfaces and optimized properties in laser remelted Al-Mg-Sc alloys remain scarce. Laser surface remelting (LSR) has been used for the final treatment of these alloys processed through additive manufacturing. However, the direct microstructure responses of the treated cast surfaces have not yet been investigated. In the present research, Al-3, 5, and 10 wt %-Mg-0.1 wt %-Sc alloys plates were processed using LSR to study the effects of local melting and rapid solidification. The morphology of the α-Al phase, microstructure coarsening, and hardness were mapped from the bottom to the top of the molten pools, varying with local Mg content and laser heat input (2.5 J/mm, 5 J/mm, and 10 J/mm). This study aimed to create a comprehensive map of the microstructures, hardness, and molten pool sizes under various conditions. The findings may help to optimize these alloys through understanding laser processing parameters. Methods used included CALPHAD computations, optical microscopy, SEM, EDS, image analysis, hardness tests, and heat flow models. The results obtained showed α-Al cell growth with bands in all alloys with hardness changes correlating with cell spacing and heat input. Higher Mg content resulted in more refined cells and a higher fraction of bands. Increased Mg content decreased the thermal diffusivity and enthalpy of melting, enlarging the molten pool size. Hardness increased with decreasing heat input and higher Mg content in the tested alloys, especially in the Al-10 wt %-Mg-0.1 wt %-Sc alloy as the heat input was varied.

Numerical Analysis of Low-Enthalpy Deep Geothermal Energy Extraction Using a Novel Gravity Heat Pipe Design

Geothermal energy, derived from the Earth’s internal heat, can be harnessed due to the geothermal gradient between the Earth’s interior and its surface. This heat, sustained by radiogenic decay, varies across regions, and is highest near volcanic areas. In 2020, 108 countries utilised geothermal energy, with an installed capacity of 15,950 MWe for electricity and 107,727 MWt for direct use in 2019. Low-enthalpy sources require binary systems for power production. Open-loop systems face issues like scaling, difficult water treatment, and potential seismicity, while closed-loop systems, using abandoned petroleum or gas wells, reduce costs and environmental impacts greatly. The novel geothermal gravity heat pipe (GGHP) design eliminates parasitic power consumption by using hydrostatic pressure for fluid circulation. Implemented in an abandoned well in north-east (NE) Slovenia, the GGHP uses a numerical finite difference method to model heat flow. The system vaporises the working fluid in the borehole, condenses it at the surface, and uses gravitational flow for circulation, maintaining efficient heat extraction. The model predicts that continuous maximum capacity extraction depletes usable heat rapidly. Future work will explore sustainable heat extraction and potential discontinuous operation for improved efficiency.

A sigmoidal model for predicting soil thermal conductivity-water content function in room temperature

Apparent thermal conductivity of soil (λ) as a function of soil water content (θ), i.e., λ(θ) is needed to determine the heat flow in soil. The function of λ(θ) can be used in heat and water flow models for simplicity. The objective of this study was to develop a sigmoidal model based on logistic equation for entire range of soil water contents and a wide range of soil textures that can be used in simulation of heat and water flow in respected modes. Further, performance of the developed sigmoidal model along with two other models in literature was evaluated. In the proposed sigmoidal model, the constants of this model are estimated based on empirical multivariate equations by using soil sand content and bulk density. The sigmoidal model was validated with good accuracy for a wide range of soil textures, as the relationship between the measured and predicted λ showed slope and intercept values of nearly 1.0 and 0.0, respectively. Comparison of the results obtained by sigmoidal model with those obtained from Johansen and Lu et al. models indicated that, the sigmoidal model was superior to the other two models in prediction of λ for a wide range of soil textures and soil water contents. Furthermore, comparison with a recently proposed model by Xiong et al. indicated that our sigmoidal model is superior. Therefore, our developed sigmoidal model can be used in heat and water flow models to predict the soil temperature and heat flow.

Structural modelling of subsurface geologic structures in Anambra and adjoining Bida Basins using aeromagnetic data: Implications for mineral explorations.

Aeromagnetic data over parts of the Northern Anambra and Southern Bida basins have been evaluated with emphasis on the subsurface geologic structures controlling mineralization and delineating the geomorphology of the basin. Four aeromagnetic data sheets (246 (Kabba), 247 (Lokoja), 266 (Auchi), and 267 (Idah)) were acquired through purchase, assemblage, analysis, and interpretation using integrated processing techniques of Euler deconvolution, analytical signal (AS), source parameter imaging (SPI), and 3D magnetic inversion model. The data are collected laterally in a series of NW-SE trajectories spaced 2 km apart, with an average flight altitude of approximately 150 m, and tielines occurring at approximately 20 km intervals. Qualitative observation of TMI and RAM shows that the research region is greatly fractured, with structures orienting in the E-W direction. The spectral analysis result shows that the sedimentary infillings range from 0.60 to 4.03 km. 3D model and curie isotherm depth subsurface maps reveal deeper curie depths at Kabba, Adavi, and Itakpe zones (28.0–40.0 km) and shallow curie depths at Auchi, Osara, Idah, and Anyigba areas (18.0–27.0 km). The result from the heat flow model shows an inverse relationship with the Curie isotherm depth. The contact depth sources calculated from Euler depth assessments show that the majority of these contact sources are trending in the E-W directions. The overall structural geomorphology of the research region conforms to the general pattern of structural trends within the research region. Most of the structures are trending NNE-SSW, NNW-SSE, and NE-SW, while the minor structures trend in the E-W direction. 3D inversion results show a discrete subsurface geologic structure that may house minerals within the region.

Regional controls on fluid flow in geothermal systems of the Taupo Volcanic Zone, New Zealand

ABSTRACT The Taupo Volcanic Zone (TVZ) provides almost 20% of New Zealand’s electricity production through seven of its 23 high-temperature geothermal systems. Fluid circulation in geothermal systems is influenced by permeable pathways, heat sources, and topography. To explore influences on systems in the TVZ, we have created regional-scale heat and fluid flow models using TOUGH2. We modelled several potential influences individually and in combination: TVZ rift boundaries, the densely faulted Taupo Fault Belt, varying thickness of volcanic cover, topographic loading, and localised heat sources inferred from magnetotelluric data. The locations of seven of the nine geothermal systems in our study area can be explained by TVZ-scale convection influenced by localised deep heat sources and hydraulic gradients due to topography. Regional-scale geology appears to be an additional influence at Rotorua and Ngatamariki, where cold recharge within the Taupo Fault Belt pushes modelled upflow towards the rift margins. At Ohaaki, modelled upflow is northwest of the system unless a dipping contact between basement and more permeable cover rocks is included. Modelled upflow zones do not correspond to geothermal systems at Te Kopia and Haroharo Volcanic Complex, which we hypothesise is related to uplift along the Paeroa Fault and shallow impermeable lava domes respectively.

Performance and modeling of freezing crystallization in desalination: Heat transfer and impurities flow

This work applies the freezing crystallization in the desalination and establishes the theoretical models with experiments, to correlate the technical parameters with product properties. We first establish the heat flow models to simulate the 3-dimensional spatiotemporal temperature distribution and deduce the two heat flows through the ice layer, in which the latent heat from the phase transition plays a dominant role in heat transfer. Next, mass transfer and boundary layer models are proposed to evaluate the desalination efficiency and demonstrate the roles of kinetic and thermodynamic factors on the impurity migration. Later, three process optimization strategies, including cooling rate regulation, air bubble-assisted, and crystal seed-induced techniques, are imported to enhance the separation efficiency. At last, sweating is incorporated to reach the secondary purification effect, in which the desalination ratio under the optimal trajectory exceeds 83 %. The established sweating models accurately predict the product properties with time and can trace the impurities flow from the ice layer and adherent inclusion. The inclusion occupies a major discharged melt at the initial sweating stage, but more ice layer begins to melt and enter into the sweating liquid, confirming that a slight increase in product purity accompanies the huge yield loss.

REHEATFUNQ (REgional HEAT-Flow Uncertainty and aNomaly Quantification) 2.0.1: a model for regional aggregate heat flow distributions and anomaly quantification

Abstract. Surface heat flow is a geophysical variable that is affected by a complex combination of various heat generation and transport processes. The processes act on different lengths scales, from tens of meters to hundreds of kilometers. In general, it is not possible to resolve all processes due to a lack of data or modeling resources, and hence the heat flow data within a region is subject to residual fluctuations. We introduce the REgional HEAT-Flow Uncertainty and aNomaly Quantification (REHEATFUNQ) model, version 2.0.1. At its core, REHEATFUNQ uses a stochastic model for heat flow within a region, considering the aggregate heat flow to be generated by a gamma-distributed random variable. Based on this assumption, REHEATFUNQ uses Bayesian inference to (i) quantify the regional aggregate heat flow distribution (RAHFD) and (ii) estimate the strength of a given heat flow anomaly, for instance as generated by a tectonically active fault. The inference uses a prior distribution conjugate to the gamma distribution for the RAHFDs, and we compute parameters for a uninformed prior distribution from the global heat flow database by Lucazeau (2019). Through the Bayesian inference, our model is the first of its kind to consistently account for the variability in regional heat flow in the inference of spatial signals in heat flow data. Interpretation of these spatial signals and in particular their interpretation in terms of fault characteristics (particularly fault strength) form a long-standing debate within the geophysical community. We describe the components of REHEATFUNQ and perform a series of goodness-of-fit tests and synthetic resilience analyses of the model. While our analysis reveals to some degree a misfit of our idealized empirical model with real-world heat flow, it simultaneously confirms the robustness of REHEATFUNQ to these model simplifications. We conclude with an application of REHEATFUNQ to the San Andreas fault in California. Our analysis finds heat flow data in the Mojave section to be sufficient for an analysis and concludes that stochastic variability can allow for a surprisingly large fault-generated heat flow anomaly to be compatible with the data. This indicates that heat flow alone may not be a suitable quantity to address fault strength of the San Andreas fault.

Numerical simulation of heat and mass transfer in laser directed energy deposition

In order to accurately predict the heat transfer and flow process of laser directed energy deposition, the numerical simulation of whole stage after the powder leaves the nozzle was carried out. First, the powder flow mode with laser body heat source in the first stage was established, and the average temperature and velocity of the powder at the convergence plane were obtained as 1300 K and 13 m/s, respectively. Based on this, the coupling model of heat flow between powder and substrate in the second stage is improved, and the influence of multiple source terms such as powder temperature, powder velocity, gas velocity, surface tension and thermal buoyancy is considered comprehensively. The evolution of the temperature and flow field of the cladding layer are analyzed in detail under the spatiotemporal variation, and it is clear that there are two convection currents with different clockwise directions in the molten pool, and the Marangoni convection occupies the main position. Finally, FeCrNiMo alloy powder was deposited on 42CrMo substrate, and its morphology, hardness and tensile strength in different directions were analyzed. The average height, width and depth deviations of the experiment and simulation are 17.3 %, 15.5 % and 22.9 %, respectively. The average strength of the tensile sample V2 prepared by the cladding layer is 11.8 % higher than that of the substrate.

Modelling Crystalline α-Mg Phase Growth in an Amorphous Alloy Mg72Zn28

A model of α-Mg grain growth in an amorphous Mg72Zn28 alloy matrix was developed together with numerical software. Its application enables tracking the growth process of the α-Mg phase in an amorphous alloy. The model was based on the diffusion-driven growth of α-Mg in an amorphous alloy under appropriate boundary conditions at an isothermal annealing temperature and taking into account the presence of a grain with an initial radius of 1 nm. The numerical model was based on a mathematical model of heat flow, described by the Fourier–Kirchhoff equation, and diffusion, described by Fick’s second law. The initial boundary conditions necessary to simulate grain growth in the amorphous phase were established. The results of the numerical simulation indicate grain growth with increasing isothermal annealing temperature and increasing isothermal annealing time.

Spherical geometry convection in a fluid with an Arrhenius thermal viscosity dependence: The impact of core size and surface temperature on the scaling of stagnant-lid thickness and internal temperature

The rock and rock-ice mixtures of the core-enveloping spherical shells comprising terrestrial body interiors have thermally determined viscosities well described by an Arrhenius dependence. Accordingly, the implied viscosity contrasts determined from the activation energies (E) characterizing such bodies can reach values exceeding 1040, for a temperature range that spans the conditions found from the lower mantle to the surface. In this study, we first explore the impact of implementing a cut-off to limit viscosity magnitude in cold regions. Using a spherical annulus geometry, we investigate the influence of core radius, surface temperature, and convective vigour on stagnant lid formation resulting from the extreme thermally induced viscosity contrasts. We demonstrate that the cut-off viscosity must be increased with decreasing curvature factor, f (=rin/rout, where rin and rout are the inner and outer radii of the annulus, respectively), if the solutions are to be not only computationally manageable but physically valid. We find that for statistically-steady systems, the mean temperature decreases with core size, and that a viscosity contrast of at least 107 is required for stagnant lid formation as f decreases below 0.5. Inverting the results from over 80 calculations featuring stagnant lids (from a total of approximately 180 calculations), we apply an energy balance model for heat flow across the thermal boundary layers and find that the non-dimensionalized temperature in the Approximately Isothermal Layer (AIL) in the convecting region under a stagnant lid is well predicted by TAIL′=12−2Tout′+γ+γ2+4γ1+Tout′ where γ is a function of E and f, and Tout′ is the non-dimensionalized surface temperature. Moreover, the normalized (i.e., non-dimensional) thickness of the stagnant lid, L′, can be obtained from a measurement of the non-dimensional surface heat flux once TAIL′ is determined. Stagnant-lid thicknesses increase from 10 to 30% of the shell thickness as f is decreased, and thick lids can overlie vigorously convecting underlying layers in small core bodies, potentially delaying secular cooling and suggesting that small objects with small cores may have developed thick elastic outer shells early in the solar system's history while maintaining vigorously convecting interiors. However, we also find that for the small number of 3-D calculations that we examined, parametrizations based on 2-D calculations overestimate the temperature of the convecting layer and the thickness of the conductive lid when f is small (less than 0.4).

Spatiotemporal variations in frost cracking measures in two dimensions: A case study for rock walls in Jotunheimen, southern Norway

The ground thermal regime has a profound impact on geomorphological processes and has been suggested to be particularly important for weathering processes in periglacial environments. Several frost-related damage indices have hitherto been developed to link climate and frost weathering potential in bedrock, although only for individual points or grid cells. Here, we model ground temperature and frost weathering potential in steep rock walls in the Jotunheimen Mountains, southern Norway, along a two-dimensional profile line for the Younger Dryas Stadial-Preboreal transition (c. 11.5 ka), the Holocene Thermal Maximum (c. 7.5 ka), the Little Ice Age (1750), and the 2010s. We use an established heat flow model and frost-cracking index based on the ice segregation theory. A central innovation of our model treatment is the implementation of ensemble simulations using distributions of automatically mapped crack radii in a rock wall, whereas previous frost damage models considered only a single characteristic crack radius. Our results allowed for the identification of sites with enhanced frost weathering. Such sites are typically found between rock walls and retreating glaciers, as well as in areas where snow depth changes abruptly, resulting in large thermal gradients. Hence, frost weathering may be highly active during glacier retreat, enhancing the damage to rock walls during deglaciation by adding to the damage from stress release. The coldest climates of the Younger Dryas Stadial-Preboreal transition and the Little Ice Age were generally most favorable for frost cracking. Such timing compares well with the knowledge about the timing of rockfall accumulations in Norway.