Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “solid dynamics”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Carbohydrate-aromatic interface and molecular architecture of lignocellulose

Abstract Plant cell walls constitute the majority of lignocellulosic biomass and serve as a renewable resource of biomaterials and biofuel. Extensive interactions between polysaccharides and the aromatic polymer lignin make lignocellulose recalcitrant to enzymatic hydrolysis, but this polymer network remains poorly understood. Here we interrogate the nanoscale assembly of lignocellulosic components in plant stems using solid-state nuclear magnetic resonance and dynamic nuclear polarization approaches. We show that the extent of glycan-aromatic association increases sequentially across grasses, hardwoods, and softwoods. Lignin principally packs with the xylan in a non-flat conformation via non-covalent interactions and partially binds the junction of flat-ribbon xylan and cellulose surface as a secondary site. All molecules are homogeneously mixed in softwoods; this unique feature enables water retention even around the hydrophobic aromatics. These findings unveil the principles of polymer interactions underlying the heterogeneous architecture of lignocellulose, which may guide the rational design of more digestible plants and more efficient biomass-conversion pathways.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Tuning the melting point of selected ionic liquids through adjustment of the cation’s dipole moment.

In prior research with thermally robust salts [Cassity et al., PCCP, 2017, 19, 31560] it was noted that an increase in the dipole moment of the cation generally led to a decrease in the melting point. Molecular dynamics simulations of the liquid state revealed that an increased dipole moment reduces cation-cation repulsions through dipole-dipole alignment. This was believed to reduce the liquid phase enthalpy, which would tend to lower the melting point of the IL. Here, we also test this principle by replacing hydrogen atoms with fluorine atoms at selected positions within the cation. This allows us to alter the electrostatics of the cation without substantially affecting the sterics. Furthermore, the strength of the dipole moment can be controlled by choosing different positions within the cation for replacement. We studied variants of four different parent cations paired with bistriflimide and determined their melting points, and enthalpies and entropies of fusion through DSC experiments. The decreases in the melting point were determined to be enthalpically driven. We found that the dipole moment of the cation, as determined by quantum chemical calculations, is inversely correlated with the melting point of the given compound. Molecular dynamics simulations of the crystalline and solid states of two isomers showed differences in their enthalpies of fusion that closely matched those seen experimentally. Moreover, this reduction in the enthalpy of fusion was determined to be caused by an increase in the enthalpy of the crystalline state. We provide evidence that dipole-dipole interactions between cations leads to the formation of cationic domains in the crystalline state. These cationic associations partially block favourable cation–anion interactions, which are recovered upon melting. If, however, the dipole- dipole interactions between cations is too strong they have a tendency to form glasses. This study provides a design rule for lowering the melting point of structurally similar ILs by altering their dipole moment.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

Incorporating the molecular-scale into a hydrodynamic description of confined aqueous systems

Hydrodynamics provides a continuum-level description of fluid motion, but its applicability at the nanoscale becomes uncertain due to the emerging importance of molecular-level effects such as spatial heterogeneity. Hydrodynamic boundary conditions that incorporate molecular details allow us to partition the system into a near-wall region and a bulk fluid region. We identify a hydrodynamic wall located inside the fluid that determines where slip begins. By extending the hydrodynamic wall with the slip length, the position of the extrapolated wall is established. This offers a unified description of both slip and stagnant flow behaviors, with wall hydrophobicity characterized by the relative location of the extrapolated wall with respect to the physical wall. Employing this concept in analyses of equilibrium molecular dynamics (MD) and non-equilibrium MD simulations of Couette and Poiseuille flows, our results demonstrate consistency between equilibrium and non-equilibrium approaches across different flow types and confinement levels. This demonstrates the robust nature of linear response theory. We then explore the effects of fluid-wall and bulk fluid interactions on the hydrodynamic properties. These findings enhance the effectiveness of molecular-based simulations for investigating complex confined systems in nanofluidics, biology, and colloidal science, offering a complementary molecular-scale perspective to traditional continuum approaches.

flow simulations↗
DOE OSTI2022

Phase diagram of the TIP4P/Ice water model by enhanced sampling simulations

In this work, we studied the phase diagram for the TIP4P/Ice water model using enhanced sampling molecular dynamics simulations. Our approach is based on the calculation of ice–liquid free energy differences from biased coexistence simulations that reversibly sample the melting and growth of layers of ice. We computed a total of 19 melting points for five different ice polymorphs, which are in excellent agreement with the melting lines obtained from the integration of the Clausius–Clapeyron equation. For proton-ordered and fully proton-disordered ice phases, the results are in very good agreement with previous calculations based on thermodynamic integration. For the partially proton-disordered ice III, we find a large increase in stability that is in line with previous observations using direct coexistence simulations for the TIP4P/2005 model. This issue highlights the robustness of the approach employed here for ice polymorphs with diverse degrees of proton disorder. Our approach is general and can be applied to the calculation of other complex phase diagrams.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Nanoparticle surfactants and structured liquids

Materials are usually classified as solids or liquids, based on their structural stability, dynamic response, and rheological properties. Structured liquid, a new state of matter, has attracted much attention in recent years. Different with either solid or liquid, structured liquid combines the desirable characteristics of fluids with the structural stability of a solid, showing a myriad of potential applications in encapsulation, biphasic reactors, and programmable liquid constructs. Here, a brief review is given, by introducing a new strategy to structure liquids based on the formation, assembly, and jamming of nanoparticle surfactants (NPSs) at liquid-liquid interfaces. Furthermore, the interfacial packing of the NPSs can be effectively manipulated using external triggers, endowing the structured liquids with adaptiveness and responsiveness to changes in their external environment.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Simulations of flaming combustion and flaming-to-smoldering transition in wildland fire spread at flame scale

Our objective in the present study is to provide basic insights into the coupling between external-gas and solid biomass vegetation processes that control the dynamics of flame spread in wildland fire problems. Here we focus on a modeling approach that resolves processes occurring at vegetation and flame scales, i.e., the formation of flammable vapors due to the thermal degradation of the solid biomass, the subsequent combustion in ambient air, the thermal feedback to the biomass through radiative and convective heat transfer, and the possible transition from flaming combustion (taking place outside of the solid biomass) to smoldering combustion (taking place inside the solid biomass). The capability uses a multiphase combustion framework and treats external-gas processes through a Large Eddy Simulation solver and solid biomass processes through a discrete particle model. The discrete particle model adopts a one-dimensional porous medium formulation, includes descriptions of drying, thermal pyrolysis, oxidative pyrolysis, and char oxidation, as well as a description of the external-gas-to-solid-biomass diffusion of oxygen mass; the discrete particle model thereby provides a treatment of in-depth oxidative processes and allows the simulation of smoldering combustion. The modeling capability is applied to the simulation of fire spread across a surrogate biomass vegetation bed corresponding to a discrete array of cylindrical-shaped, vertically-oriented, pine wood sticks, characterized by a monomodal size distribution, in horizontal flat terrain and under wind-aided conditions. The numerical results demonstrate that the model can simulate successful flaming-to-smoldering transition followed by complete biomass consumption.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2026

Site-Selective Modification of Lanthanum Oxychloride to Modulate Halide-Ion Conduction

Design principles for solid-state halide-ion conduction remain poorly defined despite the increasing importance of halide ions as charge carriers in a variety of energy storage and electrochemical computing technologies. Here, we employ a siteselective modification strategy in which aliovalent cations are preferentially introduced at the La 3+ crystallographic site of LaOCl in the 2c Wyckoff position, enabling controlled generation of chloride vacancies and modification of lattice dynamics to enhance chloride-ion conductivity. Aliovalent substitution of La 3+ with Mg 2+ , Ca 2+ , and Sr 2+ generates charge-compensating Cl vacancies while preserving the matlockite crystal structure. X-ray excited optical luminescence measurements with Dy 3+ as a reporter chromophore evidence vacancy-derived midgap electronic states and an extended energy range of radiation-less Auger electron emission corresponding to substantial modification of electronic structure and local electrostatic potentials. Ca alloying at 8−10 at. % increases the chloride-ion conductivity by three- to 4 orders of magnitude as compared to unalloyed LaOCl, whereas comparable amounts of Sr- and Mg-alloying in LaOCl imbue less pronounced conductivity enhancements. Temperature-dependent Raman spectroscopy measurements reveal that Ca- and Sr-alloying substantially soften the La−Cl sublattice and yield a more compliant crystal lattice that can deform to accommodate Cl-ion migration. Structure solutions derived from Rietveld refinements to powder Xray diffraction reveal larger O−La−Cl bond-angle deviations and enhanced out-of-plane cation displacements for Ca- and Sr-alloyed compositions as compared to Mg-alloyed LaOCl. Such local distortions enhance chloride-ion mobility by reshaping and flattening vacancy migration energy landscapes and by modulating lattice dynamics governing anion conduction. We further illustrate that coalloying of Ca with Mg and Sr induces a nonmonotonic conductivity−defect stoichiometry relationship that can be rationalized based on cooperative interactions. Together, these results establish site-selective aliovalent alloying of LaOCl as an effective route to halide-ion solid electrolytes and provide broadly generalizable design principles for site-selective modification to induce vacancy formation and lattice softening to engender facile anion transport

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

Anaerobic Dynamic Membrane Bioreactor Development to Facilitate Organic Waste Conversion to Medium-Chain Carboxylic Acids and Their Downstream Recovery

Platform chemicals such as medium-chain carboxylic acids (MCCAs) can be produced from organic waste streams via chain elongation in anaerobic mixed-culture bioreactors. A product recovery system is needed to collect MCCAs from the bioreactor effluent. Membrane-based liquid–liquid extraction, the most commonly used product recovery approach, requires suspended solids removal from the bioreactor effluent to avoid membrane fouling. An anaerobic dynamic membrane bioreactor (AnDMBR) was developed to evaluate MCCA production from brewery and prefermented food waste and to produce a permeate with low suspended solids to facilitate downstream product recovery. The AnDMBR employed an inexpensive stainless-steel mesh as the support material for the development of a biofilm or dynamic membrane, which was responsible for solids–liquid separation. The AnDMBR produced a permeate quality with an average total suspended solids (TSS) concentration of 0.12 g L –1 , while the average bioreactor TSS concentration was two orders of magnitude higher (21.6 g L –1 ). A maximum solids removal efficiency of ≥99% was achieved and good permeate quality was sustained for over 200 days without fouling control or cleaning the support material. In addition to solids–liquid separation, the dynamic membrane was responsible for a substantial fraction of the biological activity of the AnDMBR. The relative activity of Clostridiales, as determined by 16S rRNA sequencing, correlated with MCCA production and was higher in the dynamic membrane (20.0 ± 4.9%) than in the suspended biomass (5.2 ± 2.7%) in the bioreactor. This observation was consistent with MCCA production data as the permeate MCCA concentrations were significantly (p = 8.2 × 10 –5 ) higher than that in the bioreactor, suggesting that the dynamic membrane biofilm contributed to chain elongation.

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

The ultrafast Kerr effect in anisotropic and dispersive media

The ultrafast optical Kerr effect (OKE) is widely used to investigate the structural dynamics and interactions of liquids, solutions, and solids by observing their intrinsic nonlinear temporal responses through nearly collinear four-wave mixing. Non-degenerate mixing schemes allow for background free detection and can provide information on the interplay between a material’s internal degrees of freedom. Here, we show a source of temporal dynamics in the OKE signal that is not reflective of the internal degrees of freedom but arises from a group index and momentum mismatch. It is observed in two-color experiments on condensed media with sizable spectral dispersion, a common property near an optical resonance. In particular, birefringence in crystalline solids is able to entirely change the character of the OKE signal via the off-diagonal tensor elements of the nonlinear susceptibility. We develop a detailed description of the phase-mismatched ultrafast OKE and show how to extract quantitative information on the spectrally resolved birefringence and group index from time-resolved experiments in one and two dimensions.

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

Multistage nucleation pathway in LiF molten salt mirrors the crystal–melt interface structure

Despite over a century of studies, fundamental questions remain about the processes governing crystal nucleation from melts or solutions. Research over the past three decades has presented mounting evidence for kinetic pathways of crystal nucleation that are more complex than envisioned by the simplest forms of classical theory. Such observations have been presented for colloidal and elemental systems with covalent and metallic bonding. Despite the technological and geochemical importance of molten salts, similar studies for these ionically bonded systems are currently lacking. Here we develop a machine learning interatomic potential for a model ionic system: LiF. The potential features quantum-level accuracy for both liquid and multiple solid polymorphs over wide temperature and pressure ranges and accurately reproduces experimentally measured properties. Thanks to the efficiency of the potential, which enables microsecond-scale molecular dynamics simulations, induction times for nucleation of LiF solids from their melts are computed over a range of undercoolings. With the aid of a set of robust local order parameters established here, the simulations reveal that homogeneous crystal nucleation in undercooled melts preferentially initiates from liquid regions showing slow dynamics and high bond orientational order simultaneously, and the second-shell order of both precritical nuclei and the surface of postcritical nuclei is dominated by hexagonal close packing and body-centered cubic local structure, even though the nucleus core is dominated by face-centered cubic structure corresponding to the stable rocksalt crystal structure. Finally, we establish a connection between the crystallization pathway and the equilibrium crystal-melt interface structure.

Applied Physical Sciences↗
DOE OSTI2021

One-directional polarization transport in electron/nuclear spin chains with loss and gain

Understanding the joint dynamics of electron and nuclear spins is central to core concepts in solid-state magnetic resonance such as spin-lattice relaxation and dynamic nuclear polarization, but a generalization that capitalizes on competing polarization loss and gain channels is still lacking. Here, we study theoretically the non-Hermitian dynamics of hybrid electron/nuclear spin systems in the simultaneous presence of electron spin pumping and spin-lattice relaxation. Focusing on periodic, one-dimensional chains, we find that, by adjusting the electron spin pumping to a critical level, it is possible to steer the flow of nuclear polarization to create site-dependent distributions where either end of the array polarizes in opposite ways, irrespective of the initial state. By contrast, we show that ring-like patterns—where the limit nuclear polarization is uniform—exhibit a non-decaying, externally driven nuclear spin current. Interestingly, cyclic magnetic field modulation can render these processes largely robust to defects in the chain, a response featuring some interesting similarities—and differences—with recent findings in other non-Hermitian physical platforms.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Heat transfer in directly-irradiated high-temperature solid–gas flows laden with polydisperse particles

Heat transfer in directly-irradiated high-temperature solid–gas flows laden with polydisperse particles is investigated using a novel transient three-dimensional computational fluid dynamics model. The model couples particle–gas hydrodynamics of solid–gas flows laden with polydisperse particles, radiative heat transfer in non-grey absorbing, emitting and anisotropically-scattering multi-component participating media, conduction heat transfer in the gas phase, and interfacial convection heat transfer. The multiphase particle-in-cell method is used to predict high-fidelity solid–gas flow characteristics, such as the local discrete particle size distribution, with increased computational efficiency by combining the advantages of both Eulerian and Lagrangian methods. The multi-component radiative transfer model is implemented using an advanced collision-based Monte Carlo ray-tracing method. The number of the prescribed discrete particle components is found to be the key parameter affecting the computational accuracy and efficiency, which primarily depends on the size distribution of the particles. For the model particle–gas flow featuring free-falling Gamma-distributed ceramic particles exposed to concentrated solar irradiation, the particle volume fraction, radiative, fluid flow and thermal characteristics appear to converge with the increasing number of the discrete particle components. Five particle components are sufficient to obtain physically meaningful results. A further increase in the number of the particle components only slightly increases the accuracy of the numerical predictions at the expense of a rapidly increasing computational time. For five particle components, the particle vertical velocity at the receiver exit for particles with the diameter of 43.4 μm is 57% of that for the particles with the diameter of 202.8 μm. The temperatures of these two particle components increase from the initial ambient values by factors of 2 and 1.2, respectively, during the simulation time. In conclusion, the model developed allows for increased fidelity of particle–gas flow simulations with significant radiative effects.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

A hybrid smoothed-particle hydrodynamics model of oxide skins on molten aluminum

A computational model of aluminum melting is proposed which captures both the thermal fluid-solid phase transition and the mechanical effects of oxidation. The model hybridizes ideas from smoothed particle hydrodynamics and bonded particle models to simulate both hydrodynamic flows and solid elasticity. Oxidation is represented by dynamically adding and deleting spring-like bonds between surface fluid particles to represent the formation and rupture of the oxide skin. Various complex systems are simulated to demonstrate the adaptability of the method and to illustrate the significant impact of skin properties on material flow. As a result, initial comparison to experiments of a melting aluminum cantilever highlights that the computational model can reproduce key qualitative features of aluminum relocation.

Clemmer, Joel Thomas↗
DOE OSTI2021

Revealing the Configuration and Conformation of Surface Organometallic Catalysts with DNP-Enhanced NMR

Although single-site, supported organometallic catalysts were designed with homogeneity in mind, little is strictly known about their atomic-level structure and uniformity. This is in large part due to the inability of conventional characterization tools to provide structural information with adequate resolution and range for this purpose. In this work, we show that dynamic nuclear polarization (DNP)-enhanced solid-state nuclear magnetic resonance (NMR) enables the measurement of distances between the surface and individual carbons, making it possible to orient complexes and molecules on surfaces. We use this approach to determine the orientation and configuration of naturally 13 C-abundant acetate and supported organoiridium(III) pincer complex, both supported on γ-Al 2 O 3 . By combining solid-state NMR and extended X-ray absorption fine structure spectroscopy (EXAFS) experiments with periodic density functional theory (DFT) calculations, we are able to determine the three-dimensional arrangement of ligands around the metal center and thereby identify the structure of the supported complex at a resolution reminiscent of that associated with single-crystal diffractometry.

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

Evolution and Degradation Patterns of Electrochemical Cells Based on the Analysis of Interfacial Phenomena at Li Metal Anode/Electrolyte Interfaces

In this work, we report the results of a theoretical–computational analysis of the solid electrolyte interphase (SEI) growth and degradation dynamics occurring in lithium metal batteries during cycling. We use ab initio-kinetic Monte Carlo simulations to generate a synthetic data set, which is analyzed by machine learning methods. We aim to determine: (i) how modifications in interfacial interaction energies between solid electrolyte interphase (SEI) blocks and between Li ions and SEI facets impact the Coulombic efficiency (CE) of the battery and (ii) what factors, including reactions, microscopic transport, and other interfacial events, may lead to cell performance “failure” during prolonged charge and discharge cycles, signaled as a sharp decay in the CE over cycling. The demonstration of our approach is done on a cell including a Li metal surface interfacing with a previously introduced state-of-the-art electrolyte, and the idea can be applied to any electrochemical system. Outcomes include the identification of the leading chemical, physical, and structural variables causing cell failure and relating them to the electrolyte formulation, thus paving the way to future more refined analysis and electrolyte design.

batteries↗
DOE OSTI2023

Supported Electrophilic Organoruthenium Catalyst for the Hydrosilylation of Olefins

A series of supported electrophilic organoruthenium complexes has been synthesized via surface organometallic chemistry (SOMC) techniques and applied to the selective hydrosilylation of olefins. The air-sensitive 16e - complex Cp*RuMes(PCy 3 ) (1) (Cp* = pentamethylcyclopentadienyl, Mes = mesityl) was synthesized by the treatment of Cp*RuCl(PCy 3 ) with mesityl Grignard MesMgBr. This species was chemisorbed onto sulfated zirconia SO 4 /ZrO 2 , but the resulting material was inactive toward cyclohexene hydrosilylation with phenylsilane. Instead, Cp*RuMes(PCy 3 ) was treated with phenylsilane (PhSiH 3 ) to provide a ruthenium disilyl hydride complex Cp*RuH(SiH 2 Ph) 2 (PCy 3 ) (3), which was fully characterized by NMR spectroscopy and single-crystal X-ray diffraction. Grafting this species onto SO 4 /ZrO 2 resulted in the formation of phenylsilane along with the surface electrophilic species [Cp*RuH(R)(X-SiHPh)(PCy 3 )] (R = H, O 3 S-O or O 3 Zr-O; 4a, 4b, X = O 3 S-O, and O 3 Zr-O, respectively) as the major species. Material 4 was characterized via a combination of spectroscopic techniques including dynamic nuclear polarization (DNP)-enhanced solid-state NMR spectroscopy, diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS), X-ray absorption spectroscopy (XAS), and density function theory (DFT) calculations. Further, capping the remaining acid sites on 4 with Me 3 Si-SiMe 3 provides 5, which significantly reduces side reactions, such as olefin isomerization and silane redistribution. Catalyst 5 is a highly robust and selective hydrosilylation catalyst and can be recycled up to 5 times without significant diminishment of activity. Exclusive anti-Markovnikov regiochemistry, cis-addition selectivity, and the inactivity of secondary and tertiary silanes provide support for the proposed Glaser-Tilley mechanism involving cationic ruthenium silylene species analogous to homogeneous systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗