Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distribution function”

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

Roles of surface Mn(III) and MgOH + /CaOH + complexation in photochemically assisted Mn 2+ oxidation and Todorokite formation

Todorokite, a Mn oxide phase commonly found in seafloor Mn nodules, is one of the most abundant Mn oxides in nature and has wide applications in environmental engineering, including heavy metal remediation, contaminant adsorption for water purification, and catalysis for organic pollutant degradation. However, traditional abiotic formation of todorokite reported to date typically requires high temperature, pressure, and pH outside circumneutral conditions, or involves several days of phase transformations from precursor phases. In this study, we demonstrate rapid todorokite formation within 6 h via photochemically assisted oxidation of Mn 2+ (aq) in the presence of Mg 2+ or Ca 2+ , notably at ambient temperature and pressure. The concentrations of aqueous MgOH + and CaOH + , as well as the pH (initial 9 or 6), had significant collective influences on the crystalline phases of the formed Mn oxide solids. The proportions of the Mn oxidation states (IV, III, and II) and the extent of Mn(III) surface complexation with MgOH + or CaOH + controlled the crystalline phases (todorokite, feitknechtite, and birnessite) and the extent of cation incorporation into the resulting Mn oxide solids. Specifically, the oxidation state of Mn influenced the crystallization of different Mn oxide phases, and the formed Mn(III) complexed with MgOH + or CaOH + facilitated the formation of todorokite, along with cation incorporation. X-ray pair distribution function analysis revealed that incorporated Mg species were located at the inside corner sites of the 5 × 3 MnO 6 octahedral tunnel structure of todorokite, suggesting complexation between MgOH + and Mn(III) at the corners of the tunnel. Furthermore, this study provides novel insights into the elusive formation of todorokite and its potential for applications in advanced water treatment and environmental remediation.

Abiotic Mn2+ oxidation↗
DOE OSTI2020

Molecular-Level Behavior of Imidazolium-Based Ionic Liquid Mixtures

Understanding the molecular-level organization of anions and cations in a mixture of different ionic liquids (ILs) is challenging, but it can be very important for designing solvents and IL-based membrane materials. In this work, we perform all-atom molecular dynamics (MD) simulations to understand the intermolecular interactions and structuring of both neat and mixed IL compositions containing 1-benzyl-3-methylimidazolium [ZMIM + ], 1-butyl-3-methylimidazolium [BMIM + ], bistriflimide [Tf 2 N – ], and tetrafluoroborate [BF 4 – ]. We observe very strong underlying molecular-level migrations in these mixtures, and this tends to mainly involve the anion species. As compared to the systems with a single anion, the mixed anions lead to large increases/decreases in the radial distribution function peak heights, whereas the cation mixing is relatively benign. The [BF 4 – ] anions have a strong tendency to displace the [Tf 2 N – ] anions in the mixtures, and this is seen with respect to both [BMIM + ] and [ZMIM + ].

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

Laser diagnostics to characterize the in-flame growth of platinum nanoparticles manufactured by the reactive spray deposition technology

Reactive Spray Deposition Technology (RSDT) is an atmospheric pressure, flame-based advanced manufacturing method used to fabricate Membrane Electrode Assemblies (MEAs) for proton exchange fuel cells and water electrolyzers. RSDT combines the flame synthesis of catalyst nanoparticles and their deposition onto the membrane in one step. The properties of the synthesized nanoparticles, such as their Size Distribution Function (SDF), determine the manufactured electrodes’ performance, which can be evaluated when the deposition is complete via ex-situ characterization of the products. Efforts to improve RSDT for manufacturing state-of-the-art MEAs can be significantly enhanced with integrated laser diagnostics that enable in situ measurement of the synthesized catalyst nanoparticles. This paper reports the implementation of laser diagnostics in an RSDT facility and evaluates their potential to assist the manufacturing process. Laser Light Scattering (LS) and Laser-Induced Incandescence (LII) measurements are performed in the oxygen-rich zone of the flame at various distances from the flame fuel jet nozzle downstream of its luminous region. The measurements quantitatively track the evolution in size and volume fraction of platinum-based nanoparticles in two flames that yield different catalytic properties in the manufactured electrodes. The profiles of the measured nanoparticles’ volume fraction along the flame axis can be estimated a priori so that the average nanoparticle sizes can be measured in quasi-real time via LS. Nanoparticles experience an extremely slow growth rate while being convected at distances from 150 nm to 300 nm from the fuel nozzle. Complementary characterizations of the synthesized nanoparticles are performed ex-situ via High-Angle Annular Darkfield Scanning Transmission Electron Microscopy (HAADF-STEM) image analyses of samples collected on grids at a fixed distance from the fuel jet nozzle. Comparing the results from laser and microscopy techniques not only cross-validates the findings but also provides the parameters to infer the bimodal SDF in-situ and yields evidence that the coagulation efficiency of the synthesized nanoparticles has extremely low values in the investigated zone of the flames.

Rayleigh Light Scattering (LS)↗
DOE OSTI2024

Direct numerical simulations of turbulent premixed cool flames: Global and local flame dynamics analysis

The cool flame dynamics, especially in turbulent flows, is of great interest for both practical application and fundamental research. Here, in this study, a series of direct numerical simulations of turbulent premixed n-C 7 H 16 /O 2 /O 3 /N 2 cool flames are performed, with the focus on the influence of turbulence intensity (u'/S L , where S L is the laminar flame speed) on the flame structure as well as the global and local cool flame dynamics. It is found that the cool flame front is considerably wrinkled by turbulence at high u'/S L , leading to significantly thickened turbulent cool flame brush and largely altered local reactivity compared with the reference laminar flame. However, the turbulent flame structure in the temperature space is found to be insensitive to u'/S L . Besides, with increasing u'/S L , the normalized turbulent cool flame speed (S L /S L ) is monotonically increased, attributed to substantial augmentation on the flame surface area (A T /A L ), while the stretching factor (I 0 ) remains almost constant and is smaller than 1. The underlying mechanisms for such variations are revealed through local flame dynamics analysis. Specifically, the local flame displacement speed S d is found to be strongly negatively correlated with flame curvature; meanwhile, such negative correlation and the probability distribution function (PDF) of flame curvature are barely influenced by u'/S L , leading to a weak dependence of I 0 on u'/S L . In contrast, the PDF of the tangential strain rate is found to span a much wider range and shift to the positive side as u'/S L increases, suggesting that the enhanced tangential strain rate is the main cause for the increase in surface area of the turbulent premixed cool flame. Finally, the influence of equivalence ratio on above findings is found to be insignificant, indicating that although the local reactivity of turbulent premixed cool flames is altered due to the differential diffusion, the resultant flame- stretch interaction is insensitive to the equivalence ratio. This study presents some unique cool flame dynamics that are distinct from hot flames, which can help improve the understanding and modeling of turbulent cool flames.

Cool flames↗
DOE OSTI2025

Computational and experimental analysis of structural and Thermophysical properties of LiX-KX (X = chloride, iodide, or bromide) molten salts

Molten salts containing lithium (Li) and potassium (K) halides, including chlorides (Cl), bromides (Br), and iodides (I), are pivotal in advanced technological applications including nuclear reactor technologies, thermal batteries, and industrial pyrochemical processes. Despite their significance, the availability of the temperature-dependent structural and thermophysical properties of these molten salts remain scarce. Here, this study addresses this knowledge gap by predicting and measuring the key properties such as coordination number, radial distribution function, density, heat capacity, and volumetric thermal expansion of molten LiI-KI, LiBr-KBr and LiCl-KCl eutectic mixtures. Utilizing ab initio molecular dynamics (AIMD) simulations, we predict these properties of LiX-KX (X=Br,I, Cl) across various temperatures, with chlorides as a benchmark. Experimental measurements using Archimedes methods and Differential Scanning Calorimetry validate the densities and heat capacities of these molten salts. Our comprehensive analysis of the structural and thermophysical properties provides critical insights into the behavior of these molten salts at various temperatures, enhancing the understanding necessary for their application in advanced technologies.

36 - MATERIALS SCIENCE↗
DOE OSTI2023

Investigation of 3D printed lightweight hybrid composites via theoretical modeling and machine learning

Hybrid composites combine two or more different fillers to achieve multifunctional or advanced material properties, such as lightweight and enhanced mechanical properties. The properties of the composites significantly depend on their microstructures, which can be tailored via advanced 3D printing processes. Understanding the process-structure-property relationships is critical to enable the design and engineering of novel hybrid composites for applications in aerospace, automotive, and protective coatings. Here, for this work, we develop 3D printable and lightweight hybrid composites and leverage the conventional design of experiments, a theoretical hybrid model, and an image-driven machine learning (ML) method to investigate their mechanical behaviors. The hybrid composites are formulated with elastomer matrix, microfillers, and thin-shell particles, enabling a significant degree of design freedom of microstructures with densities and mechanical properties varying up to 70% and 91%, respectively. Our statistical analysis indicates that the 3D printing path direction and the microfibers fraction are dominating process parameters with contribution percentages of 45.3% and 57.7% on the specific stiffness and strength, respectively. A hybrid mechanics model is developed based on a simple Weibull distribution function and classical single-filler models to effectively capture the variations in mechanical properties, however, it overestimates the values due to its statistical constraints and idealization of experimental uncertainty. The image-driven ML model leverages the microscale images directly without losing the structural details, shows more accurate predictions with experimental data, and has 48.6% lower root mean square error than the theoretical model.

3D printing↗
DOE OSTI2020

Efficient Fourier transforms for transverse momentum dependent distributions

Hadron production at low transverse momenta in semi-inclusive deep inelastic scattering can be described by transverse momentum dependent (TMD) factorization. This formalism has also been widely used to study the Drell-Yan process and back-to-back hadron pair production in $e^+e^-$ collisions. These processes are the main ones for extractions of TMD parton distribution functions and TMD fragmentation functions, which encode important information about nucleon structure and hadronization. One of the most widely used TMD factorization formalism in phenomenology formulates TMD observables in coordinate $b_\perp$-space, the conjugate space of the transverse momentum. The Fourier transform from $b_\perp$-space back into transverse momentum space is sufficiently complicated due to oscillatory integrands that it requires a careful and computationally intensive numerical treatment in order to avoid potentially large numerical errors. Within the TMD formalism, the azimuthal angular dependence is analytically integrated and the two-dimensional $b_\perp$ integration reduces to a one-dimensional integration over the magnitude $b_\perp$. In this paper we develop a fast numerical Hankel transform algorithm for such a $b_\perp$-integration that improves the numerical accuracy of TMD calculations in all standard processes. Libraries for this algorithm are implemented in Python 2.7 and 3, C++, as well as FORTRAN77. All packages are made available open source.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Benchmark of the KGMf with a coupled Boltzmann equation solver

The Kinetic Global Model framework (KGMf) is an open-source general-purpose global model (spatially averaged) simulation code developed to explore the reaction kinetics and pathways in plasma discharge systems. It contains species continuity and electron energy balance equations, with a time-dependent evaluated electron energy distribution function (EEDF) for electron impact reactions. The EEDF is utilized to determine the rate coefficients for electron impact reactions, which can have profound impact on the temporal evolutions of plasma parameters. Previously, the EEDF was commonly assumed as Maxwellian or an analytical function of the effective electron temperature. In this work, the KGMf is coupled with a Boltzmann equation (BE) solver to self-consistently compute the EEDF. The EEDF evolution frequency is determined based on relative changes of the reduced electric field. The KGMf is benchmarked with the ZDPlasKin code based on high-pressure low-temperature argon plasma discharge cases. Additionally, the temporal evolutions of reduced electric field, electron temperature, EEDF, reaction rates, and species densities, are obtained and compared under different discharge conditions, showing good agreement between the KGMf and the ZDPlasKin simulations. The application of the KGMf for predicting breakdown times in high power microwave discharges is also presented, which shows qualitative agreement with particle-in-cell simulations. The KGMf can be further applied for more complicated plasma discharge systems, where the reaction kinetics are more intricate (e.g., plasma-assisted combustion systems).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Computing the QRPA level density with the finite amplitude method

Here, we describe a new algorithm to calculate the vibrational nuclear level density of an atomic nucleus. Fictitious perturbation operators that probe the response of the system are generated by drawing their matrix elements from some probability distribution function. We use the Finite Amplitude Method to explicitly compute the response for each such sample. With the help of the Kernel Polynomial Method, we build an estimator of the vibrational level density and provide the upper bound of the relative error in the limit of infinitely many random samples. The new algorithm can give accurate estimates of the vibrational level density. Since it is based on drawing multiple samples of perturbation operators, its computational implementation is naturally parallel and scales like the number of available processing units.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

A higher-order finite-element implementation of the nonlinear Fokker–Planck collision operator for charged particle collisions in a low density plasma

Collisions between particles in a low density plasma are described by the Fokker–Planck collision operator. In applications, this nonlinear integro-differential operator is often approximated by linearised or ad-hoc model operators due to computational cost and complexity. In this work, we present an implementation of the nonlinear Fokker–Planck collision operator written in terms of Rosenbluth potentials in the Rosenbluth–MacDonald–Judd (RMJ) form. The Rosenbluth potentials may be obtained either by direct integration or by solving partial differential equations (PDEs) similar to Poisson's equation: we optimise for performance and scalability by using sparse matrices to solve the relevant PDEs. We represent the distribution function using a tensor-product continuous-Galerkin finite-element representation and we derive and describe the implementation of the weak form of the collision operator. We present tests demonstrating a successful implementation using an explicit time integrator and we comment on the speed and accuracy of the operator. Finally, we speculate on the potential for applications in the current and next generation of kinetic plasma models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Fouling behavior of zwitterionic membranes compared to polyamide membranes

Membrane fouling remains a critical bottleneck for reverse osmosis (RO) desalination, driving energy consumption and reducing membrane lifetime. Here, we employ all-atom molecular dynamics simulations to investigate the antifouling behavior of random zwitterionic amphiphilic copolymer (r-ZAC) membranes composed of sulfobetaine methacrylate (SBMA) and allyl methacrylate (AMA), benchmarked against conventional polyamide (PA) RO membranes. Structural and dynamical analyses—including radial distribution functions, coordination numbers, tetrahedral order parameters, vector orientation, and residence-time correlation functions—reveal that r-ZAC surfaces sustain tightly bound, long-lived hydration layers with preserved tetrahedrality and anisotropic water orientation, in sharp contrast to the weak and disordered hydration of PA. Steered molecular dynamics simulations demonstrate that r-ZAC membranes impose substantial free-energy barriers to foulant approach (alginate ≈ 90 kcal/mol, sucrose ≈ 35 kcal/mol, humic acid ≈ 15 kcal/mol), whereas PA membranes exhibit negligible barriers (< 1 kcal/mol) and thermodynamically favorable adsorption. Detailed foulant–surface interaction analyses show that zwitterionic hydration and electrostatic heterogeneity in r-ZAC suppress adhesion, except in the case of amphiphilic humic acid, which exploits multiple binding modes. Together, these results establish molecular-level design principles for antifouling membranes: the combination of zwitterionic hydration, structured interfacial water, and controlled amphiphilic balance in r-ZAC membranes provides superior resistance to organic fouling relative to PA.

Cross-linked polyamide↗
DOE OSTI2021

Results from single-shock Marble experiments studying thermonuclear burn in the presence of heterogeneous mix on the National Ignition Facility

The Marble campaign on the National Ignition Facility investigates the effect of heterogeneous mix on thermonuclear burn for comparison to a probability distribution function (PDF) burn model. Marble utilizes plastic capsules filled with deuterated plastic foam and a fill gas containing tritium. As the capsules implode, the deuterium in the foam mixes with the tritium gas, and DT neutrons are produced as the shocks compress and heat the mixture. The yield of DT neutrons is dependent on the uniformity of the mix, with more heterogeneous mix producing fewer neutrons. In Marble, the heterogeneity of the mix is controlled by varying the diameter of voids introduced into the foam. The first NIF Marble campaign has been executed in which the Marble capsules were indirectly driven with a single strong shock using NIF hohlraums. The experiments produce a low-convergence, high-ion-temperature implosion. The ratio of DT to DD neutron yield is largely consistent with uniform atomic mix for fine-pore foam, and increases slightly with void diameter, contrary to 1D simulations using the PDF burn model. Recent 3D high-resolution simulations of similar experiments performed on the Omega Laser Facility suggest an explanation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Hot-melt extrudability of amorphous solid dispersions of flubendazole-copovidone: An exploratory study of the effect of drug loading and the balance of adjuvants on extrudability and dissolution

The FDA-approved anthelmintic flubendazole has shown potential to be repositioned to treat cancer and dry macular degeneration; however, its poor water solubility limits its use. Amorphous solid dispersions may overcome this challenge, but the balance of excipients may impact the preparation method and drug release. Here, we evaluate the influence of adjuvants and drug loading on the development of an amorphous solid dispersion of flubendazole-copovidone by hot-melt extrusion. The drug, copovidone, and adjuvants (magnesium stearate and hydroxypropyl cellulose) mixtures were statistically designed, and the process was performed in a twin-screw extruder. The study showed that flubendazole and copovidone mixtures were highly extrudable, except when drug loading was high (>40%). Furthermore, magnesium stearate positively impacted the extrusion and was more effective than hydroxypropyl cellulose. The extruded materials were evaluated by modulated differential scanning calorimetry and X-ray powder diffraction, obtaining positive amorphization and physical stability results. Pair distribution function analysis indicated the presence of drug rich domains with medium-range order structure and no evidence of polymer-drug interaction. All extrudates presented faster dissolution (HCl, pH 1.2) than pure flubendazole, and both adjuvants had a notable influence on the dissolution rate. In conclusion, hot-melt extrusion may be a viable option to obtain stable flubendazole: copovidone amorphous dispersions.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

A probabilistic creep model incorporating test condition, initial damage, and material property uncertainty

Uncertainty is prevalent in the creep resistance of alloys, where at elevated temperature and low pressure, rupture can range across logarithmic decades. In this study, a probabilistic continuum-damage-mechanics (CDM)-based model is derived to capture the uncertainty of creep resistance. To meet this objective, creep data for alloy 304 Stainless Steel is gathered. A constitutive model, “Sinh”, is calibrated deterministically to determine the statistical variability of the material properties. Three sources of uncertainty are injected into the model: test condition (stress and temperature), initial damage, and material properties. Probabilistic simulations are carried out by (a) calibrating probability distribution functions (pdfs) for each source of uncertainty (b) randomly sampling the pdfs using Monte Carlo methods and (c) executing simulations to replicate the uncertain creep behavior. A sensitivity analysis is performed to evaluate the relative effect of each source of uncertainty. In full probabilistic simulations, the cumulative uncertainty of creep behavior is evaluated. The probabilistic model accurately predicts the creep deformation and rupture of the available experiments. The probabilistic model is validated for interpolation but lacks extrapolation ability. Several future works are proposed to further improve the model.

36 MATERIALS SCIENCE↗
DOE OSTI2023

An anomalous structure disordering in Zr–Cu–Ag supercooled glass-forming liquids

An anomalous ordering-to-disordering transition in a Zr–Cu–Ag supercooled glass-forming liquid has been revealed by in situ high-energy x-ray diffraction. Upon heating, the supercooled liquid exhibits a structural crossover at T c = 716 K, evident by the change in statistical characteristics of the first sharp diffraction peak. When annealed in the vicinity of T c , the supercooled liquid first relaxes from its initial state to a denser and more ordered state. Subsequently, it becomes disordered at the medium-range scale with a slowing down of densification, as demonstrated by the broadening of the first sharp diffraction peak and a decrease in the peak intensities of the higher atomic shells in the reduced pair distribution function. Finally, a scenario for the evolution of less correlated locally favored clusters in the liquid–liquid phase transition is discussed.

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

Silica nanocluster binding rate coefficients from molecular dynamics trajectory calculations

Oxide nanoparticle growth from vapor phase precursors occurs in high temperature aerosol reactors first via the formation of nanoclusters, which are nanometer-scale condensed-phase species composed of 10 1 -10 2 atoms. The binding rate for nanoclusters, defined as the rate at which two nanoclusters collide with and stick to one another, is hence of interest in predicting nanoparticle size distribution functions in an aerosol. We have utilized molecular dynamics (MD) trajectory calculations to determine the homogeneous (equal-sized) and heterogeneous (disparate-sized) binding rate coefficients of SiO2 (silica) nanoclusters composed of 18, 144, and 333 atoms. MD trajectory calculations incorporated all-atom models of nanoclusters using a combined Born–Huggins–Mayer-Lennard-Jones potential model, which accounts for both short range interactions and electrostatic interactions. MD calculations were utilized to determine the binding probability as a function of both initial relative velocity and impact parameter; integration of the binding probability across impact parameter and relative velocity space yields the binding rate coefficient. MD trajectory calculations reveal that the most common type of encounters between nanoclusters leading to binding are grazing collisions, i.e. instances where collision would not occur without induced dipole potential influences. The resulting binding rate coefficients are found to be extremely weakly dependent on system temperature, which is in contrast to the use of rate coefficient models which are the product of a hard-sphere collision rate coefficient and a constant enhancement factor (leading to a rate coefficient proportional to the square root of temperature). Enhancement factors defined from MD trajectory calculations fall in the range of 3–9 as system temperature decreases from 1500 K to 300 K. Furthermore, such large values suggest that potential interactions need to be considered when modeling oxide nanocluster growth in the gas phase.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Complex magnetic ground state driving a large rotating magnetocaloric effect in Tb 3 Ni at low temperature

The rotating magnetocaloric effect (RMCE) offers a promising alternative to conventional magnetocaloric configurations by taking advantage of magnetic anisotropy to simplify device architecture and enhance refrigeration efficiency. In this study, the RMCE in a high-quality single crystal of Tb 3 Ni is investigated, a compound previously shown to exhibit significant magnetocaloric behavior along its easy axis of magnetization. By measuring the magnetization and corresponding entropy change along the three principal crystallographic axes using a discontinuous measurement protocol, we verify the anisotropic magnetic properties and derive the RMCE from rotations between hard ( a, b ) and easy ( c ) axes of magnetization. Our results show a maximum value for the rotating entropy change of 19.5 J kg −1 K −1 for µ 0 H = 7 T around 60 K, within the critical temperature window for industrial gas liquefaction applications. Neutron diffraction and magnetic Pair-Distribution Function (mPDF) analysis reveal that the origin of this large anisotropic response lies in a partially ordered incommensurate spin-density wave phase and persistent short-range ferromagnetic (FM) correlations. These complex magnetic states enable the release of a substantial amount of magnetic entropy when the field is applied along the easy c-axis, effectively driving the large RMCE. Comparison with other RMCE materials confirms Tb 3 Ni as one of the most promising candidates in this temperature regime, offering both a large magnetic entropy change and a wide operating window.

Gas liquefaction↗
DOE OSTI2021

The impact of cold electrons and cold ions in magnetospheric physics

A review of the impact of the cold-ion and cold-electron populations in the Earth’s magnetosphere is presented in this paper. The cold populations are defined by total energy less than approximately 100 eV, i.e. in the energy range which is strongly affected by spacecraft charging and that often dominates the total plasma density. We also include the warm plasma cloak in the review, since it overlaps partially with the cold energy range and is a population that is still not well understood. The known impacts of cold ions and cold electrons that are discussed are: the source of magnetospheric plasma, solar-wind/magnetosphere coupling, magnetotail reconnection and substorms, Kelvin-Helmholtz instabilities at the magnetopause, chorus, hiss, electromagnetic-ion-cyclotron and ultra-low-frequency wave-particle interactions, aurora structuring and spacecraft charging. Other possible impacts are associated with refilling on open-drift trajectories, the remnant layer and plasmapause disruption. A discussion of the difficulty of cold-plasma measurements and the need for new measurement techniques that measure the full cold-ion and cold-electron distribution functions is also presented. There remain a lot of unknowns about the cold-ion and cold-electron populations, associated with their origin, properties, drivers and impacts. These populations will need to be fully understood before the magnetosphere-ionosphere system can be fully understood.

79 ASTRONOMY AND ASTROPHYSICS↗