Search NASA⌕ Search

SEARCH · Search NASA

Results for “diffusion processes”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 793 records · Page 44

Scrambling and thermalization in a diffusive quantum many-body system

Out-of-time ordered (OTO) correlation functions describe scrambling of information in correlated quantum matter. They are of particular interest in incoherent quantum systems lacking well defined quasi-particles. Thus far, it is largely elusive how OTO correlators spread in incoherent systems with diffusive transport governed by a few globally conserved quantities. Here, we study the dynamical response of such a system using high-performance matrix-product-operator techniques. Specifically, we consider the non-integrable, one-dimensional Bose–Hubbard model in the incoherent high-temperature regime. Our system exhibits diffusive dynamics in time-ordered correlators of globally conserved quantities, whereas OTO correlators display a ballistic, light-cone spreading of quantum information. The slowest process in the global thermalization of the system is thus diffusive, yet information spreading is not inhibited by such slow dynamics. We furthermore develop an experimentally feasible protocol to overcome some challenges faced by existing proposals and to probe time-ordered and OTO correlation functions. As a result, our study opens new avenues for both the theoretical and experimental exploration of thermalization and information scrambling dynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A Probabilistic Scheme for Semilinear Nonlocal Diffusion Equations with Volume Constraints

This work presents a probabilistic scheme for solving semilinear nonlocal diffusion equations with volume constraints and integrable kernels. The nonlocal model of interest is defined by a time-dependent semilinear partial integro-differential equation (PIDE), in which the integro-differential operator consists of both local convection-diffusion and nonlocal diffusion operators. Here, our numerical scheme is based on the direct approximation of the nonlinear Feynman–Kac formula that establishes a link between nonlinear PIDEs and stochastic differential equations. The exploitation of the Feynman–Kac representation avoids solving dense linear systems arising from nonlocal operators. Compared with existing stochastic approaches, our method can achieve first-order convergence after balancing the temporal and spatial discretization errors, which is a significant improvement of existing probabilistic/stochastic methods for nonlocal diffusion problems. Error analysis of our numerical scheme is established. The effectiveness of our approach is shown in two numerical examples. The first example considers a three-dimensional nonlocal diffusion equation to numerically verify the error analysis results. The second example presents a physics problem motivated by the study of heat transport in magnetically confined fusion plasmas.

97 MATHEMATICS AND COMPUTING↗

Multiphysics and Multiscale Modeling of Coupled Transport of Chloride Ions in Concrete

Chloride ions (Cl−)-induced corrosion is one of the main degradation mechanisms in reinforced concrete (RC) structures. In most situations, the degradation initiates with the transport of Cl− from the surface of the concrete towards the reinforcing steel. The accumulation of Cl− at the steel-concrete interface could initiate reinforcement corrosion once a threshold Cl− concentration is achieved. An accurate numerical model of the Cl− transport in concrete is required to predict the corrosion initiation in RC structures. However, existing numerical models lack a representation of the heterogenous concrete microstructure resulting from the varying environmental conditions and the indirect effect of time dependent temperature and relative humidity (RH) on the water adsorption and Cl− binding isotherms. In this study, a numerical model is developed to study the coupled transport of Cl− with heat, RH and oxygen (O2) into the concrete. The modeling of the concrete microstructure is performed using the Virtual Cement and Concrete Testing Laboratory (VCCTL) code developed by the U.S. National Institute of Standards and Technology (NIST). The concept of equivalent maturation time is utilized to eliminate the limitation of simulating concrete microstructure using VCCTL in specific environmental conditions such as adiabatic. Thus, a time-dependent concrete microstructure, which depends on the hydration reactions coupled with the temperature and RH of the environment, is achieved to study the Cl− transport. Additionally, Cl− binding isotherms, which are a function of the pH of the concrete pore solution, Cl− concentration, and weight fraction of mono-sulfate aluminate (AFm) and calcium-silicate-hydrate (C-S-H), obtained from an experimental study by the same authors are utilized to account for the Cl− binding of cement hydration products. The temperature dependent RH diffusion was considered to account for the transport of Cl− with moisture transport. The temperature and RH diffusion in the concrete domain, composite theory, and Cl− binding and water adsorption isotherms are used in combination, to estimate the ensuing Cl− diffusion field within the concrete. The coupled transport process of heat, RH, Cl−, and O2 is implemented in the Multiphysics Object-Oriented Simulation Environment (MOOSE) developed by the U.S. Idaho National Laboratory (INL). The model was verified and validated using data from multiple experimental studies with different concrete mixture proportions, curing durations, and environmental conditions. Additionally, a sensitivity analysis was performed to identify that the water-to-cement (w/c) ratio, the exposure duration, the boundary conditions: temperature, RH, surface Cl− concentration, Cl− diffusion coefficient in the capillary water, and the critical RH are the important parameters that govern the Cl− transport in RC structures. In a case study, the capabilities of the developed numerical model are demonstrated by studying the complex 2D diffusion of Cl− in a RC beam located in two different climatic regions: warm and humid weather in Galveston, Texas, and cold and dry weather in North Minnesota, Minnesota, subjected to time varying temperature, RH, and surface Cl− concentrations.

composite theory↗

Dynamic Metal–Support Interaction Dictates Cu Nanoparticle Sintering on Al 2 O 3 Surfaces

Nanoparticle sintering remains a critical challenge in heterogeneous catalysis. In this work, we present a unified deep potential (DP) model based on the Perdew–Burke–Ernzerhof approximation of density functional theory for Cu nanoparticles on three Al 2 O 3 surfaces (γ-Al 2 O 3 (100), γ-Al 2 O 3 (110), and α-Al 2 O 3 (0001)). Using DP-accelerated simulations, we reveal that the nanoparticle size-mobility relationship strongly depends on the supporting surface. The diffusion of nanoparticles on the two γ-Al 2 O 3 surfaces is almost independent of the size of the nanoparticle, while the diffusion on α-Al 2 O 3 (0001) decreases rapidly with increasing size. Interestingly, nanoparticles with fewer than 55 atoms diffuse several times faster on α-Al 2 O 3 (0001) than on γ-Al 2 O 3 (100) at 800 K while expected to be more sluggish based on their larger binding energy at 0 K. The diffusion on α-Al 2 O 3 (0001) is facilitated by dynamic metal–support interaction (MSI), where Al atoms move out of the surface plane to optimize contact with the nanoparticle and relax back to the plane as the nanoparticle moves away. In contrast, the MSI on γ-Al 2 O 3 (100) and on γ-Al 2 O 3 (110) is dominated by more stable and directional Cu–O bonds, consistent with the limited diffusion observed on these surfaces. Our extended MD simulations provide insight into the sintering processes, showing that the dispersity of the nanoparticles strongly influences the coalescence driven by nanoparticle diffusion. We observed that the coalescence of Cu 13 nanoparticles on α-Al 2 O 3 (0001) can occur in a short time (10 ns) at 800 K even with an initial internanoparticle distance increased to 3 nm, while the coalescence on the two γ-Al 2 O 3 surfaces are inhibited significantly by increasing the initial internanoparticle distance. These findings demonstrate that the dynamics of the supporting surface is crucial to understanding the sintering mechanism and offer guidance for designing sinter-resistant catalysts by engineering the support morphology.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Unified Entrainment and Detrainment Closures for Extended Eddy-Diffusivity Mass-Flux Schemes

The extended Eddy Diffusivity Mass Flux (EDMF) scheme is used as a unied parameterization that simulates subgrid-scale dynamics across a range of regimes.This work presents unied entrainment and detrainment closures that allow this scheme to simulate the dry convective boundary layer, shallow convection and deep convection , while adjacent works present other necessary closures for this unied scheme and test its performance for stable, stratocumulus topped and convective boundary layers. Entrainment and detrainment are modeled as a combination of turbulent and dynamical processes. The turbulent component is taken as downgradient diffusion between subdomains. The dynamical component is based on a ratio of buoyancy and vertical velocity that ensures consistent treatment of the ratio of mass-flux at vanishing area fraction. Dynamical detrainment is controlled by the difference in relative humidity between subdomains when liquid water is present. These closures are rst evaluated with entrainment and detrainment diagnosed from Large Eddy Simulations where decaying tracers are used to identify the different subdomains and mass exchanges between them. When implemented in a single column model (SCM), these closures are tested in simulations of several observationally motivated cases of dry convective boundary layer, shallow convection and deep convection. Results show good agreement between the SCM and Large Eddy Simulations in key properties for climate models and profiles of first, second and third moments. The development of diurnal cycles and the onset of and precipitation (with simplified microphysics) are well predicted by the SCM.

54 ENVIRONMENTAL SCIENCES↗

The Application of an Evolutionary Programming Process to a Simulation of the ETEX Large-Scale Airborne Dispersion Experiment

Abstract Airborne tracer simulations are typically performed using a dispersion model driven by a high-resolution meteorological model. Besides solving the dynamic equations of momentum, heat, and moisture on the resolved model grid, mesoscale models must account for subgrid-scale fluxes and other unresolved processes. These are estimated through parameterization schemes of eddy diffusion, convection, and surface interactions, and they make use of prescribed parameters set by the user. Such “free” model parameters are often poorly constrained, and a range of plausible values exists for each. Evolutionary programming (EP) is a process to improve the selection of the parameters. A population of simulations is first run with a different set of parameter values for each member, and the member judged most accurate is selected as the “parent” of a new “generation.” After a number of iterations, the simulations should approach a configuration that is best adapted to the atmospheric conditions. We apply the EP process to simulate the first release of the 1994 European Tracer Experiment (ETEX) project, which comprised two experiments in which a tracer was released in western France and sampled by an observing network. The EP process is used to improve a simulation of the RAMS mesoscale weather model, with weather data collected during ETEX being used to “score” the individual members according to how well each simulation matches the observations. The meteorological simulations from before and after application of the EP process are each used to force a dispersion model to create a simulation of the ETEX release, and substantial improvement is observed when these are validated against sampled tracer concentrations.

54 ENVIRONMENTAL SCIENCES↗

Asynchronous domain dynamics and equilibration in layered oxide battery cathode

To improve lithium-ion battery technology, it is essential to probe and comprehend the microscopic dynamic processes that occur in a real-world composite electrode under operating conditions. The primary and secondary particles are the structural building blocks of battery cathode electrodes. Their dynamic inconsistency has profound but not well-understood impacts. In this research, we combine operando coherent multi-crystal diffraction and optical microscopy to examine the chemical dynamics in local domains of layered oxide cathode. Our results not only pinpoint the asynchronicity of the lithium (de)intercalation at the sub-particle level, but also reveal sophisticated diffusion kinetics and reaction patterns, involving various localized processes, e.g., chemical onset, reaction front propagation, domains equilibration, particle deformation and motion. These observations shed new lights onto the activation and degradation mechanisms of state-of-the-art battery cathode materials.

25 ENERGY STORAGE↗

Characteristics of Bi-metallic Interfaces Formed During Direct Energy Deposition Additive Manufacturing Processing

Various additive manufacturing processes are being evaluated to reduce the time and cost for fabrication of low volume, complex, and multifunctional assemblies. This study evaluated two direct energy deposition processes for the fabrication of large bi-metallic structures. The materials evaluated were Inconel 625 and copper alloy C18150, which are used in various high heat flux applications. Inconel was deposited onto the C18150 substrate using blown powder and wire-fed processes. Complete bonding was obtained in both processes and the resulting interfaces were evaluated using microscopy and indentation testing. Differences were observed in the interface region suggesting the kinetic energy of the blown powder process resulted in more residual stress at the interface, promoting recrystallization and enhanced diffusion. This created a broader interface in the blown powder specimens compared to a narrower mechanically mixed interface with the wire-fed process.

36 MATERIALS SCIENCE↗

Deposition pressure dependent structural and optoelectronic properties of ex-situ boron-doped poly-Si/SiO x passivating contacts based on sputtered silicon

Among common methods to form polycrystalline silicon (poly-Si) films for passivating-contact solar cells, physical vapor deposition, in particular sputtering, is the safest one as it does not require any toxic gaseous precursors. One of the critical parameters to control the properties of sputtered silicon films is their deposition pressure. In this work, structural and optoelectronic characteristics of ex-situ boron-doped poly-Si/SiO x passivating contacts, formed from sputtered intrinsic amorphous silicon (a-Si) deposited at different pressures on top of SiO x /c-Si substrates and subjected to a high-temperature boron diffusion step, are investigated. The deposition rate and density of the as-deposited a-Si films increase with reducing pressure. Low-temperature photoluminescence spectra captured from the as-deposited samples at different pressures do not show typical emissions from hydrogenated a-Si. Meanwhile, their Fourier-transform infrared absorption spectra all show Si–H stretching modes, indicating that hydrogen had been initially incorporated into the chemical SiO x layers and eventually hydrogenated the a-Si/SiO x interfaces during the sputtering process. After the high-temperature boron-diffusion step, all hydrogen-related peaks disappear. Lower pressure films (1.5 and 2.5 mTorr) show more consistent improved performance after hydrogen treatments, compared to higher pressure films (4 and 5 mTorr). Overall, the resultant passivating contacts at 2.5 mTorr achieve a low single-side recombination current density Jo of ~9 fA/cm 2 , whereas their contact resistivity is still low at 15 mO cm 2 .

14 SOLAR ENERGY↗

The effect of remnant CdSe layers on the performance of CdSeTe/CdTe photovoltaic devices

Thin film CdTe-based photovoltaic devices have achieved high efficiency above 22 %. The recent improvement in efficiency is due to Se alloying in the CdTe absorbers to form a CdSeTe/CdTe structure. The subsequent band gap grading increases the short circuit current density. The Se can be introduced by depositing a precursor thin film of either CdSe or a CdSeTe alloy and then diffusing the Se into the CdTe during the high temperature cadmium chloride activation process. Using CdSe is preferred because it is easier to control the Se concentration. However, during fabrication of the CdSeTe/CdTe devices, the CdSe thickness needs to be precisely controlled to prevent the retention of a CdSe remnant layer after the activation treatment. Retention of a remnant CdSe layer causes a dramatic reduction in device efficiency. In this work, we show that the reduction in efficiency is caused by a number of factors. The remnant CdSe layer is n-type which moves the position of the p-n junction. Also, it is widely thought that the CdSe remnants are photo-inactive. In this work, we clarify that the individual CdSe grains are actually highly photo-active. However, the grain sizes in the CdSe remnant and the adjacent CdSeTe layer are very small resulting in a high grain boundary area. Although the grain boundaries are passivated with chlorine, cathodoluminescence imaging and electrical measurements show that this is only partially effective. Also, EQE measurements show that the remnant CdSe causes parasitic absorption. Overall, the remnant CdSe layer causes a reduction in short circuit current density and device efficiency. The thickness of the CdSe precursor layer and the cadmium chloride activation process conditions must be precisely optimised to ensure that all the CdSe is consumed and inter-diffused to form the CdSeTe alloy for highest efficiency devices.

CdSe remnant↗

RF surface resistance tuning of superconducting niobium via thermal diffusion of native oxide

Recently, Nb superconducting radio frequency cavities vacuum heat treated between 300 and 400°C for a few hours have exhibited very high quality factors (~5 × 10 10 at 2.0 K). Secondary ion mass spectrometry measurements of O, N, and C show that this enhancement in RF surface conductivity is primarily associated with interstitial O alloying via dissolution and diffusion of the native oxide. We use a theory of oxide decomposition and O diffusion to quantify previously unknown parameters crucial in modeling this process. RF measurements of a vacuum heat-treated Nb superconducting radio frequency cavity confirm the minimized surface resistance (higher Q 0 ) previously expected only from 800 °C diffusive alloying with N.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A quantum encryption design featuring confusion, diffusion, and mode of operation

Abstract Quantum cryptography—the application of quantum information processing and quantum computing techniques to cryptography has been extensively investigated. Two major directions of quantum cryptography are quantum key distribution (QKD) and quantum encryption, with the former focusing on secure key distribution and the latter focusing on encryption using quantum algorithms. In contrast to the success of the QKD, the development of quantum encryption algorithms is limited to designs of mostly one-time pads (OTP) that are unsuitable for most communication needs. In this work we propose a non-OTP quantum encryption design utilizing a quantum state creation process to encrypt messages. As essentially a non-OTP quantum block cipher the method stands out against existing methods with the following features: 1. complex key-ciphertext relation (i.e. confusion) and complex plaintext-ciphertext relation (i.e. diffusion); 2. mode of operation design for practical encryption on multiple blocks. These features provide key reusability and protection against eavesdropping and standard cryptanalytic attacks.

97 MATHEMATICS AND COMPUTING↗

A model for droplet heating and evaporation of water-in-oil emulsified fuel

Water-in-oil (W/O) emulsified fuel is a promising alternative fuel by inducing flash boiling of water at high temperature that can improve the atomization of fuel spray. The heating and evaporation process of emulsified fuel droplet is affected by the diffusion and coalescence of small dispersed water droplets in oil during the heating. In this study, a model is developed for the droplet heating and evaporation of W/O emulsified fuel with these key physics considered. The diffusivity of dispersed water droplets in oil is calculated using the Stokes-Einstein equation. The deactivation temperature of surfactant is selected as a criterion for the water coalescence since the water coalescence is due to the deactivation of surfactant. The process of water coalescence is simplified that the dispersed water droplets coalesce into one single water sub-droplet at the center of the oil droplet instantly as the droplet temperature reaches the deactivation temperature of surfactant. The model is validated against experimental data of single droplet under different heating temperatures, surfactant concentrations, and sizes of dispersed water droplets. Based on the proposed model, the effects of fuel properties and heating conditions on the droplet heating and evaporation are analyzed.

Droplet coalescence↗

On the mechanism of effective chemical reactions with turbulent mixing of reactants and finite rate of molecular reactions

A generalization of the theory of chemical transformation processes under turbulent mixing of reactants and arbitrary values of the rate of molecular reactions is presented that was previously developed for the variant of an instantaneous reaction [13]. The use of the features of instantaneous reactions when considering the general case, namely, the introduction of the concept of effective reaction for the reactant volumes and writing a closing conservation equation for these volumes, became possible due to the partition of the whole amount of reactants into “active” and “passive” classes; the reactants of the first class are not mixed and react by the mechanism of instantaneous reactions, while the reactants of the second class approach each other only through molecular diffusion, and therefore their contribution to the reaction process can be neglected. The physical mechanism of reaction for the limit regime of an ideal mixing reactor (IMR) is revealed and described. Although formally the reaction rate in this regime depends on the concentration of passive fractions of the reactants, according to the theory presented, the true (hidden) mechanism of the reaction is associated only with the reaction of the active fractions of the reactants with vanishingly small concentration in the volume of the reactor. It is shown that the rate constant of fast chemical reactions can be evaluated when the mixing intensity of reactants is much less than that needed to reach the mixing conditions in an IMR.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Peak radial growth of diffuse-porous species occurs during periods of lower water availability than for ring-porous and coniferous trees

Climate models project warmer summer temperatures will increase the frequency and heat severity of droughts in temperate forests of Eastern North America. Hotter droughts are increasingly documented to affect tree growth and forest dynamics, with critical impacts on tree mortality, carbon sequestration and timber provision. The growing acknowledgement of the dominant role of drought timing on tree vulnerability to water deficit raises the issue of our limited understanding of radial growth phenology for most temperate tree species. Here, we use well-replicated dendrometer band data sampled frequently during the growing season to assess the growth phenology of 610 trees from 15 temperate species over 6 years. Patterns of diameter growth follow a typical logistic shape, with growth rates reaching a maximum in June, and then decreasing until process termination. On average, we find that diffuse-porous species take 16–18 days less than other wood-structure types to put on 50% of their annual diameter growth. However, their peak growth rate occurs almost a full month later than ring-porous and conifer species (ca. 24 ± 4 days; mean ± 95% credible interval). Unlike other species, the growth phenology of diffuse-porous species in our dataset is highly correlated with their spring foliar phenology. We also find that the later window of growth in diffuse-porous species, coinciding with peak evapotranspiration and lower water availability, exposes them to a higher water deficit of 88 ± 19 mm (mean ± SE) during their peak growth than ring-porous and coniferous species (15 ± 35 mm and 30 ± 30 mm, respectively). Given the high climatic sensitivity of wood formation, our findings highlight the importance of wood porosity as one predictor of species climatic sensitivity to the projected intensification of the drought regime in the coming decades.

59 BASIC BIOLOGICAL SCIENCES↗

Water Sorption Study of Various Polymers and Polymeric Composites

In the formulations of many composite materials, polymeric materials are often used as binder materials to enhance processability, functionality, and safety during manufacturing. For instance, poly(ester urethane), such as Estane® 5703, is used in the formulation of PBX 9501, and vinyl copolymer elastomer (VCE) is used in the formulation of VCE/filler compos ites. The stability of these binders directly impacts the overall stability and performance of their composites. Estane and VCE are prone to hydrolytic degradation, with water sorption behavior playing a significant role in their long-term stability and aging behavior. Hence, it is critical to study their water sorption behavior under various humidity and tempera ture conditions. To improve the accuracy and reduce labor-intensive work involved in water sorption experiments, an automated Dynamic Vapor Sorption (DVS) system was employed, which can generate large datasets. To efficiently process and analyze these datasets, Python scripts were developed. These scripts not only streamlined data processing, but also accu rately derived the Arrhenius parameters related to water diffusion and sorption properties for the tested materials. This approach offers a robust framework for evaluating water trans port and sorption processes in polymeric materials and adds valuable capabilities to future water sorption testing efforts.

36 MATERIALS SCIENCE↗

Hierarchical morphologies in co-sputter deposited thin films

Co-depositions of immiscible alloy films at specific processing conditions have yielded hierarchical microstructures which consist of distinct features at multiple length scales, often agglomerates and concentration modulations on the sub-micrometer-scale and fine nanoprecipitates in a matrix on the nanoscale. The present work examined a series of immiscible alloy systems: Cu-Mo, Cu-Ag, Cu-Fe, Cu-Ta, Mo-Ag, Cu-Mo-Ag, to determine the kinetic conditions favorable for hierarchical organization and the formation mechanism of such structures. Thin films of six immiscible systems were sputter co-deposited over a range of deposition rates from 0.12 to 2 nm/s and various temperatures from 400 to 800°C. The resulting microstructures indicate that hierarchical structures form with sufficient disparity in kinetic energy between the constituent atoms, one species being highly mobile (A) and the other relatively immobile (B). This condition arises typically at elevated deposition temperatures and reduced deposition rates but is also alloy dependent. The hierarchical structures form during deposition via phase separation and self-organization processes across the multiple length scales. The adatoms diffuse on the film surface with the highly mobile species swiftly agglomerating into A-rich domains within which B-rich nanoprecipitates form, often self-organizing into periodic arrays. The smallest B-rich nanoprecipitates in the A-rich domains are found to be coherent and in a metastable crystal structure (B taking the structure of A), but coarser precipitates that exhibit the equilibrium structure of B element. The A-rich domains are surrounded by a B-rich matrix that phase-separates into a concentration modulated structure. In conclusion, the observations are interpreted via a model incorporating material properties and process parameters.

36 MATERIALS SCIENCE↗

A two-dimensional numerical study of the magneto-Rayleigh–Taylor instability with FLASH: Application to the staged Z-pinch concept

Magnetically driven implosions involving a liner collapsing onto a target are inherently vulnerable to the magneto-Rayleigh–Taylor instability (MRTI). Among the various approaches proposed to achieve fusion conditions within the target, the staged Z-pinch (SZP) concept employs a high-Z liner, the advantages of which remain an active area of investigation. Consequently, ongoing design optimization efforts are essential, while critical physical processes such as magnetic field (B-field) diffusion and radiation transport increase the complexity of required simulations. In this study, we utilize the new capabilities of the FLASH code to simulate in 2D, for the first time using this code, staged Z-pinch configurations (designated SZP1 and SZP1*), focusing on their stability. First, a comparison of simulation results with theoretical predictions of MRTI growth provide new insights into the often-overlooked influences of high mesh resolution and initial perturbation seeding on instability dynamics. These findings then lead to a discussion on potential improvements for future SZP experiments. These include the use of an axial B-field stabilization and the optimization of radiation transport processes during the implosion. The results presented herein establish the framework for simulating multi-dimensional Z pinches using the FLASH code and pave the way for the development of innovative experimental configurations leveraging its advanced simulation capabilities.

Adaptive mesh refinement↗