Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical 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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Pitch angle diffusion of newborn ions due to intrinsic turbulence in the solar wind

The objective of the present study is to understand the interaction of the solar wind with newborn ions in far upstream regions of a comet where the level of intrinsic turbulence is moderately low. Based on the assumption that quasi-linear theory is adequate and applicable, the pitch angle diffusion process and the time evolution of the newborn ion distribution function are investigated. Numerical solutions to the quasi-linear diffusion equation, including the effect of resonance broadening and that of continuous creation of newborn ions due to the ionization process, are obtained under several assumptions and approximations. It is found that theoretical results are consistent with the Giotto observations recently reported by Neugebauer et al. (1989).

Ziebell, L. F.↗
DOE OSTI2023

MATBOX, an Open-Source Microstructure Analysis Toolbox for Meshing, Generation, Segmentation, and Characterization of 3D Heterogenous Volumes

Battery performance is strongly correlated with electrode microstructural properties. To account for its impact, lithium-ion battery (LIB) models either abstract the microstructural heterogeneity of composite electrodes using effective macroscopic properties (macro- or meso- scale models) or directly solve the system of equations on the microstructure geometry or mesh (microstructure-scale models). Therefore, to be adequate, both families of models require information from the microstructure geometry, which can be provided by the numerical tool presented in this work. MATBOX is a MATLAB open-source application [1] developed by NREL for performing various microstructure-related tasks including microstructure numerical generation, image filtering and microstructure segmentation, microstructure characterization and correlation, visualization, and microstructure meshing. MATBOX was originally developed for the analysis of LIB electrode microstructures; however, the algorithms provided by the toolbox are widely applicable to other heterogeneous materials. The toolbox provides a user-friendly experience thanks to a Graphical-User Interface, requires no coding by the user, and is well documented. This presentation will illustrate various MATBOX features for the characterization of a LIB electrode, including a fully automated Representative Volume Element (RVE) analysis, the numerical generation of complex 'virtual' microstructure, including dual-layer electrodes and carbon-binder additive phase, and the meshing of a complex NMC/graphite full cell microstructure suitable for 3D finite-element modeling. Other modules (segmentation, visualization, and correlation) will be briefly presented. Thanks to its modular, open-source approach, MATBOX can easily incorporate third-party algorithms to eventually build a standard in the field that will benefit the whole scientific community. Effective diffusion coefficient [2], additive phase numerical generation [3], and meshing [4] third-party algorithms have been already integrated in the toolbox with more to come.

DIRECT ENERGY CONVERSION,MATHEMATICS AND COMPUTING↗
NASA NTRS1960

Heat Diffusion in Gases, Including Effects of Chemical Reaction

The diffusion of heat through gases is treated where the coefficients of thermal conductivity and diffusivity are functions of temperature. The diffusivity is taken proportional to the integral of thermal conductivity, where the gas is ideal, and is considered constant over the temperature interval in which a chemical reaction occurs. The heat diffusion equation is then solved numerically for a semi-infinite gas medium with constant initial and boundary conditions. These solutions are in a dimensionless form applicable to gases in general, and they are used, along with measured shock velocity and heat flux through a shock reflecting surface, to evaluate the integral of thermal conductivity for air up to 5000 degrees Kelvin. This integral has the properties of a heat flux potential and replaces temperature as the dependent variable for problems of heat diffusion in media with variable coefficients. Examples are given in which the heat flux at the stagnation region of blunt hypersonic bodies is expressed in terms of this potential.

Hansen, C. Frederick↗
DOE OSTI2022

Anion diffusion in two-dimensional halide perovskites

Commercialization of halide perovskites in the semiconductor industry is hindered by their short-term stability. The instability of perovskites is closely interlinked with ionic diffusion. Historically, attempts to study diffusion in 2D perovskites mostly utilized electrical characterizations, but these characterizations pose a challenge in deconvoluting the impact of device architecture, interlayers, and ionic species. In this Perspective, we focus our attention on simple optical characterizations employed in the literature to investigate halide diffusion in 2D perovskites using lateral and vertical heterostructure platforms. We review the various synthesis techniques used for fabrication of halide perovskite heterostructures and discuss the qualitative and quantitative diffusion studies performed using these platforms. We discuss the numerical methods used to validate and supplement the experimental halide diffusion kinetics. Finally, we highlight the need to conduct further research on the impact of device operating conditions, lattice structure, and vacancy concentration on halide diffusion. Through this Perspective, we aim to emphasize the need of developing a comprehensive understanding of halide diffusion in perovskites for their successful deployment in optoelectronics.

36 MATERIALS SCIENCE↗
NASA NTRS2022

Diffusion in Melts and Magmas

Diffusion results from random motion of particles and entities. Diffusion in melts and magmas is due to thermally excited random motion of atoms, ions, and clusters, and plays a critical role in magmatic and volcanic processes. In melts and magmas, diffusion is one of the two mechanisms of mass transfer; the other being bulk flow (referred to as convection or advection). When both are present, diffusion refers to the dispersive motion relative to the mean bulk flow in a given reference frame (Richter et al. 1998). Diffusion plays critical roles in controlling magma mixing (Watson 1982; Koyaguchi 1985, 1989; Lesher 1994; Huber et al. 2009; Guo and Zhang 2020), mineral growth and dissolution rates in magmas (e.g., Watson 1982; Harrison and Watson 1983; Zhang et al. 1989; Newcombe et al. 2014; Macris et al. 2018), bubble growth and dissolution rate in magmas (Sparks 1978; Proussevitch and Sahagian 1998; Liu and Zhang 2000; Zhang 2013), and elemental and isotope fractionation during mineral 2 growth and dissolution (Jambon 1980; Richter et al. 1999, 2003; Watson and Muller 2009; Chopra et al. 2012; Watkins et al. 2014, 2017; Holycross et al. 2016, 2018). As a result, diffusion also plays an essential role in explosive volcanic eruptions and magma crystallization. Furthermore, diffusion has important applications in geospeedometry (Lasaga 1983, 1998; Zhang 1994, 2008; Trail et al. 2016; Zhang and Xu 2016). Experimental investigation of diffusion in geologically relevant silicate melts began to flourish in the 1970’s when micro-analytical measurements of diffusion profiles became available. In addition to the vast number of papers published since then, numerous books and reviews are available for diffusion in silicate melts. Hofmann et al. (1974) edited a book titled “Geochemical Transport and Kinetics” published by Carnegie Institution of Washington. This was the first landmark book summarizing the field. Lasaga and Kirkpatrick (1981) edited a book “Kinetics of Geochemical Processes” as volume 8 of the Reviews in Mineralogy (later becoming Reviews in Mineralogy and Geochemistry) series. Zhang and Cherniak (2010) edited “Diffusion in Minerals and Melts” as volume 72 of Reviews of Mineralogy and Geochemistry series, in which one chapter focused on diffusion theory, five chapters on diffusion in silicate melts (Behrens 2010; Lesher 2010; Liang 2010; Zhang and Ni 2010; Zhang et al. 2010), and other chapters were on experimental, analytical, and computational methods, and diffusion in minerals. Several textbooks covered the principles and applications of diffusion theories (Kirkaldy and Young 1987; Shewmon 1989; Cussler 1997; Lasaga 1998; Zhang 2008; Vrentas and Vrentas 2016), and two classic books covered the mathematics of diffusion (Carslaw and Jaeger 1959; Crank 1975). In preparing for this review chapter, we thought carefully about what to cover for this vast field, and decided to briefly go through the fundamentals of diffusion (more complete review can be found in Chakraborty 1995; Zhang 2008, 2010) and solutions to often-encountered 3 diffusion problems, and then focus on post-2010 diffusion studies on silicate melts and magmas. Here, melts refer to (mostly natural) silicate liquid, and magmas refer to crystal-bearing and/or bubble-bearing melts in which the continuous phase is the melt. There is a large body of work on diffusion in glasses, especially in the materials science literature, which is not covered in this review.

Youxue Zhang↗
NASA NTRS1994

Stochastic acceleration of energetic ions in Jupiter's magnetosphere

An equation governing the combined radial diffusion and stochastic acceleration of super-Alfvenic ions by magnetohydrodynamic waves in Jupiter's outer magnetosphere is derived. The formulation is based upon a total energy invariant of the adiabatic transport which applies to an isotropic distribution undergoing rapid pitch angle scattering by waves. An analytic solution to the double diffusion equation is obtained and numerical results are presented for two models of ion injection. The first model assumes S(+) and O(+) are injected throughout a broad region of space through photoionization of Jupiter's magnetospheric neutral wind and obtain an initial energy corresponding to the local corotation energy after pickup by the planetary magnetic field. The second model assumes a monoenergetic distribution of energetic protons is implanted in the middle magnetosphere by the action of field-aligned potential drops in Jupiter's auroral ionosphere. For both light and heavy ions the injection process creates a seed population of particles which are further accelerated nonadiabatically by the MHD waves and adiabatically through radial diffusion. A comparison of the theoretical results with a recent data analysis of Voyager low-energy charged particle measurements is made with very good agreement, thus providing a rigorous quantitative account of and definitive explanation for the high-energy ion component of Jupiter's magnetosphere.

Barbosa, D. D.↗
NASA NTRSDate not supplied

Experimental and Numerical Investigation of the NASA High Efficiency Centrifugal Compressor Vaned Stage Geometry and Aerodynamic Performance

Since its inception in the early 2010s, the NASA High Efficiency Centrifugal Compressor (HECC) has been enigmatic for the propulsion research community: experimental data and numerical simulations of the stage have generally not aligned in their quantifications of performance metrics. Typically, the fault for these disagreements is assigned to the numerical simulations as the simulations are models, and models are inherently incomplete representations of the experiment. This “incompleteness” may manifest in assumptions regarding roughness or heat transfer, simplifications of the flow path (i.e., neglecting bleed flows), or the oft-scapegoated turbulence model. In the case of HECC, recent work showed unexpected discrepancies between the intended impeller geometry defined in the design report (termed Design-Intent) and the manufactured impeller used in the experimental campaigns (termed As-Manufactured). That work used numerical simulations to establish that the geometric differences between the Design-Intent and As-Manufactured impellers were significant enough to result in drastically different performance predictions for the HECC vaneless diffuser configuration. The Design-Intent impeller simulations over predicted the performance relative to the experiment, whereas the As-Manufactured simulations better represented the experimental data, both in terms of one-dimensional performance metrics and spanwise flow profiles. This effort expands on that work by examining in detail the geometry and aerodynamic performance of the HECC vaned diffuser configuration. Further differences between the Design-Intent and As-Manufactured geometries have been discovered in the vaned diffuser and exit guide vanes, and these differences are documented herein. The summations of the geometric differences for all of the components were used to create two numerical models of HECC vaned diffuser configuration: the Design-Intent simulations which are generated from the original geometry definitions given in the design report and the As-Manufactured simulations which are the best available representation of the manufactured compressor hardware used in the experimental test campaigns. In congruence with the earlier vaneless diffuser work, the numerical predictions of the Design-Intent choked mass flow rate, total pressure ratio, and efficiency were notably greater than that of the As-Manufactured simulations. To increase confidence in the experimental dataset, measurements from a recent test campaign conducted in 2024 are used to validate the original experimental data acquired from 2012 to 2014 with good repeatability overall, especially considering the passage of time and differences in the data acquisition systems between the test campaigns. Both numerical simulations were then extensively evaluated against the experimental data. The As-Manufactured simulations provided better estimates of the stage performance than the Design-Intent cases in terms of most performance metrics. Nonetheless, more detailed results still show opportunities for improvement. Despite a more accurate representation of the physical hardware, characterization of the impeller work input remains a challenge even for rigorously developed numerical models.

vaned diffuser↗
DOE OSTI2019

Effect of Diffusion Length in Modeling of Equiaxed Dendritic Solidification under Buoyancy Flow in a Configuration of Hebditch–Hunt Experiment

Modeling of equiaxed solidification is vital for understanding the solidification process of metallic alloys. In this work, an extended literature review is given for the models currently used for equiaxed solidification simulations. Based on this analysis, we present a three-phase multiscale equiaxed solidification model in which some approximations regarding solute transport at microscopic scale are put together in a new way and incorporated into macroscopic transport equations. For the latter, a term relating to the momentum exchange between the two phases is revised, and a modification for the grain packing algorithm is proposed. A modernized model for equiaxed dendrite growth is tested using a case of solidification of Sn-5 wt pct Pb alloy in a parallelepiped cavity that mimics the Hebditch–Hunt experiment. The results obtained using two approaches to calculate diffusion length are presented and compared both with each other and with numerical results from elsewhere. It is demonstrated that diffusion length has a crucial effect on the final segregation pattern.

36 MATERIALS SCIENCE↗
NASA NTRS1971

Numerical calculation of turbulent boundary layers including suction or injection with binary diffusion.

The modelling of turbulent transport processes with an eddy viscosity representation of the inner, wall region for cases with suction or injection is considered in detail. It is shown how the van Driest solid wall model can be logically extended for such cases. Further, a new inductive model is developed. Finite difference calculations using the new inductive model are shown to be in good agreement with experimental data over a wide range of conditions. Lastly, a new method for calculating the wall shear within a finite difference solution is developed. This procedure greatly relaxes the requirement for fine grid spacing in the vertical direction.

Schetz, J. A.↗
NASA NTRSDate not supplied

Experimental and Numerical Investigation of the NASA High Efficiency Centrifugal Compressor Vaned Stage Geometry and Aerodynamic Performance

Since its inception in the early 2010s, the NASA High Efficiency Centrifugal Compressor (HECC) has been enigmatic for the propulsion research community: experimental data and numerical simulations of the stage have generally not aligned in their quantifications of performance metrics. Typically, the fault for these disagreements is assigned to the numerical simulations as the simulations are models, and models are inherently incomplete representations of the experiment. This “incompleteness” may manifest in assumptions regarding roughness or heat transfer, simplifications of the flow path (i.e., neglecting bleed flows), or the oft-scapegoated turbulence model. In the case of HECC, recent work showed unexpected discrepancies between the intended impeller geometry defined in the design report (termed Design-Intent) and the manufactured impeller used in the experimental campaigns (termed As-Manufactured). That work used numerical simulations to establish that the geometric differences between the Design-Intent and As-Manufactured impellers were significant enough to result in drastically different performance predictions for the HECC vaneless diffuser configuration. The Design-Intent impeller simulations over predicted the performance relative to the experiment, whereas the As-Manufactured simulations better represented the experimental data, both in terms of one-dimensional performance metrics and spanwise flow profiles. This effort expands on that work by examining in detail the geometry and aerodynamic performance of the HECC vaned diffuser configuration. Further differences between the Design-Intent and As-Manufactured geometries have been discovered in the vaned diffuser and exit guide vanes, and these differences are documented herein. The summations of the geometric differences for all of the components were used to create two numerical models of HECC vaned diffuser configuration: the Design-Intent simulations which are generated from the original geometry definitions given in the design report and the As-Manufactured simulations which are the best available representation of the manufactured compressor hardware used in the experimental test campaigns. In congruence with the earlier vaneless diffuser work, the numerical predictions of the Design-Intent choked mass flow rate, total pressure ratio, and efficiency were notably greater than that of the As-Manufactured simulations. To increase confidence in the experimental dataset, measurements from a recent test campaign conducted in 2024 are used to validate the original experimental data acquired from 2012 to 2014 with good repeatability overall, especially considering the passage of time and differences in the data acquisition systems between the test campaigns. Both numerical simulations were then extensively evaluated against the experimental data. The As-Manufactured simulations provided better estimates of the stage performance than the Design-Intent cases in terms of most performance metrics. Nonetheless, more detailed results still show opportunities for improvement. Despite a more accurate representation of the physical hardware, characterization of the impeller work input remains a challenge even for rigorously developed numerical models.

centrifugal compressor↗
NASA NTRS2001

Cool Flame Quenching

Cool flame quenching distances are generally presumed to be larger than those associated with hot flames, because the quenching distance scales with the inverse of the flame propagation speed, and cool flame propagation speeds are often times slower than those associated with hot flames. To date, this presumption has never been put to a rigorous test, because unstirred, non-isothermal cool flame studies on Earth are complicated by natural convection. Moreover, the critical Peclet number (Pe) for quenching of cool flames has never been established and may not be the same as that associated with wall quenching due to conduction heat loss in hot flames, Pe approx. = 40-60. The objectives of this ground-based study are to: (1) better understand the role of conduction heat loss and species diffusion on cool flame quenching (i.e., Lewis number effects), (2) determine cool flame quenching distances (i.e, critical Peclet number, Pe) for different experimental parameters and vessel surface pretreatments, and (3) understand the mechanisms that govern the quenching distances in premixtures that support cool flames as well as hot flames induced by spark-ignition. Objective (3) poses a unique fire safety hazard if conditions exist where cool flame quenching distances are smaller than those associated with hot flames. For example, a significant, yet unexplored risk, can occur if a multi-stage ignition (a cool flame that transitions to a hot flame) occurs in a vessel size that is smaller than that associated with the hot quenching distance. To accomplish the above objectives, a variety of hydrocarbon-air mixtures will be tested in a static reactor at elevated temperature in the laboratory (1g). In addition, reactions with chemical induction times that are sufficiently short will be tested aboard NASA's KC-135 microgravity (mu-g) aircraft. The mu-g results will be compared to a numerical model that includes species diffusion, heat conduction, and a skeletal kinetic mechanism, following the work on diffusion-controlled cool flames by Fairlie et,al., 2000.

Pearlman, Howard↗
NASA NTRS1984

Direct simulations of chemically reacting turbulent mixing layers

The report presents the results of direct numerical simulations of chemically reacting turbulent mixing layers. The work consists of two parts: (1) the development and testing of a spectral numerical computer code that treats the diffusion reaction equations; and (2) the simulation of a series of cases of chemical reactions occurring on mixing layers. The reaction considered is a binary, irreversible reaction with no heat release. The reacting species are nonpremixed. The results of the numerical tests indicate that the high accuracy of the spectral methods observed for rigid body rotation are also obtained when diffusion, reaction, and more complex flows are considered. In the simulations, the effects of vortex rollup and smaller scale turbulence on the overall reaction rates are investigated. The simulation results are found to be in approximate agreement with similarity theory. Comparisons of simulation results with certain modeling hypotheses indicate limitations in these hypotheses. The nondimensional product thickness computed from the simulations is compared with laboratory values and is found to be in reasonable agreement, especially since there are no adjustable constants in the method.

Riley, J. J.↗
NASA NTRS1983

Analytical and Numerical Solution for a Solidifying Liquid Alloy Slab

Numerical and analytical solutions are presented for the temperature and concentration distributions during the solidification of a binary liquid alloy slab. The slab is taken to be of a finite depth but infinite in the horizontal direction. The solidification process is started by withdrawing a fixed amount of heat from the lower surface of the slab. The upper surface of the slab is subjected to both radiation and convective conditions. The solution gives the concentration and temperature profiles and the interface position as a function of time. Due to the smallness of the mass diffusion coefficient in the solid, the numerical solution method breaks down whenever the ratio of the diffusivities in the solid and the liquid falls below a certain value. An analytical method is developed which gives accurate solution for any value of the diffusivity ratio.

Antar, B. N.↗
NASA NTRS2003

Elevated Temperature Creep Deformation in Solid Solution <001> NiAL-3.6Ti Single Crystals

The 1100 to 1500 K slow plastic strain rate compressive properties of <001> oriented NiAl-3.6Ti single crystals have been measured, and the results suggests that two deformation processes exist. While the intermediate temperature/faster strain rate mechanism is uncertain, plastic flow at elevated temperature/slower strain rates in NiAl-3.6Ti appears to be controlled by solute drag as described by the Cottrell-Jaswon solute drag model for gliding b = a(sub 0)<101> dislocations. While the calculated activation energy of deformation is much higher (approximately 480 kJ/mol) than the activation energy for diffusion (approximately 290 kJ/mol) used in the Cottrell-Jaswon creep model, a forced temperature compensated - power law fit using the activation energy for diffusion was able to adequately (greater than 90%) predict the observed creep properties. Thus we conclude that the rejection of a diffusion controlled mechanism can not be simply based on a large numerical difference between the activation energies for deformation and diffusion.

Whittenberger, J. Daniel↗
DOE OSTI2020

Analyzing the co-evolution of green technology diffusion and consumers’ pro-environmental attitudes: An agent-based model

Massive diffusion of green technologies is significant for building a cleaner world. However, the process of technology diffusion is usually slow and complex. An in-depth understanding of the mechanism regarding green technology diffusion is an essential precondition for effectively stimulating this process. Green technology diffusion heavily involves both social and technological changes. Although existing studies have provided rich knowledge about identifying critical drivers and barriers affecting green technology diffusion, research that considers consumers’ pro-environmental attitudes and green technology diffusion as an evolving system is still sparse. Aiming at exploring the co-evolution of consumers’ pro-environmental attitudes and green technology diffusion, this paper builds an agent-based model that integrates the relative agreement model with technology diffusion theories to conduct a sequence of controlled numerical experiments, which progressively unveil how attitudinal and technological factors impact green technology diffusion. The main findings include that (1) improving consumers’ pro-environmental attitudes is prominently beneficial to green technology diffusion and maturation; (2) technology maturity has very limited impact on consumers’ first-time purchases but significantly affects consumers’ satisfaction, which would further impact consumers’ repeat purchases; (3) consumers that do not support green technologies frequently emerge during the evolution of attitudes (despite the high technology maturity), which corresponds to the emergence of the anti-environmental groups observed in the real world; (4) active interactions between non-adopters enable their attitudes to converge, which results in the polarization of consumers’ attitudes.

ABM↗
DOE OSTI2021

Deep Learning Approaches to Surrogates for Solving the Diffusion Equation for Mechanistic Real-World Simulations

In many mechanistic medical, biological, physical, and engineered spatiotemporal dynamic models the numerical solution of partial differential equations (PDEs), especially for diffusion, fluid flow and mechanical relaxation, can make simulations impractically slow. Biological models of tissues and organs often require the simultaneous calculation of the spatial variation of concentration of dozens of diffusing chemical species. One clinical example where rapid calculation of a diffusing field is of use is the estimation of oxygen gradients in the retina, based on imaging of the retinal vasculature, to guide surgical interventions in diabetic retinopathy. Furthermore, the ability to predict blood perfusion and oxygenation may one day guide clinical interventions in diverse settings, i.e., from stent placement in treating heart disease to BOLD fMRI interpretation in evaluating cognitive function (Xie et al., 2019; Lee et al., 2020). Since the quasi-steady-state solutions required for fast-diffusing chemical species like oxygen are particularly computationally costly, we consider the use of a neural network to provide an approximate solution to the steady-state diffusion equation. Machine learning surrogates, neural networks trained to provide approximate solutions to such complicated numerical problems, can often provide speed-ups of several orders of magnitude compared to direct calculation. Surrogates of PDEs could enable use of larger and more detailed models than are possible with direct calculation and can make including such simulations in real-time or near-real time workflows practical. Creating a surrogate requires running the direct calculation tens of thousands of times to generate training data and then training the neural network, both of which are computationally expensive. Often the practical applications of such models require thousands to millions of replica simulations, for example for parameter identification and uncertainty quantification, each of which gains speed from surrogate use and rapidly recovers the up-front costs of surrogate generation. We use a Convolutional Neural Network to approximate the stationary solution to the diffusion equation in the case of two equal-diameter, circular, constant-value sources located at random positions in a two-dimensional square domain with absorbing boundary conditions. Such a configuration caricatures the chemical concentration field of a fast-diffusing species like oxygen in a tissue with two parallel blood vessels in a cross section perpendicular to the two blood vessels. To improve convergence during training, we apply a training approach that uses roll-back to reject stochastic changes to the network that increase the loss function. The trained neural network approximation is about 1000 times faster than the direct calculation for individual replicas. Because different applications will have different criteria for acceptable approximation accuracy, we discuss a variety of loss functions and accuracy estimators that can help select the best network for a particular application. We briefly discuss some of the issues we encountered with overfitting, mismapping of the field values and the geometrical conditions that lead to large absolute and relative errors in the approximate solution.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Partitioned exponential methods for coupled multiphysics systems

Multiphysics problems involving two or more coupled physical phenomena are ubiquitous in science and engineering. This work develops a new partitioned exponential approach for the time integration of multiphysics problems. After a possible semi-discretization in space, the class of problems under consideration is modeled by a system of ordinary differential equations where the right-hand side is a summation of two component functions, each corresponding to a given set of physical processes. The partitioned-exponential methods proposed herein evolve each component of the system via an exponential integrator, and information between partitions is exchanged via coupling terms. Here, the traditional approach to constructing exponential methods, based on the variation-of-constants formula, is not directly applicable to partitioned systems. Rather, our approach to developing new partitioned-exponential families is based on a general-structure additive formulation of the schemes. Two method formulations are considered, one based on a linear-nonlinear splitting of the right hand component functions, and another based on approximate Jacobians. The paper develops classical (non-stiff) order conditions theory for partitioned exponential schemes based on particular families of T-trees and B-series theory. Several practical methods of third order are constructed that extend the Rosenbrock-type and EPIRK families of exponential integrators. Several implementation optimizations specific to the application of these methods to reaction-diffusion systems are also discussed. Numerical experiments reveal that the new partitioned-exponential methods can perform better than traditional unpartitioned exponential methods on some problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗