Search NASA⌕ Search

SEARCH · Search NASA

Results for “Diffuse interface methods”

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

A higher-order finite element reactive transport model for unstructured and fractured grids

Abstract This work presents a new reactive transport framework that combines a powerful geochemistry engine with advanced numerical methods for flow and transport in subsurface fractured porous media. Specifically, the PhreeqcRM interface (developed by the USGS) is used to take advantage of a large library of equilibrium and kinetic aqueous and fluid-rock reactions, which has been validated by numerous experiments and benchmark studies. Fluid flow is modeled by the Mixed Hybrid Finite Element (FE) method, which provides smooth velocity fields even in highly heterogenous formations with discrete fractures. A multilinear Discontinuous Galerkin FE method is used to solve the multicomponent transport problem. This method is locally mass conserving and its second order convergence significantly reduces numerical dispersion. In terms of thermodynamics, the aqueous phase is considered as a compressible fluid and its properties are derived from a Cubic Plus Association (CPA) equation of state. The new simulator is validated against several benchmark problems (involving, e.g., Fickian and Nernst-Planck diffusion, isotope fractionation, advection-dispersion transport, and rock-fluid reactions) before demonstrating the expanded capabilities offered by the underlying FE foundation, such as high computational efficiency, parallelizability, low numerical dispersion, unstructured 3D gridding, and discrete fraction modeling.

58 GEOSCIENCES↗
DOE OSTI2025

Inference of phase field fracture models

The phase field approach to modeling fracture uses a diffuse damage field to represent cracks. This representation mollifies singularities that arise in computations with sharp interface models and some of the resultant difficulties in the mathematical and numerical treatment of fracture. Phase field fracture models have proven effective at representing crack propagation, branching, and merging. Specific formulations, beginning with brittle fracture, have also been shown to converge to classical solutions. Extensions to cover the range of material failure, including ductile and cohesive fracture, lead to an array of possible models. There exists a large body of literature focusing on this class of models and on the impact of model form on the predicted crack evolution. However, there have not been systematic studies into how optimal models may be chosen. Here, we take a first step in this direction by developing formal methods for identification of the best parsimonious model of phase field fracture given full-field data on the damage and deformation fields. We consider some of the main models that have been used for the degradation of elastic response due to damage and its propagation. Our approach builds upon Variational System Identification (VSI), a weak form variant of the Sparse Identification of Nonlinear Dynamics (SINDy). Furthermore, in this first communication we focus on synthetically generated data but we also consider central issues associated with the use of experimental full-field data, such as data sparsity and noise.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Tailoring the Physicochemical Properties of Nb Thin Films via Surface Engineering Methods

The modification of surface oxide layers formed on niobium (Nb) thin films via chemical mechanical planarization (CMP) and accelerated neutral atom beam (ANAB) processing provides a promising route toward tailoring their emergent properties and performance when used as superconducting qubits. Here, in this study, we show that CMP- and ANAB-formed Nb oxides are significantly thinner and smoother than the native oxide, as revealed by transmission electron microscopy (TEM) and atomic force microscopy. Scanning TEM and energy-dispersive X-ray spectroscopy along with X-ray photoelectron spectroscopy identified an oxidation gradient within the native and surface-engineered oxides. The topside layer is dominated by Nb 5+ (Nb 2 O 5 ), with various Nb suboxides present closer to the oxide/metal interface. Time-of-flight secondary ion mass spectrometry (ToF-SIMS) depth profiling confirmed the presence of an oxygen content gradient and demonstrated the enhanced resistance of the CMP- and ANAB-formed oxides to oxygen surface exchange and subsequent diffusion via 18 O 2 isotopic labeling experiments. ToF-SIMS also identified an interfacial layer containing trapped hydrogen (H)-containing species at the Nb oxide/metal interface. In situ ToF-SIMS and TEM revealed migration of the H/OH interfacial layer coinciding with decomposition of the surface oxide. Furthermore, our density functional theory calculations indicated that both H from moisture present in ambient air and bulk H in Nb films tend to segregate at the interface. These findings underscore the importance of understanding surface oxidation mechanisms, hydrogen incorporation, and their impact on the designed functionalities of Nb-based devices.

36 MATERIALS SCIENCE↗
NASA NTRS1992

Study of the oxygen transport through Ag (110), Ag (poly), and Ag 2.0 Zr

The transport of oxygen through high-purity membranes of Ag (110), Ag (poly), Ag (nano), and Ag 2.0 Zr has been studied by an ultrahigh vacuum permeation method over the temperature range of 400-800 C. The data show that there are substantial deviations from ordinary diffusion-controlled transport. A surface limitation has been confirmed by glow-discharge studies where the upstream O2 supply has been partially converted to atoms, which, for the same temperature and pressure, gave rise to over an order of magnitude increase in transport flux. Further, the addition of 2.0 wt percent Zr to the Ag has provided increased dissociative adsorption rates, which, in turn, increased the transport flux by a factor of 2. It was also observed that below a temperature of 630 C, the diffusivity exhibits an increase in activation energy of over 4 kcal/mol, which has been attributed to trapping of the atomic oxygen and/or kinetic barriers at the surface and subsurface of the vacuum interface. Above 630 C, the activation barrier decreases to the accepted value of about 11 kcal/mol for Ag (poly), consistent with zero concentration at the vacuum interface.

Outlaw, R. A.↗
DOE OSTI2023

SA508 low alloy steel to 316L stainless steel dissimilar metal joint made by powder metallurgy hot isostatic pressing

Joining ferritic SA508 low alloy steel (LAS) and austenitic 316 L stainless steel (SS) via powder metallurgy hot isostatic pressing (PM-HIP) was evaluated as an alternative method to welding. This study investigated the mechanical and microstructural evolutions of the bimetallic interface under different joint designs and heat treatments. Here, the direct joining of dissimilar metal alloys by PM-HIP method resulted in two designs: 1) powder SA508 to solid bar 316 L (P508–B316L) and 2) powder 316 L solid bar SA508 (P316L-B508). In both cases, P508–B316L and P316L-demonstrated satisfactory tensile strength, however, high hardness and severe depreciation in toughness were located on the bimetallic interface. The mechanisms responsible for the detrimental mechanical properties were verified. Large oxides were observed in P508–B316L due to the prior powder boundary (PPB) oxides present in SA508 powder. The intense sensitization occurred from the formation of M 23 C 6 carbides, consequently from the slow cooling after PM-HIP in P316L-B508. Post-HIP heat treatments were explored to reduce the distance of carbide formation; however, the heat treatment could not eliminate the carbides. The experimental results were compared to the diffusion couple simulation as a function of carbide formation with distance. The analysis also showed the high hardness at the bimetallic interface was primarily contributed by solid solution strengthening. In conclusion, the direct joining of P316L-B508 and P508–B316L via PM-HIP was deemed to be unfeasible, and a transitional material is necessary to impede the diffusion of carbon.

316L stainless Steel↗
NASA NTRS1995

Method for fabricating semiconductor devices

A process for fabricating gold/gallium arsenide structures, in situ, on molecular beam epitaxially grown gallium arsenide. The resulting interface proves to be Ohmic, an unexpected result which is interpreted in terms of increased electrode interdiffusion. More importantly, the present invention surprisingly permits the fabrication of Ohmic contacts in a III-V semiconductor material at room temperature. Although it may be desireable to heat the Ohmic contact to a temperature of, for example, 200 degrees Centigrade if one wishes to further decrease the resistance of the contact, such low temperature annealing is much less likely to have any deleterious affect on the underlying substrate. The use of the term in situ herein, contemplates continuously maintaining an ultra-high vacuum, that is a vacuum which is at least 10.sup.-8 Torr, until after the metallization has been completed. An alternative embodiment of the present invention comprising an additional step, namely the termination of the gallium arsenide by a two monolayer thickness of epitaxial aluminum arsenide as a diffusion barrier, enables the recovery of Schottky barrier behavior, namely a rectified I-V characteristic. The present invention provides a significant breakthrough in the fabrication of III-V semiconductor devices wherein excellent Ohmic contact and Schottky barrier interfaces to such devices can be achieved simply and inexpensively and without requiring the high temperature processing of the prior art and also without requiring the use of exotic high temperature refractory materials as substitutes for those preferred contact metals such as gold, aluminum and the like.

Kaiser, William J.↗
DOE OSTI2024

Low-energy radiative backgrounds in CCD-based dark-matter detectors

The reach of sub-GeV dark-matter detectors is at present severely affected by low-energy events from various origins. We present the theoretical methods to compute the single- and few-electron events that arise from secondary radiation emitted by high-energy particles as they pass through detector materials and perform a detailed simulation to quantify them at (Skipper) CCD-based experiments, focusing on the SENSEI data collected at Fermilab near the MINOS cavern. The simulations account for the generation of secondaries from Cherenkov and luminescent recombination radiation; photo-absorption in the bulk, backside layer, pitch adapter, and epoxy; the photon reflection and refraction at interfaces; thin-film interference; the roughness of the interfaces; the dynamics of charges produced in the highly doped CCD-backside-layers; and the partial charge collection on the CCD backside. We consider several systematic uncertainties, notably those stemming from the backside modeling, which we estimate with a “fiducial” and an “extreme” charge-diffusion model, with the former model being preferred due to better agreement with partial-charge collection data. We find that Cherenkov photons constitute about 30% of the observed single-electron events for both diffusion models; radiative recombination contributes negligibly to the event rate for the fiducial model, although it can dominate over Cherenkov for the extreme model. We also estimate the fraction of 2-electron events that arise from 1-electron event coincidences in the same pixel, finding that the entire 2-electron rate can be explained by coincidences of radiative events and spurious charge. Accounting for both radiative and non-radiative backgrounds, we project the sensitivity of future Skipper-CCD-based experiments to different dark-matter models. For light-mediator models with dark-matter masses of 1, 5, and 10 MeV, we find that future experiments with 10-kg-year exposures and successful background mitigation could have a sensitivity that is larger by 9, 3, and 2 orders of magnitude, respectively, when compared to an experiment without background improvements.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1991

Hydrogen no-vent fill testing in a 1.2 cubic foot (34 liter) tank

Experimental results of no-vent fill testing with liquid hydrogen in a 1.2 cubic foot (34 liter) stainless steel tank are presented. More than 40 tests were performed with various liquid inlet temperatures, inlet flowrates, initial tank wall temperatures, and liquid injection techniques. Fill levels equal to or exceeding 90 percent by volume were achieved in 40 percent of the tests with the tank pressure limited to a maximum of 30 psia. Three liquid injection techniques were employed; top spray, upward pipe discharge, and bottom diffuser. Effects of each of the varied parameters on the tank pressure history and final fill level are evaluated. The final fill level is found to be indirectly proportional to the initial wall and inlet liquid temperatures and directly proportional to the inlet liquid flowrate. Furthermore, the top spray is the most efficient no-vent fill method of the three configurations examined. The success of this injection method is primarily due to condensation of the ullage vapor onto the incoming liquid droplets. Ullage condensation counteracts the tank pressure rise resulting from energy exchange between the fluid and the warmer tank walls, and ullage compression. Upward pipe discharge from the tank bottom is the next most efficient method. Fluid circulation induced by this fill configuration tends to diminish thermal stratification in the bulk liquid, thus enhancing condensation at the liquid gas interface.

Moran, Matthew E.↗
DOE OSTI2023

Boosting the Low-Temperature Performance of Graphite Anodes by Creating an Electrochemically Active Interface

Graphite is the major anode material used in commercial lithium-ion batteries (LIBs). However, the sluggish ion-transfer kinetics associated with graphite anodes significantly restrict the operation of LIBs over a wide temperature. This is primarily due to their low reversible capacity and the substantial overpotential exhibited under low-temperature conditions. To address this limitation, we demonstrate herein an approach that involves grafting an electrochemically active lithium benzenesulfonate layer onto a graphite surface through a typical reduction reaction of diazonium cations, followed by ion exchange process. This surface modification reduces the charge transfer resistance of graphite anodes, leading to an excellent reversible capacity of ~150 mAh g –1 at low-temperatures (-20 °C, 0.1C). Electrochemical impedance spectroscopy indicates that both desolvation of the lithium ions outside the graphite, and lithium diffusion within the solid electrolyte interphase and graphite lattice are two crucial rate-limiting steps during the Li (de)lithiation, with the latter dominating during the low-temperature operation. In conclusion, these findings demonstrate a facile method for enhancing the low-temperature performance of graphite through surface modification and provide valuable insights into fundamental understandings that can guide the future design of better -low-temperature graphite anodes.

25 ENERGY STORAGE↗
NASA NTRS2021

SNPP VIIRS Day Night Band: Ten Years of On-Orbit Calibration and Performance

Aboard the polar-orbiting SNPP satellite, the VIIRS instrument has been in operation since launch in October 2011. It is a visible and infrared radiometer with a unique panchromatic channel capability designated as a day-night band (DNB). This channel covers wavelengths from 0.5 to 0.9 µm and is designed with a near-constant spatial resolution for Earth observations 24 h a day. The DNB operates at 3 gain stages (low, middle, and high) to cover a large dynamic range. An onboard solar diffuser (SD) is used for calibration in the low gain stage, and to enable the derivation of gain ratios between the different stages. In this paper, we present the SNPP VIIRS DNB calibration performed by the NASA VIIRS characterization support team (VCST). The DNB calibration algorithms are described to generate the calibration coefficient look up tables (LUTs) for the latest NASA Level 1B Collection 2 products. We provide an evaluation of DNB on-orbit calibration performance. This activity supports the NASA Earth science community by delivering consistent VIIRS sensor data products via the Land Science Investigator-led Processing Systems, including the SD degradation applied for DNB calibrations in detector gain and gain ratio trending. The DNB stray light contamination and its correction are highlighted. Performance validations are presented using comparisons to the calibration methods employed by NOAA’s operational Interface Data Processing Segment. Further work on stray light corrections is also discussed.

VIIRS↗
NASA NTRS1992

Hydrogen leak detection in the Space Shuttle

This study focuses on a helium gas jet flowing into room air. Measurements of helium concentration and velocity in the jet-air mixture are reported. The objective is to learn about jet characteristics so that dynamically similar hydrogen leaks may be located in the Space Shuttle. The hazardous gas detection system (HGDS) in the mobile launch pad uses mass spectrometers to monitor the shuttle environment for leaks. The mass spectrometers are fed by long sample tubes which draw gas from the payload bay, mid body, aft engine compartment and external tank. The overall purpose of this study is to improve the HGDS especially in its potential for locating hydrogen leaks. A rapid-response leak detection experiment was designed, built, and tested, following on the work done in this program last summer. The apparatus included a Perkin Elmer MGA-1200 mass spectrometer and air velocity transducer, both monitored by a Macintosh IIFX computer using LabVIEW software. A jet of helium flowing into the lab air simulated a gas leak. Steady helium or hydrogen-nitrogen jets were logged for concentration and velocity, and the power spectral density of each was computed. Last year, large eddies and vortices were visually seen with Schlieren imaging, and they were detected in the time plots of the various instruments. The response time of the MGA-1200 was found in the range of 0.05 to 0.1 sec. Pulsed concentration waves were clearly detected at 25 cycles per sec by spectral analysis of MGA data. No peaks were detected in the power spectrum, so in the present study, 10 Hz bandwidth-averaged power levels were examined at regular frequency intervals. The practical consequences of last year's study are as follows: sampling frequency should be increased above the present rate of 1 sample per second so that transients could be observed and analyzed with frequency response methods. Many more experiments and conditions were observed in this second summer, including the effects of orifice diameter, jet velocity, sample tube design, radial effects, vertical flow, and low hydrogen concentrations (1 percent). A frequent observation was that the power spectrum, calculated from the Fourier transform of concentration fluctuations, gives a separate piece of information from concentration. Many of the tests suggest that power is high where mixing occurs at the helium-air interface. This fact is apparently independent of the concentration level, which could be high or low, but depends on the sample location relative to the jet (leak) origin, whereas high concentration may be due to a strong leak far away or a small leak close to the sample tube. If the power is low for any concentration level, this would signify helium is arriving at the sample tube by diffusion, not chaotic mixing caused by the jet interaction with air. The practical result is to propose a modification of the HGDL mass spectrometer data sampling and software so that sampling rates could be capable of observing at least 25 Hz fluctuations.

Barile, Ronald G↗
DOE OSTI2023

Neural Network Enhanced RKPM for Electrochemical-Mechanical Coupled Damage Modeling of Energy Storage Materials

Energy storage materials undergo significant charge cycling, which makes understanding their reliability and durability fundamental in predicting performance and service life. Strong electrochemical-mechanical coupling and highly anisotropic material properties contribute to the formation and propagation of micro-cracking, largely along material interfaces and grain boundaries. For Li-ion batteries, for example, lithium moving between electrodes during charging and discharging process causes expansion and contraction of grains, and the strongly anisotropic and nonlinearly [Li]-dependent grain material properties can cause grains to expand into and contract away from each other, leading to chemo-mechanical cracking. In the first part of this work, a RKPM based computational framework for solving the coupled solid-phase lithium conservation with Fickian diffusion and the lithium concentration dependent anisotropic mechanical problem subjected to a highly nonlinear Butler-Volmer boundary condition is introduced. The choice of RKPM completeness conditions for lithium concentration and mechanical deformation fields, and the variational consistency condition for the domain integration of the coupled problem is first determined. In the second part of this work, a neural network-enhanced reproducing kernel particle method (NN-RKPM) [1] is leveraged to accurately capture damage and crack propagation throughout the material, by learning the location, orientation, and sharpness of discontinuity while allowing for a coarser nodal distribution than that is necessary for capturing sharp solution transitions using traditional mesh-based methods. NN-RKPM is used to inform how crack opening and closure in turn affect the coupled chemical equations and material microstructure.

damage modeling↗
NASA NTRS1999

Step Bunching: Influence of Impurities and Solution Flow

Step bunching results in striations even at relatively early stages of its development and in inclusions of mother liquor at the later stages. Therefore, eliminating step bunching is crucial for high crystal perfection. At least 5 major effects causing and influencing step bunching are known: (1) Basic morphological instability of stepped interfaces. It is caused by concentration gradient in the solution normal to the face and by the redistribution of solute tangentially to the interface which redistribution enhances occasional perturbations in step density due to various types of noise; (2) Aggravation of the above basic instability by solution flowing tangentially to the face in the same directions as the steps or stabilization of equidistant step train if these flows are antiparallel; (3) Enhanced bunching at supersaturation where step velocity v increases with relative supersaturation s much faster than linear. This v(s) dependence is believed to be associated with impurities. The impurities of which adsorption time is comparable with the time needed to deposit one lattice layer may also be responsible for bunching; (4) Very intensive solution flow stabilizes growing interface even at parallel solution and step flows; (5) Macrosteps were observed to nucleate at crystal corners and edges. Numerical simulation, assuming step-step interactions via surface diffusion also show that step bunching may be induced by random step nucleation at the facet edge and by discontinuity in the step density (a ridge) somewhere in the middle of a face. The corresponding bunching patterns produce the ones observed in experiment. The nature of step bunching generated at the corners and edges and by dislocation step sources, as well as the also relative importance and interrelations between mechanisms 1-5 is not clear, both from experimental and theoretical standpoints. Furthermore, several laws controlling the evolution of existing step bunches have been suggested, though unambiguous conclusions are still missing. Addressing these issues is the major goal of the present project. The theory addressing the above problem, experimental methods, several figures which include: (1) the spatial wave numbers at which the system is neutrally stable as a function of growth velocity for linear kinetics and supersaturation for nonlinear kinetics; (2) a schematic of the experiment of lysozyme crystal growing under conditions of natural convection; (3) fluctuations in time, t, of the normal growth rate, R(t), vicinal slope, p(t) and Fourier Spectra of R(t), discussions and conclusions are presented.

Chernov, A. A.↗
DOE OSTI2021

Analysis of a new multispecies tumor growth model coupling 3D phase-fields with a 1D vascular network

In this study, we present and analyze a mathematical model for tumor growth incorporating ECM erosion, interstitial flow, and the effect of vascular flow and nutrient transport. The model is of phase-field or diffused-interface type in which multiple phases of cell species and other constituents are separated by smooth evolving interfaces. The model involves a mesoscale version of Darcy’s law to capture the flow mechanism in the tissue matrix. Modeling flow and transport processes in the vasculature supplying the healthy and cancerous tissue, one-dimensional (1D) equations are considered. Since the models governing the transport and flow processes are defined together with cell species models on a three-dimensional (3D) domain, we obtain a 3D–1D coupled model.

3D–1D coupled blood flow models↗
NASA NTRS2001

2-D Modeling of Nanoscale MOSFETs: Non-Equilibrium Green's Function Approach

We have developed physical approximations and computer code capable of realistically simulating 2-D nanoscale transistors, using the non-equilibrium Green's function (NEGF) method. This is the most accurate full quantum model yet applied to 2-D device simulation. Open boundary conditions and oxide tunneling are treated on an equal footing. Electrons in the ellipsoids of the conduction band are treated within the anisotropic effective mass approximation. Electron-electron interaction is treated within Hartree approximation by solving NEGF and Poisson equations self-consistently. For the calculations presented here, parallelization is performed by distributing the solution of NEGF equations to various processors, energy wise. We present simulation of the "benchmark" MIT 25nm and 90nm MOSFETs and compare our results to those from the drift-diffusion simulator and the quantum-corrected results available. In the 25nm MOSFET, the channel length is less than ten times the electron wavelength, and the electron scattering time is comparable to its transit time. Our main results are: (1) Simulated drain subthreshold current characteristics are shown, where the potential profiles are calculated self-consistently by the corresponding simulation methods. The current predicted by our quantum simulation has smaller subthreshold slope of the Vg dependence which results in higher threshold voltage. (2) When gate oxide thickness is less than 2 nm, gate oxide leakage is a primary factor which determines off-current of a MOSFET (3) Using our 2-D NEGF simulator, we found several ways to drastically decrease oxide leakage current without compromising drive current. (4) Quantum mechanically calculated electron density is much smaller than the background doping density in the poly silicon gate region near oxide interface. This creates an additional effective gate voltage. Different ways to. include this effect approximately will be discussed.

Svizhenko, Alexei↗
NASA NTRS1998

Annual Research Briefs, 1998

The topics contained in this progress report are direct numerical simulation of turbulent non-premixed combustion with realistic chemistry; LES of non-premixed turbulent reacting flows with conditional source term estimation; measurements of the three-dimensional scalar dissipation rate in gas-phase planar turbulent jets; direct simulation of a jet diffusion flame; on the use of interpolating wavelets in the direct numerical simulation of combustion; on the use of a dynamically adaptive wavelet collocation algorithm in DNS (direct numerical simulation) of non-premixed turbulent combustion; 2D simulations of Hall thrusters; computation of trailing-edge noise at low mach number using LES and acoustic analogy; weakly nonlinear modeling of the early stages of bypass transition; interactions between freestream turbulence and boundary layers; interfaces at the outer boundaries of turbulent motions; largest scales of turbulent wall flows; the instability of streaks in near-wall turbulence; an implementation of the v(sup 2) - f model with application to transonic flows; heat transfer predictions in cavities; a structure-based model with stropholysis effects; modeling a confined swirling coaxial jet; subgrid-scale models based on incremental unknowns for large eddy simulations; subgrid scale modeling taking the numerical error into consideration; towards a near-wall model for LES of a separated diffuser flow; on the feasibility of merging LES with RANS (Reynolds Averaging Numerical simulation) for the near-wall region of attached turbulent flows; large-eddy simulation of a separated boundary layer; numerical study of a channel flow with variable properties; on the construction of high order finite difference schemes on non-uniform meshes with good conservation properties; development of immersed boundary methods for complex geometries; and particle methods for micro and macroscale flow simulations.

Spinks, Debra↗
DOE OSTI2021

Toward predictive permeabilities: Experimental measurements and multiscale simulation of methanol transport in Nafion

A polymer membrane's permeability to solutes determines its suitability for various applications: a permeability value is essential for predicting performance in diverse contexts. Using aqueous methanol permeation through Nafion as an example, we describe a methodology for determining membrane permeability that accounts for boundary layer effects and the possibility of swelling. For the materials and apparatus used herein, analysis of a permeance measurement and computational fluid dynamics simulations show that the mass transfer boundary layer is on the order of ones to tens of microns. Additionally, the data are used to develop and validate a multiscale model describing solute permeation through a hydrated membrane as a series of physical mechanistic steps: reversible adsorption from solution at the membrane interface, diffusion driven by a concentration gradient within the membrane, and reversible desorption into solution at the opposite membrane interface. The validated model is used to predict methanol transport across a solar-driven CO 2 reduction device and to assess the impact of polymer changes on the measured value. The approach of combining experimental data, computational fluid dynamics, and the mechanistic multiscale model is expected to provide more accurate analysis of membrane permeation data in cases with polymer swelling or unusual device geometries, among others.

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

Regulating hydrogel mechanical properties with an electric field

Stimuli-responsive polymeric materials have attracted significant attention due to their ability to change properties in response to various external stimuli. Using an electric field as the stimulus is of particular interest as it possesses the potential for seamless integration of materials with electronic systems. While many materials with electric field responsive actuation have an associated mechanical property change, it is beneficial to develop materials that exhibit mechanical property changes without accompanying significant shape deformation. To address this challenge, here we designed a semi-interpenetrating polymer network (semi-IPN) hydrogel system containing both polyelectrolytes and salt ions, which enables electric field induced changes in mechanical properties while minimizing actuation. We first successfully verified the viability of our design by removing salt ions through a diffusion-only method where we witnessed the stiffness increased to 4.5 times the initial value while still being highly deformable. After this, we applied an electric field to transport the salt ions out of the hydrogel, as shown by both Raman spectroscopy and scanning electron microscopy. We were able to show a time-dependent stiffness increase, the maximum of which was 5 times the original stiffness. We quantified ion transport and water-splitting in the hydrogel by both experiments and simulations. Following this, we showed functional system reversibility by reversing the direction of the current to reinject salt ions into the semi-IPN hydrogel and reducing its stiffness to below the initial value. It's worth noting that our simulations enable us to understand the governing mechanisms behind ion generation and salt transport that leads to mechanical property changes. Finally, we were able to fabricate a spatially variable stiffness haptic interface with our hydrogel, with demonstrated reversibility and cyclability. This research can possibly find applications in soft robotics and haptics and also inspire the development of bio-compatible electronics related devices.

36 MATERIALS SCIENCE↗