Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multicomponent”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

A tool and a methodology to use macros for abstracting variations in code for different computational demands

Scientific software used on high-performance computing platforms is in a phase of transformation because of the combined increase in the heterogeneity and complexity of models and hardware platforms. Having separate implementations for different platforms can easily lead to combinatorial explosions; therefore, the computational science community has been looking for mechanisms to express code through abstractions that can be specialized for different platforms. Most existing approaches use template meta-programming in C++, and are, therefore language specific. Here, we have developed a tool that uses customized expansion of macros to mimic some of C++ behavior in other languages. It enables unification of any code variants that may be necessary to run efficiently on different target architectures and different computational environments through use of macros with multiple alternative definitions and ability to arbitrate on definition selection for expansion. Combined with two other tools, a custom runtime, and a user specified recipe translator, our custom macroprocessor becomes a part of an overall performance portability solution that does not depend on any specific programming language. We also use macros as code-shorthand that lets code snippets become building blocks that allow variations in control flow to explore performance options. We demonstrate use of macros in Flash-X, a multiphysics multicomponent code with many Fortran legacy components derived from an earlier community code FLASH.

Heterogenous computing↗
DOE OSTI2021

TOUGH-RFPA: Coupled thermal-hydraulic-mechanical Rock Failure Process Analysis with application to deep geothermal wells

This paper presents a coupled thermal-hydraulic-mechanical (THM) model for simulating failure processes in geological media. The coupling between stress/strain and heat transfer describing how thermal expansion affects stress is based on the theory of thermo-elasticity. Additionally, the coupling between stress/strain and fluid flow describing how pore-pressure affects stress is governed by Biot's consolidation theory. Stress/strain-dependent porosity/permeability/capillary-pressure defines the influence of stress on fluid flow and thermal fields. The THM model is implemented in the solver module that links the multiphase, multicomponent fluid flow simulator TOUGH2 with the geomechanical simulator RFPA (Rock Failure Process Analysis), which is based on the concept that heterogeneity leads to non-linearity and progressive failure behavior. Capability of TOUGH-RFPA simulator was verified by two benchmark problems related to hydraulic-mechanical and thermal-mechanical response in a hollow cylinder. Then the applicability of TOUGH-RFPA was demonstrated for modeling cooling-induced damage around a very deep geothermal well with mesoscopic heterogeneous properties, which showed complex damage evolution and fracturing during well cooling and fluid pressurization.

42 ENGINEERING↗
DOE OSTI2022

Targeting the Y-box Binding Protein-1 Axis to Overcome Radiochemotherapy Resistance in Solid Tumors

Multifunctional Y-box binding protein-1 (YB-1) is highly expressed in different human solid tumors and is involved in various cellular processes. DNA damage is the major mechanism by which radiochemotherapy (RCT) induces cell death. On induction of DNA damage, a multicomponent signal transduction network, known as the DNA damage response, is activated to induce cell cycle arrest and initiate DNA repair, which protects cells against damage. YB-1 regulates nearly all cancer hallmarks described to date by participating in DNA damage response, gene transcription, mRNA splicing, translation, and tumor stemness. YB-1 lacks kinase activity, and p90 ribosomal S6 kinase and AKT are the key kinases within the RAS/mitogen-activated protein kinase and phosphoinositide 3-kinase pathways that directly activate YB-1. Thus, the molecular targeting of ribosomal S6 kinase and AKT is thought to be the most effective strategy for blocking the cellular function of YB-1 in human solid tumors. In this review, after describing the prosurvival effect of YB-1 with a focus on DNA damage repair and cancer cell stemness, clinical evidence will be provided indicating an inverse correlation between YB-1 expression and the treatment outcome of solid tumors after RCT. In the interest of being concise, YB-1 signaling cascades will be briefly discussed and the current literature on YB-1 posttranslational modifications will be summarized. Finally, the current status of targeting the YB-1 axis, especially in combination with RCT, will be highlighted.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2023

Uncertainty quantification of a deep learning fuel property prediction model

Deep learning models are being widely used in the field of combustion. Given the black-box nature of typical neural network based models, uncertainty quantification (UQ) is critical to ensure the reliability of predictions as well as the training datasets, and for a principled quantification of noise and its various sources. Deep learning surrogate models for predicting properties of chemical compounds and mixtures have been recently shown to be promising for enabling data-driven fuel design and optimization, with the ultimate goal of improving efficiency and lowering emissions from combustion engines. In this study, UQ is performed for a multi-task deep learning model that simultaneously predicts the research octane number (RON), Motor Octane Number (MON), and Yield Sooting Index (YSI) of pure components and multicomponent blends. The deep learning model is comprised of three smaller networks: Extractor 1, Extractor 2, and Predictor, and a mixing operator. The molecular fingerprints of individual components are encoded via Extractor 1 and Extractor 2, the mixing operator generates fingerprints for mixtures/blends based on linear mixing operation, and the predictor maps the fingerprint to the target properties. Two different classes of UQ methods, Monte Carlo ensemble methods and Bayesian neural networks (BNNs), are employed for quantifying the epistemic uncertainty. Combinations of Bernoulli and Gaussian distributions with DropConnect and DropOut techniques are explored as ensemble methods. All the DropConnect, DropOut and Bayesian layers are applied to the predictor network. Aleatoric uncertainty is modeled by assuming that each data point has an independent uncertainty associated with it. The results of the UQ study are further analyzed to compare the performance of BNN and ensemble methods. Although this study is confined to UQ of fuel property prediction, the methodologies are applicable to other deep learning frameworks that are being widely used in the combustion community.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2020

Solid-solid phase transition via the liquid in a Pd 43 Cu 27 Ni 10 P 20 bulk metallic glass under conventional conditions

The microscopic mechanism of solid-solid phase transitions is a long-standing fundamental issue in materials science. Here we directly report the experimental evidence of the existence of an intermediate liquid during the solid-solid phase transition in a bulk Pd 43 Cu 27 Ni 10 P 20 alloy by in situ high temperature high energy X-ray diffraction, combined with the differential scanning calorimetry at slow heating rates down to 5 K/min. We elucidate that this intermediate liquid is attributed to the melting of the low-melting-point crystalline phases, rather than induced by the elastic stress energy along the interface proposed in literatures. This could also be a general feature in multicomponent alloys if a low-melting-point crystalline phase is formed together with other crystalline phases during heating. Furthermore, all results obtained here are scientifically sound and provide new insight into the phase transition theory, which will trigger more studies on the microscopic mechanism of solid-solid phase transition via an intermediate liquid.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Revisit the VEC criterion in high entropy alloys (HEAs) with high-throughput ab initio calculations: A case study with Al-Co-Cr-Fe-Ni system

Valence electron concentration (VEC) was treated as a useful parameter to predict the stability of solid solution phases. However, the available experimental data to support this criterion is far from enough. In the current study, the high-throughput ab initio modeling is applied to investigate the relative stability of FCC and BCC single crystals of the Al-Co-Cr-Fe-Ni high entropy alloys (HEAs) by using the special quasi-random structure (SQS) approach. Furthermore, the predictions start with pure elements of the Al-Co-Cr-Fe-Ni system and are continued with binaries, ternaries, and quaternary compositions, which come up with 180 compositions (360 structures). After that, the reliability of the VEC criterion is testified. The results show that the VEC criterion not only works for the stable structure but also works effectively for metastable structure when both FCC and BCC are not thermodynamic stable. However, it is found that the old VEC criterion proposed by Guo et al. fails to work effectively for compositions containing high concentrations of light-weight metals such as Al at VEC<5. To solve this problem, the present work proposes a new VEC rule to define the stability of FCC and BCC structures at the ground state. With the implementation of the new VEC rule, the effectiveness of the VEC rule (EVEC) of both FCC and BCC structures is enhanced, especially for pure elements and binary compositions, indicating that this rule does not only work effectively for multicomponent systems but also works for low-order systems.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Heterogeneous phase transformation pathways in additively manufactured Al-Ce-Mn alloys

Heat treatment of additively manufactured Al-Ce based multicomponent alloys leads to complex microstructure evolution. In this research, the ability to extend the phase transformation theories involving nucleation of a product phase from a heterogeneous multi-phase microstructure typical to that of additively manufactured samples is explored. The Al-10Ce-8Mn (wt%) was used as a model alloy system. Under additive manufacturing conditions different solidification microstructures were obtained due to spatial and temporal variations of thermal gradients (G) and liquid-solid interface velocities (R) within a given melt pool. Near the melt pool boundary (high G and low R, referred as MPB region), initially, Al 20 Mn 2 Ce forms from the liquid followed by a eutectic of FCC Al and Al 11 Ce 3 . In the melt pool interiors (low G and high R referred as ES region) a eutectic structure between FCC Al and Al 20 Mn 2 Ce is observed. During subsequent heat treatments, the MPB and ES regions transform into different sets of microstructures. In the MPB region, a fine globular microstructure containing FCC Al, Al 11 Ce 3 , Al 6 Mn, and Al 12 Mn results from the decomposition of Al 20 Mn 2 Ce. In the ES region a faceted Al 51 Mn 7 Ce 4 plate phase results from the decomposition of Al 20 Mn 2 Ce. The formation of the Al 51 Mn 7 Ce 4 phase within the eutectic microstructure at the boundaries of FCC Al and Al 20 Mn 2 Ce has not been reported in the literature. Further, these two distinct phase transformation pathways are rationalized based on the role of driving force on the nucleation of (Al 6 Mn) and/or metastable intermetallic (Al 51 Mn 7 Ce 4 ) phases at the interface of aluminum (FCC) and the non-equilibrium intermetallic (Al 20 Mn 2 Ce) phases.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Assessment of diffuse-interface methods for compressible multiphase fluid flows and elastic-plastic deformation in solids

This work describes three diffuse-interface methods for the simulation of immiscible, compressible multiphase fluid flows and elastic-plastic deformation in solids. The first method is the localized-artificial-diffusivity approach of Cook, Subramaniam et al., and Adler and Lele, in which artificial diffusion terms are added to the individual phase mass fraction transport equations and are coupled with the other conservation equations. The second method is the gradient-form approach that is based on the quasi-conservative method of Shukla et al., in which the diffusion and sharpening terms (together called regularization terms) are added to the individual phase volume fraction transport equations and are coupled with the other conservation equations. The third approach is the divergence-form approach that is based on the fully conservative method of Jain et al., in which the regularization terms are added to the individual phase volume fraction transport equations and are coupled with the other conservation equations. In the present study, all three diffuse-interface methods are used in conjunction with a four-equation, multicomponent mixture model, in which pressure and temperature equilibria are assumed among the various phases. The primary objective of this work is to compare these three methods in terms of their ability to: maintain constant interface thickness throughout the simulation; conserve mass, momentum, and energy; and maintain accurate interface shape for long-time integration. The second objective of this work is to consistently extend these methods to model interfaces between solid materials with strength. To assess and compare the methods, they are used to simulate a wide variety of problems, including (1) advection of an air bubble in water, (2) shock interaction with a helium bubble in air, (3) shock interaction and the collapse of an air bubble in water, and (4) Richtmyer–Meshkov instability of a copper–aluminum interface. The current work focuses on comparing these methods in the limit of relatively coarse grid resolution, which illustrates the true performance of these methods. In conclusion, this is because it is rarely practical to use hundreds of grid points to resolve a single bubble or drop in large-scale simulations of engineering interest.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

A new re-redistribution scheme for weighted state redistribution with adaptive mesh refinement

State redistribution (SRD) is a recently developed technique for stabilizing cut cells that result from finite-volume embedded boundary methods. SRD has been successfully applied to a variety of compressible and incompressible flow problems. When used in conjunction with adaptive mesh refinement (AMR), additional steps are needed to preserve the accuracy and conservation properties of the solution if the embedded boundary is not restricted to a single level of the mesh hierarchy. In this work, we extend the weighted state redistribution algorithm to cases where cut cells live at or near a coarse-fine interface within the domain. Here, we present numerical results that demonstrate that the algorithm is conservative when the coarse-fine interface intersects the embedded boundary. Additionally we compare the numerical solution of the Sod shock tube problem in an inclined cylinder with the analytic solution, and we compare the simulation of a shock hitting a cylindrical obstacle with experimental data. Finally we demonstrate the methodology for simulation of the multicomponent compressible Navier-Stokes equations in a piston-bowl geometry, and discuss the computational efficiency gained by not requiring the entire embedded boundary to be defined at the finest level.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

GrainGNN: A dynamic graph neural network for predicting 3D grain microstructure

We propose GrainGNN, a surrogate model for the evolution of polycrystalline grain structure under rapid solidification conditions in metal additive manufacturing. High fidelity simulations of solidification microstructures are typically performed using multicomponent partial differential equations (PDEs) with moving interfaces. The inherent randomness of the PDE initial conditions (grain seeds) necessitates ensemble simulations to predict microstructure statistics, e.g., grain size, aspect ratio, and crystallographic orientation. Here, currently such ensemble simulations are prohibitively expensive and surrogates are necessary.In GrainGNN, we use a dynamic graph to represent interface motion and topological changes due to grain coarsening. We use a reduced representation of the microstructure using hand-crafted features; we combine pattern finding and altering graph algorithms with two neural networks, a classifier (for topological changes) and a regressor (for interface motion). Both networks have an encoder-decoder architecture; the encoder has a multi-layer transformer long-short-term-memory architecture; the decoder is a single layer perceptron.We evaluate GrainGNN by comparing it to high-fidelity phase field simulations for in-distribution and out-of-distribution grain configurations for solidification under laser power bed fusion conditions. GrainGNN results in 80%–90% pointwise accuracy; and nearly identical distributions of scalar quantities of interest (QoI) between phase field and GrainGNN simulations compared using Kolmogorov-Smirnov test. GrainGNN's inference speedup (PyTorch on single x86 CPU) over a high-fidelity phase field simulation (CUDA on a single NVIDIA A100 GPU) is 150×–2000× for 100-initial grain problem. Further, using GrainGNN, we model the formation of 11,600 grains in 220 seconds on a single CPU core.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Phosphate adsorption kinetics and equilibria on natural iron and manganese oxide composites

We report that it is well known that phosphate retention in soils and sediments is strongly influenced by binding to secondary iron oxides, there have been relatively few studies examining its adsorption/desorption behavior on multicomponent particles of realistic natural complexity. In this study, natural Mn-rich limonite (LM), was used to prepare naturally complex Fe- and Mn-oxide composite materials to examine phosphate adsorption/desorption. To clarify the role of the Mn-oxides, results for the LM sample were compared to those for an acid treated version (LAT), in which the acid-extractable Mn-oxide fraction has been selectively eliminated while leaving the Fe-oxide fraction intact. The saturated adsorption capacity on LAT was almost double that on LM, suggesting that phosphate adsorption to the iron oxides is strongly occluded by the Mn-oxide fraction. This result is reinforced by the comparing the pH dependence and fits to adsorption isotherms, and by desorption experiments and STEM-EDS mapping showing that phosphate loading on Mn-oxides was limited. Hence, although the collective results confirm that phosphate uptake and strong binding is selectively controlled by the Fe-oxide fraction, our study reveals that the Mn-oxide fraction strongly interferes with this process. Therefore, phosphate uptake behavior on metal oxides cannot be predicted solely on the basis of the Fe-oxide fraction present, but instead must take into account the deleterious impacts of other intimately associated phases. For co-diagenetic Fe/Mn-oxide composites in particular, Mn-oxides appear to severely limit phosphate uptake on the Fe-oxide fraction, either by hindering access to binding sites on the Fe-oxide or by lowering their affinity for P.

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

In situ measurement of three-dimensional intergranular stress localizations and grain yielding under elastoplastic axial-torsional loading

The three-dimensional grain-averaged response of solid bar samples under non-proportional (NP) elastoplastic axial-torsional loading was investigated using in situ high energy diffraction microscopy (HEDM) and companion crystal plasticity finite element (CPFE) modeling. Important stress metrics including applied shear (σ θZ ) and axial (σ ZZ ) stress tensor components, stress and stress deviator tensor invariants (I 1 , J 2 , and J 3 ), von Mises equivalent stress (σ$^{grain}_{VM}$), maximum resolved shear stress (mRSS), stress triaxiality (η), and lode angle parameter ($\barθ$) values were tracked for ~300 grains under two different loading conditions: (1) Torsion-dominated loading (low NP) and (2) Tension-torsion loading (high NP) in equiatomic NiCoCr, a representative multicomponent face-centered cubic (FCC) superalloy. Overall, significant stress localizations existed within both samples as evidenced by the radial dependence of grain-resolved σ θZ , σ$^{grain}_{VM}$, and J 2 ; by comparison, I 1 , J 3 , η, and $\barθ$ metrics did not show discernible trends within the volume. These stress localizations reveal a complex interplay between axial and shear stress components (e.g., stress coupling) resulting in grain yielding near the sample surface largely driven by shear stress, whereas internal grain yielding was largely accommodated by axial stress. Grain-resolved stress localization trends were described well by the CPFE model, although some discrepancies in magnitude occurred, particularly for volumetric stress metrics (I 1 and η) due to initial type II residual stress distributions. The superposition of initial residual stress states onto CPFE grain-resolved data significantly improved model accuracy for η. This suggests that residual stresses more strongly influence the simulation of volumetric rather than deviatoric (yield) stress metrics.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Augmented Adam-Gibbs model for glass melt viscosity and configuration entropy as functions of temperature and composition

As the temperature of glass melt increases, its structure approaches the state of a simple liquid while the configuration entropy approaches a maximum value. We describe this gradual change using a power law function of inverse temperature. The Adam-Gibbs model for glass viscosity as a function of temperature and glass composition augmented in this way is greatly simplified when applied to common glass families occupying moderate composition regions, such as float glass or nuclear waste glasses, on which properties can be approximated as linear functions of composition. The parsimonious model thus obtained is preferable for use in optimizing glass formulation and mathematical modeling of glass melting and forming. For multicomponent glasses with N viscosity-affecting components, the augmented Adam-Gibbs model requires 2N + 3 adjustable parameters. The model efficacy is demonstrated by fitting the model to a viscosity-temperature-composition dataset for low-activity nuclear waste glasses.

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

Examining phase separation and crystallization in glasses with X-ray nano-computed tomography

X-ray nano-computed tomography (nano-CT) is a powerful technique to characterize and visualize 3 dimensional (3D) phenomena in complex glasses at the nanoscale. This technique can offer a unique opportunity to explore the intricate morphology of multicomponent glass and glass-ceramic samples, due to its low X-ray energy and high spatial resolution (down to 50 nm). In the current demonstration paper, nano-CT provided insight into crystallization and phase separation, as well as a variety of other phenomena, including corrosion, fracture, and porosity. At scales ranging from 100 to 20 μm, the microstructures of phase separation in borosilicate, complex silicate, and chalcogenide glass compositions were examined. In addition, crystallites present in simulant nuclear waste glasses, volcanic glasses, and other related samples were explored. Nano-CT analysis can add to the understanding of (1) the formation processes, (2) distributions, interactions, and compositions of phases, and (3) the mechanisms, structures, and growth of crystallization and phase separation. Nano-CT offers wide potential in the field of glass science, especially when (1) common characterization methods are insufficient, (2) simple sample preparation is required, (3) 3D rendering is needed, or (4) compelling images of small, complex features are desired.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Viscosity anomaly of a metallic glass-forming liquid under high pressure

Viscosity, as a critical property closely associated with the glass-forming ability of a liquid, has been extensively studied with varying temperatures. However, its pressure dependence has not been well explored yet due to experimental difficulties. Here, we measured the viscosity of a metallic glass-forming liquid, Zr 46 Cu 37.6 Ag 8.4 Al 8 , at pressures up to 6.1 GPa above the melting points by falling-sphere viscometry with ultrafast synchrotron x-ray imaging. Overall, the viscosity increases with pressure, while surprisingly, there is an abrupt drop between 3.2 GPa and 3.7 GPa, indicating the possible existence of a pressure-induced liquid-to-liquid transition. Pressure could change the short- and medium-range orders in the multicomponent glass-forming liquid, as suggested by the different crystalline outcomes after cooling to room temperature at high- and low-pressure ranges. In conclusion, our work extended the viscosity investigation of metallic glass-forming liquids to the high-pressure regime, which will expand our understanding of liquid-liquid transitions and metallic glass formation.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Separation of volatile organic contaminants from water using a direct-contact dehumidifier: An experimental study and modeling

A packed-bed, direct-contact dehumidifier (PB-DCD) has been developed to selectively condense water from a gaseous mixture of organic contaminants. In this process, hot air is humidified by mixing with wastewater vapor in a static mixer and is dehumidified by subcooled fresh water passing through the PB-DCD in a counter-current configuration. The goal of this study is to understand the relationship between various operating parameters and how they affect separation of isopropyl alcohol (IPA) from the gaseous mixture with the purpose of expanding the investigation to other contaminants in the future. The process has been successfully used to separate IPA from wastewater. Different parameters including the concentration of IPA in water, feed temperature, air-to-vapor ratio, and cooling-water-to-vapor mass ratio, as well as residence time were examined to determine their influence on contaminant separation and clean water recovery rate. The Buckingham Pi theorem was applied to reduce the number of studying variables by generating non-dimensional groups. Experiments were performed to determine the relationship between the non-dimensional parameters and model mass transfer in the system. Mass transfer coefficients for multicomponent gaseous mixture condensation in a packed-bed direct contact condenser column are calculated and a modified correlation of mass transfer is developed. The parameters of the mass transfer model are obtained from operating conditions of the PB-DCD, stream conditions and thermodynamic properties. A modified Sherwood correlation is developed to predict IPA separation from a gaseous mixture consisting of air, water, and IPA. The modified Sherwood correlation is applied to successfully predict the Sherwood number with a mean absolute error equal to 6.6 %. According to the experimental results, the maximum IPA separation achieved was 76.3 %, which corresponds to a water recovery rate of 58.7 %.

42 ENGINEERING↗
DOE OSTI2023

Scalable aqueous-phase fabrication of reduced graphene oxide nanofiltration membranes by an integrated roll-to-roll (R2R) process

Economic fractionation of multicomponent biomass-derived streams is a key challenge in biobased fuels and chemicals production. Kraft black liquor (BL) is generated at ~1 billion tons/yr globally from biomass pulping and contains ~15wt% total solids including lignin, hemicellulose fragments, and inorganics. Membrane-based BL concentration is attractive but challenged by low solute rejections and poor stability in BL, which combines alkaline pH (~13), high dissolved solids content (15+ wt%) and high temperature (70–85°C). Here we discuss the scaled fabrication and characterization of reduced graphene oxide (rGO) nanofiltration membranes by slot die coating on a roll-to-roll (R2R) with integrated vacuum filtration. Here, we obtain high-quality 90×30 cm rGO membranes (referred to as “R2R-rGO membranes” in this work) with ~100 nm average thickness supported on porous poly(ethersulfone) (PES) sheets, and an effective area of 2700cm 2 . We present characterizations of their microstructure and uniformity as mapped over large areas by scanning electron microscopy (SEM), ellipsometry, adhesion test and X-ray photoelectron spectroscopy (XPS). The constructed ellipsometry model for the rGO membrane/PES substrate composite was successfully applied to 48 different locations of the R2R-rGO membranes. The R2R-rGO membranes showed >98% lignin rejection with a steady state flux of ~12 LMH at 50bar in a 15.7wt% total solids (TS) kraft BL stream at 72°C and 2–3 gal/min crossflow rates. Slot die coating on a R2R platform with integrated vacuum filtration can enable rapid fabrication of high-quality GO membranes at scale without organic solvents or volatile organic compounds.

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

Amphiphile conformation impacts aggregate morphology and solution structure across multiple lengthscales

Although the self-assembly of amphiphiles is well-studied in aqueous solutions, much less is understood about the fundamental driving forces and structure property relationships in non-polar media. In recent work [Journal of Physical Chemistry B, 2020, 124, 10822.] the authors have studied a series of malonomide-based amphiphiles that are relevant to liquid-liquid extraction. That work demonstrated that aggregation is largely driven by local dipole-dipole interactions between molecules. Here, we build upon this observation to develop a more detailed understanding of how the balance of dipole-dipole interactions (controlled by conformation) and molecular architecture influences the morphology of the aggregates across lengthscales. Using constrained molecular dynamics about key degrees of freedom, we demonstrate that the conformation of N,N-dimethyl,N,N-dioctylhexylethoxy malonamide (DMDOHEMA) and N,N-dimethyl,N,N-dibutyltetradecyl malonamide (DMDBTDMA) has a significant impact upon self-association - where appropriate conformational sampling is essential. To quantify the aggregate morphology, several graph theoretic and persistent homology based properties are determined. The former examines the patterns of intermolecular interactions within clusters, while the latter examines the 3-dimensional spatial distribution. Based upon these analyses, we find that the morphology of aggregates, particularly at higher malonamide concentration, depends on a balance of dipole alignment and alkyl tail sterics. Dipole alignment encourages linear patterns of the intermolecular interactions within aggregates, while the alkyl tail steric interactions between the malonamide result in noticeably less linear aggregates for DMDOHEMA than DMDBTDMA. This is reflected in the spatial distribution, where more holes or voids exist between extractants within the DMDOHEMA that distribute within the solution in more of a "swiss cheese" arrangement as opposed to the more filamentous distribution of DMBDTDMA. This study links conformation and molecular structure to the morphology of amphiphile assemblies, and serves as a basis for ongoing study of multicomponent amphiphile solutions with polar and other solutes, and how these impact aggregation phenomena.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗