SEARCH · Engineering Papers

Results for “complex mixtures”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 667 records · Page 37

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Synthesis and Characterization of Divalent Samarium and Thulium N , N -Dimethylaminodiboranates

Here, the syntheses and molecular structures of new Sm II and Tm II N,N-dimethylaminodiboranate (DMADB) complexes are described. Treating SmI 2 (THF) 2 with Na(H 3 BNMe 2 BH 3 ) in THF results in the formation of Sm(H 3 BNMe 2 BH 3 ) 2 (THF) 3 (1), which can be readily converted to Sm(H 3 BNMe 2 BH 3 ) 2 (DME) 2 (DME = 1,2-dimethoxyethane) or Sm(H 3 BNMe 2 BH 3 ) 2 (diglyme) by exchange with the corresponding ether. We also show that Sm(H 3 BNMe 2 BH 3 ) 2 (THF) 3 can be prepared by reduction of the SmIII compound Sm(H 3 BNMe 2 BH 3 ) 3 (THF) with KC 8 and that addition of 18-crown-6 to this reaction mixture results in the formation of the Sm II compound Sm(H 3 BNMe 2 BH 3 ) 2 (18-crown-6). In a similar fashion, two new Tm II complexes have been synthesized: treatment of TmI 2 in THF with Na(H 3 BNMe 2 BH 3 ) results in the formation of Tm(H 3 BNMe 2 BH 3 ) 2 (THF) 2 and Tm(H 3 BNMe 2 BH 3 ) 2 (THF) 3 , which form a cocrystal. IR data and elemental analyses are reported for all the new compounds, as are their crystal structures. 1 H and 11 B NMR data are provided where available.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Chemical specificity in polyzwitterion-polyelectrolyte coacervates: polycations vs polyanions

Aqueous solutions containing polyzwitterions and polyelectrolytes were studied to probe the effects of chemical specificity on the complexation among different types of chains and resulting liquid–liquid phase separation (coacervation). Two kinds of blends were studied. A polyzwitterion, poly(sulfobetaine methacrylate) (PSBMA), was blended with either (1) a polycation, poly(diallydimethylammonium chloride) (PDADMAC), or (2) a polyanion, sodium poly(styrene sulphonate) (NaPSS). Coacervation was observed in both blends after equilibration, while the freshly prepared blend solutions exhibited distinct behavior dependent on the mixing protocols. Mixing solutions of the polyzwitterions and the polycations led to coacervation almost instantaneously. In contrast, the blends containing the same polyzwitterions and the polyanions exhibited no such coacervation by following the same mixing protocol. However, blends prepared after dissolving a solid mixture containing the same polyzwitterions and the polyanions in water exhibited coacervation. Based on the small-angle X-ray scattering (SAXS), cryo-electron microscopy, and rheology measurements, complexation and coacervation in the polyzwitterion-polyelectrolyte blends is rationalized in terms of the entropy change associated with reorganization of solvent (water), which appears to be larger in the polyzwitterion-polycation case in comparison with the polyzwitterion-polyanion blends. These results highlight the importance of chemical specificity and mixing protocols leading to metastable structures in designing a new class of materials based on macromolecular complexation of polyzwitterions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1993

Dynamics of flames near the rich-flammability limit of hydrogen-air mixtures

Flames near the rich-flammability limit of hydrogen-air mixtures are studied using a detailed, time-dependent, one-dimensional Lagrangian model. Results from the numerical simulations indicate that a steady burning velocity is not obtained for very rich hydrogen-air mixtures. As the amount of hydrogen is increased, a damped oscillation is observed in the flame and burning velocities, and then, with further increase in the amount of hydrogen, an undamped oscillation with a complex set of frequencies is observed. Simulations with a simplified one-step irreversible chemical reaction do not show these oscillations, suggesting that chemical kinetics plays a strong role in inducing these oscillations. Further analysis shows that the oscillations are due to a competition for H atoms between chain branching and chain-terminating reactions. Simulations of spherically expanding flames suggest that stretch effects (due to curvature) will cause the oscillations to occur in less rich mixtures than that observed for planar flames. The implications of these oscillations on the rich-flammability limit as well as the role of chemical kinetics in creating a fundamental flammability limit is discussed.

Kailasanath, K.↗

DOE OSTI2024

Molecular dynamics simulations of uranyl and plutonyl cations in a task-specific ionic liquid

Ionic liquids (ILs) are a unique class of solvents with potential applications in advanced separation technologies relevant to the nuclear industry. ILs are salts with low melting points and a wide range of tunable physical properties, such as viscosity, hydrophobiciy, conductivity, and liquidus range. ILs have negligible vapor pressure, are often non-flammable, and can have high thermal stability and a wide electrochemical window, making them attractive for use in separations processes relevant to the nuclear industry. Metal salts generally have a low solubility in ILs; however, by incorporating new functional groups into the IL cation or anion that promote complexation with the metal, the solubility can be greatly increased. One such task-specific ionic liquid (TSIL) is 1-carboxy-N, N, N-trimethylglycine bis(trifluoromethylsulfonyl)imide ([Hbet][Tf 2 N]). Water, which is detrimental for electrochemical separations, is a common impurity in ILs and can coordinate with actinyl cations, particularly in ILs containing only weakly coordinating components. Understanding the behavior of actinides in TSIL/water mixtures on a molecular level is vital for designing improved separations processes. Classical molecular dynamics simulations of uranyl(VI) and plutonyl(VI) in 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ([EMIM][Tf 2 N]) with deprotonated Hbet (betaine) and water have been performed to understand the coordination and dynamics of the actinyl cations. We find that betaine is a much stronger ligand than water and prefers to coordinate the metal in a bidentate manner. Potential of mean force simulations yield a relative free energy for betaine coordination of approximately -120 to -90 kJ/mol in mixtures with water. As the amount of betaine coordinated to the actinide increases, the diffusion coefficient of the actinyl cation decreases. Moreover, the betaine ligand is able to bridge between two metal centers, resulting in dimeric complexes with actinide–actinide distances of ~5 Å. Potential of mean force simulations show that these structures are stable, with relative free energies of up to -40 kJ/mol. The crystal structure for [(UO 2 ) 2 (bet) 6 (H 2 O) 2 ][Tf 2 N] 4 shows that the betaine bridges between two uranium atoms to form dimeric complexes similar to those found in our simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Verification of a specialized hydrodynamic simulation code for modeling deflagration and detonation of high explosives

A specialized hydrodynamic simulation code has been developed and verified for the simulation of one-dimensional unsteady problems involving the detonation and deflagration of high explosives. To model all the relevant physical processes in these problems, a code is required to simulate compressible hydrodynamics, unsteady thermal conduction, and chemical reactions with complex rate laws. Several verification exercises are presented which test the implementation of these capabilities. The code also requires models for physics processes such as equations of state and conductivity for pure materials and mixtures as well as rate laws for chemical reactions. Additional verification tests are required to ensure that these models are implemented correctly. Though this code is limited in the types of problems it can simulate, its computationally efficient formulation allows it to be used in calibration studies for reactive burn models for high explosives. Furthermore, this study demonstrates how a series of verification tests can be used to ensure that the various physics processes needed to simulate complex phenomenon can be tested to ensure that they are correctly implemented.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Comprehensive Synthesis and Structural Trends in Tetramethyl Diglycolamide (TMDGA) Nitrate Complexes with Lanthanides and Americium

Complexes of N,N,N',N'-tetramethyl diglycolamide (TMDGA), a hydrophilic diglycolamide (DGA) proposed as an aqueous phase holdback reagent, have been crystallized for the majority of the lanthanide series (excluding promethium), yttrium, and americium to deepen our structural understanding of trivalent metal ion (M 3+ ) DGA coordination compounds in the presence of nitrate counter-anions. The presented collection of 16 complexes with accompanying single-crystal structures, taking formulas [M(TMDGA) 3 ][M(NO 3 ) 6 ] (M = La, Ce, Pr, Nd, Sm, Am), [M(TMDGA) 3 ][M(NO 3 ) 5 (H 2 O)] 1–x [M(NO 3 ) 4 (H 2 O) 2 ] x (NO 3 ) 1+x (M = Eu, Gd, Tb, Dy, Ho, Er, Tm, Yb), [M(TMDGA) 3 ] 2 [M(NO 3 ) 4 (H 2 O) 2 ] 0.75 [M(NO 3 ) 5 (H 2 O)] 1.25 (NO 3 ) 2.75 ·H 2 O (M = Lu), and [M(TMDGA) 3 ][M(NO 3 ) 5 (H 2 O)](NO 3 )·CH 3 OH (M = Y) were all synthesized via solvent diffusion of reaction mixtures containing the metal nitrate M(NO 3 ) 3 ·nH 2 O and TMDGA. Single-crystal X-ray diffraction analyses of these new structures show that each TMDGA complex comprises of three TMDGA ligands coordinating the metal ion via carbonyl and etheric oxygen atoms forming [M(TMDGA) 3 ] 3+ cations. Spectroscopy measurements under high pressure displayed notable differences in the f → f transition shifting between that of Nd(III) and Am(III). Shifting of transitions by 2.0 nm were observed in Nd up to 8.50 ± 0.09 GPa, while Am saw shifting between 11.0 to 13.5 nm at pressures up to 20.06 ± 1.90 GPa. The local geometry in these complexes is a distorted spherical capped square antiprism (CSAPR-9) except for of the yttrium complex, which exhibits a distorted spherical tricapped trigonal prismatic (TCTPR-9) geometry. Further, the anions that form concomitantly with the TMDGA complexes are composed of hexanitrato species for the early lanthanide ions (lanthanum to samarium); whereas the remaining smaller lanthanides did not possess sufficiently large ionic radii to coordinate six bidentate nitrate anions, instead, one or two nitrate anions are situated in the outer sphere. The systematic progression of changes in the anionic environments of these complexes outlines the changing coordination habits afforded by the lanthanide contraction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Computational analysis of flame initiation, quenching, and re-ignition in a prechamber natural gas engine under varying EGR-dilution levels

The on-road natural-gas (NG) fueled transportation relies on stoichiometric spark-ignition engines for the advantages of simple after-treatment system despite the efficiency penalty relative to lean-burn combustion strategies. Exhaust gas recirculation (EGR) has the potential to reduce this efficiency gap at low to moderate loads without the need for complex lean-exhaust aftertreatment systems. However, EGR dilution leads to reduced combustion stability and increased cycle-to-cycle variability. A promising technology that has the potential to achieve reliable operation under diluted conditions is the prechamber ignition (or turbulent jet ignition) which uses chemically active turbulent jets generated from combustion inside a prechamber to initiate, stabilize and accelerate combustion of the mixture inside the main chamber. The present work focusses on developing a RANS-based CFD approach to accurately reproduce in-cylinder phenomena in a stoichiometric NG prechamber-assisted heavy-duty engine without relying on complex combustion models that account for turbulence-chemistry interactions. This is necessary because reactive prechamber jets at high EGR dilution tend to extinguish while emerging into the main chamber, which is followed by a phase of re-ignition — a phenomenon that conventional G-equation or well-stirred reactor combustion models cannot reproduce. With addition of a damping multiplier to the well-stirred reactor model, the predictions are seen to show good agreement with experimental pressure evolution and combustion images acquired from a single cylinder Cummins N-14 optical diesel engine retrofitted with a prechamber ignition system. Model predictions of local heat release in the flame and temperature evolution inside the flame are used to investigate combustion dynamics in the prechamber and the main chamber. It is seen that the well-stirred reactor model with the inclusion of damping is able to reproduce the temporary reduction in heat release within the flame, which can be considered equivalent to quenching of jets, and the subsequent re-ignition of the flame inside the main chamber. The delay between quenching and re-ignition depends on the amount of dilution, as explained by an illustration of flame evolution in a Borghi diagram.

Prechamber ignition↗

DOE OSTI2023

Direct Transformation of SiH 4 to a Molecular L(H) 2 Co=Si=Co(H) 2 L Silicide Complex

The synthesis of bimetallic molecular silicide complexes is reported, based on the use of multiple Si–H bond activations in SiH 4 at the metal centers of 14-electron LCo I fragments (L = Tp", HB(3,5-diisopropylpyrazolyl) 3 – ; [BP 2 tBu Pz], PhB(CH 2 P t Bu 2 ) 2 (pyrazolyl)). Upon exposure of (Tp"Co) 2 (μ-N 2 ) (1) to SiH 4 , a mixture of (Tp"Co) 2 (μ-H) (2) and (Tp"Co) 2 (μ-H) 2 (3) was formed and no evidence for Si–H oxidative addition products was observed. In contrast, [BP 2 tBu Pz]-supported Co complexes led to Si–H oxidative additions with the generation of silylene and silicide complexes as products. Notably, the reaction of ([BP 2 tBu Pz]Co) 2 (μ-N 2 ) (5) with SiH 4 gave the dicobalt silicide complex [BP 2 tBu Pz](H) 2 Co=Si=Co(H) 2 [BP 2 tBu Pz] (8) in high yield, representing the first direct route to a symmetrical bimetallic silicide. Here, the effect of the [BP 2 tBu Pz] ligand on Co–Si bonding in 7 and 8 was explored by analysis of solid-state molecular structures and density functional theory (DFT) investigations. Upon exposure to CO or DMAP (DMAP = 4-dimethylaminopyridine), 8 converted to the corresponding [BP 2 tBu Pz]Co(L) x adducts (L = CO, x = 2; L = DMAP, x = 1) with concomitant loss of SiH 4 , despite the lack of significant Si–H interactions in the starting complex. On heating to 60 °C, 8 underwent reaction with MeCl to produce small quantities of Me x SiH 4–x (x = 1–3), demonstrating functionalization of the μ-silicon atom in a molecular silicide to form organosilanes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Strategic Lifetime Tuning of Visible-Light Absorbing Two-Coordinate Metal Complexes

This paper highlights how the singlet and triplet amide ligand centered ( 1,3 LC) and interligand charge transfer ( 1,3 ICT) states in (carbene)M(amide), M = Cu, Au (cMa) complexes influence the excited state properties when the triplet states are close in energy. To that end we prepared a set of five cMa complexes, in which the amide (i.e., 5H-benzo[b]carbazole, H-BnCz) was kept constant giving a 3 LC energy of ca. 2.15 eV. Four different carbene ligands were selected to develop M BnCz Carbene complexes which have energies for the ICT state that vary from being either markedly higher, lower, or close to the energy of the 3 LC state on the amide ligand. Steady-state and time-resolved spectroscopic studies show that the emission spectrum of the cMa complex mirrors the phosphorescence of the parent H-BnCz amide and has a luminescence decay lifetime in the millisecond regime when the lowest energy excited state is 3 LC. When the lowest energy states are 1,3 ICT in nature, the emission band is broad and featureless, giving a lifetime in the nanosecond regime. As the 3 LC and 1,3 ICT states have comparable energies, dynamic equilibrium between the states is observed with the luminescence consisting of a mixture of 3 LC and 1,3 ICT transitions. The measured lifetime of this equilibrating system is between 45 and 350 μs depending on the solvent, for both copper- and gold-based cMa complexes. Furthermore, these complexes demonstrate the ability to “park” the excited state population in the 3 LC state, allowing it to act as a reservoir for thermally activated emission, while still maintaining a very rapid equilibrium between LC and ICT states.

14 SOLAR ENERGY↗

DOE OSTI2023

Mosaic Cu I −Cu II −In III 2D Perovskites: Pressure‐Dependence of the Intervalence Charge Transfer and a Mechanochemical Alloying Method

Abstract The perovskite (BA) 4 [Cu II (Cu I In III ) 0.5 ]Cl 8 ( 1 BA ; BA + =butylammonium) allows us to study the high‐pressure structural, optical, and transport properties of a mixed‐valence 2D perovskite. Compressing 1 BA reduces the onset energy of Cu I/II intervalence charge transfer from 1.2 eV at ambient pressure to 0.2 eV at 21 GPa. The electronic conductivity of 1 BA increases by 4 orders of magnitude upon compression to 20 GPa, when the activation energy for conduction decreases to 0.16 eV. In contrast, Cu II perovskites achieve similar conductivity at ≈50 GPa. The solution‐state synthesis of these perovskites is complicated, with more undesirable side products likely from the precursor mixtures containing three different metal ions. To circumvent this problem, we demonstrate an efficient mechanochemical synthesis to expand this family of halide perovskites with complex composition by simply pulverizing together powders of 2D Cu II single perovskites and Cu I In III double perovskites.

Li, Jiayi↗

DOE OSTI2026

Synthesis and Reactivity of Heteroleptic U 4+ Alkyl, Benzyl, and Hydride Imidophosphorane Complexes

A series of heteroleptic U 4+ benzyl, neopentyl, and methyl complexes supported by the imidophosphorane ligand, [N = P(N,N′-ditert-butylethylenediamide)(diethylamide)] 1− (NP*), were synthesized from the monoiodide precursor, [UI(NP*) 3 ]. These heteroleptic complexes were synthesized through the selective formation of [UI(NP*) 3 ] under transmetalation conditions in the reaction between [UI 4 (1,4-dioxane) 2 ] and K[NP*]. Formation of the homoleptic complex [U(NP*) 4 ] was not observed even in the presence of excess K[NP*]. The oxidation and hydrogenolysis reactivity of the neopentyl complex, [U(Npt)- (NP*) 3 ] (Npt = neopentyl) was explored. While cyclic voltammetry indicates a potentially isolable U5+ alkyl cation, chemical oxidation of the neopentyl complex results in the isolation of a cationic U 4+ complex with a bound diethyl ether in the primary coordination sphere, [U4+(NP*)) 3 (Et 2 O)][(BArF 24 )] (BArF 24 = tetrakis(3,5-bis(trifluoromethyl)phenyl)borate). Notably, hydrogenolysis of [U(Npt)(NP*) 3 ] with H2 gas at −20 °C results in the formation of a terminal hydride intermediate confirmed by in situ NMR spectroscopy and deuterium labeling with D 2 . The connectivity and structural parameters of this hydride intermediate, [UH(NP*) 3 ], which rapidly thermally decomposes to the homoleptic complex, [U(NP*) 4 ], can be confirmed by single-crystal X-ray diffraction studies of a crystal grown by chilling the reaction mixture. The identity of [U(NP*) 4 ] was confirmed by its direct, bulk synthesis from [U(Me)(NP*) 3 ] and HNP* in a protonolysis reaction.

Alkyls↗

DOE OSTI2023

Epitaxial growth of hexagonal BaNi O 3–δ thin films on SrTi O 3 (111) substrates

Transition metal oxides containing nickel species in oxidation states higher than 3+ often exhibit high catalytic activity, which makes them promising for applications as advanced electrocatalysts for water splitting and fuel cells. Here, we examine the structure and properties of BaNiO 3 (BNO) thin films, containing formally Ni 4+ , grown on SrTiO 3 (111) substrates using oxygen-plasma–assisted molecular beam epitaxy. X-ray diffraction and scanning transmission electron microscopy measurements reveal that BNO films have a hexagonal structure with $c$- and $a$-axes of mixed textures showing epitaxial relationships BNO (0001) ∥ SrTiO 3 (111) and BNO (10$\overline{1}$0) ∥ SrTi O 3 (111), respectively. The formation of the $a$-axis texture is dominant due to the smaller lattice mismatch with the substrate. Here, density functional theory calculations confirm that the hexagonal BNO film with the $a$-axis texture is energetically more favorable than the competing $c$-axis texture. Detailed spectroscopy data analysis indicates that hexagonal BNO films contain mixtures of Ni 2+ , Ni 3+ , and Ni 4+ species, with Ni 4+ being dominant. Our study provides insights into stabilizing Ni 4+ in complex oxides, which is important for further exploration of the potential of materials containing Ni 4+ .

36 MATERIALS SCIENCE↗

DOE OSTI2025

FuelLib (Fuel Library) [SWR-25-26]

FuelLib is a library that utilizes the group contribution method (GCM) for calculating thermodynamic properties of hydro-carbon jet fuels. FuelLib utilizes the tables and functions of the GCM as proposed by Constantinou and Gani (1994) and Constantinou, Gani and O'Connel (1995), with additional physical properties discussed in Govindaraju & Ihme (2016). The code is based on Pavan B. Govindaraju's Matlab implementation of the GCM, and has been expanded to include additional thermodynamic properties and mixture properties. The fuel library contains gas chromatography (GC x GC) data for a variety of fuels ranging from simple single component fuels to complex jet fuels. The GC x GC data for POSF jet fuels comes from Edwards (2020).

Montgomery, David [National Renewable Energy Labor↗

DOE OSTI2023

Fuel Injection Dynamics and Composition Effects on RDE Performance

Rotating detonation engines (RDEs) provide a promising route to substantially increasing cycle efficiency in stationary gas turbines. Much of this increase relies on the ability to achieve consistent pressure gain within the combustor. In particular, the design of injectors that feed fuel and air into the detonation channel plays a crucial role. Such injectors have to ensure proper mixing of fuel and oxidizer, while minimizing backflow of detonation products into the feed plenums, and reduce susceptibility to the complex wave structures that exist within the combustor. From a practical perspective, such RDEs also need to operate with variable fuel composition. When fuel mixtures with components that possess vastly different oxidation pathways and time-scales are used, there could be additional losses through deflagrative burning instead of detonation-driven heat release. Such sensitivity to the complete flow path is akin to the physics of thermoacoustic instabilities in conventional gas turbines. In this sense, RDEs pose a unique research challenge: the performance of the device relies on the small-scale heat release process, which is highly dependent on the flow interactions within the full-scale system. As a result, canonical flow configurations, instrumented with detailed diagnostics or modeled using high-fidelity tools, but only focus on the small-scale processes will not contain the key system-level interactions. At the same time, macroscopic measurements and models that only capture system-level performance will not provide insight into the key sources of pressure losses. These couplings and sensitivities provide a formidable challenge to both experimental and simulation studies of the effects Thus, a joint experimental/computational program designed specifically to address these challenges was undertaken in this program. The focus of this program was on two key topics: a) the interaction between injector flow and the overall wave dynamics within the combustor, and b) the deflagration/detonation structure in multi-component fuels that are of practical interest. Both topics involve interaction of small-scale heat release processes with the geometry-dependent wave structure. Studies focused on the study of full-scale RDE systems, based on a 6-inch conventional annular geometry. Experimentally RDEs were studied using a combination of diagnostics. A combination of optical diagnostics and aero-thermo-acoustic analysis based on a combination of spectral and mode decomposition analysis was used to identify the dynamics of the detonation wave and other secondary waves that exist in the system. These studies have helped the identification and investigation of injector and detonation dynamics arising from coupling, and how they affect RDE mixing, detonation structure, operability and performance. Performance of RDEs was investigated through thrust stand measurements, which was used to evaluate the effective pressure gain generated by the system through the concept of equivalent available pressure. Optical diagnostics were developed and implemented to investigate the distribution of heat release, across the detonation wave. Novel optical diagnostics of NIR imaging was also developed and applied to investigate the high temperature / high pressure distribution across the detonation wave. In order to complement the experiments, the computational tools were geared to simulate the full experimental setup. GPU-based acceleration of the models and computations were developed to enable rapid simulation of the full system. In addition, the use of adaptive mesh refinement, and unstructured grid formulation, enabled the investigation of realistic geometries studied in the laboratory. The simulations produced a wealth of detail on the structure of the detonation wave under different operating conditions. Emphasis was placed on quantifying mixture pre-burning and the impact on wave propagation and structure.

03 NATURAL GAS↗

NASA NTRS1998

Interaction of Vascular Smooth Muscle Cells Under Low Shear Stress

The blood vessel wall consists of three cellular layers, an outer adventitial, a middle medial and an inner intimal layer. When the blood vessel forms in the embryo it begins as a tube composed of a single cell type called endothelial cells. Over time, other cells are recruited from the surrounding tissue to form additional layers on the outer surface of the endothelial tube. The cells that are recruited are called mesenchymal cells. Mesenchymal cells are responsible for the production of connective tissue that holds the blood vessel together and for developing into vascular smooth muscle cells that are responsible for regulating the diameter of the vessel (1) and therefore, blood flow. In a fully developed blood vessel, the endothelial cells make- up the majority of cells in the intimal layer while the mesenchymal cells make-up the majority of cells in the medial and adventitial layers. Within the medial layer of a mature vessel, cells are organized into multiple circular layers of alternating bands of connective tissue and cells. The cell layer is composed of a mixture of mesenchymal cells that have not developed into smooth muscle cells and fully developed smooth muscle cells (2). The assembly and organization of complex tissues is directed in part by a signaling system composed of proteins on the cell surface called adhesion molecules. Adhesion molecules enable cells to recognize each other as well as the composition of the connective tissue in which they reside (3). It was hypothesized that the different cell types that compose the vascular wall possess different adhesion molecules that enable them to recognize each other and through this recognition system, form the complex layered organization of the vascular wall. In other words, the layered organization is an intrinsic property of the cells. If this hypothesis is correct then the different cells that make up the vessel wall, when mixed together, should organize themselves into a layered structure resembling an intact blood vessel. Experiments described below were designed to test this hypothesis.

Seidel, Charles L.↗

DOE OSTI2023

Development of Binder Jet Compatible Silicon Carbide Mixtures for Increased Sintered Density

Binder jet additive manufacturing can process a variety of ceramic powder feedstocks to manufacture near net shape green parts that are subsequently sintered to achieve densification. In particular, silicon carbide (SiC) is a promising ceramic material for heat exchangers in energy storage applications and binder jet is a well-suited process for manufacturing these heat exchangers due to the often geometrically complex designs. However, the need for good powder flowability in binder jet printing often limits the feedstock powder to be larger particles rather than that of ultra fine submicron powders used in conventional ceramic processing, which are more readily densified during sintering. The present work explores the use of various powder mixtures and sintering settings to identify strategies for achieving high density sintered parts in binder jet compatible systems without the use of infiltration. Liquid-phase sintering additives and bimodal powder mixtures are explored alongside different sintering temperatures, rates, atmospheres, and hold times. Two mixtures were found to be good candidates for further exploration: 1) 90 wt% 0.6 µm SiC + 6.25 wt% Al 2 O 3 + 3.75 wt% Y 2 O 3 mixture achieved the highest sintered relative density of 91.2%; and 2) a bimodal mixture comprised of 45 wt% 10 µm SiC + 45 wt% 0.6 µm SiC + 6.25 wt% Al 2 O 3 + 3.75 wt% Y 2 O 3 had the best powder flowability of the powders measured and a 75.1% sintered relative density.

25 ENERGY STORAGE↗

DOE OSTI2026

From molecular to macroscopic: predicting liquid–liquid phase equilibria and small-angle scattering of mixtures of organic liquids from atomistic simulation using Kirkwood–Buff theory

Macroscopic phase equilibria between solutions define the functionality of many biological and industrial processes, yet they are challenging to predict due to the inherent complexity of liquids containing large molecules. This work introduces an approach for the purely predictive calculation of such phase equilibria in temperature-composition space from molecular dynamics (MD) simulations at one temperature in the single-phase region. We use an approach developed previously to obtain the entropic and enthalpic contributions to the free energy of mixing from the atomic-scale information given by MD simulations via Kirkwood–Buff theory. This allows us to accurately estimate the free energy of mixing as a function of temperature, and thus obtain liquid–liquid phase equilibria, including liquid–liquid critical points, associated binodal and spinodal lines, and composition fluctuations across a region of temperature and composition. Results for binary malonamide–alkane systems are validated by comparison to a direct experimental probe of the fluctuations: the small angle X-ray scattering intensity near zero wavenumber. The MDKB → Phase method demonstrated here provides a significant improvement in predicting liquid–liquid equilibria and free energy as a function of temperature for our systems of interest compared to conventional thermodynamic models. The accurate performance of this purely predictive approach lies in its preservation of atomistic details when determining thermodynamic properties. Furthermore, its inherent extensibility to multi-component systems will likely make the MDKB → Phase approach a valuable general tool for connecting molecular interactions to macroscopic phase equilibria and for the computational screening of materials for targeted thermodynamic behavior.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Mathematical Models and Numerical Methods for High-Fidelity Simulation of Ignition of Reactive Mixtures by Nanosecond Plasma Discharges in Realistic Configurations

We present a newly developed framework for the numerical simulation of ignition of reactive mixtures using single or repeated nanosecond discharge pulses. The framework builds upon the AMReX library, using the existing compressible solver PeleC and low-Mach solver PeleLMeX and allowing for adaptive mesh refinement, complex geometries, and execution on next-generation high-performance computing (HPC) systems. High-fidelity elementary models are adopted for weakly-ionised plasma discharges with significant energy deposition, consistent with nanosecond discharge pulses, and then implemented in the solver. The treatment of non-thermal electrons and charged species, thermodynamics of non-equilbrium species, plasma kinetics, limiting time scales, and boundary conditions for charged species are discussed and addressed for computational efficiency. The framework is demonstrated for three relevant applications: single and multi-pulse discharges in air, single pulse ignition of an ethylene/air mixture, and a three-dimensional plasma discharge in air with temperature stratification. The successful application of the framework demonstrates the feasibility of high-fidelity simulation of ignition of air/hydrocarbon mixtures in three-dimensions with multiple discharge pulses.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗