Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Microstructures, and Modelling”

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,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Focused Ion Beam Tomography of Alloy 617 Corroded in Molten Chloride Salt

Materials qualification of reactor structural materials is a critical step in rapid implementation of advanced nuclear reactor technologies, particularly to assess the corrosion performance in these designs. Accelerated qualification of reactor structural materials requires incorporating powerful computational toolsets, such as phase field modelling in the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework, to predict the evolution of structural materials due to corrosion. Accordingly, computational toolsets will require experimental data generated at appropriate length scales to validate accuracy. Focused ion beam (FIB) provides a high degree of control over manipulation of materials for analytical purposes, including capturing data on the evolution in the microstructure and elemental composition of materials at the mesoscale, an appropriate length scale for phase field modelling of intergranular diffusion phenomena using the MOOSE framework. For instance, the FEI Helios G4 UX dual beam plasma FIB microscope at the Irradiated Materials Characterization Laboratory (IMCL) is capable of backscatter diffraction (EBSD) and energy-dispersive x-ray spectroscopy (EDS) documenting the evolution in the microstructure and elemental composition, respectively. The Helios can perform EDS and EBSD three-dimensionally (3D) using tomography, which is then combined using different software packages to visualize 3D volumes correlating elemental composition to microstructural data. The purpose of this investigation was to develop a streamlined characterization and data processing workflow for 3D tomography studies on the FEI Helios G4 plasma FIB. The investigation is segmented into three parts: 1) Optimizing the data collection workflow, 2) identifying appropriate data processing and visualization software (i.e. DREAM.3D, MIPAR, and VGStudioMax), and 3) establishing an infrastructure for public release. The optimization of the data collection workflow is in collaboration with members of the U220 department to setup formal training on the tomography operation of the G4, through ThermoFisher Scientific, and exploring DREAM.3D, MIPAR, and VGStudioMax data processing/visualization software packages. VGStudioMax currently demonstrates the most promise for future use. Optimization of the data collection and processing workflow is still ongoing. A collaboration with INL High Performance Computing (HPC) established an open-source license for expediting the public release of FIB tomography datasets through HPC. FIB tomography data generated by the G4 will provide comprehensive data for validating 3D phase field mesoscale modelling tools within the MOOSE framework for accelerated qualification of reactor structural materials.

Copeland-Johnson, Trishelle↗
NASA NTRS2010

Microscopic Scale Simulation of the Ablation of Fibrous Materials

This slide presentation reviews the ablation by oxidation of carbon-fiber preforms impregnated in carbonized phenolic-formaldehyde matrix is modeled at microscopic scale. Direct numerical simulations show that the matrix ablates in volume leaving the carbon fibers exposed. This is due to the fact that the reactivity of carbonized phenolics is higher than the reactivity of carbon fibers. After the matrix is depleted, the fibers ablate showing progressive reduction of their diameter. The overall material recession occurs when the fibers are consumed. Two materials with the same carbon-fiber preform, density and chemical composition, but with different matrix distributions are studied. These studies show that at moderate temperatures (< 1000K) the microstructure of the material influences its recession rate; a fact that is not captured by current models that are based on chemical composition only. Surprisingly, the response of these impregnated-fiber materials is weakly dependent on the microstructure at very high temperatures (e.g., Stardust peak heating conditions: 3360K).

Lachaud, Jean↗
DOE OSTI2020

An Accelerated Creep Testing Program For Nickel Based Superalloys

In order to improve the energy production within the United States, new materials are developed with the goal of increased functionality. These materials are required to perform at extreme pressures and temperatures (up to and beyond 200MPa and 650â??) as well as be sustainable beyond 105 hours. Before developing these materials into components for power generation, mechanical testing is required to better understand the behavior of these materials at various loads and temperatures. The longest of these tests is the conventional creep test which is a real time test that can beyond 105 hours depending on service conditions. Hundreds of specimens must be tested before being approved for use in power generation, but it is extremely difficult to rapidly qualify when each material is required to exceed 105 hours making testing through conventional means unfeasible. While others have designed accelerated testing methods to help qualify materials, many of these experiments use an empirical approach that requires test data to fit already existing creep data. In this study, an accelerated creep test called the stepped isostress method (SSM) is then applied to wrought Inconel 718 (IN 718) specimens. A test matrix is designed using deformation mechanism maps, time-temperature-transformation diagrams, and time-temperature-precipitation diagrams, in order to ensure that microstructurally, the accelerated creep tests fail similarly to the conventional creep tests. The data from the stepped isostress method is then applied to the McVetty and Sin-hyperbolic (Sinh) constitutive models in order to generate a creep master curve. The results from the model fits are then compared to the conventional creep data done from experiments as well as from the NIMS database for Inconel 718. Optical microscopy and scanning electron microscopy are used to compare the evolution of the microstructure from the conventional creep and accelerated creep tests. It is found that using the Sinh model on the accelerated creep data provides reasonable creep rupture predictions when compared to the scatter band of creep tests. It was also observed that the microstructural evolution between the SSM experiments and creep tests were very similar in terms of the volume fraction of precipitates and phases.

Mach, Robert↗
NASA NTRSDate not supplied

Overview of Ablative TPS Modeling at NASA Ames

Over the past decade, NASA has invested in efforts to build predictive thermal protection system (TPS) material models from the micro-scale to the macro-scale. To complement the mission design cycle process and reduce the need for extensive testing, NASA is developing modeling and simulation tools that enable characterizing material properties and response to hot plasma experienced during atmospheric entry. Traditional material response and ablation modeling tools, such as the heritage code FIAT, and its multidimensional siblings, TITAN and 3dFIAT, are being complemented with newly developed software such as Icarus and PATO. Both of these programs are three-dimensional, finite-volume solvers that use unstructured meshes and 21st century programming paradigms to allow for efficient parallel simulations. FIAT and Icarus are also used for TPS sizing purposes. Today, these traditional tools are being supplemented with computational materials models at the atomistic level. The scales of interest range from computational chemistry (Density Functional Theory [DFT]), to atomistic simulations (Molecular Dynamics [MD]), to the microscale with the Porous Microstructure Analysis (PuMA) software that was recently awarded the 2022 NASA Software of the Year award. Finally, thermo-structural modeling is also of interest to the TPS Materials branch and done using commercial tools such as MSC MARC, MENTAT, NASTRAN and PATRAN. The present talk will also link the use of these computational tools to current NASA missions and projects associated with challenging and complex vehicles entries/reentries.

materials modeling↗
NASA NTRS2010

Microscopic Scale Simulation of the Ablation of Fibrous Materials

Ablation by oxidation of carbon-fiber preforms impregnated in carbonized phenolic matrix is modeled at microscopic scale. Direct numerical simulations show that the carbonized phenolic matrix ablates in volume leaving the carbon fibers exposed. This is due to the fact that the reactivity of carbonized phenolic is higher than the reactivity of carbon fibers. After the matrix is depleted, the fibers ablate showing progressive reduction of their diameter. The overall material recession occurs when the fibers are consumed. Two materials with the same carbon-fiber preform, density and chemical composition, but with different matrix distributions are studied. These studies show that at moderate temperatures (< 1000 K) the microstructure of the material influences its recession rate; a fact that is not captured by current models that are based on chemical composition only. Surprisingly, the response of these impregnated-fiber materials is weakly dependent on the microstructure at very high temperatures (e.g., Stardust peak heating conditions: 3360K).

Lachaud, Jean Romain↗
DOE OSTI2023

Mechanisms of Zn removal from water by amorphous geopolymer: Molecular-level insights from X-ray absorption spectroscopy, isotope fractionation, and surface complexation modeling

Porous geopolymers have attracted widespread attention as promising heavy metal adsorbents that can be synthesized from aluminosilicate solid wastes. However, the precise microstructural evidence for adsorbed heavy metals on geopolymers remains unclear due to the insensitivity of conventional characterization techniques on minerals with amorphous structure and surface disorder. Batch adsorption and column experiments coupled with X-ray absorption spectroscopy (XAS), Zn stable isotope, and surface complexation model (SCM) were employed to reveal the Zn removal mechanisms with coal fly ash porous geopolymer (CFAPG) at a molecular scale. The macroscopic kinetic and isothermal adsorption of Zn on CFAPG were well described by the pseudo-second-order model and Bi_Langmuir equation, respectively, indicating the presence of abundant heterogeneous active sites on the CFAPG surface. Further, two types of active sites on the CFAPG surface were identified by XAS coupled with Zn isotopes in batch experiments at pH <= 6.0: one is pH-dependent and associated with tetrahedral zinc coordi-nation, and the other is pH-insensitive and associated with octahedral zinc coordination; these sites were confirmed by the bidentate SCM as the variable charge site (surface complexation, >S-OH) and the permanent negative charge site (cation exchange, >X - ), respectively. Furthermore, the important contribution of surface co -precipitation besides surface complexation and cation exchange to the Zn adsorption on CFAPG was identified by XAS coupled with SCM in a flow-through column experiment at pH >6.0. These investigations provide a systemic understanding of the Zn adsorption mechanisms on CFAPG and an SCM reference for the application and prediction of geopolymers in heavy metal-contaminated water remediations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1991

Property and microstructural nonuniformity in the yttrium-barium-copper-oxide superconductor determined from electrical, magnetic, and ultrasonic measurements

The purpose of this dissertation was the following: (1) to characterize the effect of pore fraction on a comprehensive set of electrical and magnetic properties for the yttrium-barium-copper-oxide (YBCO) high temperature ceramic superconductor; and (2) to determine the viability of using a room-temperature, nondestructive characterization method to aid in the prediction of superconducting (cryogenic) properties. The latter involved correlating ultrasonic velocity measurements at room temperature with property-affecting pore fraction and oxygen content variations. The use of ultrasonic velocity for estimating pore fraction in YBCO is presented, and other polycrystalline materials are reviewed, modeled, and statistically analyzed. This provides the basis for using ultrasonic velocity to interrogate microstructure. The effect of pore fraction (0.10-0.25) on superconductor properties of YBCO samples was characterized. Spatial (within-sample) variations in microstructure and superconductor properties were investigated, and the effect of oxygen content on elastic behavior was examined. Experimental methods used included a.c. susceptibility, electrical, and ultrasonic velocity measurements. Superconductor properties measured included transition temperature, magnetic transition width, transport and magnetic critical current density, magnetic shielding, a.c. loss, and sharpness of the voltage-current characteristics. An ultrasonic velocity image constructed from measurements at 1mm increments across a YBCO sample revealed microstructural variations that correlated with variations in magnetic shielding and a.c. loss behavior. Destructive examination using quantitative image analysis revealed pore fraction to be the varying microstructural feature.

Roth, Don J.↗
NASA NTRS2014

Achieving ICME with Multiscale Modeling: The Effects of Constituent Properties and Processing on the Performance of Laminated Polymer Matrix Composite Structures

Integrated computational materials engineering (ICME) is a useful approach for tailoring the performance of a material. For fiber-reinforced composites, not only do the properties of the constituents of the composite affect the performance, but so does the architecture (or microstructure) of the constituents. The generalized method of cells is demonstrated to be a viable micromechanics tool for determining the effects of the microstructure on the performance of laminates. The micromechanics is used to predict the inputs for a macroscale model for a variety of different fiber volume fractions, and fiber architectures. Using this technique, the material performance can be tailored for specific applications by judicious selection of constituents, volume fraction, and architectural arrangement given a particular manufacturing scenario

computational mechanics↗
DOE OSTI2024

From Fundamental Interfacial Reaction Kinetics to Macroscopic Current–Voltage Characteristics: Case Study of Solid Acid Fuel Cell Limitations and Possibilities

The unique properties of solid acid electrolytes, in particular CsH 2 PO 4 , are in many ways ideal for fuel cell operation. However, the technology is constrained by high cathode overpotentials. Here a simplified cathode geometry is employed to obtain the fundamental electrochemical parameters (exchange current density and charge transfer coefficient) describing the oxygen reduction reaction (ORR) at the CsH 2 PO 4 -Pt-gas interface. The parameters are incorporated into a 1D model of the voltage–current characteristics of realistic SAFC cathodes, which reproduced the measured polarization behavior of such cathodes without recourse to fitting adjustable parameters. Following this validation, the model is utilized to evaluate the impact of changes to cathode properties, microstructure, and operating conditions. Of these, the charge transfer coefficient, measured to have a value of ≈0.6 for ORR on Pt in the SAFC cathode environment, is found to have the greatest impact on power output. Nevertheless, even without material modifications, a combination of microstructural and operational modifications are identified with projected performance metrics meeting Department of Energy targets (0.8 V at 300 mA cm –2 , and peak power density of 1 W cm –2 ), albeit at high Pt loadings. However, the analysis indicates that truly meaningful advances will likely necessitate the discovery of alternative ORR catalysts.

36 MATERIALS SCIENCE↗
DOE OSTI2024

HTO and selenate diffusion through compacted Na-, Na–Ca-, and Ca-montmorillonite

Radionuclide transport in smectite clay barrier systems used for nuclear waste disposal is controlled by diffusion, with adsorption significantly retarding transport rates. While a relatively minor component of spent nuclear fuel, 79 Se is a major driver of the safety case for spent fuel disposal due to its long half-life (3.3×10 5 yr) and its low adsorption to clay (K D < 10 L/kg), thus a thorough understanding of Se diffusion through clay is critical for understanding the long-term safety of spent fuel disposal systems. Through-diffusion experiments with tritiated water (HTO, conservative tracer) and Se(VI) were conducted with a well-characterized, purified montmorillonite source clay (SWy-2) under a constant ionic strength (0.1 M) and three different electrolyte compositions: Na + , Ca 2+ , and a Na + -Ca 2+ mixture at pH 6.5 in order to probe the effects of electrolyte composition and interlayer cation composition on clay microstructure, Se(VI) aqueous speciation, and ultimately diffusion. Further, the results were modeled using a reactive transport modeling approach to determine values of porosity (ε), D e (effective diffusion coefficient), and K D (distribution coefficient for adsorption). HTO diffusive flux was higher in Ca-montmorillonite (D e =1.68×10 -10 m 2 s -1 ) compared to Na-montmorillonite (De=7.83×10 -11 m 2 s -1 ). This increase in flux is likely due to a greater degree of clay layer stacking in the presence of Ca 2+ compared to Na + , which leads to larger inter-particle pores. Overall, the Se(VI) flux was much lower than the HTO flux due to anion exclusion, with Se(VI) flux following the order Ca (D e = 1.03×10 -11 m 2 s -1 ) > Na–Ca (D e = 2.12×10 -12 m 2 s -1 ) > Na (D e = 1.28×10 -12 m 2 s -1 ). These differences in Se(VI) flux are due to a combination of factors, including (1) larger accessible porosity in Ca-montmorillonite due to clay layer stacking and smaller electrostatic effects compared to Na-montmorillonite, (2) larger accessible porosity for neutral-charge CaSeO4 species which makes up 32% of aqueous Se(VI) in the pure Ca system, and (3) possibly higher Se(VI) adsorption for Ca-montmorillonite. Through a combination of experimental and modeling work, this study highlights the compounding effects that electrolyte and counterion compositions can have on radionuclide transport through clay. Diffusion models that neglect these effects are not transferable from laboratory experimental conditions to in situ repository conditions.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Effect of non-solvent from the phase inversion method on the morphology and performance of the anode supported microtubular solid oxide fuel cells

The microstructure of the anode in anode-supported solid oxide fuel cells has significant influence on the cell performance. Here, microtubular Ni-yttria stabilized zircona (Zr 0.8 Y 0.2 O 2 , YSZ) anode support has been prepared by the phase inversion method. Different compositions of non-solvent have been used for the fabrication of the Ni-YSZ anode support, and the correlation between non-solvent composition and characteristics of the microstructure of the anode support has been investigated. The presence of ethanol or isopropanol in the non-solvent can inhibit the growth of the finger-like pores in the anode support. With the increase of the concentration of ethanol or isopropanol in the non-solvent, a thin dense layer can be observed on the top of the prepared tubular anode support. In addition, the mechanism of pore formation is explained based on the inter-diffusivity between the solvent and the non-solvent. Lastly, the prepared microtubular solid oxide fuel cells (MT-SOFCs) have been tested, and the influence of the anode microstructure on the cell electrochemical performance is analyzed based on a polarization model.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Correlating dislocation structures to basal and prismatic slip bands near grain boundaries in tensile-tested Mg–4Al using a multiscale electron microscopy approach

Dislocation structures where basal and prismatic slip bands meet grain boundaries in tensile-tested Mg with 4 wt% Al (Mg–4Al) were studied using a multiscale electron microscopy approach to explore the assumptions made for the dislocation pileup theory: a) slip bands consist of dislocation arrays piling up near grain boundaries, and b) stress concentration in the adjacent grain is solely caused by dislocation pileup. After post-testing electron backscatter diffraction (EBSD) scans, focused-ion-beam (FIB) lift-out specimens were prepared from regions of interest so that the specimen plane is parallel to the bulk sample surface to allow for direct correlation between mesoscale and microscale characterization. High dislocation density and the dislocation arrays in the slip bands corroborated with the first assumption. Evidence of plastic deformation in all grains, however, showed that the second assumption was false. Additionally, dislocation structures including crystal rotation boundaries and dislocations dissociating to the basal plane provide additional stress relief mechanisms near grain boundaries that are not accounted for in many crystal plasticity models, which requires auxiliary models to more accurately predict the stress and strain distribution in a microstructure.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Finite-Temperature Structures of Supported Subnanometer Catalysts Inferred via Statistical Learning and Genetic Algorithm-Based Optimization

Single-atom catalysts (SACs) minimize noble metal utilization and can alter the activity and selectivity of supported metal nanoparticles. However, the morphology of active centers, including single atoms and subnanometer clusters of a few atoms, remains elusive due to experimental challenges. The computational cost to describe numerous cluster shapes and sizes makes direct first-principles calculations impractical. We present a computational framework to enable structure determination for single-atom and subnanometer cluster catalysts. As a case study, we obtained the low energy structures of Pd n (n = 1-21) clusters supported on CeO 2 (111), which are critical components of automobile three-way catalysts. Trained on density functional theory data, a three-dimensional cluster expansion is established using statistical learning to describe the Hamiltonian and predict energies of supported Pdn clusters of any structure. Low energy stable and metastable structures are identified using a Metropolis Monte Carlo-based genetic algorithm in the canonical ensemble at 300 K. We observe that supported single atoms sinter to form bilayer clusters and large cluster isomers share similarities in both shape and energy, and elucidate the significance of the support and microstructure on cluster stability. We discovered a simple surrogate structure-energy model, where the energy per atom scales with the square root of the average first coordination number, which can be used to estimate energies and compare the stability of clusters. Our framework, applicable to any metal/support system, fills an important methodological gap to predict the stability of supported metal catalysts in the subnanometer regime.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Novel method to control explosive shock sensitivity: A mesoscale study to understand the effect of thermally expandable microsphere (TEM) inclusions in high explosives (HE) microstructure

When the void content and/or void structure of a high explosive (HE) is altered by some means (i.e., bulk heating or mechanical damage), the shock initiation behavior of the material changes. The ability to precisely predict the change in shock sensitivity after an HE has undergone microstructural changes is a crucial capability in multi-scale reactive flow models. Here, we utilize thermally expandable microspheres (TEMs) as a dopant in a polymer bonded explosive (PBX) matrix to alter the shock initiation properties in a controlled fashion. Here, using a mesoscale modeling approach, we evaluated how a single TEM (before and after thermal expansion) behaves under shock compression, as well as how the matrix PBX in the direct vicinity of the TEM is affected. We first examined the effect of an unexpanded TEM in the explosive matrix and found that its presence does not significantly perturb the bulk flow and by extension will not affect bulk sensitivity. Next, we examined the effect of an expanded TEM and found that its presence significantly perturbs the flow via hydrodynamic jetting, which causes a secondary shock wave with a strength that exceeds that of the incident wave. Finally, we showed that this secondary shock interacts with the downstream porosity to ignite a larger fraction of the overall pore volume, commensurate with the secondary shock strength and the affected volume, increasing the global (bulk) shock sensitivity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Direct Measurements of Ionic Transport Behavior of Dual-Layer Porous Electrodes

To improve power and cycling performance of lithium-ion batteries, dual-layer or porosity-gradient electrodes have been proposed. By engineering a higher porosity close to the separator, the intention is to improve ion transport where it is most impactful. In this research, MacMullin numbers of two dual-layer anode samples are tested using an impedance measurement technique developed previously. To characterize the microstructure of each layer independently, we developed an improved transmission-line model that accounts for each layer's properties.Virtual experiments using COMSOL Multiphysics to simulate impedance measurements are used to examine and improve the accuracy of the numerical inversion procedure. The results for the two dual-layer anodes studied show that MacMullin numbers follow expected trends, though the anodes are quite different from each other.

Dual-Layer Electrodes↗
DOE OSTI2022

Investigating the Influence of Ni, ZrO 2 , and Y 2 O 3 from SOFC Anodes on Siloxane Deposition

Siloxanes, as a type of impurity in biogas, can poison the Ni-YSZ anode of SOFCs. However, the influence of individual components of the anode, such as Ni, ZrO 2 , and Y 2 O 3 , on the siloxane deposition process has not been investigated extensively. In this study, Ni, ZrO 2 , and Y 2 O 3 pellets were exposed to H 2 + N 2 + H 2 O + D4 (octamethylcyclotetrasiloxane, 2.5 ppmv) and H 2 + N 2 + D4 (2.5 ppmv siloxane) gas mixtures at 750 °C to investigate their affinity and tolerance for siloxane degradation. Surface morphology analysis and electrochemical analysis including electrochemical impedance spectroscopy (EIS), related distribution of relaxation times (DRT) analysis and equivalent circuit modeling with complex nonlinear least square (CNLS) fitting were conducted. Here, a microstructure parameter—tortuosity factor to porosity ratio $\tau /\varepsilon $ calculated by diffusion polarization resistance was utilized for siloxane deposition evaluation. After comparing pellets surface morphology changes before and after experiments and $\tau /\varepsilon $ change following the contamination, Ni is considered as a major factor in siloxane deposition reactions in Ni-YSZ anode.

25 ENERGY STORAGE↗
DOE OSTI2022

MEUMAPPS (C++ Version)

Many materials, metal alloys in particular, have features on the on micrometer or nanometer scale that have a large impact on the properties of the material. These features are known as the microstructure of the material. Understanding why and how the microstructure forms in a material is of fundamental scientific interest as well as of significant technological interest. The capability to predict microstructure evolution in a material allows the intentional design of microstructures and hence the intentional design of material properties. The phase-field method is one of the leading methods for predicting microstructure evolution. One of the most significant problems for phase-field models is their computational expense. Even limited phase-field simulations can easily require thousands of CPU core-hours to complete, which significantly limits their use. This code provides both a general framework for creating scalable, GPU-accelerated phase-field model applications as well as several applications themselves. The code is capable of using hundreds of GPUs efficiently, which greatly reduces the time required to perform simulations. The code is written with an emphasis on performance portability, that is the ability for the code to run efficiently on a number of different computing architectures without modification of the source code. The performance portability of this code is primarily enabled through the use of two libraries, Kokkos (performance portable data structures and execution patterns) and heFFTe (performance portable distributed 3D fast Fourier transforms). The code consists of a core library, applications, and tests. The core library includes shared functionality between applications. This includes interfaces with fast Fourier transform (FFT) libraries such as heFFTe, data structures based on Kokkos, file input and output capabilities, and a solver for infinitesimal strain mechanical equilibrium problems. Five applications are included in the code. The flagship application is the MEUMAPPS-SS application, which implements the Kim-Kim-Suzuki phase-field model for precipitation for an arbitrary number of phases and components in a metal alloy. Five simpler applications are also included that solve the Eshelby inclusion problem, Allen-Cahn equation, the coupled Allen-Cahn and diffusion equations, and the Cahn-Hilliard equation. The code includes two applications to solve the Cahn-Hilliard equation, one with constant-step-size first-order time integration and the second with adaptive high-order time integration.

DeWitt, Stephen [Oak Ridge National Lab. (ORNL), O↗