Search NASA⌕ Search

SEARCH · Search NASA

Results for “Spatial Generative Model”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
NASA NTRS1983

Satellite-derived sea surface temperature: Introduction

Satellites now play an increasing role in systematic monitoring of the global oceans. Measurements of sea surface temperature (SST) are of primary importance in understanding heat storage and transport within the ocean and cross the ocean-atmosphere boundary. In some regions, local changes in SST of only 1 to 2 C have major effects on global climate and weather patterns. The satellite measurements provide a data base complementary to the (sometimes) accurate but sparsely-distributed point measurements available from ships and buoys. The demands placed on satellite sensors are stringent. Accuracies of better than 1 C are required and are often desired to a few tenths of a degree. Furthermore, measurement accuracies must be stable spatially and temporally in order for satellite data to be used with confidence in models of air-sea interaction and climate. There now exists a need to evaluate objectively the performance of the latest generation of sensors under a sufficient variety of environmental conditions to indicate present accuracies, deficiencies, and potential for improvement.

Njoku, E. G.↗
NASA NTRS1992

Inverse kinematics problem in robotics using neural networks

In this paper, Multilayer Feedforward Networks are applied to the robot inverse kinematic problem. The networks are trained with endeffector position and joint angles. After training, performance is measured by having the network generate joint angles for arbitrary endeffector trajectories. A 3-degree-of-freedom (DOF) spatial manipulator is used for the study. It is found that neural networks provide a simple and effective way to both model the manipulator inverse kinematics and circumvent the problems associated with algorithmic solution methods.

Choi, Benjamin B.↗
DOE OSTI2018

Estimating consistent biased source probabilities over arbitrary bins - 25621

This paper describes computational algorithms for estimating consistent biased source probabilities for the CADIS and FW-CADIS methods that are compatible with most of the general fixed-source capabilities of the MCNP continuous-energy transport code. The basic algorithm is a novel Monte Carlo sampling technique that enables the straightforward estimation of biased source probabilities over arbitrarily complex distribution trees. In recent work, a method has been developed for generating correct biased spatial source probabilities in the presence of rejection from an arbitrary number of geometry cells. Two numerical test problems that demonstrate the correctness of this new algorithm are presented, including a simplified model of a cooling water system with an N-17 neutron source. (authors)

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2024

Autocalibration of the E3SM Version 2 Atmosphere Model Using a PCA-Based Surrogate for Spatial Fields

Global Climate Model tuning (calibration) is a tedious and time-consuming process, with high-dimensional input and output fields. Experts typically tune by iteratively running climate simulations with hand-picked values of tuning parameters. Many, in both the statistical and climate literature, have proposed alternative calibration methods, but most are impractical or difficult to implement. We present a practical, robust, and rigorous calibration approach on the atmosphere-only model of the Department of Energy's Energy Exascale Earth System Model (E3SM) version 2. Our approach can be summarized into two main parts: (a) the training of a surrogate that predicts E3SM output in a fraction of the time compared to running E3SM, and (b) gradient-based parameter optimization. To train the surrogate, we generate a set of designed ensemble runs that span our input parameter space and use polynomial chaos expansions on a reduced output space to fit the E3SM output. We use this surrogate in an optimization scheme to identify values of the input parameters for which our model best matches gridded spatial fields of climate observations. To validate our choice of parameters, we run E3SMv2 with the optimal parameter values and compare prediction results to expertly-tuned simulations across 45 different output fields. This flexible, robust, and automated approach is straightforward to implement, and we demonstrate that the resulting model output matches present day climate observations as well or better than the corresponding output from expert tuned parameter values, while considering high-dimensional output and operating in a fraction of the time.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

HTTR 3 D Cross Section Generation with Serpent and MAMMOTH

The High Temperature Test Reactor (HTTR) is a graphite moderated thermal reactor with a nominal power output of 30 MWth. The purpose of this work is to prepare a set of neutron cross sections for use in the transient analyses of three Loss of Forced Cooling (LOFC) events planned for this reactor. The cross sections are generated with the Serpent 2 3 D Monte Carlo code to include the axial heterogeneity and the strong axial coupling of this core. The approach yields good agreement of results when compared to previously used 2 D lattice models, which require a fine spatial discretization and, surprisingly, comes at a higher computational cost. In addition, the temperature distribution must be considered to obtain reasonable cross sections and the axial leakage component dominates the spectral effects in the core. Consequently, a temperature model as a function of the fuel and moderator temperatures is adopted, and the core is depleted to 390 effective full-power days (EFPDs). A full tabulation of cross sections is prepared at this burnup point for a variety of fuel and moderator temperatures. Next, the multiphysics reactor application MAMMOTH is used to evaluate the quality of these cross sections using a diffusion solver. As expected, the homogenization error in the cross sections is significant and the Super Homogenization (SPH) correction from MAMMOTH is necessary to preserve key reaction rates. The SPH corrected MAMMOTH results are in excellent agreement with the Serpent results and reproduce both the reference power profile and temperature coefficients at each tabulation point. This fact confirms the accuracy of the SPH corrected cross sections and demonstrates the maturity of using Serpent 2 and MAMMOTH for 3 D cross section generation, even for reactors as complex as the HTTR. Finally, the authors recommend that the cross sections should be parametrized with the local burnup and that a full core burnup calculation with coupled thermal fluids would provide a better estimate of the initial condition for any subsequent transient analysis. Unfortunately, this was not within the scope of this work, and direct analysis with these tabulations should be considered as a first order approximation, at best.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

MACHINE LEARNING BASED CHEMICAL EXPLOSIVE MODE ANALYSIS (ML-CEMA)

The software consists of a Machine Learning based Chemical Explosive Mode Analysis (ML-CEMA) tool for advanced computational flame diagnostics. CEMA, originally based on eigen-analysis of the local thermochemical system, is capable of identifying reaction fronts and limit phenomena such as auto-ignition and extinction in practical combustion systems, such as internal combustion engines and gas turbine combustors. However, the original CEMA is computationally expensive for large reaction mechanisms that are typically needed to describe fuel chemistry of practical large-hydrocarbon fuels. This novel ML-CEMA tool employs a ML technique to accelerate the eigen-analysis of the basic CEMA approach by orders of magnitude, thus making it suitable for practical fuels. In ML-CEMA,zero-dimensional (0D) reactors and one-dimensional (1D) premixed flames are first used to generate a large number of data points for neural network based ML training. The trained ML model is then used to perform CEMA prediction. This ML-CEMA tool has been demonstrated in canonical 0D and 1D configurations as well as highly-transient three-dimensional spray flames exhibiting multi-mode turbulent combustion, showing promising results. ML-CEMA, as a standalone tool, can be used for computationally-efficient diagnostics of massive datasets generated from both experiments and simulations. For example, based on spatially resolved measurements of a small set of reactive scalars(such as temperature, hydroxyl radical and formaldehyde), ML-CEMA can effectively identify flame fronts and rare events. ML-CEMA also provides a robust online or offline flame feature detection tool. When used for on-the-fly simulations, ML-CEMA further enables zone-adaptive combustion modeling, in which the predicted eigenvalue is used as a robust mode indicator for judicious assignment of locally-valid combustion models. This ML-CEMA based zone-adaptive model can lead to substantial computational cost savings when used for large-scale simulations of multi-mode combustion systems. Third Party Code Web Page to Download Code Web Page Location of Third Party License

Xu, Chao↗
NASA NTRS2023

Shape from spectra

We introduce a new unified atmospheric–topographic correction approach that estimates surface geometry directly from the radiance measurement. Surface topography influences the at-sensor radiance measurement, making precise topography modeling critical in applications like vegetation or snow studies in mountainous terrain. Currently, elevation maps are used to derive topographic variables such as the slope and sky-view factor. This process is error-prone since static global digital elevation models do not generally achieve the accuracy required, and even minor mismatches in spatial resolution can introduce significant artifacts in downstream processing. Here we demonstrate that it is possible to estimate topographic parameters directly from spectral data, ensuring perfect physical consistency, temporal coincidence, and spatial alignment. We present experiments estimating topographic slope in two scenes in Southern California, with data from NASA’s Next Generation Airborne Visible/Near Infrared Imaging Spectrometer (AVIRIS-NG). We compared our radiance-based estimates against high-resolution lidar datasets. Our initial validation result showed a correlation of R 2 = 0.864 (n = 160) over the homogeneous surface of Beckman Auditorium’s cone-shaped roof on the Caltech campus in Pasadena, California. We then validate the model over a larger study site near Santa Clarita, California, finding R 2 = 0.923 (n = 40,000) in a 350 x 350m area. The accuracy of our model estimates, combined with its systematic advantages over the alternative, show the potential of the approach for use in both airborne campaigns and orbital missions.

Nimrod Carmon↗
NASA NTRS2008

Stochastic Transport Theory for Investigating the Three-Dimensional Canopy Structure from Space Measurements

Radiation reflected from vegetation canopies exhibits high spatial variation. Satellite-borne sensors measure the mean intensities emanating from heterogeneous vegetated pixels. The theory of radiative transfer in stochastic media provides the most logical linkage between satellite observations and the three-dimensional canopy structure through a closed system of simple equations which contains the mean intensity and higher statistical moments directly as its unknowns. Although this theory has been a highly active research field in recent years, its potential for satellite remote sensing of vegetated surfaces has not been fully realized because of the lack of models of a canopy pair-correlation function that the stochastic radiative transfer equations require. The pair correlation function is defined as the probability of finding simultaneously phytoelements at two points. This paper presents analytical and Monte Carlo generated pair correlation functions. Theoretical and numerical analyses show that the spatial correlation between phytoelements is primarily responsible for the effects of the three-dimensional canopy structure on canopy reflective and absorptive properties. The pair correlation function, therefore, is the most natural and physically meaningful measure of the canopy structure over a wide range of scales. The stochastic radiative transfer equations naturally admit this measure and thus provide a powerful means to investigate the three-dimensional canopy structure from space. Canopy reflectances predicted by the stochastic equations are assessed by comparisons with the PARABOLA measurements from coniferous and broadleaf forest stands in the BOREAS Southern Study Areas. The pair correlation functions are derived from data on tree structural parameters collected during field campaigns conducted at these sites. The simulated canopy reflectances compare well with the PARABOLA data.

Huang, Dong↗
DOE OSTI2023

Current data are consistent with flat spatial hypersurfaces in the Λ CDM cosmological model but favor more lensing than the model predicts

Here, we study the performance of three pairs of tilted, and a pair of untilited, ΛCDM cosmological models, with three of these four pairs allowing for non-flat spatial hypersurfaces, against cosmic microwave background (CMB) temperature and polarization power spectrum data (P18), measurements of the Planck 2018 lensing potential power spectrum (lensing), and a large compilation of non-CMB data (non-CMB). For the eight models, we measure cosmological parameters and study whether or not pairs of the data sets (as well as subsets of them) are mutually consistent in these models. Half of these models allow the lensing consistency parameter A L , which re-scales the gravitational potential power spectrum, to be an additional free parameter to be determined from data, while the other three have A L = 1 which is the theoretically expected value. The pair of untilted non-flat ΛCDM models are incompatible with P18 data. The tilted spatially-flat models assume the usual primordial spatial inhomogeneity power spectrum that is a power law in wave number. The tilted non-flat models assume either the primordial power spectrum used in the Planck group anal yses [Planck P(q)], that has recently been numerically shown to be a good approximation to what is quantum-mechanically generated from a particular choice of closed inflation model initial conditions, or a recently computed power spectrum [new P(q)] that quantum-mechanically follows from a different set of non-flat inflation model initial conditions. In the tilted non-flat models with A L = 1 we find differences between P18 data and non-CMB data cosmological parameter constraints, which are large enough to rule out the Planck P(q) model at 3σ but not the new P(q) model. No significant differences are found when cosmological parameter constraints obtained with two different data sets are compared within the standard tilted flat ΛCDM model. While both P18 data and non-CMB data separately favor a closed geometry, with spatial curvature density parameter Ω k < 0, when P18+non-CMB data are jointly analyzed the evidence in favor of non-flat hypersurfaces subsides. Differences between P18 data and non-CMB data cosmological constraints subside when A L is allowed to vary. From the most restrictive P18+lensing+non-CMB data combination we get almost model-independent constraints on the cosmological parameters and find that the A L > 1 option is preferred over the Ω k < 0 one, with the A L parameter, for all models, being larger than unity by ~ 2.5σ. According to the deviance information criterion, in the P18+lensing+non-CMB analysis, the varying A L option is on the verge of being strongly favored over the A L = 1 one, which could indicate a problem for the standard tilted flat ΛCDM model. These data are consistent with flat spatial hypersurfaces but more and better data could improve the constraints on Ω k and might alter this conclusion. Error bars on some cosmological parameters are significantly reduced when non-CMB data are used jointly with P18+lensing data. For example, in the tilted flat ΛCDM model for P18+lensing+non-CMB data the Hubble constant H 0 = 68.09 ± 0.38 km s -1 Mpc -1 , which is consistent with that from a median statistics analysis of a large compilation of H 0 measurements as well as with a number of local measurements of the cosmological expansion rate. This H 0 error bar is 31% smaller than that from P18+lensing data alone.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2017

Creating Aerosol Types Form CHemistry (CATCH): a New Algorithm to Extend the Link Between Remote Sensing and Models

Current remote sensing methods can identify aerosol types within an atmospheric column, presenting an opportunity to incrementally bridge the gap between remote sensing and models. Here a new algorithm was designed for Creating Aerosol Types from CHemistry (CATCH). CATCH-derived aerosol types—dusty mix, maritime, urban, smoke, and fresh smoke—are based on first-generation airborne High Spectral Resolution Lidar (HSRL-1) retrievals during the Ship-Aircraft Bio-Optical Research (SABOR) campaign, July/August 2014. CATCH is designed to derive aerosol types from model output of chemical composition. CATCH-derived aerosol types are determined by multivariate clustering of model-calculated variables that have been trained using retrievals of aerosol types from HSRL-1. CATCH-derived aerosol types (with the exception of smoke) compare well with HSRL-1 retrievals during SABOR with an average difference in aerosol optical depth (AOD) <0.03. Data analysis shows that episodic free tropospheric transport of smoke is underpredicted by the Goddard Earth Observing System- with Chemistry (GEOS-Chem) model. Spatial distributions of CATCH-derived aerosol types for the North American model domain during July/August 2014 show that aerosol type-specific AOD values occurred over representative locations: urban over areas with large population, maritime over oceans, smoke, and fresh smoke over typical biomass burning regions. This study demonstrates that model-generated information on aerosol chemical composition can be translated into aerosol types analogous to those retrieved from remote sensing methods. In the future, spaceborne HSRL-1 and CATCH can be used to gain insight into chemical composition of aerosol types, reducing uncertainties in estimates of aerosol radiative forcing.

K W Dawson↗
DOE OSTI2025

Spectral line identification from a photoionised silicon plasma in emission

Next-generation X-ray satellite telescopes such as XRISM, NewAthena and Lynx will enable observations of exotic astrophysical sources at unprecedented spectral and spatial resolution. Proper interpretation of these data demands that the accuracy of the models is at least within the uncertainty of the observations. One set of quantities that might not currently meet this requirement is transition energies of various astrophysically relevant ions. Current databases are populated with many untested theoretical calculations. Accurate laboratory benchmarks are required to better understand the coming data. We obtained laboratory spectra of X-ray lines from a silicon plasma at an average spectral resolving power of ∼7500 with a spherically bent crystal spectrometer on the Z facility at Sandia National Laboratories. Many of the lines in the data are measured here for the first time. We report measurements of 53 transitions originating from the K-shells of He-like to B-like silicon in the energy range between ∼1795 and 1880 eV (6.6–6.9 Å). The lines were identified by qualitative comparison against a full synthetic spectrum calculated with ATOMIC. The average fractional uncertainty (uncertainty/energy) for all reported lines is ∼5.4 × 10 −5 . We compare the measured quantities against transition energies calculated with RATS and FAC as well as those reported in the NIST ASD and XSTAR’s uaDB. Average absolute differences relative to experimentally measured values are 0.20, 0.32, 0.17 and 0.38 eV, respectively. All calculations/databases show good agreement with the experimental values; NIST ASD shows the closest match overall.

astrophysical plasmas↗
DOE OSTI2018

A point process model for generating biofilms with realistic microstructure and rheology

Biofilms are communities of bacteria that exhibit a multitude of multiscale biomechanical behaviours. Recent experimental advances have led to characterisations of these behaviours in terms of measurements of the viscoelastic moduli of biofilms grown in bioreactors and the fracture and fragmentation properties of biofilms. These properties are macroscale features of biofilms; however, a previous work by our group has shown that heterogeneous microscale features are critical in predicting biofilm rheology. In this paper, we use tools from statistical physics to develop a generative statistical model of the positions of bacteria in biofilms. Specifically, the model is a type of pairwise interaction model (PIM). We show through simulation that the macroscopic mechanical properties of biofilms depend on the choice of microscale spatial model. A key finding is that uniform and non-uniform sets of points lead to differing mechanical properties. This distinction appears not to have been previously considered in mathematical biofilm literature. We also found that realisations of a biologically informed PIM have realistic in silico mechanical properties, and have statistical properties that closely match experimental data. We also note that a Poisson spatial point process of suitable number density also yields realistic mechanical properties, but that the spatial distribution of points does not reflect those occurring in our experimentally observed biofilm.

Mathematics↗
DOE OSTI2019

Valuation of ecosystem services in alternative bioenergy landscape scenarios

Abstract Agricultural land in the Midwest is a source of food and fuel, as well as biodiversity. It is also a cause of excess nutrients that make their way to the Mississippi River and the Gulf of Mexico. To address unsustainable changes to biogeochemical cycles and ecosystem functions, a multidisciplinary approach involving social science, natural science, and engineering is often effective. Given the potential of second‐generation biofuels, and capitalizing on the deep‐rooted perennial bioenergy crops capable of thriving in poor soils, we demonstrated an integrated socio‐environmental analysis of the impacts of growing switchgrass within row‐crop landscapes in Illinois. In this study, we model land use scenarios that incorporate switchgrass as a biofuel crop in a Midwest corn‐belt watershed using the Soil Water Assessment Tool coupled with an economic analysis for the Vermilion Basin in Illinois. We estimated the values of ecosystem services under an alternative bioenergy landscape, including commodity and bioenergy crops, changes in biogeochemistry, and recreational services. The estimated annual values of nitrate and sediment reduction attributed to bioenergy crops range from $38 million to $97 million and $16,000 to $197,000, respectively. The annual value of carbon dioxide emission reduction ranges from $1.8 million to $6.1 million based on the initial crop rotation pattern. Estimated average annual values for wildlife viewing, water‐based recreation, and pheasant hunting are $1.24 million, $0.17 million, and $0.3 million, respectively. To our knowledge, this study represents the first effort to comprehensively quantify ecosystem services using a process‐based model, and estimate their value in an alternative bioenergy landscape. The information we generate could aid in understanding the potential for biomass production from marginal land and the total economic value of the landscape at various spatial scales. The framework is useful in fostering alternative bioenergy landscapes with synergies in a food, energy, and conservation nexus.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Power Grid Behavioral Patterns and Risks of Generalization in Applied Machine Learning

Recent years have seen a rich literature of data-driven approaches designed for power grid applications. However, insufficient consideration of domain knowledge can impose a high risk to the practicality of the methods. Specifically, ignoring the grid-specific spatiotemporal patterns (in load, generation, and topology, etc.) can lead to outputting infeasible, unrealizable, or completely meaningless predictions on new inputs. To address this concern, this paper investigates real-world operational data to provide insights into power grid behavioral patterns, including the time-varying topology, load, and generation, as well as the spatial differences (in peak hours, diverse styles) between individual loads and generations. Then based on these observations, we evaluate the generalization risks in some existing ML works caused by ignoring these grid-specific patterns in model design and training.

Li, Shimiao↗
DOE OSTI2021

A HPC Theory-Guided Machine Learning Cyberinfrastructure for Communicating Hydrometeorological Data Across Scales

High-resolution predictions of hydrometeorological variables are critical for supporting hydropower generation decisions and flood control at hydroelectric power plants. Traditional climate and hydrologic models rely on the numerical simulation of detailed physical processes. Therefore, running these simulations is time-, labor-, and computation-intensive. Improving the spatial and temporal resolution in these modeling outputs could lead to cubic increases in both the simulation time and computational demands, rendering high-resolution hydrometeorological predictions expensive and impractical. Many past studies apply the super resolution (SR) technique to downscale climate models using deep learners. However, deep learners are deemed “black-boxes,” as their derivation processes from low-resolution outputs to high-resolution outputs are often hidden. Their results are difficult for domain scientists to interpret and validate. Thus, there is a need for an exploratory machine learning approach that can partially integrate domain-specific theory and knowledge into the data-driven mapping process between simulation outputs of different spatial scales. The domain-specific theory and knowledge can be incorporated into the data model through an inductive approach in which process-related environmental variables are used and analyzed as key drivers (i.e., environmental surrogates) to reflect the complex physical processes. Many of these variables, such as land use land cover, soil types, topography, digital elevation, air temperature, and various watershed characteristics, can be directly measured through sensors or remote sensing techniques. Additionally, SR applications that can downscale hydrological and hydrodynamics models to efficiently produce high-resolution (1 m) flood depth grids are still rare. Since the flood depth grid can be used to support critical decisions for flood control operation at hydroelectric power plants, it is crucial to enable an SR-based capability for interpolating high-resolution flood inundation maps.

13 HYDRO ENERGY↗
NASA NTRS1987

Interseismic deformation for normal fault earthquakes

Comparison between the coseismic vertical deformation associated with normal fault earthquakes and permanent deformation indicated by geologic structure indicates that large deformation must occur during the period between earthquakes. When sufficient geodetic and geologic information is available it is possible to estimate the spatial character of this interseismic deformation. A case in point is the 1983 Borah Peak, Idaho earthquake. This normal fault earthquake produced roughly 5 times as much basin subsidence as it did uplift of the adjacent mountain ranges. In contrast, geophysical and geological observations show that the basin is roughly as deep as the bounding range is high. A similar normal fault event in the Basin and Range, the 1959 Hebgen Lake earthquake, was also accompanied by substantially larger subsidence than uplift. In this case, postseismic geodetic measurements show broad regional uplift with a spatial pattern which is roughly consistent with that proposed for the Borah Peak event. Modeling suggests that postseismic viscoelastic relaxation and strain accumulation in an elastic lithosphere overlying a viscoelastic asthenosphere are possible physical mechanisms to generate interseismic uplift. These mechanisms may be the contemporary expression of those processes responsible for the high elevation of the Basin and Range Province.

Reilinger, Robert↗
NASA NTRS1991

Homogeneity of spatial correlation statistics of tropical oceanic rainfall

The possibility of uniform horizontal correlation scales for tropical oceanic rainfall has been examined by a study of satellite-observed microwave data as a proxy measure of rain rates. From the brightness temperatures from the electrically scanning microwave radiometer on Nimbus 5 near nadir during the year 1974, the mean spatial autocorrelation function as a function of simultaneous pixel separation was calculated in each 5 by 5 deg grid box over the tropical Pacific and Atlantic for each season. The equal-time spatially lagged correlations were compared for geographical dependence to investigate the hypothesis of homogeneity. A simple model of the spatial statistics of the microwave brightness temperatures was used, consisting of a mixture of uncorrelated spatial white noise incoherently superimposed on a spatially correlated field (spatial red noise). The red noise signals are presumed to be generated by convective activity in the tropical atmosphere. The parameters of the red noise are consistent with this scheme over the tropical oceans, yielding a uniform spatial scale of about 50 km throughout.

Shin, Kyung-Sup↗
DOE OSTI2019

A space-time adaptive finite element method with exponential time integrator for the phase field model of pitting corrosion

In this paper we propose a space-time adaptive finite element method for the phase field model of pitting corrosion, which is a parabolic partial differential equation system consisting of a phase variable and a concentration variable. A major challenge in solving this phase field model is that the problem is very stiff, which makes the time step size extremely small for standard temporal discretizations. Another difficulty is that a high spatial resolution is required to capture the steep gradients within the diffused interface, which results in very large number of degrees of freedom for uniform meshes. To overcome the stiffness of this model, we combine the Rosenbrock–Euler exponential integrator with Crank–Nicolson scheme for the temporal discretization. Moreover, by exploiting the fact that the speed of the corroding interface decreases with time, we derive an adaptive time stepping formula. For the spatial approximation, we propose a simple and efficient strategy to generate adaptive meshes that reduces the computational cost significantly. Thus, the proposed method utilizes local adaptivity and mesh refinement for efficient simulation of the corrosive dissolution over long times in heterogeneous media with complex microstructures. We also present an extensive set of numerical experiments in both two and three dimensional spaces to demonstrate efficiency and robustness of the proposed method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗