Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computational Geometry”

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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Rarefied xenon flow in orificed hollow cathodes

A parametric study is conducted to quantify the effect of the keeper electrode geometry on the xenon neutral flow quantities within orificed hollow cathodes, prior to ignition. The keeper impinges directly on the flow out of the cathode orifice and its geometry influences the product between the pressure in the orifice–keeper region and the cathode-to-keeper distance. A representative cathode is simulated using the Direct Simulation Monte Carlo method. The numerical model is first validated with computational results from the literature. A parametric study is then conducted. Parameters include the cathode pressure–diameter in the range of 1–5 Torr cm and the following geometric ratios (and ranges): cathode orifice-to-inner radii (0.1–0.7), keeper orifice-to-cathode orifice radii (1–5), and keeper distance-to-cathode-orifice diameter (0.5–10). It is found that, if both keeper and cathode have identical orifice radii, the flow remains subsonic in the orifice-to-keeper region. In most cases, however, the flow becomes underexpanded and supersonic, and the static pressure within the orifice-to-keeper region is, on average, 4% that of the upstream pressure value. The orifice–keeper region pressure increases with either a decrease in the keeper orifice diameter or an increase in the distance between cathode and keeper, in agreement with literature data. Both trends are explained through conservation laws. A statistical study of numerical results reveals that the ratio of ignition-to-nominal mass flow rates has a most probable value of 50, which suggests that heaterless cathode ignition at a minimum DC voltage may be achieved by increasing the input mass flow rate by a factor of 50.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Calibrating uncertain parameters in melt pool simulations of additive manufacturing

Melt pool scale numerical modeling of additive manufacturing (AM) processes can provide predictive capabilities and theoretical insight into the process-property-structure-performance relationships for AM parts. Despite capabilities of numerical models to solve complex multi-physics problems, it is often important to consider a tradeoff between detailed physics and computational cost. Therefore, sources of uncertainty in both experimental conditions and the parameters needed for modeling require models to be validated against empirical evidence. Here, a method is proposed to calibrate uncertain parameters used in continuum-scale melt pool models for powder bed fusion (PBF) AM. Both a simplified heat transfer model and a heat transfer and fluid flow model were investigated. A surrogate model and Markov chain-based optimization algorithm calibrated melt pool geometry for models within experimental variation of the target melt pool width and depth from the NIST AM-Bench 2018-02 dataset. The melt pool temperature distributions, solidification parameters, and simulated multi-layer solidification microstructures were compared between the two models. Similar results from both models indicate that calibrated, lower fidelity numerical models may be used in place of higher fidelity models to generate melt pool solidification data. Finally, these calibrated models therefore enable lower computational cost melt pool simulations without a noticeable decrease in simulation accuracy for grain-scale microstructure simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

3D ion gyro-orbit heat load predictions for NSTX-U

Abstract High power tokamaks operate with divertor heat loads capable of destroying the plasma facing components (PFCs). High fidelity heat load predictions are necessary to ascertain the PFC state for design and during operation. Typical heat flux calculations are 2D, time invariant, and assume that power flows directly along the magnetic field lines (the optical approximation). These assumptions neglect the complex 3D geometries employed to protect the PFCs, the time varying nature of the plasma and PFC thermal state, and the helical trajectories of ions with finite Larmor radii (the gyro-orbit approximation). An integrated software framework, the heat flux engineering analysis toolkit (HEAT), was developed to generate time varying optical heat loads applied to real engineering computer aided design (CAD) (Looby et al 2022 Fusion Sci. Technol. 78 10–27). Recently, an ion-gyro orbit module has been added to HEAT. This module calculates the helical trajectories of ions as they gyrate about the magnetic field lines using kinetic theory macro-particles to accelerate the calculation. First, the new gyro-orbit module will be presented. Next, a comparison to existing research is performed. Finally, an analysis of the gyro-orbit heat loads for NSTX-U is presented for diverted discharges using the engineering CAD models utilized for PFC fabrication. Including these gyro-orbit effects can enhance the PFC performance by ‘smearing’ out the magnetic shadows associated with the castellated fish-scaled geometry. Simultaneously, the helical trajectories can degrade performance when they load narrow regions on edges and corners with high heat fluxes. Analysis of the trade-offs between these competing effects is included, and regions for further investigation are identified.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Characterizing the Geometry and Quantifying the Impact of Nanoscopic Electrocatalyst/Semiconductor Interfaces under Solar Water Splitting Conditions

The materials that are receiving the most attention in photoelectrochemical water splitting are metallic nanoparticle electrocatalysts (np-EC) attached to the surface of a semiconductor (SC) light absorber. In these multicomponent systems, the interface between the semiconductor and electrocatalysts critically affects performance. However, the np-EC/SC interface remains poorly understood as it is complex on atomic scales, dynamic under reaction conditions, and inaccessible to direct experimental probes. This contribution sheds light on how the electrocatalyst/semiconductor interface evolves under reaction conditions by investigating the behavior of nickel electrocatalysts (as nanoparticles and films) deposited on silicon semiconductors. Rigorous electrochemical experiments, interfacial atomistic characterization, and computational modeling are combined to demonstrate critical links between the atomistic features of the interface and the overall performance. It is shown that electrolyte-induced atomistic changes to the interface lead to (1) modulation of the charge carrier fluxes and a dramatic decrease in the electron/hole recombination rates and (2) a change in the barrier height of the interface. Additionally, the critical roles of nonidealities and electrocatalyst coverage due to interfacial geometry are explored. Each of these factors must be considered to optimize the design of metal/semiconductor interfaces which are broadly applicable to photoelectrocatalysis and photovoltaic research.

36 MATERIALS SCIENCE↗
DOE OSTI2025

OpenSn: A massively parallel, open-source simulation environment for discrete ordinates radiation transport

OpenSn is an open-source, massively parallel deterministic radiation transport code for solving the discrete-ordinates ( S N ) form of the Boltzmann transport equation on unstructured, arbitrary polyhedral meshes. It supports high-fidelity simulations involving steady-state, eigenvalue, and adjoint problems for neutral particles (e.g., neutrons, photons, multi-particles), using the multigroup approximation in energy. OpenSn combines angular discretization via discrete ordinates with a discontinuous Galerkin finite element method (DGFEM) in space, enabling accurate resolution of transport physics on arbitrary polyhedral cells, included locally refined spatial grids. It includes multiple angular quadrature types, including locally refined angular quadratures. Written in modern C++ with a Python API, OpenSn runs efficiently on platforms ranging from laptops to supercomputers. The transport sweep algorithm is implemented using a task-based, directed-acyclic-graph (DAG) approach for each angle and supports asynchronous parallelism across thousands of MPI ranks. Group-set aggregation improves compute intensity, and synthetic acceleration techniques (e.g., diffusion synthetic acceleration, second-moment method) enhance solver convergence. OpenSn has been verified on reactor physics problems and demonstrated excellent weak and strong scaling performance on more than 32,768 processes, making it a versatile and robust platform for large-scale transport simulations in complex geometries.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

tomoCAM : fast model-based iterative reconstruction via GPU acceleration and non-uniform fast Fourier transforms

X-ray-based computed tomography is a well established technique for determining the three-dimensional structure of an object from its two-dimensional projections. In the past few decades, there have been significant advancements in the brightness and detector technology of tomography instruments at synchrotron sources. These advancements have led to the emergence of new observations and discoveries, with improved capabilities such as faster frame rates, larger fields of view, higher resolution and higher dimensionality. These advancements have enabled the material science community to expand the scope of tomographic measurements towards increasingly in situ and in operando measurements. In these new experiments, samples can be rapidly evolving, have complex geometries and restrictions on the field of view, limiting the number of projections that can be collected. In such cases, standard filtered back-projection often results in poor quality reconstructions. Iterative reconstruction algorithms, such as model-based iterative reconstructions (MBIR), have demonstrated considerable success in producing high-quality reconstructions under such restrictions, but typically require high-performance computing resources with hundreds of compute nodes to solve the problem in a reasonable time. Here, tomoCAM , is introduced, a new GPU-accelerated implementation of model-based iterative reconstruction that leverages non-uniform fast Fourier transforms to efficiently compute Radon and back-projection operators and asynchronous memory transfers to maximize the throughput to the GPU memory. The resulting code is significantly faster than traditional MBIR codes and delivers the reconstructive improvement offered by MBIR with affordable computing time and resources. tomoCAM has a Python front-end, allowing access from Jupyter -based frameworks, providing straightforward integration into existing workflows at synchrotron facilities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Improved natural convection heat transfer correlations for reactor cavity cooling systems of high-temperature gas-cooled reactors: From computational fluid dynamics to Pronghorn

The Reactor Cavity Cooling System (RCCS) is a common reactor safety system in High Temperature Gas Cooled Reactors (HTGR) that removes heat from the Reactor Pressure Vessel (RPV) by radiation ($\sim 80\%$) and natural convection ($\sim 20\%$). For simulation of accident scenarios of HTGRs, intermediate fidelity and system codes models must be employed for limiting the models' execution time. While accurate quantification of the radiative heat transfer is available in these models, quantification of natural convection must rely on correlations of questionable accuracy for the Nusselt number. Commonly used correlations are based in experiments performed at low Rayleigh numbers and/or using isothermal walls in simplified geometries. Here, this work improves on the accuracy of natural convection heat transfer correlations in support for HTGR designs. These correlations include both local and average Nusselt numbers as a function of the global Rayleigh number, the local Rayleigh number, and the temperature profile at the hot wall of the RCCS. In the absence of dedicated experiments and the difficulty of performing high-fidelity simulations at realistic Rayleigh numbers, the data to fit the correlations are generated with Computational Fluid Dynamics (CFD) using Reynolds Averaged Navier-Stokes (RANS) models. First, a careful selection of the RANS turbulence model is performed by comparing the results obtained with different RANS turbulence models against high fidelity simulations of natural convection at $Ra \ 1 \times 10^{11}$ in a rectangular cavity. Next, the selected model is used to perform simulations of an HTGR cavity at different high Rayleigh numbers $\in [6.1 \times 10^{11},2.9 \times 10^{13}]$ to encompass several HTGR designs, assuming an isothermal RPV wall. The results obtained are used to fit a correlation for the average and space-varying Nusselt number as a function of the global and local Rayleigh numbers via a sparsity-promoting least-squares method. The selected RANS model is then used to perform simulations of a PBMR 400 HTGR cavity with the temperature profiles at the RPV wall obtained during a PLOFC transient. We use the results obtained to fit a temperature-dependent correction to the space-varying Nusselt number with the sparsity-promoting least-squares method. The results obtained in this work, enable system-level codes, such as Pronghorn, to perform higher-fidelity simulations of the heat exchange process in the RCCS while still maintaining a low computational cost.

42 ENGINEERING↗
DOE OSTI2020

An Alternative Casting Technique to Improve the Creep Resistance of Cast INCONEL Alloy 740H

The increasing performance requirements of power plant designs, such as advanced-ultra supercritical (A-USC), require the use of Ni-based superalloys to replace high-strength, ferritic-martensitic steels for components subjected to temperatures above 898 K (625 °C) and for austenitic stainless steels components at temperatures above 973 K (700 °C). To date, commercial Ni-based superalloy INCONEL 740H has been shown to be appropriate for use in A-USC power plants as boiler components in a wrought product. However, large complex components in boilers as well as other casings in the turbine and valve chest require castings of a thick-wall nature. Using the alloy in its cast form would be significantly valuable in terms of range of component size, geometry and complexity. Previous investigations revealed short creep lives from cast INCONEL alloy 740H. Here in this investigation, an alternative casting route that utilized a melt procedure resulting in a fine-grain casting, and in conjunction with a computationally optimized homogenization heat treatment, not only controlled the grain size and grain boundary structure but minimized chemistry variability and segregation. A primarily equiaxed and homogenous grain size distribution was obtained from this approach with better repartition of M 23 C 6 carbides along the grain boundaries. Furthermore, better than 38 pct increase was obtained for this material in comparison to the creep life obtained from the best performing conventionally cast material. More importantly, the fine-grain homogenized (FGH) casting route resulted in the Larson–Miller plot for this material that coincided with that of wrought alloy 740. At low creep stresses (with a test still in progress), the FGH casting is resulting in higher values of the LMP than the wrought alloy.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Impact of anisotropy on TRISO fuel performance

Manufacturing of tristructural isotropic (TRISO) particles involves the deposition of pyrolytic carbon (PyC) and silicon carbide (SiC) layers using the fluidized bed chemical vapor deposition (CVD) process. The CVD process is known to generate polycrystalline layers with crystallographic textures, which imparts anisotropic thermophysical properties to the layers. Past studies have shown the risk for particle failure increases with an increase in anisotropy. The limit beyond which the anisotropy of PyC layers becomes unacceptable due to failure risk has been identified as a high-priority knowledge gap. This work presents a first systematic study on the effects of anisotropic thermal and mechanical properties on TRISO fuel performance. This computational study, performed using the fuel performance code BISON, investigates how the anisotropy in elasticity and thermal properties affect the stresses, temperature, and failure of a TRISO particle. The influence of other factors, such as operating temperature and particle geometry on the anisotropy effects, also has been analyzed. The studies utilize the recently published anisotropic elasticity and thermal behavior models for TRISO PyC and SiC layers implemented using tensors with full anisotropic capability. The spherical TRISO particles with anisotropic properties were found to have greater maximum tensile stress and significantly higher failure probability than the spherical particles with isotropic properties. In conclusion, the fuel performance predicted using these recently developed models was found to be comparable with the performance obtained using the historical models.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Libration of hydroxyl groups in layered aluminum (oxy)hydroxides and other material analogs: insights from inelastic neutron scattering and theory

We analyzed the hydroxyl librational signatures of five structurally related aluminum (oxy)hydroxides, using inelastic neutron scattering (INS) and plane-wave lattice dynamics simulations. A clear trend across these aluminum-containing phases illustrates the relationship between hydrogen bonding, local atomic structure, and the spectral location and profile of the librational bands. The INS spectra have been compared to previous optical spectroscopy and computational studies, highlighting the complementary nature of the INS technique. Taking into account other structurally or chemically related material analogs, we have identified a correlation between a blueshift (to higher energy) of the upper librational band edge and the geometry of the hydrogen bond interactions, mirroring (with opposite correlation) the well-known redshift in the intramolecular O–H stretching energy with increasing hydrogen bond strength. For hydroxyl groups that do not participate in hydrogen bonding effectively, the bending librations occur at lower energies and hybridize with metal–oxygen lattice modes. Standard density functional theory approximations, including dispersion corrections, struggle to correctly predict vibrational frequencies of motions dominated by H but perform well for metal–oxygen modes, allowing us to make detailed mode assignments in several cases, including a demonstration of how layer-to-layer disorder in boehmite hydrogen bond orientations is reflected in the sharp but minor low energy peaks (at ∼70–80 meV) of the INS spectrum.

Wang, Hsiu-Wen [Oak Ridge National Laboratory (ORN↗
DOE OSTI2026

Molecular concentration field design using closed-form steady-state solutions

Control over spatial concentration fields represents a fundamental challenge in designing synthetic biological systems and programmable soft materials. While nature creates morphogen gradients that orchestrate complex developmental processes, synthetic approaches have largely relied on empirical optimization and computationally intensive simulations. Here, we present an analytical framework for steady-state concentration fields generated by finite-sized localized sources in diffusion–degradation systems and derive closed-form solutions for one-, two-, and three-dimensional geometries. By expressing these solutions in dimensionless form, we show that gradient steepness and spatial structure are organized by the Thiele modulus, which captures the competition between diffusion and degradation length scales. The analysis reveals distinct design regimes: in degradation-dominated systems, gradient shape is governed by exponential decay and becomes dimension-independent, whereas in diffusion-dominated systems, gradient magnitude and extent follow dimension-dependent power-law scaling. Building on these results, we introduce a quantitative design strategy that uses threshold-based criteria to program concentration ranges by tuning physically accessible parameters, most directly the production rate, while holding transport and degradation properties fixed. Comparisons with numerical solutions and reported experimental systems demonstrate consistency with the predicted scaling behavior. Together, this work provides a generalizable and physically transparent framework for designing steady-state concentration fields in synthetic biological and soft matter systems, enabling predictive control of gradient-mediated organization without reliance on extensive numerical optimization.

Kim, Dong Woo [Johns Hopkins University, Baltimore↗
DOE OSTI2021

Reconstructing magnetic deflections from sets of proton images using differential evolution

Proton imaging is a powerful technique for imaging electromagnetic fields within an experimental volume, in which spatial variations in proton fluence are a result of deflections to proton trajectories due to interaction with the fields. When deflections are large, proton trajectories can overlap, and this nonlinearity creates regions of greatly increased proton fluence on the image, known as caustics. The formation of caustics has been a persistent barrier to reconstructing the underlying fields from proton images. We have developed a new method for reconstructing the path-integrated magnetic fields, which begins to address the problem posed by caustics. Our method uses multiple proton images of the same object, each image at a different energy, to fill in the information gaps and provide some uniqueness when reconstructing caustic features. We use a differential evolution algorithm to iteratively estimate the underlying deflection function, which accurately reproduces the observed proton fluence at multiple proton energies simultaneously. We test this reconstruction method using synthetic proton images generated for three different, cylindrically symmetric field geometries at various field amplitudes and levels of proton statistics and present reconstruction results from a set of experimental images. Here, the method we propose requires no assumption of deflection linearity and can reliably solve for fields underlying linear, nonlinear, and caustic proton image features for the selected geometries and is shown to be fairly robust to noise in the input proton intensity.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Modeling of convective cells, turbulence, and transport induced by a radio-frequency antenna in the tokamak boundary plasma

The edge turbulence model Hermes (Dudson et al 2017 Plasma Phys. Control. Fusion 59 05401) is set up for plasma boundary simulations with an radiofrequency (RF) antenna, using parameters characteristic of a tokamak edge. Cartesian slab geometry is used with thin plate limiters representing the ion cyclotron range of frequency (ICRF) antenna side-wall limiters. Ad-hoc DC electric biasing of the limiters, motivated by calculations with VSim (Nieter et al 2004 J. Comput. Phys. 196 448), represents an induced RF sheath rectified potential in the plasma turbulence model. Flux-driven turbulence simulations demonstrate a realistic distribution of plasma profiles and fluctuations. There is a clear effect of the antenna sheath voltage leading to formation of convective cells; bias-induced convective transport flattens the scrape-off layer density profile and fluctuations penetrate into the shadow region of the limiters as the bias voltage increases. Turbulent transport for impurity ions is inferred by following ion trajectories in the simulated plasma turbulence fields, showing Bohm-like effective diffusion rates. All in all, the model elucidates the key physical phenomena governing the effects of ICRF-induced antenna biasing on the tokamak boundary plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Development of a neural network model for peeling–ballooning stability analysis in the KSTAR tokamak pedestals

The neural network model, MISHKA-NN is developed to mitigate the computational burden associated with the linear ideal magnetohydrodynamic (MHD) stability analysis of the pedestal based on the peeling–ballooning (P–B) model. By utilizing both 1D plasma profiles (current density, pressure gradient, and safety factor) and 0D parameters (plasma geometry, total current, and toroidal mode number), the model predicts linear growth rate of edge-localized ideal MHD instability in a given equilibrium state. By enabling the prediction of each instability within a second, the model reduces the time required for plotting a pedestal P–B stability diagram (j - α diagram) from approximately 100 CPU hours to a few CPU minutes. Notably, even with the utilization of parametric pressure and current profiles and plasma boundary shapes for the training dataset, the model shows a satisfactory level of performance in benchmarking the j - α diagram for the reconstructed equilibrium from a KSTAR tokamak experiment. We anticipate the model to serve as a versatile alternative to 2D linear MHD stability codes, alleviating numerical costs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A novel approach for tetrahedral-element-based finite element simulations of anisotropic hyperelastic intervertebral disc behavior

Intervertebral discs are microstructurally complex spinal tissues that add greatly to the flexibility and mechanical strength of the human spine. Attempting to provide an adjustable basis for capturing a wide range of mechanical characteristics and to better address known challenges of numerical modeling of the disc, we present a robust finite-element-based model formulation for spinal segments in a hyperelastic framework using tetrahedral elements. We evaluate the model stability and accuracy using numerical simulations, with particular attention to the degenerated intervertebral discs and their likely skewed and narrowed geometry. To this end, 1) annulus fibrosus is modeled as a fiber-reinforced Mooney-Rivlin type solid for numerical analysis. 2) An adaptive state-variable dependent explicit time step is proposed and utilized here as a computationally efficient alternative to theoretical estimates. 3) Tetrahedral-element-based FE models for spinal segments under various loading conditions are evaluated for their use in robust numerical simulations. For flexion, extension, lateral bending, and axial rotation load cases, numerical simulations reveal that a suitable framework based on tetrahedral elements can provide greater stability and flexibility concerning geometrical meshing over commonly employed hexahedral-element-based ones for representation and study of spinal segments in various stages of degeneration.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

First principles reaction discovery: from the Schrodinger equation to experimental prediction for methane pyrolysis

Our recent success in exploiting graphical processing units (GPUs) to accelerate quantum chemistry computations led to the development of the ab initio nanoreactor, a computational framework for automatic reaction discovery and kinetic model construction. In this work, we apply the ab initio nanoreactor to methane pyrolysis, from automatic reaction discovery to path refinement and kinetic modeling. Elementary reactions occurring during methane pyrolysis are revealed using GPU-accelerated ab initio molecular dynamics simulations. Subsequently, these reaction paths are refined at a higher level of theory with optimized reactant, product, and transition state geometries. Reaction rate coefficients are calculated by transition state theory based on the optimized reaction paths. The discovered reactions lead to a kinetic model with 53 species and 134 reactions, which is validated against experimental data and simulations using literature kinetic models. We highlight the advantage of leveraging local brute force and Monte Carlo sensitivity analysis approaches for efficient identification of important reactions. Both sensitivity approaches can further improve the accuracy of the methane pyrolysis kinetic model. The results in this work demonstrate the power of the ab initio nanoreactor framework for computationally affordable systematic reaction discovery and accurate kinetic modeling.

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

Assessment of Dynamic Fracture in Ultra-High Performance Concrete Using Synchrotron X-ray Source

Ultra-high-performance concretes (UHPCs) retain higher ultimate (compressive) strengths and greater durability relative to traditional concrete materials. As a result, they are enabling construction materials for blast resistant structures, where a higher capacity for energy absorption is required. In past investigations, the strength of UHPCs was shown to retain a dependency on strain rate. Still, the failure mechanisms which govern this rate dependency have yet to be completely elucidated. In the current study, the dynamic fracture behavior of Cor-tuf UHPC was assessed, in-situ, using a Synchrotron X-ray source in phase contrast imaging (PCI) mode. Prior to testing, the material’s morphology was characterized using X-ray computed tomography (CT) and scanning electron microscopy (SEM). A polychromatic pulsed X-ray source, high-speed camera and Kolsky bar apparatus were used in tandem, to observe the transient failure modes. A pin was placed at the end of the incident bar to control the location of damage initiation, and tests were performed at a constant velocity of ~6 m/s. These test conditions were observed to yield cracks with median and/or cone like geometries. Despite the variation in the crack paths, a significant difference in peak load was not observed. The failure mode was also observed to depend on the microstructure of the concrete samples.

Kolsky bar↗
DOE OSTI2020

On positive geometry and scattering forms for matter particles

We initiate the study of positive geometry and scattering forms for tree-level amplitudes with matter particles in the (anti-)fundamental representation of the color/flavor group. As a toy example, we study the bi-color scalar theory, which supplements the bi-adjoint theory with scalars in the (anti-)fundamental representations of both groups. Using a recursive construction we obtain a class of unbounded polytopes called open associahedra (or associahedra with certain facets at infinity) whose canonical form computes amplitudes in bi-color theory, for arbitrary number of legs and flavor assignments. In addition, we discuss the duality between color factors and wedge products, or "color is kinematics", for amplitudes with matter particles as well.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗