Search NASA⌕ Search

SEARCH · Search NASA

Results for “parallel processing”

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 1,315 records · Page 73
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Reducing software mass through behavior control

Attention is given to the tradeoff between communication and computation as regards a planetary rover (both these subsystems are very power-intensive, and both can be the major driver of the rover's power subsystem, and therefore the minimum mass and size of the rover). Software techniques that can be used to reduce the requirements on both communciation and computation, allowing the overall robot mass to be greatly reduced, are discussed. Novel approaches to autonomous control, called behavior control, employ an entirely different approach, and for many tasks will yield a similar or superior level of autonomy to traditional control techniques, while greatly reducing the computational demand. Traditional systems have several expensive processes that operate serially, while behavior techniques employ robot capabilities that run in parallel. Traditional systems make extensive world models, while behavior control systems use minimal world models or none at all.

Miller, David P.↗
NASA NTRS1992

Particle acceleration at comets

This paper compares calculated and measured energy spectra of implanted H(+) and O(+) ions on the assumption that the pickup geometry is quasi-parallel and about 1 percent of the waves generated by the cometary pickup process propagates backward (towards the comet). The model provides a good description of the implanted O(+) and H(+) energy distribution near the pickup energies.

Gombosi, Tamas I.↗
NASA NTRS2001

A Debugger for Computational Grid Applications

This viewgraph presentation gives an overview of a debugger for computational grid applications. Details are given on NAS parallel tools groups (including parallelization support tools, evaluation of various parallelization strategies, and distributed and aggregated computing), debugger dependencies, scalability, initial implementation, the process grid, and information on Globus.

Hood, Robert↗
NASA NTRS1994

Dermal Aged and Fetal Fibroblasts Realign in Response to Mechanical Strain

Integrins specifically recognize and bind extracellular matrix components, providing physical anchor points and functional setpoints. Focal adhesion complexes, containing integrin and cytoskeletal proteins, are potential mechanoreceptors, poised to distribute applied forces through the cytoskeleton. Pursuing the hypothesis that cells both perceive and respond to external force, we applied a stretch/relaxation regimen to normal human fetal and aged dermal fibroblast monolayers cultured on flexible membranes. The frequency and magnitude of the applied force is precisely controlled by the Flexercell Unit(Trademark). A protocol of stretch (20% elongation of the monolayer) at a frequency of 6 cycles/min caused a progressive change from a randomly distributed pattern of cells to a symmetric, radial distribution with cells aligned parallel to the applied force. We have coined the term 'orienteering' as the process of active alignment of cells in response to applied force. Cytochalasin D was added in graded doses to investigate the role of the actin cytoskeleton in force perception and transmission. A clear dose response was found; at high concentrations orienteering was abolished; and the drug's impact was reversible. The two cell strains used were similar in their alignment behavior and in their responses to cytochalasin D. Orienteering was influenced by cell density, and the cell strains studied differed in this respect. Fetal cells, unlike their aged counterparts, failed to orient at high cell density. In both cell strains, mid-density cultures aligned rapidly and sparse cultures lagged. These results indicate that both cell-cell adhesion and cytoskeleton integrity are critical in mediating the orienteering response. Differences between these two cell strains may relate to their expression of extracellular matrix molecules (fibronectin, collagen type 1) integrins and their relative binding affinities.

Sawyer, Christine↗
NASA NTRSDate not supplied

Particle Interaction Physics Model Formulation for Plume-Surface Interaction Erosion and Cratering

As part of the Game Changing Development (GCD) Program, funded by NASA’s Space Technology Mission Directorate (STMD), the development of simulation capability for the prediction of extra-terrestrial Plume Surface Interaction (PSI) environments has been undertaken by the Fluid Dynamics Branch at NASA/MSFC. The Predictive Simulation Capability (PSC) Element is focused on creating simulation capability for the reliable and accurate prediction of PSI in Martian (~650 Pa) and Lunar (vacuum) ambient environments. In addition to the predictive simulation capability, the GCD Program also contains a companion Ground Testing Element for development of focused datasets for validation of predictive capability as well as a Flight-focused Instrumentation Element. This paper will present the status of implementing and maturing particle-particle interaction constituent physics models essential in simulating the landing surface granular material flow under PSI effects. This gas-particle multi-phase interaction modeling of plume impingement flow on the extra-terrestrial soil material is performed with the Gas-Granular Flow Solver (GGFS) addressed in a companion paper. The response of regolith particle flow induced by lander PSI requires accurate representation of the regolith granular material fluidic behavior and gas-granular interactions. The lunar regolith, as the extreme example, is poorly sorted with broad particle size distributions and large fines content. It has significant cohesion, due to interlocking particle shapes for the very jagged particles. The combination of particle shape and size distribution has been identified as major drivers in the complex particle flow response and resulting crater shape characteristics of extraterrestrial granular material. Constituent models for spherical particles can be formulated directly from particle kinetics theory. Complex particle shapes can be modeled by gluing together elemental spherical shapes into composite particles, requiring a Discrete Element Model (DEM) particle kinetics modeling approach to extract data and formulate constituent models. Mixture constituent models for poly-disperse mixtures (i.e, containing distribution of particle sizes) have recently been developed. The required non-spherical particle mixture granular material response closure models are then obtained through small-scale unit physics DEM simulations for the range of particle shapes, mixtures and packing densities. The granular material response closure models are then implemented in the Eulerian granular flow formulation. This DEM-based constituent model extraction process and formulation of poly-disperse particle mixtures has been successfully developed by small business and academic partners in the development of the Gas-Granular Flow Solver (GGFS) simulation program simulation framework. The currently implemented capabilities have reached the capability level of modeling bi-disperse, non-spherical particle mixtures is being continuously extended towards computational modeling of full range irregular particle mixtures. Under the GCD project, this technology is being further developed, transferred to NASA analysts, and matured towards application readiness. The predictive simulation capability team under the GCD project has acquired the modeling tools and processes of the DEM based constituent model formulation from the GGFS development team and is developing the capability to replicate the existing process. This is the first important step towards the ability of the NASA team to independently perform such model development in a production setting. Further efforts are underway to migrate the DEM based model simulation process performed with the academic based tools to more capable Open Source, highly parallelized simulation tools for efficient operation on NASA HPC assets. Evaluation of the currently implemented (such as mono-disperse and bi-disperse spherical and irregular shape particle constituent model applications) and continuously evolving full-range particle physics models in the GGFS tool is performed by the NASA team to advance application readiness of the simulations. Application testing for complex PSI erosions and cratering scenarios such as the Apollo LM is performed for axi-symmetric and full 3D simulations to aid the tool developers in achieving practical application readiness for NASA projects. Important validation and application testing will further be performed against experimental data generated under the GCD PSI project experimental component.

Peter A Liever↗
NASA NTRS2012

Radiation-Hard SpaceWire/Gigabit Ethernet-Compatible Transponder

A radiation-hard transponder was developed utilizing submicron/nanotechnology from IBM. The device consumes low power and has a low fabrication cost. This device utilizes a Plug-and-Play concept, and can be integrated into intra-satellite networks, supporting SpaceWire and Gigabit Ethernet I/O. A space-qualified, 100-pin package also was developed, allowing space-qualified (class K) transponders to be delivered within a six-month time frame. The novel, optical, radiation-tolerant transponder was implemented as a standalone board, containing the transponder ASIC (application specific integrated circuit) and optical module, with an FPGA (field-programmable gate array) friendly parallel interface. It features improved radiation tolerance; high-data-rate, low-power consumption; and advanced functionality. The transponder utilizes a patented current mode logic library of radiation-hardened-by-architecture cells. The transponder was developed, fabricated, and radhard tested up to 1 MRad. It was fabricated using 90-nm CMOS (complementary metal oxide semiconductor) 9 SF process from IBM, and incorporates full BIT circuitry, allowing a loop back test. The low-speed parallel LVCMOS (lowvoltage complementary metal oxide semiconductor) bus is compatible with Actel FPGA. The output LVDS (low-voltage differential signaling) interface operates up to 1.5 Gb/s. Built-in CDR (clock-data recovery) circuitry provides robust synchronization and incorporates two alarm signals such as synch loss and signal loss. The ultra-linear peak detector scheme allows on-line control of the amplitude of the input signal. Power consumption is less than 300 mW. The developed transponder with a 1.25 Gb/s serial data rate incorporates a 10-to-1 serializer with an internal clock multiplication unit and a 10-1 deserializer with internal clock and data recovery block, which can operate with 8B10B encoded signals. Three loop-back test modes are provided to facilitate the built-in-test functionality. The design is based on a proprietary library of differential current switching logic cells implemented in the standard 90-nm CMOS 9SF technology from IBM. The proprietary low-power LVDS physical interface is fully compatible with the SpaceWire standard, and can be directly connected to the SFP MSA (small form factor pluggable Multiple Source Agreement) optical transponder. The low-speed parallel interfaces are fully compatible with the standard 1.8 V CMOS input/output devices. The utilized proprietary annular CMOS layout structures provide TID tolerance above 1.2 MRad. The complete chip consumes less than 150 mW of power from a single 1.8-V positive supply source.

Katzman, Vladimir↗
NASA NTRS1997

Novel Highly Parallel and Systolic Architectures Using Quantum Dot-Based Hardware

VLSI technology has made possible the integration of massive number of components (processors, memory, etc.) into a single chip. In VLSI design, memory and processing power are relatively cheap and the main emphasis of the design is on reducing the overall interconnection complexity since data routing costs dominate the power, time, and area required to implement a computation. Communication is costly because wires occupy the most space on a circuit and it can also degrade clock time. In fact, much of the complexity (and hence the cost) of VLSI design results from minimization of data routing. The main difficulty in VLSI routing is due to the fact that crossing of the lines carrying data, instruction, control, etc. is not possible in a plane. Thus, in order to meet this constraint, the VLSI design aims at keeping the architecture highly regular with local and short interconnection. As a result, while the high level of integration has opened the way for massively parallel computation, practical and full exploitation of such a capability in many applications of interest has been hindered by the constraints on interconnection pattern. More precisely. the use of only localized communication significantly simplifies the design of interconnection architecture but at the expense of somewhat restricted class of applications. For example, there are currently commercially available products integrating; hundreds of simple processor elements within a single chip. However, the lack of adequate interconnection pattern among these processing elements make them inefficient for exploiting a large degree of parallelism in many applications.

Fijany, Amir↗
NASA NTRS2002

Automated CFD Parameter Studies on Distributed Parallel Computers

The objective of the current work is to build a prototype software system which will automated the process of running CFD jobs on Information Power Grid (IPG) resources. This system should remove the need for user monitoring and intervention of every single CFD job. It should enable the use of many different computers to populate a massive run matrix in the shortest time possible. Such a software system has been developed, and is known as the AeroDB script system. The approach taken for the development of AeroDB was to build several discrete modules. These include a database, a job-launcher module, a run-manager module to monitor each individual job, and a web-based user portal for monitoring of the progress of the parameter study. The details of the design of AeroDB are presented in the following section. The following section provides the results of a parameter study which was performed using AeroDB for the analysis of a reusable launch vehicle (RLV). The paper concludes with a section on the lessons learned in this effort, and ideas for future work in this area.

Rogers, Stuart E.↗
DOE OSTI2025

Multi-physics Topology OPtimization and Additive Manufacturing for High-temperature Heat Exchangers

This research significantly advances the understanding of high-temperature heat exchanger design through an integrated approach that combines topology optimization (TO), triply periodic minimal surface (TPMS) structures, additive manufacturing (AM) and thermohydraulic testing. Each of these components contributes uniquely to a unified, high-performance design, fabrication and testing workflow. Topology optimization serves as the foundation of the design methodology by providing a systematic way to determine the most effective material layout for separating hot and cold fluids while maximizing thermal performance. The researchers introduced a novel three-material optimization framework using two density fields to represent hot fluid, cold fluid, and solid domains. This approach enables automated discovery of optimal shapes and flow paths that cannot be intuitively designed, especially under constraints imposed by manufacturing technologies. Furthermore, constraints such as minimal wall thickness and overhang angles were embedded into the optimization process, ensuring that resulting designs are not only thermally efficient but also manufacturable using modern additive techniques. In parallel, the study delves into the use of Gyroid-based TPMS geometries for constructing the core of the heat exchanger. TPMS structures are known for their high surface area, excellent fluid mixing capabilities, and minimal pressure drop characteristics. The researchers applied a data-driven modeling framework using Heteroscedastic Sparse Gaussian Process Regression (HSGPR) combined with genetic algorithms. This allowed for the rapid evaluation and optimization of key geometric parameters such as frequency, iso-value, and phase shift. The result was a set of Gyroid structures tailored for high heat transfer and low flow resistance, demonstrating clear improvements over conventional straight-channel designs. After the designing process, additive manufacturing played a critical role by turning these highly complex, optimized geometries into physical components. Utilizing Laser Powder Bed Fusion (LPBF) with Haynes 282, the study demonstrated the feasibility of fabricating these heat exchangers at high precision. Post-processing methods, including dilation-erosion operations, were applied to ensure local features adhered to self-supporting constraints. The fabricated structures were then subjected to thermohydraulic testing under conditions representative of supercritical CO 2 Brayton cycles, validating the predicted performance and confirming the viability of the full design-to-fabrication pipeline. Finally, thermohydraulic testing across the above studies served as a crucial experimental validation of advanced heat exchanger. Under consistent high-temperature and high-pressure conditions using supercritical CO 2 , the testing demonstrated that both TO and Gyroid-based TPMS designs significantly outperformed conventional straight-channel HXs. The TO design achieved a 115% increase in UA and NTU and a 27.6% boost in gravimetric power density, while the data-driven optimized Gyroid design delivered a 166% increase in UA and NTU and improved effectiveness from 68.7% to 86.1%. These results validate the simulation models, confirm the manufacturability of complex geometries under AM constraints, and provide key insights into design-performance trade-offs, thereby advancing the development of high-efficiency, compact heat exchangers for extreme environments.

36 MATERIALS SCIENCE↗
NASA NTRS1997

Description and Analysis of Core Samples: The Lunar Experience

Although no samples yet have been returned from a comet, extensive experience from sampling another solar system body, the Moon, does exist. While, in overall structure, composition, and physical properties the Moon bears little resemblance to what is expected for a comet, sampling the Moon has provided some basic lessons in how to do things which may be equally applicable to cometary samples. In particular, an extensive series of core samples has been taken on the Moon, and coring is the best way to sample a comet in three dimensions. Data from cores taken at 24 Apollo collection stations and 3 Luna sites have been used to provide insight into the evolution of the lunar regolith. It is now well understood that this regolith is very complex and reflects gardening (stirring of grains by micrometeorites), erosion (from impacts and solar wind sputtering), maturation (exposure on the bare lunar surface to solar winds ions and micrometeorite impacts) and comminution of coarse grains into finer grains, blanket deposition of coarse-grained layers, and other processes. All of these processes have been documented in cores. While a cometary regolith should not be expected to parallel in detail the lunar regolith, it is possible that the upper part of a cometary regolith may include textural, mineralogical, and chemical features which reflect the original accretion of the comet, including a form of gardening. Differences in relative velocities and gravitational attraction no doubt made this accretionary gardening qualitatively much different than the lunar version. Furthermore, at least some comets, depending on their orbits, have been subjected to impacts of the uppermost surface by small projectiles at some time in their history. Consequently, a more recent post-accretional gardening may have occurred. Finally, for comets which approach the sun, large scale erosion may have occurred driven by gas loss. The uppermost material of these comets may reflect some of the features of this erosional process, such as crust formation, and variations with depth might be expected. Overall, the upper few meters of a comet may be as complex in their own way as the upper few meters of the lunar regolith have proven to be, and by analogy, detailed studies of core samples containing this depth information will be needed to understand these processes and the details of the accretional history and the subsequent alteration history of comets.

McKay, David S.↗
NASA NTRS1992

Two dimensional PIC simulations of plasma heating by the dissipation of Alfven waves

Two dimensional plasma particle simulations of the evolution of large amplitude circularly polarized Alfven waves propagating parallel to the magnetic field show that the waves decay via both one- and two- dimensional parametric decay instabilities. For parameters studied, one-dimensional processes dominate the simulations, but two-dimensional decay processes, including the recently predicted filamentation instability are also observed. The daughter waves generated by the parametric decay are primarily damped by the ions, leading to ion heating. The parametric decay processes efficiently convert the ordered fluid ion motion in the Alfven wave into ion thermal energy. These processes may be important for the dissipation of Alfven waves in the solar wind, the corona and other space plasma environments. The computations were performed on the Intel Touchstone parallel supercomputer.

Liewer, P. C.↗
DOE OSTI2025

Mechanical properties, strain hardening, and fracture behavior of ultrasonic additively manufactured Zircaloy-4 after low-temperature neutron irradiation

Ultrasonic additive manufacturing (UAM) is a solid-state, layer-by-layer advanced manufacturing process that has the potential to create custom spatially controlled composites with embedded wires and sensors for nuclear component manufacture. For this work, to assess the feasibility of using UAM for nuclear-relevant materials research, the technique was used to produce a 3.5-mm-thick Zircaloy-4 plate for irradiation testing. The UAM Zircaloy-4 specimens were irradiated in the High Flux Isotope Reactor at a target irradiation temperature of 117 °C to 2.9 displacements per atom (dpa) to assess differences in irradiation-hardening behavior as a function of alloy processing path. The UAM and reference baseplate (BP) materials increased in yield strength by 372±27 MPa and 346±21 MPa, respectively, and both suffered significant reductions in uniform and total elongation attributed to irradiation hardening at low-temperature. Although the materials had similar nanoscale defect structures, including nanoscale black dot/loop features and strain-induced dislocation channels, the UAM material’s processing-related defects resulted in accelerated strain localization and failure as demonstrated by lower post-irradiation uniform elongation of UAM specimens (0.5 %) compared to BP (1.5 %) material. The UAM material also showed considerable anisotropy in mechanical response due to crack propagation along weld boundaries, resulting in differences in strength & ductility when tested parallel and perpendicular to the prior UAM build orientation. Therefore, although the fundamental irradiation response of UAM-processed Zircaloy-4 was phenomenologically comparable to that of BP reference material, additional optimization of the UAM processing is needed to produce irradiation-resistant and nuclear-relevant materials.

Digital image correlation↗
NASA NTRS1999

An Adaptive Kalman Filter Using a Simple Residual Tuning Method

One difficulty in using Kalman filters in real world situations is the selection of the correct process noise, measurement noise, and initial state estimate and covariance. These parameters are commonly referred to as tuning parameters. Multiple methods have been developed to estimate these parameters. Most of those methods such as maximum likelihood, subspace, and observer Kalman Identification require extensive offline processing and are not suitable for real time processing. One technique, which is suitable for real time processing, is the residual tuning method. Any mismodeling of the filter tuning parameters will result in a non-white sequence for the filter measurement residuals. The residual tuning technique uses this information to estimate corrections to those tuning parameters. The actual implementation results in a set of sequential equations that run in parallel with the Kalman filter. A. H. Jazwinski developed a specialized version of this technique for estimation of process noise. Equations for the estimation of the measurement noise have also been developed. These algorithms are used to estimate the process noise and measurement noise for the Wide Field Infrared Explorer star tracker and gyro.

Harman, Richard R.↗
NASA NTRS1994

Society of the plastic industry process emission initiatives

At first view, plastics process emissions research may not seem to have much bearing on outgassing considerations relative to advanced composite materials; however, several parallel issues and cross currents are of mutual interest. The following topics are discussed: relevance of plastics industry research to aerospace composites; impact of clean air act amendment requirements; scope of the Society of the Plastics Industry, Inc. activities in thermoplastic process emissions and reinforced plastics/composites process emissions; and utility of SPI research for advanced polymer composites audiences.

Mcdermott, Joseph↗
NASA NTRS1995

Fossilization Processes in Thermal Springs

To create a comparative framework for the study of ancient examples, we have been carrying out parallel studies of the microbial biosedimentology, taphonomy and geochemistry of modem and sub-Recent thermal spring deposits. One goal of the research is the development of integrated litho- and taphofacies models for siliceous and travertline sinters. Thermal springs are regarded as important environments for the origin and early evolution of life on Earth, and we seek to utilize information from the fossil record to reconstruct the evolution of high temperature ecosystems. Microbial contributions to the fabric of thermal spring sinters occur when population growth rates keep pace with, or exceed rates of inorganic precipitation, allowing for the development of continuous biofilms or mats. In siliceous thermal springs, microorganisms are typically entombed while viable. Modes of preservation reflect the balance between rates of organic matter degradation, silica precipitation and secondary infilling. Subaerial sinters are initially quite porous and permeable and at temperatures higher than about 20 C, organic materials are usually degraded prior to secondary infilling of sinter frameworks. Thus, organically-preserved microfossils are rare and fossil information consists of characteristic biofabrics formed by the encrustation and underplating of microbial mat surfaces. This probably accounts for the typically low total organic carbon values observed in thermal spring deposits. In mid-temperature, (approx. 35 - 59 C) ponds and outflows, the surface morphology of tufted Phormidium mats is preserved through mat underplating by thin siliceous: crusts. Microbial taxes lead to clumping of ceils and/or preferred filament orientations that together define higher order composite fabrics in thermal spring stromatolites (e.g. network, coniform, and palisade). At lower temperatures (less than 35 C), Calothrix mats cover shallow terracette pools forming flat carpets or pustular surfaces that produce palisade and "shrub" fabrics, respectively. At finer scales, composite fabrics are seen to consist distinctive associations of microstructures formed by the encrustation of individual cells and filaments. Composite fabrics survive the diagenetic transitions from primary opaline silica to quartz and are known from subaerial thermal spring deposits as old as Lower Carboniferous. However, fossil microorganisms tend to be rare in older deposits, and are usually preserved only where cells or sheaths have been stained by iron oxides. In subaqueous mineralizing springs at lower temperatures, early infilling leads to a more rapid and complete reduction in porosity and permeability. This process, along with the slower rates of microbial degradation at lower temperatures, creates a more favorable situation for organic matter preservation. Application of this taphonomic model to the Rhynie Chert, previously interpreted as subaerial, suggest it was probably deposited in a subaqueous spring setting at lower temperatures.

Farmer, Jack D.↗
NASA NTRS1991

Texture analysis of radiometric signatures of new sea ice forming in Arctic leads

Analysis of 33.6-GHz, high-resolution, passive microwave images suggests that new sea ice accumulating in open leads is characterized by a unique textural signature which can be used to discriminate new ice forming in this environment from adjacent surfaces of similar radiometric temperature. Ten training areas were selected from the data set, three of which consisted entirely of first-year ice, four entirely of multilayer ice, and three of new ice in open leads in the process of freezing. A simple gradient operator was used to characterize the radiometric texture in each training region in terms of the degree to which radiometric gradients are oriented. New ice in leads has a sufficiently high proportion of well-oriented features to distinguish it uniquely from first-year ice and multiyear ice. The predominance of well-oriented features probably reflects physical processes by which new ice accumulates in open leads. Banded structures, which are evident in aerial photographs of new ice, apparently give rise to the radiometric signature observed, in which the trend of brightness temperature gradients is aligned parallel to lead trends. First-year ice and multiyear ice, which have been subjected to a more random growth and process history, lack this banded structure and therefore are characterized by signatures in which well-aligned elements are less dominant.

Eppler, Duane T.↗
NASA NTRS1975

Feasibility study: Liquid hydrogen plant, 30 tons per day

The design considerations of the plant are discussed in detail along with management planning, objective schedules, and cost estimates. The processing scheme is aimed at ultimate use of coal as the basic raw material. For back-up, and to provide assurance of a dependable and steady supply of hydrogen, a parallel and redundant facility for gasifying heavy residual oil will be installed. Both the coal and residual oil gasifiers will use the partial oxidation process.

Source record↗
DOE OSTI2024

A meshless stochastic method for Poisson–Nernst–Planck equations

A plethora of biological, physical, and chemical phenomena involve transport of charged particles (ions). Its continuum-scale description relies on the Poisson–Nernst–Planck (PNP) system, which encapsulates the conservation of mass and charge. The numerical solution of these coupled partial differential equations is challenging and suffers from both the curse of dimensionality and difficulty in efficiently parallelizing. We present a novel particle-based framework to solve the full PNP system by simulating a drift–diffusion process with time- and space-varying drift. We leverage Green’s functions, kernel-independent fast multipole methods, and kernel density estimation to solve the PNP system in a meshless manner, capable of handling discontinuous initial states. The method is embarrassingly parallel, and the computational cost scales linearly with the number of particles and dimension. We use a series of numerical experiments to demonstrate both the method’s convergence with respect to the number of particles and computational cost vis-à-vis a traditional partial differential equation solver.

Chemistry↗