Search NASA⌕ Search

SEARCH · Search NASA

Results for “bulk viscosity”

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

Atmospheric clusters to nanoparticles: Recent progress and challenges in closing the gap in chemical composition

Aerosol nanoclusters (AN), defined here as molecular aggregates suspended in a gas with dimensions between 2 and 10 nm, are the link between substances that we think of as molecules, or the “gas phase,” and those that we consider as particles, or the “condensed phase.” The ability to measure and model the physical and chemical properties of size-resolved AN, which at present is rudimentary at best, is crucial for understanding how particles form and evolve in a number of environments that are natural or influenced by human activities. This review describes the current state-of-the-art for measuring and modeling the size-resolved composition of atmospheric AN. We focus specifically on instruments, many relying on mass spectrometry, that show promise for closing the measurement gap under atmospherically relevant conditions by increasing the size of measurable gas-phase clusters (bottom-up approaches) and by decreasing the size of measurable nanoparticles (top-down approaches). Theoretical methods for predicting AN composition have similarly relied on bottom-up approaches that extend the accuracy of quantum chemistry calculations to larger molecular systems, as well as top-down approaches that correct bulk composition aerosol models for size-dependent properties such as viscosity and volatility. Finally, current measurement and modeling challenges that must be overcome in order to close the gap are discussed.

42 ENGINEERING↗
DOE OSTI2020

Shear-Induced Gelation of Self-Yielding Active Networks

An enticing feature of active materials is the possibility of controlling macroscale rheological properties through the activity of the microscopic constituents. Using a unique combination of microscopy and rheology we study three dimensional microtubule-based active materials whose autonomous flows are powered by a continually rearranging connected network. We quantify the relationship between the microscopic dynamics and the bulk mechanical properties of these nonequilibrium networks. Experiments reveal a surprising nonmonotonic viscosity that strongly depends on the relative magnitude of the rate of internally generated activity and the externally applied shear. Here, a simple two-state mechanical model that accounts for both the solidlike and yielded fluidlike elements of the network accurately describes the rheological measurements.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Transport Phenomena in Low Temperature Lithium-Ion Battery Electrolytes

Lithium-ion batteries face low temperature performance issues, limiting the adoption of technologies ranging from electric vehicles to stationary grid storage. This problem is thought to be exacerbated by slow transport within the electrolyte, which in turn may be influenced by ion association, solvent viscosity, and cation transference number. How these factors collectively impact low temperature transport phenomena, however, remains poorly understood. Here we show using all-atom classical molecular dynamics (MD) simulations that the dominant factor influencing low temperature transport in LP57 (1 M LiPF 6 in 3:7 ethylene carbonate (EC)/ethyl methyl carbonate (EMC)) is solvent viscosity, rather than ion aggregation or cation transference number. We find that ion association decreases with decreasing temperature, while the cation transference number is positive and roughly independent of temperature. In an effort to improve low temperature performance, we introduce γ -butyrolactone (GBL) as a low viscosity co-solvent to explore two alternative formulations: 1 M LiPF 6 in 15:15:70 EC/GBL/EMC and 3:7 GBL/EMC. While GBL reduces solution viscosity, its low dielectric constant results in increased ion pairing, yielding neither improved bulk ionic conductivity nor appreciably altered ion transport mechanisms. We expect that these results will enhance understanding of low temperature transport and inform the development of superior electrolytes.

25 ENERGY STORAGE↗
DOE OSTI2021

Incorporation and Assembly of a Light-Emitting Enzymatic Reaction into Model Protein Condensates

Eukaryotic cells partition enzymes and other cellular components into distinct subcellular compartments to generate specialized biochemical niches. A subclass of these compartments form in the absence of lipid membranes, via liquid-liquid phase separation of proteins to form biomolecular condensates or “membraneless organelles” such as nucleoli, stress granules, and P-bodies. Because of their ability to form compartments from simple starting materials, membraneless organelles are an attractive target for engineering new functionalities in both living cells and protocells. In this work, we demonstrate incorporation of novel enzymatic activity in protein coacervates with a light-generating enzyme, NanoLuc, to produce bioluminescence. Using condensates comprised of the disordered RGG domain of C. elegans LAF-1, we show functionalization of condensates with enzymatic activity in vitro and that localization to protein coacervates enhances the assembly and activity of split enzymes. To build condensates that function as light emitting reactors, we designed a NanoLuc enzyme flanked by RGG domains. Furthermore, the resulting condensates concentrated NanoLuc by 10-fold over bulk solution and display significantly increased net reaction rates. We further show that condensate viscosity impacts light emission due to diffusion-limited behavior. By splitting NanoLuc enzyme into its constituent components, we demonstrate that NanoLuc activity can be reconstituted via co-condensation. Further, we demonstrate control of the spatial localization of enzyme within condensates by targettng NanoLuc to the surface of in vitro condensates. Collectively, this work demonstrates that membraneless organelles can be endowed with localized enzymatic activity, and that this activity can be spatially and temporally controlled via enzyme reconstitution and design of protein surfactants.

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

Optimized purification methods for metallic contaminant removal from directly recycled Li-ion battery cathodes

Metallic contaminants pose a significant challenge to the viability of directly recycling Li-ion batteries. To date, few strategies exist to selectively remove metallic impurities from mixtures of shredded end-of-life material (black mass; BM) without concurrently damaging the structure and electrochemical performance of the target active material. We herein present tailored methods to selectively ionize two major contaminants—Al and Cu—while retaining a representative cathode (LiNi 0.33 Mn 0.33 Co 0.33 O 2 ; NMC-111) intact. This BM purification process is conducted at moderate temperatures in a KOH-based solution matrix. We rationally evaluate approaches to increase both the kinetic corrosion rate and the thermodynamic solubility of Al 0 and Cu 0 , and evaluate the impact of these treatment conditions on the structure, chemistry, and electrochemical performance of NMC. Specifically, we explore the impacts of chloride-based salts, a strong chelating agent, elevated temperature, and sonication on the rate and extent of contaminant corrosion, while concurrently evaluating the effects on NMC. The reported BM purification process is then demonstrated on samples of “simulated BM” containing a practically relevant 1 wt% concentration of Al or Cu. Increasing the kinetic energy of the purifying solution matrix through elevated temperature and sonication accelerates the corrosion of metallic Al and Cu, such that ∼100% corrosion of 75 μm Al and Cu particles is achieved within 2.5 hr. Further, we determine that effective mass transport of ionized species critically impacts the efficacy of Cu corrosion, and that saturated Cl – hinders rather than accelerates Cu corrosion by increasing solution viscosity and introducing competitive pathways for Cu surface passivation. The purification conditions do not induce bulk structural damage to NMC, and electrochemical capacity is maintained in half-cell format. Testing in full cells suggests that a limited quantity of residual surface species are present after treatment, which initially disrupt electrochemical behavior at the graphite anode but are subsequently consumed. Process demonstration on simulated BM suggests that contaminated samples—which prior to treatment show catastrophic electrochemical performance—can be recovered to pristine electrochemical capacity. The reported BM purification method offers a compelling and commercially viable solution to address contamination, particularly in the “fine” fraction of BM where contaminant sizes are on the same order of magnitude as NMC and where traditional separation approaches are unfeasible. Thus, this optimized BM purification technique offers a pathway towards viable direct recycling of BM feedstocks that would otherwise be unusable.

25 ENERGY STORAGE↗
DOE OSTI2025

Predicting multi-nodal in-nozzle particle interactions in high-viscosity fluid mediums for acoustophoretic direct-ink writing of line-patterned composites

Patterned functional materials offer improved properties (electrical, thermal, etc.) over their bulk counterparts in many applications, including energy storage, flexible electronics, and sensors. However, manufacturing approaches for patterning materials over large areas with features on the order of hundreds of microns or less are limited. Acoustophoresis, which uses acoustic forces to control particle arrangement in a fluid medium, is a pathway to address this challenge. This process is dependent on particle and fluid properties and enables patterning of a broad range of materials. Herein, a model with experimental validation is presented to demonstrate that acoustophoresis can be combined with direct-ink writing (DIW) to fabricate line patterns over large cm-scale areas. An in-nozzle particle interaction model was developed to investigate the impact of processing conditions on multi-nodal acoustophoretic DIW. The model predicts patterned line widths within a factor of two relative to experimental results for a high viscosity case study. Here, the model was used to investigate the impact of frequency, particle loading, particle radius, and acoustic pressure on line width and patterning time, providing critical feedback regarding the processing conditions suitable for a target application. Model results illustrate that frequency has the greatest impact on line patterns: increasing from 1 to 3 MHz resulted in a greater than 65% reduction in line width and a greater than 85% reduction in patterning time. Additionally, experiments were conducted with an alumina-epoxy ink and a ~21 cm 2 area pattern was rastered in ~5.5 minutes, demonstrating a path towards large-area line-patterned composite fabrication.

25 ENERGY STORAGE↗
DOE OSTI2022

Random Walks and Sticky Surfaces: Single-Molecule Measurements of Solute Diffusion in Ethanol/Water-Filled Anodic Alumina Nanopores

Nanoporous anodic aluminum oxide (AAO) membranes are now being explored for use in advanced chemical separations, including in the dehydration of biofuels such as ethanol. Optimization of membrane performance requires an in-depth understanding of how solvent mixtures and solutes behave under nanoconfinement. In this work, the diffusion of rhodamine B (RhB) dye through 10 nm and 20 nm AAO nanopores filled with a series of ethanol/water mixtures (0 - 33% water) is explored by fluorescence correlation spectroscopy (FCS). RhB was found to diffuse through the pores by two distinct mechanisms with mean diffusion coefficients, $D_f$ and $D_s$, reflecting fast and slow diffusive motions, respectively, with values that differ by nearly two orders of magnitude. Further, both $D_f$ and $D_s$ increased with pore size and were significantly smaller than $D_b$, the RhB diffusion coefficient in bulk liquid. Mean $D_f$ values follow a composition-dependent trend that closely mimics the viscosity dependence of $D_b$. Additional slowing of fast RhB diffusion is attributed to both hydrodynamic drag and electrostatic interactions with the nanopore surface. The mean $D_s$ values exhibit a different trend with increasing water content, revealing an increase in $D_s$ and a decrease in the contributions of slow diffusion to the observed dynamics. The fluorescence time transient data used in the analysis show that the slow diffusion process is strongly hindered and likely involves frequent adsorption of RhB to the pore surfaces. These results provide new insights into the detailed molecular-level mechanisms of mass transport in nanoporous AAO membranes.

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

Enhancing Adaptive Physics Refinement Simulations Through the Addition of Realistic Red Blood Cell Counts

Simulations of cancer cell transport require accurately modeling mm-scale and longer trajectories through a circulatory system containing trillions of deformable red blood cells, whose intercellular interactions require submicron fidelity. Using a hybrid CPU-GPU approach, we extend the advanced physics refinement (APR) method to couple a finely-resolved region of explicitly-modeled red blood cells to a coarsely-resolved bulk fluid domain. We further develop algorithms that: capture the dynamics at the interface of differing viscosities, maintain hematocrit within the cell-filled volume, and move the finely-resolved region and encapsulated cells while tracking an individual cancer cell. Comparison to a fully-resolved fluid-structure interaction model is presented for verification. Finally, we use the advanced APR method to simulate cancer cell transport over a mm-scale distance while maintaining a local region of RBCs, using a fraction of the computational power required to run a fully-resolved model.

Roychowdhury, Sayan↗
DOE OSTI2020

Why Is the Range of Timescale So Wide in Glass-Forming Liquid?

The viscosity and the relaxation time of a glass-forming liquid vary over 15 orders of magnitude before the liquid freezes into a glass. The rate of the change with temperature is characterized by liquid fragility. The mechanism of such a spectacular behavior and the origin of fragility have long been discussed, but it remains unresolved because of the difficulty of carrying out experiments and constructing theories that bridge over a wide timescale from atomic (ps) to bulk (minutes). Through the x-ray diffraction measurement and molecular dynamics simulation for metallic liquids we suggest that large changes in viscosity can be caused by relatively small changes in the structural coherence which characterizes the medium-range order. Here the structural coherence does not imply that of atomic-scale structure, but it relates to the coarse-grained density fluctuations represented by the peaks in the pair-distribution function (PDF) beyond the nearest neighbors. The coherence length is related to fragility and increases with decreasing temperature, and it diverges only at a negative temperature. This analysis is compared with several current theories which predict a phase transition near the glass transition temperature.

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

Pseudo-viscous modeling of transport in dense granular flows for thermal energy storage applications

Dense, granular flows were examined to effectively capture and model bulk viscous properties in thin packed beds. A modified Couette cell with particle image velocimetry was used to experimentally determine pseudo-viscosity properties of four particulate media with varying morphologies: (1) iron oxide-coated SiO 2 particles, (2) CARBOBEAD CP30-60 particles, (3) CARBOBEAD CP40-100 particles, and (4) Al 2 O 3 beads. The pseudo-viscosity functions were fitted using a power law to correlate the measured shear stress as a function of measured shear rate. The pseudo-viscous functions were used as inputs to computation fluid dynamics models for a single-phase viscous fluid to predict granular flow profiles. Steady-state free surface velocity profiles at angular velocities <7 rad/s predicted by the model were in good agreement with the experimental particle image velocimetry measurements, resulting in Pearson correlation coefficients of 0.97 for iron-oxide coated SiO 2 particles and 0.95 for CP30-60 particles. As a result, this alternative approach to measuring pseudo-viscous properties under shearing and modeling bulk transport behavior of granular flow using computation fluid dynamics model offered significant reduction in computational load compared to discrete element methods.

14 SOLAR ENERGY↗
DOE OSTI2024

Laser‐Sintered Silver Metallization for Silicon Heterojunction Photovoltaic Cells

Herein, a novel metallization technique is reported for crystalline silicon heterojunction (SHJ) solar cells in which silver (Ag) fingers are printed on the SHJ substrates by dispensing Ag nanoparticle‐based inks through a needle and then sintered with a continuous‐wave carbon dioxide (CO 2 ) laser. The impact of the Ag ink viscosity on the line quality and the line resistance is investigated on three Ag inks with different viscosities. Increasing ink viscosity yields higher Ag contact heights, larger aspect ratios, and lower line resistance values. The Ag line height increases from less than a micrometer to ≈18.62 ± 3.48 μm with the increasing viscosity. Photoluminescence imaging shows that the low‐resistance Ag metal contacts obtained do not result in any passivation damage of the SHJ substrate. This is because the wavelength of light emitted from the CO 2 laser (i.e., 10.6 μm) leads to optical absorption in the Ag, but this light is effectively transparent to the transparent conductive oxide film, amorphous silicon films, and crystalline silicon substrate. Bulk resistivity values as low as 6.5 μΩ cm are obtained for the laser‐sintered Ag contact and printed using the Ag ink with the highest viscosity in this work.

Mousumi, Jannatul Ferdous↗
DOE OSTI2022

Bioremediation of heavy oily sludge: a microcosms study

Oily sludge is a residue from the petroleum industry composed of a mixture of sand, water, metals, and high content of hydrocarbons (HCs). The heavy oily sludge used in this study originated from Colombian crude oil with high density and low American Petroleum Institute (API) gravity. The residual waste from heavy oil processing was subject to thermal and centrifugal extraction, resulting in heavy oily sludge with very high density and viscosity. Biodegradation of the total petroleum hydrocarbons (TPH) was tested in microcosms using several bioremediation approaches, including: biostimulation with bulking agents and nutrients, the surfactant Tween 80, and bioaugmentation. Select HC degrading bacteria were isolated based on their ability to grow and produce clear zones on different HCs. Degradation of TPH in the microcosms was monitored gravimetrically and with gas chromatography (GC). The TPH removal in all treatments ranged between 2 and 67%, regardless of the addition of microbial consortiums, amendments, or surfactants within the tested parameters. The results of this study demonstrated that bioremediation of heavy oily sludge presents greater challenges to achieve regulatory requirements. Additional physicochemical treatments analysis to remediate this recalcitrant material may be required to achieve a desirable degradation rate.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Connecting particle interactions to agglomerate morphology and rheology of boehmite nanocrystal suspensions

Rheology imposes significant challenges on processing of complex suspensions such as nuclear waste slurries at the Hanford and Savannah River sites. Understanding rheology connecting to microstructures and underlying particle interactions in complex slurries is therefore important for both fundamental knowledge and practical applications. Here, we use suspensions of aluminum oxyhydroxide minerals in the form of boehmite as an analog of the radioactive waste slurry to gain physical insights on the correlation between particle interactions, microstructures, and slurry rheology. Specifically, we use a combination of Couette rheometry and small-angle scattering techniques (independently and simultaneously) to understand how the slurry microstructure changes under flow and how these structural changes manifest themselves in the bulk rheology of the suspensions. Our experiments show that the boehmite slurries are thixotropic, with the rheology and structure of the suspensions changing with increasing exposure to flow. In the slurries, particle aggregates begin as loose, system-spanning clusters, but exposure to moderate shear rates causes the aggregates to irreversibly consolidate into denser clusters of finite size. The microstructural changes directly influence the rheological properties of the slurries such as viscosity and viscoelasticity. More importantly, our study shows that solution pH affects the amount of structural rearrangement and the kinetics of the rearrangement process, with an increase in pH leading to faster and more dramatic changes in the bulk rheology. Such dynamic microstructural changes and resultant rheology are understood via correlations between particle interactions and strength of particle network, coupled with surface chemistry and anisotropic nature of particle interactions. Nearly identical structural changes are also observed in Poiseuille flow geometries, implying that the observed changes are relevant in the pipe flow conditions present during waste processing.

Weston, Javen S.↗
DOE OSTI2021

Ultrasound-Assisted Nonthermal, Nonequilibrium Separation of Organic Molecules from Their Binary Aqueous Solutions: Effect of Solute Properties on Separation

The separation of organic molecules in downstream processes, such as distillation, plays a vital role in biobased chemical production. A comprehensive study was performed here using ultrasonic energy to separate 11 organic molecules from their binary aqueous solutions under nonthermal, nonequilibrium, and phase-change free conditions. A batch sonoseparator with a three-stage mist collection unit allowed direct analysis of ultrasound-generated mist with an HPLC/GC and differentiation and quantification of the vapor and mist. The results showed that the alcohols (1-butanol, ethanol, 1-propanol, and methanol) and acetone were enriched in the ultrasound-generated mist with enrichment ratios of 3.2–5. On the other hand, ethylene glycol, glycerol, γ-valerolactone (GVL), glucose, and sucrose were diluted with dilution ratios in the range of 2.04–15.7 in the mist and vapor generated by ultrasound (or concentrated in the bulk solution). No enrichment nor dilution of acetic acid was observed. The role of various physicochemical parameters such as dynamic viscosity, surface tension, Henry’s law of solubility constant, vapor pressure, and octanol–water partition coefficient in the enrichment was examined. Furthermore, the hydrophobicity of the organic molecules represented by the octanol–water partition coefficient was found to play a crucial role in determining the separation characteristics of the molecules in ultrasound-generated mist and vapor.

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

A modified many-body dissipative particle dynamics model for mesoscopic fluid simulation: methodology, calibration, and application for hydrocarbon and water

The many-body dissipative particle dynamics (mDPD) is a prominent mesoscopic multiphase model for fluid transport in mesoconfinement. However, it has been a long-standing challenge for mDPD (and other multiphase-enabled DPD models) to accurately predict real-fluid static and dynamic properties simultaneously. We have developed a modified mDPD model to overcome the issue and a rigorous calibration approach that uses reference data, including experimental and/or molecular dynamics (MD) simulations to parameterise the modified mDPD for real fluids. We choose heptane as a representative example of hydrocarbon in source rocks to demonstrate the model's capability to accurately predict the equation of state (EOS), free surface tension, diffusivity, and viscosity. Our timing test shows that the modified mDPD is 400–500 times faster than its MD counterpart for simulating bulk heptane in equivalent volumes. To further demonstrate the robustness of the model, we revisited the benchmark problem of mesoscopic modelling of water, in which all the previous works on DPD reported only a limited portion of the water properties. Here, we show that the modified mDPD can provide accurate modelling of water static and dynamic properties and an EOS that matches the experimental data to a large range of confinement pressure.

42 ENGINEERING↗
DOE OSTI2022

Extremely Slow Diffusion of Gold Nanoparticles under Confinement in Mesoporous Silica

Here, we report an X-ray photon correlation spectroscopy (XPCS) study of the mobility of colloidal gold nanoparticles with a diameter of approximately 4 nm in glycerol under confinement within the channels of SBA-15 mesoporous silica with pore diameters of approximately 6 and 12 nm. The XPCS correlation functions result from an effective heterodyne signal due to the mixing of the coherent scattering from the nanoparticles and the mesoporous silica. Over the range of wave vectors and hence length scales probed, the nanoparticle dynamics are well described by one-dimensional diffusion. The nanoparticle diffusivity varies with temperature in a manner expected based on the temperature-dependent viscosity of glycerol; however, the magnitudes of the diffusion coefficients are several orders of magnitude smaller than those of the nanoparticles in bulk glycerol. We consider mechanisms that might contribute to this reduction in diffusivity including enhancement in hydrodynamic drag under confinement, effects of nanoparticle adsorption to the pore walls, and slowing structural dynamics in glycerol due to the combination of the nanometer-scale confinement and the presence of the nanoparticles.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Anionic Effects on Concentrated Aqueous Lithium Ion Dynamics

In this article, the structural dynamics, chemical reactivity, anisotropy, diffusivity, viscosity, and density were measured for concentrated lithium salt solutions, including lithium chloride (LiCl), lithium bromide (LiBr), lithium nitrite (LiNO 2 ), and lithium nitrate (LiNO 3 ), with methyl thiocyanate as an infrared vibrational probe molecule, using two-dimensional infrared spectroscopy (2D IR), nuclear magnetic resonance (NMR) spectroscopy, and viscometry. 2D IR, NMR, and viscosity results show that LiNO 2 exhibits longer correlation times, lower diffusivity, and nearly four times greater viscosity when compared to the other lithium salt solutions of the same concentration, suggesting that nitrite anions may strongly facilitate structure formation via strengthening water-ion interactions, directly impacting bulk solution properties at sufficiently high concentrations. Additionally, the LiNO 2 and LiNO 3 solutions show significantly reduced chemical reactivity with respect to lithium cations coordinating with the methyl thiocyanate when compared to the lithium halide salts.

2D IR↗
DOE OSTI2024

Evanescent and inertial-like waves in rigidly rotating odd viscous liquids

Three-dimensional non-rotating odd viscous liquids give rise to Taylor columns and support axisymmetric inertial-like waves ( J. Fluid Mech ., vol. 973, 2023, A30). When an odd viscous liquid is subjected to rigid-body rotation however, there arise in addition a plethora of other phenomena that need to be clarified. In this paper, we show that three-dimensional incompressible or two-dimensional compressible odd viscous liquids, rotating rigidly with angular velocity 𝛺, give rise to both oscillatory and evanescent inertial-like waves or a combination thereof (which we call of mixed type) that can be non-axisymmetric. By evanescent, we mean that along the radial direction, typically when moving away from a solid boundary, the velocity field decreases exponentially. These waves precess in a prograde or retrograde manner with respect to the rotating frame. The oscillatory and evanescent waves resemble respectively the body and wall-modes observed in (non-odd) rotating Rayleigh–Bénard convection ( J. Fluid Mech ., vol. 248, 1993, pp. 583–604). We show that the three types of waves (wall, body or mixed) can be classified with respect to pairs of planar wavenumbers 𝜅 which are complex, real or a combination, respectively. Experimentally, by observing the precession rate of the patterns, it would be possible to determine the largely unknown values of the odd viscosity coefficients. This formulation recovers as special cases recent studies of equatorial or topological waves in two-dimensional odd viscous liquids which provided examples of the bulk–interface correspondence at frequencies 𝜔 < 2⁢𝛺. We finally point out that the two- and three-dimensional problems are formally equivalent. Their difference then lies in the way data propagate along characteristic rays in three dimensions, which we demonstrate by classifying the resulting Poincaré–Cartan equations.

Mechanics↗