SEARCH · Engineering Papers

Results for “Shape Matching”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 667 records · Page 37

All sources NASA NTRS DOE OSTI

NASA NTRS2022

Rotation of Low-mass Stars in Upper Centaurus-Lupus and Lower Centaurus-Crux with TESS

We present stellar rotation rates derived from Transiting Exoplanet Survey Satellite (TESS) light curves for stars in Upper Centaurus–Lupus (UCL; ∼136 pc, ∼16 Myr) and Lower Centaurus–Crux (LCC; ∼115 pc, ∼17 Myr). We find spot-modulated periods (P) for ∼90% of members. The range of light-curve and periodogram shapes echoes that found for other clusters with K2, but fewer multiperiod stars may be an indication of the different noise characteristics of TESS, or a result of the source selection methods here. The distribution of P as a function of color as a proxy for mass fits nicely in between that for both older and younger clusters observed by K2, with fast rotators being found among both the highest and lowest masses probed here, and a well-organized distribution of M-star rotation rates. About 13% of the stars have an infrared excess, suggesting a circumstellar disk; this is well matched to expectations, given the age of the stars. There is an obvious pileup of disked M stars at P ∼ 2 days, and the pileup may move to shorter P as the mass decreases. There is also a strong concentration of disk-free M stars at P ∼ 2 days, hinting that perhaps these stars have recently freed themselves from their disks. Exploring the rotation rates of stars in UCL/LCC has the potential to help us understand the beginning of the end of the influence of disks on rotation, and the timescale on which stars respond to unlocking.

L M Rebull↗

DOE OSTI2023

Optics Studies for Multipass Energy Recovery at CEBAF: ER@CEBAF

Energy recovery linacs (ERLs), focus on recycling the kinetic energy of electron beam for the purpose of accelerating a newly injected beam within the same accelerating structure. The rising developments in the super conducting radio frequency technology, ERL technology has achieved several noteworthy milestones over the past few decades. In year 2003, Jefferson Lab has successfully demonstrated a single pass energy recovery at the CEBAF accelerator. Furthermore, they conducted successful experiments with IR-FEL demo and upgrades, as well as the UV FEL driver. This multi-pass, multi-GeV range energy recovery demonstration proposed to be carried out at CEBAF accelerator at Jefferson Lab focuses on demonstrating highest energy recovery in super conducting linac in the low-current range. Continuous electron beam accelerate up to 7.5 GeV within 5-passes and decelerate in the next 5-passes recovering RF energy and dumps at a low energy dump. The beamline optics design for recirculating linacs require special attention to avoid beam instabilities due to RF wakefields. Usually, multi-pass linac beam lines require stronger focusing at lower energies as that is necessary to avoid beam breakup (BBU) instabilities, even with this small beam current. The CEBAF linac optics optimization is focused on balancing over-focusing at higher energies and beta excursions at lower energies. The race-track-shaped geometry of CEBAF accelerator allows its linacs to accommodate multiple energy beams simultaneously, while individual recirculating arcs transporting one beam energy, are shared between accelerating/decelerating beams. For the linac optics optimization process, an extended strategy is used that is originally used in 6-pass Recirculating Linac design of the LHeC, to represent the ten passes through a single linac. Using proper mathematical expressions, linac optics optimization can be achieved with evolutionary genetic algorithms, with Multi-Objective optimization. This thesis introduces a CEBAF optics redesign tailored to accommodates the ER@CEBAF multi-pass ER scheme. The isochronous arcs were retuned to match into optics solutions for optimized 10-pass linacs. Within this work, a single bunch particle tracking analysis presented here focuses on the further improvements of the beamline and beam transportation.

Neththikumara, Isurumali↗

DOE OSTI2022

CLARION2-TRINITY: A Compton-suppressed HPGe and GAGG:Ce-Si-Si array for absolute cross-section measurements with heavy ions

The design and performance of a new Compton-suppressed HPGe and charged-particle array, CLARION2-TRINITY, are described. The TRINITY charged-particle array is comprised of 64 Cerium-doped Gadolinium Aluminum Gallium Garnet (GAGG:Ce) crystals configured into five rings spanning 7–54 degrees, and two annular silicon detectors that can shadow or extend the angular coverage to backward angles with minimal $\gamma$-ray attenuation. GAGG:Ce is a non-hygroscopic, bright, and relatively fast scintillator with a light distribution well matched to SiPMs. Count rates up to 40 kHz per crystal are sustainable. Fundamental characteristics of GAGG:Ce are measured and presented, including light- and heavy-ion particle identification (PID) capability, pulse-height defects, radiation hardness, and emission spectra. The CLARION2 array consists of up to 16 Compton-suppressed HPGe Clover detectors (efficiency at 1 MeV) configured into four rings (eight HPGe crystal rings) using a non-Archimedean geometry that suppresses back-to-back coincident 511-keV gamma rays. The entire array is instrumented with 100- and 500-MHz (14 bit) waveform digitizers which enable triggerless operation, pulse-shape discrimination, fast timing, and pileup correction. Lastly, two examples of experimental data taken during the commissioning of the CLARION2-TRINITY system are given: a PID spectrum from 16 O + 18 O fusion-evaporation, and PID and Doppler-corrected -ray spectra from 48 Ti + 12 C Coulomb excitation.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2023

Skin-integrated stretchable actuators toward skin-compatible haptic feedback and closed-loop human-machine interactions

Abstract Skin-integrated haptic interfaces that can relay a wealth of information from the machine to the human are of great interest. However, existing haptic devices are not yet able to produce haptic cues that are compatible with the skin. In this work, we present the stretchable soft actuators for haptic feedback, which can match the perception range, spatial resolution, and stretchability of the skin. Pressure-amplification structures are fabricated using a scalable self-assembly process to ensure an output pressure beyond the skin perception threshold. Due to the minimized device size, the actuator array can be fabricated with a sufficiently high spatial resolution, which makes the haptic device applicable for skin locations with the highest spatial acuity. A haptic feedback system is demonstrated by employing the developed soft actuators and highly sensitive pressure sensors. Two proof-of-concept applications are developed to illustrate the capability of transferring information related to surface textures and object shapes acquired at the robot side to the user side.

Chen, Si↗

NASA NTRS1994

The 8-13 micron spectra of comets and the composition of silicate grains

We have analyzed the existing spectra of seven comets which show an emission feature at 7.8-13 micrometers. Most have been converted to a common calibration, taking into account the SiO feature in late-type standard stars. The spectra are compared with spectra of the Trapezium, interplanetary dust particles (IDPs), laboratory mineral samples, and small particle emission models. The emission spectra show a variety of shapes; there is no unique 'cometary silicate'. A peak at 11.20-11.25 micrometers, indicative of small crystalline olivine particles, is seen in only three comets of this sample, P/Halley, Bradfield 1987 XXIX, and Levy 1990 XX. The widths of the emission features range from 2.6 to 4.1 micrometers (FWHM). To explain the differing widths and the broad 9.8 micrometers maximum, glassy silicate particles, including both pyroxene and olivine compositions, are the most plausible candidates. Calculations of emission models confirm that small grains of glassy silicate well mixed with carbonaceous material are plausible cometary constituents. No single class of chondritic aggregate IDPs exhibits spectra closely matching the comet spectra. A mixture of IDP spectra, particularly the glass-rich aggregates, approximately matches the spectra of comets P/Halley, Levy, and Bradfield 1987 XXIX. Yet, if comets are simply a mix of IDP types, it is puzzling that the classes of IDPs are so distinct. None of the comet spectra match the spectrum of the Trapezium. Thus, the mineralogy of the cometary silicates is not the same as that of the interstellar medium. The presence of a component of crystalline silicates in comets may be evidence of mixing between high- and low-temperature regions in the solar nebula.

Hanner, Martha S.↗

DOE OSTI2025

Evaluating pulse-shaping capabilities of next-generation pulsed power architectures

This project evaluated the pulse shaping capabilities of next-generation pulsed power (NGPP) architectures. NGPP architectures share several common attributes including multiple independent pulse-generation lines, a radial water-insulated impedance transformer, and a central vacuum insulated load region. A multi-module circuit model was developed, incorporating independent pulse-generation lines and a 2-D transmission line mesh of the radial impedance transformer to assess the effects of azimuthal asymmetry in pulse-shaped experiments. Circuit model simulations demonstrated that NGPP architectures are able to produce the the desired current pulse shapes for exemplar NGPP experiments. Additionally, the project explored automated methods for experiment design, including derivative -ree optimization and machine learning. Pulse-shaped experiments require designers to determine machine parameters that reliably produce the desired current pulse at the load, a process that typically relies on expert knowledge and iterative adjustments using the Z circuit model. Given the increased complexity of NGPP systems, this manual approach may be impractical. While the evaluated methods do not eliminate the need for manual iteration, they can reduce the time required for experiment design. Derivative-free optimization automates much of the trial-and-error process, providing a close starting point for manual adjustments or making small modifications to near-final designs. Meanwhile, deep neural network methods can generate a good qualitative match to the desired current pulse in under one second without requiring circuit model simulations.

42 ENGINEERING↗

NASA NTRS1991

Sensor fusion III: 3-D perception and recognition; Proceedings of the Meeting, Boston, MA, Nov. 5-8, 1990

The volume on data fusion from multiple sources discusses fusing multiple views, temporal analysis and 3D motion interpretation, sensor fusion and eye-to-hand coordination, and integration in human shape perception. Attention is given to surface reconstruction, statistical methods in sensor fusion, fusing sensor data with environmental knowledge, computational models for sensor fusion, and evaluation and selection of sensor fusion techniques. Topics addressed include the structure of a scene from two and three projections, optical flow techniques for moving target detection, tactical sensor-based exploration in a robotic environment, and the fusion of human and machine skills for remote robotic operations. Also discussed are K-nearest-neighbor concepts for sensor fusion, surface reconstruction with discontinuities, a sensor-knowledge-command fusion paradigm for man-machine systems, coordinating sensing and local navigation, and terrain map matching using multisensing techniques for applications to autonomous vehicle navigation.

Schenker, Paul S.↗

NASA NTRS2005

Conservative Patch Algorithm and Mesh Sequencing for PAB3D

A mesh-sequencing algorithm and a conservative patched-grid-interface algorithm (hereafter Patch Algorithm ) have been incorporated into the PAB3D code, which is a computer program that solves the Navier-Stokes equations for the simulation of subsonic, transonic, or supersonic flows surrounding an aircraft or other complex aerodynamic shapes. These algorithms are efficient, flexible, and have added tremendously to the capabilities of PAB3D. The mesh-sequencing algorithm makes it possible to perform preliminary computations using only a fraction of the grid cells (provided the original cell count is divisible by an integer) along any grid coordinate axis, independently of the other axes. The patch algorithm addresses another critical need in multi-block grid situation where the cell faces of adjacent grid blocks may not coincide, leading to errors in calculating fluxes of conserved physical quantities across interfaces between the blocks. The patch algorithm, based on the Stokes integral formulation of the applicable conservation laws, effectively matches each of the interfacial cells on one side of the block interface to the corresponding fractional cell area pieces on the other side. This approach is comprehensive and unified such that all interface topology is automatically processed without user intervention. This algorithm is implemented in a preprocessing code that creates a cell-by-cell database that will maintain flux conservation at any level of full or reduced grid density as the user may choose by way of the mesh-sequencing algorithm. These two algorithms have enhanced the numerical accuracy of the code, reduced the time and effort for grid preprocessing, and provided users with the flexibility of performing computations at any desired full or reduced grid resolution to suit their specific computational requirements.

Pao, S. P.↗

DOE OSTI2023

Taxonomic and functional restoration of tallgrass prairie soil microbial communities in comparison to remnant and agricultural soils

Abstract Restoring ecosystems requires the re-establishment of diverse soil microbial communities that drive critical ecosystem functions. In grasslands, restoration and management require the application of disturbances like fire and grazing. Disturbances can shape microbial taxonomic composition and potentially functional composition as well. We characterized taxonomic and functional gene composition of soil communities using whole genome shotgun metagenomic sequencing to determine how restored soil communities differed from pre-restoration agricultural soils and original remnant soils, how management affects soil microbes, and whether restoration and management affect the number of microbial genes associated with carbohydrate degradation. We found distinct differences in both taxonomic and functional diversity and composition among restored, remnant, and agricultural soils. Remnant soils had low taxonomic and functional richness and diversity, as well as distinct composition, indicating that restoration of agricultural soils does not re-create soil microbial communities that match remnants. Prescribed fire management increased functional diversity, which also was higher in more recently planted restorations. Finally, restored and post-fire soils included high abundances of genes encoding cellulose-degrading enzymes, so restorations and their ongoing management can potentially support functions important in carbon cycling.

Mason, Cayla N. (ORCID:0000000345983824)↗

DOE OSTI2024

Effects of drive pulse shape on graded metal pushered single shell capsule implosions on the National Ignition Facility

Graded metal pushered single shells (PSS) are a viable alternative to low-Z capsules (Z is the atomic number) for indirect drive inertial confinement fusion implosions due to enhanced core tamping and radiation trapping, but they can be compromised by the pusher mixing with the fuel. We compare 2-shock and 3-shock laser pulses for Be/Cr PSS capsules filled with deuterium–tritium gas fuel at 6 mg/cc density. 1D radiation-hydrodynamic simulations predict higher core compression and, hence, ∼2× higher fusion yield for the 3-shock drive than for 2-shock. Nevertheless, we observe similar core ion temperatures and fusion yields for both drives. The implosion burn duration is 25% shorter and the core volume is ∼2.5× smaller for the 3-shock drive than for 2-shock, consistent with a higher compression. 1D LASNEX mix simulations using a buoyancy-drag model matching the measured yields also agree with the observed core sizes and burn durations and suggest ∼40% and ∼70% yield degradations for 2-shock and 3-shock drives due to hydrodynamic instabilities and atomic mix at the pusher–fuel interface. At the same time, 2D HYDRA simulations show that mid-mode (2–250) instability degradations are negligible for the 2-shock implosion (9%) and significant (45%) for 3-shock. Subtracting these from the 1D mix simulations, we infer similar degradations from high-mode instabilities and atomic mix for both drives. Due to its robustness to mid-mode instabilities, future pusher–gas mix studies will use the 2-shock drive.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Velocity profile shapes in Alcator C-Mod plasmas

Abstract Toroidal rotation velocity spatial profiles ( r / a < 0.8) have been obtained from C-Mod over a wide range of operational conditions, including H-mode, I-mode, ICRF-heated L-mode and Ohmic L-mode (LOC and SOC), and in plasmas with ITBs, LH wave injection and MCFD. Peaked, flat and hollow rotation profiles have been observed. In H- and I-mode plasmas, generally with co-current peaked profiles, the peaking is correlated with temperature profile peaking, and both increase with toroidal magnetic field (decrease with ρ ∗ ). Any dependence on density peaking is unclear. For Ohmic L-mode discharges, with LOC, the velocity profiles are usually flat and most often directed co-current, while with SOC the profiles are hollow, mostly co-current at the edge and counter-current in the core. Both of these Ohmic rotation states exist with matched density and temperature profiles (and gradients), indicating that neither gradient is relevant during rotation reversals. For plasmas with LH wave injection and discharges with ITBs, the velocity profiles are hollow while the density and temperature profiles exhibit substantial peaking. Broadly speaking for all operational regimes, there is no unifying ordering of the velocity gradient with plasma parameters.

Rice, J. E. (ORCID:0000000183195971)↗

NASA NTRSDate not supplied

Flow Characterization of the NASA Langley Unitary Plan Wind Tunnel, Test Section 2: Computational Results

This is an abstract for an invited paper at the AIAA Aviation Conference, June 2021. The work described here is part of an effort of coordinated experiments in the Unitary Plan Wind Tunnel (UPWT) facility at the NASA Langley Research Center (LaRC) and matching CFD simulations. The primary goal of the work is to assess the productivity and true predictive accuracy of CFD, absent any guidance from experiment, in the high supersonic speed range as compared to experiments performed at the NASA LaRC’s UPWT facility. This report concerns CFD simulation of the primary flow-path in the empty wind tunnel, including the settling chamber, nozzle, test section, and some of the tunnel downstream of the test section. The Mach number in the test section ranges from M~2.4 to M~4.6, and the required area ratio variation is achieved by translation of a nozzle block which constricts the area of a loosely S-shaped throat. Flow past protuberances in the settling chamber and into this S-bend throat are predicted by CFD to generate streamwise vorticity that may, or may not, persist through the throat and into the test section as coherent vortices. Some flow conditions are notably unsteady at frequencies well below those of turbulence, due to unsteady separated flow ahead of the nozzle block. The bulk flow moves at velocities ranging from 'walking speed' in the settling chamber to M~4.6 in the test section. Heat transfer to the settling chamber walls and buoyancy are significant at high Mach number. Subtle variations in surface curvature in the nozzle generate Mach waves that propagate into the test section. CFD of the empty tunnel serves two purposes. Firstly, the full-tunnel simulations are used to provide upstream boundary conditions for CFD of vehicle aerodynamics which are generally performed in a domain confined to the wind tunnel test section; these companion studies are addressed in other papers. Secondly, it is a challenging test for CFD to resolve all of the empty tunnel flow phenomena relevant to flow in the test section. It requires a more complete definition of geometry than was originally anticipated. In addition, it requires good spatial and temporal accuracy, and turbulence modeling that performs well on specific phenomena such as corner flows. The boundary conditions and solution algorithms must perform well from incompressible to almost hypersonic speeds. The final state of the pre-test CFD was a product of an iterative self-improvement process. The initial simulations of the empty tunnel were deficient in many respects, but hints to those deficiencies were recognized in the solutions, and remedies were implemented. Possible further improvements will be studied in the post-test phase when comparisons with experimental data are possible. The CFD was performed by five separate collaborative teams using four different flow solvers: FUN3D, Overflow, Star-CCM+ and USM3D. The level of effort of these teams varied significantly, but each made important contributions to the goals of the work. A concerted effort to use uncertainty quantification methods (UQ) in CFD is also a goal of this work. To this end, variations in CFD results due to grid refinement, turbulence modeling, and boundary conditions have been characterized. Code-to-code variation is another means of assessing CFD uncertainty. All CFD solvers predict similarity among the primary flow features; these include the variation of Mach number due to changes in Reynolds number, and the bulk flow angularity due to tunnel-wall curvature. All CFD solvers also predict similar trends in secondary flows, such as the downwash in the side-wall boundary layers. Three of the high-spatial resolution simulations give similar predictions of a complex secondary flow phenomena, streamwise vortices generated in the S-bend throat that persist into the side-wall boundary layers of the test section. Two of the highest-resolution simulations, run with the same turbulence model in different CFD solvers, gave encouragingly similar predictions of a complex tertiary flow phenomenon, small transient "sprites" of upwash flows, resulting from vortices that presumably originate in the separated flow near the leading edge of the nozzle block. While the CFD was done in a "blind pre-test" mode, requests from the experimental team for CFD results pertaining to unsteadiness and total temperature variations in the test section resulted in CFD runs that included a cooled wall in settling chamber. This then led to a change in the standard practice for running the Overflow results, which would not have occurred without this "release" of this experimental information. CFD was also used to guide some measurements. The paper will focus on establishing the consensus among CFD results and understanding differences among those results. Initial findings from the efforts to characterize CFD uncertainty have been done and will be included in the paper.

Robert Edward Childs↗

NASA NTRS1988

Unsteady heat transfer in turbine blade ducts: Focus on combustor sources

Thermal waves generated by either turbine rotor blades cutting through nonuniform combustor temperature fields or unsteady burning could lead to thermal fatigue cracking in the blades. To determine the magnitude of the thermal oscillation in blades with complex shapes and material compositions, a finite element Galerkin formulation has been developed to study combustor generated thermal wave propagation in a model two-dimensional duct with a uniform plug flow profile. The reflection and transmission of the thermal waves at the entrance and exit boundaries are determined by coupling the finite element solutions at the entrance and exit to the eigenfunctions of an infinitely long adiabatic duct. Example solutions are presented. In general, thermal wave propagation from an air passage into a metallic blade wall is small and not a problem. However, if a thermal barrier coating is applied to a metallic surface under conditions of a high heat transfer, a good impedance match is obtained and a significant portion of the thermal wave can pass into the blade material.

Baumeister, Kenneth J.↗

NASA NTRS1988

Unsteady heat transfer in turbine blade ducts - Focus on combustor sources

Baumeister, K. J.↗

NASA NTRS2006

The Zero-Degree Detector System

We will report on a detector system used for accelerator measurement of nuclear fragmentation cross sections. This system consists of two detector planes, each carrying a ring of 8 detectors. Each detector has 64 pads. These two detector planes are arranged facing each other so that the matching detector pads on each plane form a two element charged particle telescope. Each of these telescopes is capable of determining the elemental identity of nuclear fragments passing through it. The system is used to measure light fragment production in the presence of heavier fragments. We will present a detailed discussion of the 64-pad detector design, the substrate design. The front-end electronics used to read out the signals is based on a custom VLSI chip developed for the Advanced Thin Ionization Calorimeter experiment which has been flown successfully twice in Antarctica. Each of these chips has 16 channels and each channel consists of a charge-sensitive preamplifier followed by a shaping amplifier and a track-and-hold circuit. The track-and-hold circuits are connected via a multiplexer to an output line driver. This allows the held signals to be presented, one-by-one via a common data line to a analog-to-digital converter. Because the output line driver can be placed in a high input impedance state when not in use, it is possible to daisy-change many chips on the same common data line. The front-end electronics and data readout scheme will be discussed in detail. The Zero Degree Detector has been used in several accelerator experiments conducted at the NASA Space Radiation Laboratory and the Alternating Gradient Synchrotron at Brookhaven National Laboratory as well as at the HIMAC accelerator in Japan. We will show examples of data taken at these accelerator runs to demonstrate how the system works.

Adams, James H.↗

DOE OSTI2025

Rapid quantification of whole seed fatty acid amount, composition, and shape phenotypes from diverse oilseed species with large differences in seed size

Seed oils are widely used in the food, biofuel, and industrial feedstock industries, with their utility and value determined by total oil content and fatty acid composition. Current high throughput seed oil analysis methods either lack accuracy in total fatty acid profiling or require extensive labor for lipid extraction prior to derivatization to fatty acid methyl esters (FAME) and quantification by gas chromatography (GC). Alternatively, direct whole seed FAME production methods have been developed for the very small seeds in the model species Arabidopsis thaliana but these have generally not been adapted to larger seeds of most oilseed crops. High-throughput direct whole seed FAME production methods were optimized for seeds up to 5 mg each utilizing acid-catalyzed esterification. For the oilseed species Camelina sativa, Thlaspi avernse (pennycress), Cuphea viscosissima, and Brassica napus (var. Canola), the total seed fatty acid content and composition from direct seed esterification to FAME matched that of lipid extract derivatization demonstrating the accuracy of the methods. In combination with seed phenotyping using GridFree, this approach enabled the development of a rapid pipeline for simultaneous seed weight, count, size/shape phenotyping, and oil analysis. For the larger and tougher seeds produced by Limnanthes alba (Meadowfoam) and Cannabis sativa L. (hemp) the whole seed acid-based method proved insufficient, and prior laborious homogenization of seeds was required. Therefore, a rapid one-tube bead homogenization and base catalyzed-esterification method was developed. Base-derived fatty acid esterification cannot derivatize free fatty acids leading to slightly lower total seed fatty acid than acid-catalyzed methods, however the seed oil content and fatty acid composition that is valuable for screening large numbers of samples in research populations was accurately measured. New rapid whole seed fatty acid esterification and phenotyping protocols were developed to accurately assess oilseed lipid content. These methods are particularly valuable in oilseed research, breeding, and engineering applications where efficient analysis of large numbers of samples and accurate oil fatty acid profiling is essential. While having been developed for current and emerging oilseed crops, these methods also provide a foundation from which protocols might be established for new and emerging crop species.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2025

Coupling of Charge Regulation and Geometry in Soft Ionizable Molecular Assemblies

The size, shape, and charge of structures, such as proteins and amphiphile assemblies, respond in an interconnected manner to solution ionic conditions. Here, we analyze assemblies of an amphiphile (C 16 K 2 ), with two ionizable amino acids [lysine (K)] coupled to a 16-carbon alkyl tail, via small-angle X-ray scattering (SAXS), nonlinear Poisson–Boltzmann theory (nl-PB), and hybrid Monte Carlo-molecular dynamics (MC-MD) simulations. SAXS revealed structural transitions from spherical micelles to cylindrical micelles to bilayers with increasing pH. By combining SAXS-determined structural information and nl-PB, we derived the molecular degree of ionization as a function of pH. The back-calculated titration curves matched the experimental data over an extended pH range, without adjustable parameters. Similarly, the SAXS data on the evolution of spherical micelle structure with ionic strength were combined with nl-PB and MC-MD to derive the bare and effective charges. MC-MD, which considered finite ion sizes, showed that bare and effective charges saturate quickly with increasing salt concentration. Furthermore, the calculated effective charges closely matched results from Zeta-potential measurements. The presented approach has advantages over customary methods for charge regulation, such as the Henderson–Hasselbalch (HH) or Hill models, where molecular ionization/deionization in assemblies is described by effective pKs that are distinct from the pK for isolated molecules. However, these models lack a physical explanation for these pK shifts. By contrast, our approach of combining structural details with an electrostatic model and simulations provides a more intuitive understanding of structure-charge coupling and a framework for understanding charge regulation in many synthetic and biological systems.

ionization↗

DOE OSTI2026

Antarctic ice sheet model comparison with uncurated geological constraints shows that higher spatial resolution improves deglacial reconstructions

Accurately reconstructing past changes to the shape and volume of the Antarctic ice sheet relies on the use of physically based and thus internally consistent ice sheet modeling, benchmarked against spatially limited geologic data. The challenge in model benchmarking against geologic data is diagnosing whether model-data misfits are the result of an inadequate model, inherently noisy or biased geologic data, and/or incorrect association between modeled quantities and geologic observations. In this work we address this challenge by (i) the development and use of a new model-data evaluation framework applied to an uncurated data set of geologic constraints, and (ii) nested high-spatial-resolution modeling designed to test the hypothesis that model resolution is an important limitation in matching geologic data. While previous approaches to model benchmarking employed highly curated datasets, our approach applies an automated screening and quality control algorithm to an uncurated public dataset of geochronological observations (specifically, cosmogenic-nuclide exposure-age measurements from glacial deposits in ice-free areas). This optimizes data utilization by including more geological constraints, reduces potential interpretive bias, and allows unsupervised assimilation of new data as they are collected. We also incorporate a nested model framework in which high-resolution domains are downscaled from a continent-wide ice sheet model. We highlight the application of this framework by applying these methods to a small ensemble of deglacial ice-sheet model simulations, and demonstrate that the nested approach improves the ability of model simulations to match exposure age data collected from areas of complex topography and ice flow. We develop a range of diagnostic model-data comparison metrics to provide more insight into model performance than possible from a single-valued misfit statistic, showing that different metrics capture different aspects of ice sheet deflation.

Geosciences↗