Search NASA⌕ Search

SEARCH · Search NASA

Results for “CONJUNCTION”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Autoignition behavior of gasoline/ethanol blends at engine-relevant conditions

Ethanol is an attractive oxygenate increasingly used for blending with petroleum-derived gasoline yielding beneficial combustion and emissions behavior for a range of internal combustion engine schemes, including stoichiometric spark-ignition and low temperature combustion (LTC). As such, it is important to fundamentally understand the autoignition behavior of gasoline/ethanol blends. This work utilizes a rapid compression machine (RCM) and a homogeneous charge compression ignition (HCCI) engine to experimentally quantify changes in fuel reactivity, through ignition delay times and preliminary heat release, for blends of 0 to 30% vol./vol. into a full boiling range research gasoline (FACE-F). Diluted/stoichiometric and undiluted/fuel-lean conditions are explored covering a wide range of compressed temperatures and pressures relevant to conventional and advanced, gasoline combustion engines. Detailed chemical kinetic modeling is undertaken using a recently updated gasoline surrogate model in conjunction with a five-component surrogate to model the RCM experiments and provide chemical insight into the perturbative effects of ethanol on the autoignition process. The diluted/stoichiometric RCM measurements reveal that within the low-temperature regime ethanol retards first-stage and main ignition delay times, and suppresses both the rates and extents of low-temperature heat release (LTHR), while within the intermediate-temperature regime ethanol only causes slight changes. Good agreement of ignition delay time and preliminary heat release prediction is found between model and experimental results. Sensitivity and flux analyses further show that ethanol blending effects are dominated by the competition between the H-atom abstraction from ethanol and other fuel components by OH radical at low temperatures and by HO 2 radical at intermediate temperatures. These findings are consistent across both fuel loading conditions explored in this study. In addition, when HCCI engine experiments are mapped onto undiluted/lean RCM measurements under a constant combustion phasing scenario, good correspondence between the two apparatuses is observed for LTHR and start of high-temperature heat release. Finally, the current study highlights the importance of characterizing LTHR in predicting fuel behaviors in high-boost/low-temperature engines, and demonstrates that RCM experiments can provide an alternative, and more-efficient avenue for such characterization.

02 PETROLEUM↗
DOE OSTI2021

Probing intermediate temperature heat release in autoignition of C3-C4 iso-alcohol/gasoline blends

In this work, we report an experimental and modeling study on the blending effects of C3–C4 iso-alcohols, namely iso-propanol and iso-butanol, on the characteristics of intermediate temperature heat release (ITHR) for a research-grade gasoline (FACE-F) in a rapid compression machine operating at diluted/stoichiometric fuel loading, compressed pressure of 43 bar, and compressed temperatures from 700 to 970 K. Changes in ITHR behavior are characterized through ITHR extent and evolution at 0 to 30 vol% iso-alcohol blending levels. Experimental observations reveal the strong promoting and decelerating effects on ITHR extent and evolution, respectively, within the low-temperature regime for both iso-alcohols, with stronger effects seen for iso-propanol; within the intermediate-temperature regime, the influence of both iso-alcohols diminishes greatly. Comprehensive chemical kinetic modeling is undertaken using a recently developed gasoline surrogate/alcohol model in conjunction with a five-component gasoline surrogate (FGF-LLNL), with good agreement obtained with the experiments for all the blends. Sensitivity, rate of production and flux analyses highlight the importance of chemical kinetic interactions via both fuel-specific and non-fuel-specific reactions. Disabling all the chemical kinetic interactions between iso-alcohol and FGF-LLNL in carbonated/non-fuel-specific species leads to somewhat extended ITHR duration and reduced ignition reactivity. Direct comparison between FGF-LLNL/iso-propanol and FGF-LLNL/iso-butanol blends is also made, where iso-propanol and iso-butanol are found to influence ITHR characteristics via significantly different chemical kinetic interactions.

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

O 2 -Dependence of reactions of 1,2-dimethoxyethanyl and 1,2-dimethoxyethanylperoxy isomers

Reaction mechanisms of Ṙ and ROȮ radicals derived from low-temperature oxidation of 1,2-dimethoxyethane (CH 3 O(CH 2 ) 2 OCH 3 ) were investigated using speciation from multiplexed photoionization mass spectrometry (MPIMS) measurements via Cl-initiated oxidation, in conjunction with electronic structure calculations. The experiments were conducted at 5 bar, from 450 K – 650 K, and O 2 concentrations from 1 · 10 14 cm –3 – 6 · 10 18 cm –3 to probe the effects on competing reaction channels of 1,2-dimethoxyethanyl (Ṙ) and 1,2-dimethoxyethanylperoxy (ROȮ) isomers. Several species were detected with photoionization spectral fitting – ethene, formaldehyde, methyl vinyl ether, and 2-methoxyacetaldehyde – and, as determined by electronic structure calculations, may form via unimolecular decomposition of 1,2-dimethoxyethanyl or 1,2-dimethoxyethanylperoxy. O 2 -dependent yield ratios show that the formation pathways for all species undergo a competition between O 2 -addition and unimolecular decomposition. Here, adiabatic ionization energies were also calculated and utilized along with exact mass determinations to infer contributions for other species derived exclusively from first- and second-O 2 -addition, including 1,2-dimethoxyethene, cyclic ethers, and dicarbonyls.

1,2-dimethoxyethane↗
DOE OSTI2023

Simulation toolkit for digital material characterization of large image-based microstructures

In this paper, an efficient image-based simulation toolkit for material characterization is presented, which is scalable to work from personal computers to workstations. The effective thermal conductivity, elasticity, and permeability are evaluated employing a computational homogenization framework based on the Finite Element Method (FEM). Two complementary open-source packages are presented: one developed in Python, which can convert digital images into voxel meshes (pyTomoviewer); the other developed in Julia, that can run numerical simulations to compute effective material properties (chpack). Also, a CUDA C version of chpack is provided (chfem_gpu). They were designed to deal with large multi-phase models, so strategies were devised to minimize their memory footprint, while avoiding a high toll on execution time. The voxel-based approach significantly simplifies the FEM meshes and allows efficient matrix-free implementations. In that sense, to handle large linear systems of equations, the element-by-element (EBE) technique is adopted, in conjunction with a low-memory implementation of the Preconditioned Conjugate Gradient (PCG) method. Finally, the code was thoroughly tested on an artificial geometry made of a square array of cylinders, for which analytical solutions exist, as well as on a real micro-tomographic reconstruction of FiberForm TM , a carbon preform commonly used in thermal protection systems.

36 MATERIALS SCIENCE↗
DOE OSTI2020

An Eulerian multimaterial framework for simulating high-explosive aquarium tests

Aquarium tests of cylindrical high-explosive charges provide optical data of the detonation front velocity and shape, propagation of the shock wave in the surrounding water, and expansion rates of the detonation products behind the front. Data from aquarium experiments is often used for calibration of reactive burn models based on phenomenological equations of state (EOS) and reaction rate laws. This paper presents a multimaterial numerical modeling framework to solve the 2D axisymmetric reactive Euler equations for high-explosive aquarium tests, in particular for ammonium nitrate - fuel oil (ANFO) explosives. An extension of the Ghost Fluid Method (GFM) is used to handle the dynamic material interfaces for the ANFO explosion products, the charge-confining material (polymethyl methacrylate PMMA), and the surrounding water. This study analyzes the sensitivity of calculations (both computational efficiency and numerical accuracy) to different algorithms for the material interface models including the original GFM versus Riemann solver-based strategies. A novel method for defining the left and right states in the interfacial Riemann problem eliminates the need for sorting or nodal interpolation during the projection along the material interface. Numerical tests indicate that populating the interface node values using the Riemann solution mitigate the overheating error observed in steady-state calculations. Solution convergence and computational efficiency are explored as a function of the spatial and temporal order of the schemes. Results from the computational model with analytical equations of state and fitted reaction rate parameters show very good quantitative agreement with experimentally observed detonation front velocity, reaction products expansion, and shock wave propagation in the surrounding water for a cylindrical ANFO charge encased in PMMA. Finally, the proposed modeling framework, in conjunction with experimental tests, provides a reliable tool to assess equations of state and reaction rate expressions for reactive burn models of confined high explosives.

42 ENGINEERING↗
DOE OSTI2025

Influence of gradation on shear behavior and pore water response of sand

Here, this paper examines the influence of particle size distribution (PSD) on the shearing behavior of sand using axisymmetric triaxial compression (ATC) testing in conjunction with 3D imaging using synchrotron micro-computed tomography (SMT). Angular sand was sieved to obtain a narrow gradation (uniform) and a wide gradation (non-uniform) with uniformity coefficients (C u ) of 1.2 and 3.0, respectively. Experiments were conducted under drained conditions and a constant effective cell pressure (σ$^{'}_3$) of 50 kPa at low and high back pressures (BP). Results of paired ATC-SMT were analyzed in terms of global constitutive behavior, strain localization, and localized flow characteristics. A new particle tracking algorithm is developed to expose strain localization. Relative particle translation gradient (RPTG) results revealed that gradation affects how early the shear band develops and a tendency for narrow gradation to exhibit varying degrees of diffused shear band at high axial strain (ε 1 ) even though a clear single shear band was observed at the peak strength. Geometric tortuosity (τ) is used as a quantitative parameter to assess flow changes and the efficiency of the flow path. Representative elementary volumes (REVs) inside and outside the shear band were analyzed and the results showed a variation in flow path characteristics where wide gradation specimens had preferential flow paths.

Computed tomography↗
DOE OSTI2020

Innovative large overbraided carbon fiber composite beams in historical Friendship Bell Pavilion

The International Friendship Bell Pavilion located in the City of Oak Ridge, Tennessee symbolizes the spirit of the peace and friendship shared by the city of Oak Ridge and city of Naka, Japan. The Bell Pavilion features an 8000-pound bronze cast bell along with futuristic materials and advanced manufacturing technologies. The team has developed innovative overbraided carbon fiber arched beams in conjunction with the bronze bell. Carbon fiber was utilized in these beams because of its lightweight, strong, stiff, and durable characteristics. A sandwich design used a high-performance foam core, overbraided with triaxial carbon fibers. This is the world’s largest overbraided composite structure for parks and recreation. For the first time in such an application, 17 free floating, 33-foot (10 m) long arched carbon fiber composite beams cantilever from a central concrete-steel superstructure over the bronze bell. The arc radius of the bottom of each beam is uniquely different giving an appearance of an open dome when standing underneath the structure.

Parks and Recreation↗
DOE OSTI2025

History, present status, and future directions of vanadium alloys for fusion reactors

Since the need for low-activation fusion reactor materials was recognized in the mid-1980s, vanadium alloys have been researched and developed as promising candidates for blanket structural materials. Vanadium alloys are non-ferromagnetic and ductile materials and thus have advantages different from other candidates. However, since vanadium alloys are still at a premature industrial stage, research and development have been carried out not only on issues specific to fusion reactors but also on those related to industrial materials. Blankets using vanadium alloys as structural materials and liquid Li as a tritium breeder and coolant (self-cooled V/Li blankets) have unique characteristics, and research and development of issues specific to V/Li blankets have been conducted in conjunction with the development of vanadium alloys. In this paper, the progress of this research and development is reviewed, and the remaining issues are pointed out. As a result, the future direction of vanadium alloy research and development is discussed, considering recent changes in the environment surrounding fusion energy.

Liquid Li blanket↗
DOE OSTI2021

A multi-center quadrature scheme for the molecular continuum

We report a common way to evaluate electronic integrals for polyatomic molecules is to use Becke’s partitioning scheme (Becke and Chem, 1988) in conjunction with overlapping grids centered at each atomic site. The Becke scheme was designed for integrands that fall off rapidly at large distances, such as those approximating bound electronic states. When applied to states in the electronic continuum, however, Becke scheme exhibits slow convergence and it is highly redundant. Here, we present a modified version of Becke scheme that is applicable to functions of the electronic continuum, such as those involved in molecular photoionization and electron–molecule scattering, and which ensures convergence and efficiency comparable to those realized in the calculation of bound states. In this modified scheme, the atomic weights already present in Becke’s partition are smoothly switched off within a range of few bond lengths from their respective nuclei, and complemented by an asymptotically unitary weight. The atomic integrals are evaluated on small spherical grids, centered on each atom, with size commensurate to the support of the corresponding atomic weight. The residual integral of the interstitial and long-range region is evaluated with a central master grid. The accuracy of the method is demonstrated by evaluating integrals involving integrands containing Gaussian Type Orbitals and Yukawa potentials, on the atomic sites, as well as spherical Bessel functions centered on the master grid. These functions are representative of those encountered in realistic electron-scattering and photoionization calculations in polyatomic molecules.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

The hazard components of representative key risks. The physical climate perspective

The framework of Representative Key Risks (RKRs) has been adopted by the Intergovernmental Panel on Climate Change Working Group II in order to categorize, assess and communicate a wide range of regional and sectoral key risks, expected to become severe as a consequence of the potentially detrimental convergence of changing climate conditions, exposure and vulnerability of human and natural systems, and responses through (mal)adaptation measures. Other papers in this special issue treat each of eight RKRs by holistically assessing their current status and future evolution as a result of this convergence. However, their treatment cannot always organize such assessment according to a systematic gradation of climatic changes, often times having to extrapolate the big-picture evolution of risk from either qualitative effects of “low” medium” and “high” warming, or limited/focused analysis of the consequences of particular mitigation choices (e.g., benefits of limiting warming to 1.5 or 2C). In this study we offer as systematic a representation as possible – given current literature and assessments – of the future evolution of the hazard components of RKRs. We draw on the assessment of changes in climatic impact-drivers relevant to RKRs described in the 6th Assessment Report by Working Group 1 (WGI) supplemented when needed by more recent literature. We identify the relevant hazards for each RKR, based upon the WGII authors’ assessment, and we report on their current state and expected future changes in magnitude, intensity and/or frequency, linking these changes to Global Warming Levels to the extent possible. For some of these quantities - like regional trends in oceanic and atmospheric temperature and precipitation, some heat and precipitation extremes, permafrost thaw and Northern Hemisphere snow cover - a strong and quantitative relationship with global average temperature change has been assessed. For others - like frequency and intensity of tropical cyclones and extra-tropical storms, and fire weather - the link to increasing global temperatures can only be described qualitatively. For some processes - like the behavior of ice-sheets, or changes in circulation dynamics - large uncertainties about the effects of different levels of warming remain, and for a few others - like ocean pH and air pollution - the composition of the scenario of anthropogenic emissions is most relevant, rather than the warming reached. In almost all cases, however, the basic message remains that every small increment in CO2 concentration in the atmosphere and associated warming will bring changes in climate phenomena that will contribute to increasing risk of impacts on human and natural systems when these changes happen in conjunction with these systems’ conditions of exposure and vulnerability, and in lack of effective adaptation. Our picture of the evolution of RKR-relevant climatic impact-drivers complements and enriches the treatment of RKRs in the other papers in at least two ways: by filling in an often only cursory or limited representation of the physical climate aspects driving impacts, and by providing a fuller representation of their future possible evolution which represents an important component – if never the only one – of the potential future evolution of risk severity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Energy impact of heating electrification in mid-rise multifamily buildings in mixed-humid climates

Decarbonizing the electric grid in conjunction with electrifying residential heating is a critical step to combat climate change. Heating in multifamily buildings with the existing natural gas-fired central boiler is a complex process that not only leads to overheating in some apartment units but also results in energy waste and high gas bills. In this study, we consider a multifamily building in New York City, USA, to evaluate the performance of five different heating systems, which represent a step-by-step transition from the conventional to a fully electrified heating system, and determine their impact on the site energy consumption and source CO 2 emissions. Results indicate that overheating in a multifamily building can raise the indoor temperature by as much as 8°C above comfortable limits. Transitioning from conventional steam radiators to cold climate heat pumps can reduce annual site heating energy by up to 70% and source CO 2 emissions by up to 21%.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

Operation of helium sub-atmospheric multistage cryogenic centrifugal compressor trains: Part 2 – Transient modeling and pump-down path selection

Low-pressure conditions required for operation of helium cryogenic systems below the normal boiling point (i.e. 4.2 K) are established through a transient process, commonly referred to as ‘pump-down’. This process is defined as the transition from pressures above atmospheric conditions to the saturation pressure which corresponds to a specified operational temperature. The FRIB 2 K system consists of five cryogenic centrifugal compressors which are operated in series. Historically, the pump-down process path has been established through empirical methods and system operator experience. Investigation into the pump-down process at FRIB aimed to develop a pump-down methodology which relies on theoretical model predictions rather than empirically developed process paths. Ensuring stable operation during the pump-down process involved application of a centrifugal compressor performance prediction model, which is described in Part 1 of this paper. Compressor performance maps can be directly used to evaluate the stability of a selected pump-down path and anticipate the overall reliability of the selected path. In conjunction with the compressor performance maps, a system pressure model was developed to estimate the transient pressure response during the pump-down process. Lastly, an explicit equation was developed to establish a mass flow rate profile for the pump-down process. Implementation of the presented methodology (including the developed models) allows for the system operator to determine a continuous pump-down path which maintains compressor stability while conforming to overall system capabilities. Altogether, the methodology presented has resulted in simplification of transient pump-down operations and increased the reliability, stability and efficiency of the pump-down process.

Compressor train control↗
DOE OSTI2023

Synergy between cosmological and laboratory searches in neutrino physics

The intersection of the cosmic and neutrino frontiers is a rich field where much discovery space still remains. Neutrinos play a pivotal role in the hot big bang cosmology, influencing the dynamics of the universe over numerous decades in cosmological history. Recent studies have made tremendous progress in understanding some properties of cosmological neutrinos, primarily their energy density. Upcoming cosmological probes will measure the energy density of relativistic particles with higher precision, but could also start probing other properties of the neutrino spectra. When convolved with results from terrestrial experiments, cosmology can become even more acute at probing new physics related to neutrinos or even Beyond the Standard Model (BSM). Any discordance between laboratory and cosmological data sets may reveal new BSM physics and/or suggest alternative models of cosmology. Here we give examples of the intersection between terrestrial and cosmological probes in the neutrino sector, and briefly discuss the possibilities of what different laboratory experiments may see in conjunction with cosmological observatories.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Digitalization of an experimental electrochemical reactor via the smart manufacturing innovation platform

The exponential increase in data produced over the last two decades has revolutionized the way we collect, store, process, analyze, model, and interpret information to improve profitability. Manufacturing is no exception. How- ever, Smart Manufacturing, the digital practice, organization, workforce, and infrastructure transformation for collection and deployment of data and models at scale and at all levels of manufacturing, is a complex, costly, and labor-intensive journey that is still seeing slow adoption. The Clean Energy Smart Manufacturing Innovation Institute (CESMII), a national Manufacturing USA public-private partnership sponsored by the Department of Energy, is addressing this scaled use of data and modeling in manufacturing. CESMII has focused on how to col- lect and use operating data for numerous applications that improve productivity, precision, and performance of manufacturing operations from factory floor to supply chain using process simulation, predictive analytics, mon- itoring and control, and real-time optimization. Because contextualized data are key, CESMII has developed the Smart Manufacturing Innovation Platform (SMIP) to lower the barriers to the data that are needed to accelerate data-based model building, improve data visualization, and more quickly gain insights. Reusable, standards-based ways of doing data collection, ingestion, and contextualization are particularly important for scaling access and use of data. The SMIP uses a standards-based definition and construct for reusable information models called an SM Profile. When an SM Profile is used in conjunction with the SMIP, the SMIP ensures the availability of contextualized, operational data for model building. The present work demonstrates Smart Manufacturing and the application of the SMIP for building several data-centered models for the operation and control of an ex- perimental electrochemical reactor that reduces carbon dioxide (CO 2 ) gas to valuable liquid and gas chemicals, such as alcohols, olefins, and syngas. We describe how the SMIP plays a central role in more effective model building and we demonstrate how the electochemical reactor can be controlled and optimized for the desired products. Use of the SMIP involves the transmission of real-time sensor measurements to a cloud resource so that the operating data are available to all model building experts. The data collection and transmission process is fully automated to greatly reduce the need for manual manipulation of the data. Data-driven machine learning models are used for advanced real-time state estimation, real-time optimization, and model-based feedback control for the reactor. The application models are implemented as a system to monitor the data flow and control the electrochemical reactor with a single visualization interface. SM Profiles are used to demonstrate reusability of the information models for the reactor and the instrumentation. The application packages, algorithms, and user interfaces developed are cast as Docker images in a library to facilitate reusability of the application models.

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

Atomic structure and grain shape evolution of nanodiamond during annealing in oxidizing atmosphere from neutron diffraction and MD simulations

In this work, real space pair-distribution function (PDF) analysis in conjunction with molecular dynamics (MD) simulations was applied to characterize crystal structure, shape and surface structure of diamond nanoparticles annealed in air up to 880 °C. MD simulated models of individual nanodiamond particles of various shapes and sizes were used to calculate theoretical radial distribution functions G(r) for powders with different grain size distributions. The calculated curves, showing distinct differences for the explored parameter space, were compared to the experimental G(r) functions obtained from neutron scattering for annealed detonated nanodiamond samples in a search for a model best describing the actual powders. It was found that as-synthesized nanodiamond grains are terminated by equally abundant (100), (110) and (111) crystallographic faces and their structure does not change upon annealing in air up to 280 °C. The shape and size of grains annealed above 480 °C changes due to surface etching. The etching rate differs between crystals faces and depends on temperature, resulting in the observed differences in grain shape. Below 680 °C the share of the (110) facets gradually increases and at 680 °C the grains become rhombic dodecahedra terminated solely by (110) surfaces. Further increase of the annealing temperature promotes formation of (111) surfaces and at 880 °C the preferred grain shape becomes an octahedron terminated by a variant of (111) surfaces possessing three dangling sp3 bonds per atom. At lower temperatures all grains become smaller and the size distribution moves towards the lower values. Furthermore, as the annealing temperature increases, the smallest grains disappear entirely and the average grain size grows back.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Towards high data-rate diffusive molecular communications: A review on performance enhancement strategies

Diffusive molecular communications (DiMC) have recently gained attention as a candidate for nano- to micro- and macro-scale communications due to its simplicity and energy efficiency. As signal propagation is solely enabled by Brownian motion mechanics, DiMC faces severe inter-symbol interference (ISI), which limits reliable and high data-rate communications. Herein, recent literature on DiMC performance enhancement strategies is surveyed; key research directions are identified. Signaling design and associated design constraints are presented. Here, studies on fundamental information theoretic limits of DiMC channel are reviewed. Classical and novel transceiver designs are discussed with an emphasis on methods for ISI mitigation and performance-complexity tradeoffs. Key parameter estimation strategies such as synchronization and channel estimation are considered in conjunction with asynchronous and timing error robust receiver methods. Finally, source and channel coding in the context of DiMC is presented.

42 ENGINEERING↗
DOE OSTI2023

Beneficial effect of copper on pitting resistance of Ni-Cr-Fe alloys

This study examines the effect of copper alloying on pitting resistance in a model solid solution FCC Ni-13%Cr-10%Fe alloy through potentiodynamic and potentiostatic polarization in 0.1 M NaCl in conjunction with an analysis using first-principles competitive electro-chemisorption modeling. The pitting potential increased with increasing Cu content in the alloy. Furthermore, the extent of metastable pit growth was suppressed and the incubation time for metastable to stable pit transition increased with Cu content. The first-principles competitive adsorption calculations suggested that Cu alloying suppresses chloride ion adsorption on the alloy surface in a simulated pit environment, which inhibits active dissolution at the pit bottom, enhances proton adsorption, and thereby increases the local pH at the pit bottom. In conclusion, we propose that these two effects of Cu in solid solution combine to reduce pit stability and may act in addition to the enrichment of Cu on the corroding pit surface.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A reaction–diffusion model for grayscale digital light processing 3D printing

We report that digital light processing (DLP) 3D printing is an additive manufacturing process that utilizes light patterns to photopolymerize a liquid resin into a solid. Due to the accuracy of modern digital micromirror devices (DMD) and recent advances in resin chemistry, it is now possible to create functionally graded structures using different light intensity values, also known as grayscale DLP (g-DLP). Different intensities of light lead to differences in the polymer crosslinking density after curing, which ultimately produces a part with gradients of material properties. However, g-DLP is a complicated process. First, the DLP printing is a highly coupled chemical and physical process that involves light propagation, chemical reactions, species diffusion, heat transfer, volume shrinkage, and changes in mechanical behaviors of the curing resin. Second, in g-DLP, light gradients create strong in plane gradients of chemical species concentrations in the curing liquid resin due to the strong dependence of light intensity on the rate of monomer crosslinking. Furthermore, light gradients through the depth create concentration gradients due to the degree of cure dependent light absorption and the use of photoabsorbers. These complex physical features of the printing process must be understood in order to properly control printing parameters such as light exposure time, printing speed, and grayscale variations to achieve accurate mechanical properties. In this paper, a photopolymerization reaction–diffusion model is developed and used in conjunction with experiments to investigate the coupled effects of light propagation, chemical reaction rates, and species diffusion during g-DLP 3D printing. The model is implemented numerically utilizing the finite difference method and simulation results are compared to experimental findings of simple printed structures. The agreement between experimental and model predictions of simple quantities of interest, such as geometric feature sizes, shows that the model can capture the overcure due to free-radical and other species diffusion during printing when grayscale patterns are employed. This model lays the groundwork for future extensions that can incorporate more complex coupled physics such as heat transfer, volume shrinkage, and material property evolution, which are critically important in utilizing g-DLP 3D printing for the fabrication of high-performance parts which excellent geometric and material property tolerances.

36 MATERIALS SCIENCE↗