Search NASA⌕ Search

SEARCH · Search NASA

Results for “matrix diffusion”

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 541 records · Page 30

Randomized Sketching Algorithms for Low-Memory Dynamic Optimization

This paper develops a novel limited-memory method to solve dynamic optimization problems. The memory requirements for such problems often present a major obstacle, particularly for problems with PDE constraints such as optimal flow control, full waveform inversion, and optical tomography. In these problems, PDE constraints uniquely determine the state of a physical system for a given control; the goal is to find the value of the control that minimizes an objective. While the control is often low dimensional, the state is typically more expensive to store. This paper suggests using randomized matrix approximation to compress the state as it is generated and shows how to use the compressed state to reliably solve the original dynamic optimization problem. Concretely, the compressed state is used to compute approximate gradients and to apply the Hessian to vectors. The approximation error in these quantities is controlled by the target rank of the sketch. This approximate first- and second-order information can readily be used in any optimization algorithm. As an example, we develop a sketched trust-region method that adaptively chooses the target rank using a posteriori error information and provably converges to a stationary point of the original problem. Numerical experiments with the sketched trust-region method show promising performance on challenging problems such as the optimal control of an advection-reaction-diffusion equation and the optimal control of fluid flow past a cylinder.

97 MATHEMATICS AND COMPUTING↗

The role of irradiation-enhanced interstitial diffusion in over-pressurizing fission gas bubbles in UO 2

Fission gas bubbles in UO 2 nuclear fuel have been observed to exhibit pressures in excess of the equilibrium bubble pressure; however, the cause of bubble over-pressurization has not yet been demonstrated. The mechanical interaction between a bubble and the surrounding matrix or grain boundary depends on the internal pressure of the bubble and local stress state, such that over-pressurized bubbles are thought to be responsible for fragmentation and pulverization, when exposed to a temperature ramp. Here, in this work, we investigate the role of U interstitials, produced through irradiation, in over-pressurizing bubbles by using a combined molecular dynamics (MD) and cluster dynamics approach. Firstly, the energies for the capture of interstitials and vacancies by bubbles have been determined from MD as a function of the ratio of gas atoms to vacancies that make up the bubble. Secondly, these reaction energies have been implemented in the cluster dynamics code Centipede to predict bubble over-pressurization as a function of temperature for typical fission rates. It was found that there is a transition from low pressure bubbles (at high temperatures) to high pressure bubbles (at lower temperatures). The cause of this behavior was shown to be the creation of irradiation-induced interstitials that are highly mobile relative to vacancies at low temperature; whereas, vacancies are sufficiently mobile at high temperatures to limit bubble pressures. This result supports the hypothesis that over-pressurized bubbles form during steady-state operation and that this behavior is highly sensitive to the local pellet temperature.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Transport phenomena in electrolyte solutions: Nonequilibrium thermodynamics and statistical mechanics

The theory of transport phenomena in multicomponent electrolyte solutions is presented here through the integration of continuum mechanics, electromagnetism, and nonequilibrium thermodynamics. The governing equations of irreversible thermodynamics, including balance laws, Maxwell's equations, internal entropy production, and linear laws relating the thermodynamic forces and fluxes, are derived. Green–Kubo relations for the transport coefficients connecting electrochemical potential gradients and diffusive fluxes are obtained in terms of the flux–flux time correlations. Furthermore, the relationship between the derived transport coefficients and those of the Stefan–Maxwell and infinitely dilute frameworks are presented, and the connection between the transport matrix and experimentally measurable quantities is described. To exemplify the application of the derived Green–Kubo relations in molecular simulations, the matrix of transport coefficients for lithium and chloride ions in dimethyl sulfoxide is computed using classical molecular dynamics and compared with experimental measurements.

42 ENGINEERING↗

Super-Resolution Study of Polymer Mobility Fluctuations near c *

Nanoscale dynamic heterogeneities in synthetic polymer solutions are detected using super-resolution optical microscopy. To this end, we map concentration fluctuations in polystyrene–toluene solutions with spatial resolution below the diffraction limit, focusing on critical fluctuations near the polymer overlap concentration, c*. Two-photon super-resolution microscopy was adapted to be applicable in an organic solvent, and a home-built STED-FCS system with stimulated emission depletion (STED) was used to perform fluorescence correlation spectroscopy (FCS). The polystyrene serving as the tracer probe (670 kg mol –1 , radius of gyration R G ≈ 35 nm, end-labeled with a bodipy derivative chromophore) was dissolved in toluene at room temperature (good solvent) and mixed with matrix polystyrene (3,840 kg mol –1 , R G ≈ 97 nm, M w /M n = 1.04) whose concentration was varied from dilute to more than 10c*. Whereas for dilute solutions the intensity–intensity correlation function follows a single diffusion process, it splits starting at c* to imply an additional relaxation process provided that the experimental focal area does not greatly exceed the polymer blob size. We identify the slower mode as self-diffusion and the increasingly rapid mode as correlated segment fluctuations that reflect the cooperative diffusion coefficient, D coop . These real-space measurements find quantitative agreement between correlation lengths inferred from dynamic measurements and those from determining the limit below which diffusion coefficients are independent of spot size. Furthermore, this study is considered to illustrate the potential of importing into polymer science the techniques of super-resolution imaging.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Vacuum ultra-violet damage and damage mitigation for plasma processing of highly porous organosilicate glass dielectrics

Porous organosilicate glass thin films, with k-value 2.0, were exposed to 147 nm vacuum ultra-violet (VUV) photons emitted in a Xenon capacitive coupled plasma discharge. Strong methyl bond depletion was observed, concomitant with a significant increase of the bulk dielectric constant. This indicates that, besides reactive radical diffusion, photons emitted during plasma processing do impede dielectric properties and therefore need to be tackled appropriately during patterning and integration. The detrimental effect of VUV irradiation can be partly suppressed by stuffing the low-k porous matrix with proper sacrificial polymers showing high VUV absorption together with good thermal and VUV stability. In addition, the choice of an appropriate hard-mask, showing high VUV absorption, can minimize VUV damage. Particular processing conditions allow to minimize the fluence of photons to the substrate and lead to negligible VUV damage. For patterned structures, in order to reduce VUV damage in the bulk and on feature sidewalls, the combination of both pore stuffing/material densification and absorbing hard-mask is recommended, and/or the use of low VUV-emitting plasma discharge.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A two-phase model that unifies and extends the classical models of membrane transport

Two models describe solvent transport through swollen, nonporous membranes. The pore-flow model, based on fluid mechanics, works for porous membranes, whereas the solution-diffusion model invokes molecular diffusion to treat nonporous membranes. Both approaches make valid arguments for swollen polymer membranes, but they disagree in their predictions of intramembrane pressure and concentration profiles. Using a fluid-solid model that treats the solvent and membrane matrix as separate phases, we show both classical models to be valid, to represent complementary approaches to the same phenomenon, and to make identical predictions. The fluid-solid model clarifies recent reverse osmosis measurements; provides a predictive and mechanistic basis for empirical high-pressure limiting flux phenomena, in quantitative agreement with classic measurements; and gives a framework to treat nonporous but mechanically heterogeneous membrane materials.

Science & Technology - Other Topics↗

Investigation of interfacial strength in nacre-mimicking tungsten heavy alloys for nuclear fusion applications

Abstract Tungsten heavy alloys have been proposed as plasma facing material components in nuclear fusion reactors and require experimental investigation in their confirmation. For this purpose, a 90W–7Ni–3Fe alloy has been selected and microstructurally manipulated to present a multiphase brick-and-mortar structure of W-phase ‘bricks’ surrounded by a ductile ‘mortar’. This work draws inspiration from nature to artificially imitate the extraordinary combination of strength and stiffness exhibited by mollusks and produce a nacre-mimicking metal matrix composite capable of withstanding the extremely hostile environment of the reactor interior and maintaining structural integrity. The underlying mechanisms behind this integrity have been probed through high-resolution structural and chemical characterization techniques and have revealed chemically diffuse phase boundaries exhibiting unexpected lattice coherency. These features have been attributed to an increase in the energy required for interfacial decohesion in these systems and the simultaneous expression of high strength and toughness in tungsten heavy alloys.

36 MATERIALS SCIENCE↗

Calculating Alpha Eigenvalues and Eigenfunctions with a Markov Transition Rate Matrix Monte Carlo Method

For a nuclear system in which the entire α-eigenvalue spectrum is known, eigenfunction expansion yields the time-dependent flux response to any arbitrary source. Applications in which this response is of interest include pulsed-neutron experiments, accelerator-driven subcritical systems, and fast burst reactors, where a steady-state assumption used in neutron transport is invalid for characterizing the time-dependent flux. To obtain the α-eigenvalue spectrum, the transition rate matrix method (TRMM) tallies transition rates describing neutron behavior in a discretized position-direction-energy phase space using Monte Carlo. Interpretation of the resulting Markov process transition rate matrix as the operator in the adjoint α-eigenvalue problem provides an avenue for determining a large finite set of α eigenvalues and eigenfunctions of a nuclear system. Results from the TRMM are verified using analytic solutions, time-dependent Monte Carlo simulations, and modal expansion from diffusion theory. For simplified infinite-medium and one-dimensional geometries, the TRMM accurately calculates eigenvalues, eigenfunctions, and eigenfunction expansion solutions. Applications and comparisons to measurements are made for the small fast burst reactor CALIBAN and the Fort St. Vrain high-temperature gas-cooled reactor. For large three-dimensional geometries, discretization of the large position-energy-direction phase space limits the accuracy of eigenfunction expansion solutions using the TRMM, but it can still generate a fair estimate of the fundamental eigenvalue and eigenfunction. These results show that the TRMM generates an accurate estimate of a large number of α eigenvalues. This is not possible with existing Monte Carlo–based methods.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

AGR-3/4 Irradiation Test Train Disassembly and Component Metrology First Look Report

The AGR-3/4 experiment was designed to study fission product transport within graphitic matrix material and nuclear-grade graphite. To this end, this experiment consisted of 12 capsules, each fueled with 4 compacts containing UCO TRISO particles as driver fuel and 20 UCO designed-to-fail (DTF) fuel particles in each compact. The DTF fuel was fabricated with a thin pyrocarbon layer which was intended to fail during irradiation and provide a source of fission products. These fission products could then migrate through the compact and into the surrounding concentric rings of graphitic matrix material and/or nuclear graphite. Through post-irradiation examination (PIE) of the rings (including physical sampling and gamma scanning) fission product concentration profiles within the rings can be determined. These data can be used to elucidate fission product transport parameters (e.g. diffusion coefficients within the test materials) which will be used to inform and refine models of fission product transport. After irradiation in the Advanced Test Reactor (ATR) had been completed in April 2014, the AGR-3/4 experiment was shipped to the Hot Fuel Examination Facility (HFEF) at the Materials and Fuels Complex (MFC) for inspection, disassembly, and metrology. The AGR-3/4 test train was received at MFC in two separate shipments between February and April 2015. Visual examinations of the test train exterior did not indicate dimensional distortion, and only two small discolored areas were observed at the bottom of Capsules 8 and 9. No corresponding discoloration was found on the inside of these capsules, however. Prior to disassembly, the two test train sections were subject to analysis via the Precision Gamma Scanner (PGS), which did not indicate that any gross fuel relocation had occurred. A series of specialized tools (including clamps, cutters, and drills) had been designed and fabricated in order to carry out test train disassembly and recovery of capsule components (graphite rings and fuel compacts). This equipment performed well for separating each capsule in the test train and extracting the capsule components. Only a few problems were encountered. In one case, the outermost ring (the sink ring) was cracked during removal of the capsule through tubes. Although the sink ring will be analyzed in order to obtain a mass balance of fission products in the experiment, these cracks do not pose a major concern because the sink ring will not be analyzed in detail to obtain the spatial distribution of fission products. In Capsules 4 and 5, the compacts could not be removed from the inner rings. Strategies for removing the compacts are being evaluated. Sampling the inner rings with the compacts in-place is also an option. Dimensional measurements were made on the compacts, inner rings, outer rings, and sink rings. The diameters of all compacts decreased by 0.5 to 2.0 %. Generally, the extent of diametric shrinkage increased linearly with increasing neutron fluence. Most compact lengths also decreased. Compact lengths decreased with increasing fluence, reaching maximum shrinkage of about 0.9 % at a fast fluence of 4.0x1025 n/m2 E > 0.18 MeV. Above this fluence, the extent of length shrinkage appeared to decrease with fluence, and two compacts from Capsule 7 were found to have slightly increased in length (< 0.1 %) after a fluence of 5.2x1025 n/m2.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Migration of inclusions in a matrix due to a spatially varying interface energy

An interfacial energy can be a function of a bulk field such as temperature or electric field. Here we find that in a system with a gradient in temperature or electric potential, the resulting variation in interfacial energy can induce a particle to migrate by either surface or bulk diffusion. For a circular particle under a constant unidirectional gradient in the bulk field, the field dependence of the interfacial energy induces a circular particle to move as a circle at a constant velocity along the direction from higher to lower interfacial energy. A linear stability analysis of this steady state migration suggests that perturbations will damp out as time evolves, and thus under these conditions a migrating circular particle is morphologically stable. Other spatial distributions of interface energy can lead to the distortion of an initially circular shaped particle during migration. A phase field model is developed that captures these distortions and verifies the theoretical results mentioned above.

36 MATERIALS SCIENCE↗

Electrostatic Barriers to Nanoparticle Accessibility of a Porous Matrix

Translocation from one cavity to another through a narrow constriction (i.e., a “hole”) represents the fundamental elementary process underlying hindered mass transport of nanoparticles and macromolecules within many natural and synthetic porous materials, including intracellular environments. This process is complex and may be influenced by long-range (e.g., electrostatic) particle–wall interactions, transient adsorption/desorption, surface diffusion, and hydrodynamic effects. In this study, we used a three-dimensional (3D) tracking method to explicitly visualize the process of nanoparticle diffusion within periodic porous nanostructures, where electrostatic interactions were mediated via ionic strength. The effects of electrostatic interactions on nanoparticle transport were surprisingly large. For example, an increase in the Debye length of only a few nanometers (in a material with a hole diameter of ~100 nm) increased the mean cavity escape time 3-fold. A combination of computational and experimental analyses indicated that this hindered cavity escape was due to an electrostatic energy barrier in the region of the hole, which was quantitatively explained using DLVO theory. These findings explicitly demonstrate that the cavity escape process was barrier-limited and dominated by electrostatic effects.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An Alternative Low-Cost Process for Deposition of MCrAlY Bond Coats for Advanced Syngas/Hydrogen Turbine Applications

The objective of this project was to develop and optimize MCrAlY bond coats for syngas/hydrogen turbine applications using a low-cost electrolytic codeposition process. Prealloyed CrAlY-based powders were codeposited into a metal matrix of Ni, Co or Ni-Co during the electroplating process, and a subsequent post-deposition heat treatment converted it to the MCrAlY coating. Our research efforts focused on: (1) investigation of the effects of electro-codeposition configuration and parameters on the CrAlY particle incorporation in the NiCo-CrAlY composite coatings; (2) development of the post-deposition heat treating procedure; (3) characterization of coating properties and evaluation of coating oxidation performance; (4) exploration of a sulfurfree electroplating solution; (5) cost analysis of the present electrolytic codeposition process. Different electro-codeposition configurations were investigated, and the rotating barrel system demonstrated the capability of depositing NiCo-CrAlY composite coatings uniformly on the entire specimen surface, with the CrAlY particle incorporation in the range 37-42 vol.%. Post-deposition heat treatment at 1000-1200 °C promoted interdiffusion between the CrAlY particles and the Ni-Co metal matrix, resulting in β/γ’/γ or β/γ’ phases in the heat-treated coatings. The results also indicate that the post-deposition heat treatment should be conducted at temperatures ≤1100 °C to minimize Cr evaporation and outward diffusion of Ti. The electro-codeposited NiCrAlY coatings in general showed lower hardness and surface roughness than thermal spray MCrAlY coatings. Coating oxidation performance was evaluated at 1000-1100 °C in dry and wet air environments. The initial electro-codeposited NiCoCrAlY coatings containing relatively high sulfur did not show good oxidation resistance. After modifications of the coating process, the cleaner NiCoCrAlY coating exhibited good oxidation performance at 1000 °C during the 2,000 1-h cyclic oxidation test. A sulfur-free fluoboratebased plating solution was explored in order to reduce the sulfur level in the electro-codeposited MCrAlY coatings. However, chemical reactions occurred between the CrAlY powder and the fluoborate plating solution, resulting in dark powdery coatings. The fluoborate-based plating bath can thus only be used to codeposit more inert particles (e.g., oxides or carbides) instead of the present Al-containing metallic particles.

36 MATERIALS SCIENCE↗

AGR-1, AGR-2 and AGR-3/4 Release-to-Birth Ratio Data Analysis

A series of Advanced Gas Reactor (AGR) irradiation tests is being conducted in the Advanced Test Reactor at Idaho National Laboratory in support of development and qualification of tristructural isotropic (TRISO) low enriched fuel used in the High Temperature Gas-cooled Reactor (HTGR). Each AGR test consists of multiple independently controlled and monitored capsules containing fuel compacts placed in a graphite cylinder shrouded by a steel shell. These capsules are instrumented with thermocouples embedded in the graphite enabling temperature control. AGR configuration and irradiation conditions are based on prismatic HTGR technology distinguished primarily by the use of helium coolant, a low-power-density ceramic core capable of withstanding very high temperatures, and TRISO coated particle fuel. Thus, these tests provide valuable irradiation performance data to support fuel process development, qualify fuel for normal operating conditions, and support development and validation of fuel performance and fission product transport models and codes. The release-rate-to-birth-rate ratio (R/B) for each of fission product isotopes (i.e., krypton and xenon) is calculated from release rates in the sweep gas flow measured by the germanium detectors used in the AGR Fission Product Monitoring (FPM) System installed downstream from each irradiated capsule. Birth rates are calculated based on the fission power in the experiment and fission product generation models. Thus, this R/B is a measure of the ability of fuel kernel, particle coating layers, and compact matrix to retain fission gas atoms preventing their release into the sweep gas flow, especially in the event of particle coating failures that occurred during AGR-2 and AGR-3/4 irradiations. The major factors that govern gaseous radioactive decay, diffusion, and release processes are found to be material diffusion coefficient, temperature, and isotopic decay constant. For each of all AGR capsules, ABAQUS-based three-dimensional finite-element thermal models are created to predict daily averages of fuel compact temperatures for the entire irradiation period, which are used in establishing the R/B correlation with temperature and decay constant. This correlation can be used by reactor designers to estimate fission gas release from postulated failed fuel particles in HTGR cores, which is the key safety factor for fuel performance assessment.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

AGR-2 and AGR-3/4 Release-to-Birth Ratio Data Analysis

A series of Advanced Gas Reactor (AGR) irradiation tests is being conducted in the Advanced Test Reactor at Idaho National Laboratory in support of development and qualification of tristructural isotropic (TRISO) low enriched fuel used in the High Temperature Gas-cooled Reactor (HTGR). Each AGR test consists of multiple independently controlled and monitored capsules containing fuel compacts placed in a graphite cylinder shrouded by a steel shell. These capsules are instrumented with thermocouples embedded in the graphite enabling temperature control. AGR configuration and irradiation conditions are based on prismatic HTGR technology distinguished primarily by the use of helium coolant, a low-power-density ceramic core capable of withstanding very high temperatures, and TRISO coated particle fuel. Thus, these tests provide valuable irradiation performance data to support fuel process development, qualify fuel for normal operating conditions, and support development and validation of fuel performance and fission product transport models and codes. The release-rate-to-birth-rate ratio (R/B) for each of fission product isotopes (i.e., krypton and xenon) is calculated from release rates in the sweep gas flow measured by the germanium detectors used in the AGR Fission Product Monitoring (FPM) System installed downstream from each irradiated capsule. Birth rates are calculated based on the fission power in the experiment and fission product generation models. Thus, this R/B is a measure of the ability of fuel kernel, particle coating layers, and compact matrix to retain fission gas atoms preventing their release into the sweep gas flow, especially in the event of particle coating failures that occurred during AGR-2 and AGR-3/4 irradiations. The major factors that govern gaseous radioactive decay, diffusion, and release processes are found to be material diffusion coefficient, temperature, and isotopic decay constant. For each of all AGR capsules, ABAQUS-based three-dimensional finite-element thermal models are created to predict daily averages of fuel compact temperatures for the entire irradiation period, which are used in establishing the R/B correlation with temperature and decay constant. This correlation can be used by reactor designers to estimate fission gas release from postulated failed fuel particles in HTGR cores, which is the key safety factor for fuel performance assessment.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Microstructural Characteristics, Mechanical Properties, and Corrosion Behavior of Field-Assisted Diffusion Welded Alloy 617

This study investigated the microstructural characteristics, mechanical response, and corrosion resistance of Alloy 617 diffusion welded by electric field-assisted sintering (EFAS). The microstructure exhibited good grain boundary (GB) migration across the interface and flat interface in the samples bonded at different conditions. Graded grain size was observed, with grain coarsening being more significant near sample surface due to EFAS-induced temperature gradient. Tensile testing revealed that the specimens fractured in the matrix away from the interface. The peak strength of 807 MPa was obtained in the samples welded at 1000°C. The materials obtained at 1100°C exhibited reduced tensile strength but improved ductility. Strain maps revealed by digital image correlation showed the flat interfaces were less ductile compared to the matrix. The hardness reduced near sample surfaces due to enlarged grains induced by temperature gradient. The interfaces with GB migration exhibited improved corrosion resistance in molten salt compared to the flat interfaces.

36 - MATERIALS SCIENCE↗

Characterization of thermal destruction behavior of hybrid composites based on polyoxymethylene, ethylene-octene copolymer impact modifier and ZnO nanofiller

Hybrid polymer nanocomposites, composed of polyoxymethylene (POM), ethylene octene copolymer (EOC) and plasma synthesized tetrapod shaped zinc oxide (ZnO), were prepared by using melt compounding. The content of EOC in the POM based composites was varied between 10 and 50 mass %, while the content of ZnO was constant (2 mass %). Thermal behaviour of POM based systems was studied by using thermogravimetric analysis coupled with Fourier transform infrared spectroscopy. The influence of the elastomer content and/or ZnO addition on the thermal stability of POM based systems was evaluated. The influence of the α-octene content in the elastomer on the thermal decomposition behaviour of POM and its nanocomposites with ZnO was also evaluated. Results of thermogravimetric analysis showed that, by rising either the elastomer or ZnO content, thermal stability of the investigated POM composites was increased. The modifying effect of EOC17 in respect of thermal resistance was somewhat larger than that of EOC38 because of the smaller amount of tertiary carbon atoms in the macromolecular structure of the former elastomer. Improved thermal resistance of ZnO containing POM based composites was because of impermeable structure the inorganic nanofiller allowing decrease gas exchange rate and facilitating non-combustible gases, such as CO{sub 2}, stay in the zone of burning. Addition of ZnO have a potential to influence structure of the polymer blend matrix itself by improving its barrier characteristics.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Effects of high-dose neutron irradiation at light-water reactor relevant temperature on the mechanical properties of SiC/SiC composites

For this study, the neutron dose-dependent evolutions and the underlying mechanisms of properties of SiC fiber-reinforced SiC matrix (SiC/SiC) composites at a temperature relevant to light-water reactors (∼600 K) were investigated and analyzed. Chemical vapor infiltrated (CVI) SiC/SiC composites reinforced with Hi-Nicalon Type S or Tyranno SA3 fiber were neutron-irradiated to doses up to 30 dpa. The irradiated composites retained their flexural strengths. The thermal diffusivity and dimensional changes were mostly retained from 2.0 to 30.2 dpa. Discrepancies in irradiation responses among CVI SiC/SiC composites from different sources were found. Additional microstructural analysis using Raman spectroscopy and numerical analysis on irradiation effect on residual stress were used to explain how the microstructural variables, especially of carbon interphases, affect the mechanical properties in the 30–40 dpa dose range.

Continuous fiber-reinforced ceramic matrix composi↗

A fractional model for anomalous diffusion with increased variability: Analysis, algorithms and applications to interface problems

Fractional equations have become the model of choice in several applications where heterogeneities at the microstructure result in anomalous diffusive behavior at the macroscale. Here, we introduce a new fractional operator characterized by a doubly-variable fractional order and possibly truncated interactions. Under certain conditions on the model parameters and on the regularity of the fractional order we show that the corresponding Poisson problem is well-posed. Additionally, we introduce a finite element discretization and describe an efficient implementation of the finite-element matrix assembly in the case of piecewise constant fractional order. Through several numerical tests, we illustrate the improved descriptive power of this new operator across media interfaces. Furthermore, we present one-dimensional and two-dimensional h-convergence results that show that the variable-order model has the same convergence behavior as the constant-order model.

97 MATHEMATICS AND COMPUTING↗