SEARCH · Engineering Papers

Results for “enhanced kinetics”

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 1,459 records · Page 81

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Design and Analysis of Liquid Lithium Plasma Facing Components

In this study, liquid metal (LM) on plasma-facing components (PFCs) creates a renewable protective cover enhancing power exhaust and protecting the solid substrate. Additionally, energy confinement improvement through particle pumping can be achieved in the case of liquid lithium (LL). A new variant of PFCs was recently introduced at Princeton Plasma Physics Laboratory (PPPL) where a porous wall is used to stabilize the LM surface, while using magnetohydrodynamic (MHD) drive to push the LM flow inside the component. This arrangement allows efficient heat exhaust, and enhanced control of the LM surface temperature, leading to spatial control of evaporation and condensation of LL on the plasma interface. This feature is particularly attractive when vapor shielding is introduced to allow heat flux redistribution. This system has the advantage that as the heat flux increases, the evaporation rate will increase, while decreasing heat flux will decrease the evaporation rate of lithium, ideally creating a feedback effect which could self-regulate the amount of lithium evaporated. Analytical and numerical models for LL PFCs were developed in PPPL. To calculate the target temperature distribution for the case of evaporation from the divertor, we apply an iterative process allowing two-way coupling between the fluid-kinetic analysis of plasma using SOLPS-ITER code and the flow and heat transfer analysis of the PFC using an analytical model. In the last stage, the results are validated using computational fluid dynamics (CFD) analysis with customized version of the CFX code from ANSYS. CFX was modified at PPPL to allow MHD analysis at the high magnetic field typical for fusion applications. Results of the analysis for NSTX-U tokamak conditions will be presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1996

Progress Toward an Efficient and General CFD Tool for Propulsion Design/Analysis

The simulation of propulsive flows inherently involves chemical activity. Recent years have seen substantial strides made in the development of numerical schemes for reacting flowfields, in particular those involving finite-rate chemistry. However, finite-rate calculations are computationally intensive and require knowledge of the actual kinetics, which are not always known with sufficient accuracy. Alternatively, flow simulations based on the assumption of local chemical equilibrium are capable of obtaining physically reasonable results at far less computational cost. The present study summarizes the development of efficient numerical techniques for the simulation of flows in local chemical equilibrium, whereby a 'Black Box' chemical equilibrium solver is coupled to the usual gasdynamic equations. The generalization of the methods enables the modelling of any arbitrary mixture of thermally perfect gases, including air, combustion mixtures and plasmas. As demonstration of the potential of the methodologies, several solutions, involving reacting and perfect gas flows, will be presented. Included is a preliminary simulation of the SSME startup transient. Future enhancements to the proposed techniques will be discussed, including more efficient finite-rate and hybrid (partial equilibrium) schemes. The algorithms that have been developed and are being optimized provide for an efficient and general tool for the design and analysis of propulsion systems.

Cox, C. F.↗

NASA NTRS2000

NCC: A Physics-Based Design and Analysis Tool for Combustion Systems

The National Combustion Code (NCC) is an integrated system of computer codes for physics-based design and analysis of combustion systems. It uses unstructured meshes and runs on parallel computing platforms. The NCC is composed of a set of distinct yet closely related modules. They are: (1) a gaseous flow module solving 3-D Navier-Stokes equations; (2) a turbulence module containing the non-linear k-epsilon models; (3) a chemistry module using either the conventional reduced kinetics approach of solving species equations or the Intrinsic Low Dimensional Manifold (ILDM) kinetics approach of table looking up in conjunction with solving the equations of the progressive variables; (4) a turbulence-chemistry interaction module including the option of solving the joint probability density function (PDF) for species and enthalpy; and (5) a spray module for solving the liquid phase equations. In early 1995, an industry-government team was formed to develop the NCC. In July 1998, the baseline beta version was completed and presented in two NCC sessions at the 34th AIAA/ASME/SAE/ASEE Joint Propulsion Conference & Exhibit, July 1998. An overview of this baseline beta version was presented at the NASA HPCCP/CAS Workshop 98, August 1998. Since then, the effort has been focused on the streamlining, validation, and enhancement of the th baseline beta version. The progress is presented in two NCC sessions at the AIAA 38 Aerospace Sciences Meeting & Exhibit, January 2000. At this NASA HPCCP/CAS Workshop 2000, an overview of the NCC papers presented at the AIAA 38 th Aerospace Sciences Meeting & Exhibit is presented, with emphasis on the reduction of analysis time of simulating the (gaseous) reacting flows in full combustors. In addition, results of NCC simulation of a modern turbofan combustor will also be reported.

Liu, Nan-Suey↗

DOE OSTI2025

Synergistic torrefaction of plastic polymers and biomass

This study explores the synergistic torrefaction of biomass and plastics, aimed at enhancing bioenergy production and promoting a circular economy. By leveraging the unique properties of both materials, we investigated the thermochemical transformations occurring during the torrefaction process, from material preparation to the final characteristics of the torrefied product. The biomass used included corn stover (CS) and loblolly pine (LP), while various plastics were categorized from #1 to #7. Torrefaction was conducted at temperatures of 200, 225, and 250 °C, with subsequent extrusion of the torrefied materials and raw materials to produce composite filaments. The results show a consistent decrease in mass yield with increased torrefaction temperature, with notable variations among different biomass-plastic combinations. Co-torrefaction of biomass with polyvinyl chloride and polypropylene resulted in accelerated reaction kinetics, with an observed mass loss rate increase of 15 % at 250 °C compared to the expected rates for individual components. This synergy was quantified, indicating a 20.3 % increase in mass loss for the loblolly pine-polypropylene combination and 23.9 % for corn stover-polypropylene. In contrast, other plastics, including polyethylene terephthalate, high-density polyethylene, low-density polyethylene, polystyrene, and polycarbonate, did not exhibit significant synergistic effects. Mechanical testing indicated that the torrefaction process alters the strength and brittleness of the resulting materials, with implications for their application in bioenergy production and bio-renewable materials. Overall, this research highlights the potential of synergistic torrefaction as a viable strategy for co-processing biomass and plastics, paving the way for innovative solutions in waste management and renewable energy resource development.

09 - BIOMASS FUELS↗

DOE OSTI2025

Three-dimensional core-collapse supernova models with phenomenological treatment of neutrino flavor conversions

Abstract We perform three-dimensional supernova simulations with a phenomenological treatment of neutrino flavor conversions. We show that the explosion energy can increase to as high as $\sim 10^{51}$ erg depending on the critical density for the onset of flavor conversions, due to a significant enhancement of the mean energy of electron antineutrinos. Our results confirm previous studies showing such energetic explosions, but for the first time in three-dimensional configurations. In addition, we predict neutrino and gravitational wave (GW) signals from a nearby supernova explosion aided by flavor conversions. We find that the neutrino event number decreases because of the reduced flux of heavy-lepton neutrinos. In order to detect GWs, next-generation GW telescopes such as Cosmic Explorer and the Einstein Telescope are needed even if the supernova event is located at the Galactic Center. These findings show that the neutrino flavor conversions can significantly change supernova dynamics and highlight the importance of further studies on the quantum kinetic equations to determine the conditions of the conversions and their asymptotic states.

Mori, Kanji↗

DOE OSTI2024

Experimental and Computational Study of Pyrogenic Carbonaceous Matter Facilitated Hydrolysis of 2,4-Dinitroanisole (DNAN)

This study investigated the reaction pathway of 2,4-dinitroanisole (DNAN) on the pyrogenic carbonaceous matter (PCM) to assess the scope and mechanism of PCM-facilitated surface hydrolysis. DNAN degradation was observed at pH 11.5 and 25 °C with a model PCM, graphite, whereas no significant decay occurred without graphite. Experiments were performed at pH 11.5 due to the lack of DNAN decay at pH below 11.0, which was consistent with previous studies. Graphite exhibited a 1.78-fold enhancement toward DNAN decay at 65 °C and pH 11.5 relative to homogeneous solution by lowering the activation energy for DNAN hydrolysis by 54.3 ± 3.9%. This is supported by our results from the computational modeling using Car–Parrinello simulations by ab initio molecular dynamics/molecular mechanics (AIMD/MM) and DFT free energy simulations, which suggest that PCM effectively lowered the reaction barriers by approximately 8 kcal mol –1 compared to a homogeneous solution. Quaternary ammonium (QA)-modified activated carbon performed the best among several PCMs by reducing DNAN half-life from 185 to 2.5 days at pH 11.5 and 25 °C while maintaining its reactivity over 10 consecutive additions of DNAN. We propose that PCM can affect the thermodynamics and kinetics of hydrolysis reactions by confining the reaction species near PCM surfaces, thus making them less accessible to solvent molecules and creating an environment with a weaker dielectric constant that favors nucleophilic substitution reactions. Nitrite formation during DNAN decay confirmed a denitration pathway, whereas demethylation, the preferred pathway in homogeneous solution, produces 2,4-dinitrophenol (DNP). Denitration catalyzed by PCM is advantageous to demethylation because nitrite is less toxic than DNAN and DNP. These findings provide critical insights for reactive adsorbent design that has broad implications for catalyst design and pollutant abatement.

2,4-dinitroanisole (DNAN)↗

DOE OSTI2024

Active doping controls the mode of failure in dense colloidal gels

Mechanical properties of disordered materials are governed by their underlying free energy landscape. In contrast to external fields, embedding a small fraction of active particles within a disordered material generates nonequilibrium internal fields, which can help to circumvent kinetic barriers and modulate the free energy landscape. In this work, we investigate through computer simulations how the activity of active particles alters the mechanical response of deeply annealed polydisperse colloidal gels. We show that the “swim force” generated by the embedded active particles is responsible for determining the mode of mechanical failure, i.e., brittle vs. ductile. We find, and theoretically justify, that at a critical swim force the mechanical properties of the gel decrease abruptly, signaling a change in the mode of mechanical failure. The weakening of the elastic modulus above the critical swim force results from the change in gel porosity and distribution of attractive forces among gel particles, while below the critical swim force, the ductility enhancement is caused by an increase of gel structural disorder. Above the critical swim force, the gel develops a pronounced heterogeneous structure characterized by multiple pore spaces, and the mechanical response is controlled by dynamical heterogeneities. We contrast these results with those of a simulated monodisperse gel that exhibits a nonmonotonic trend of ductility modulation with increasing swim force, revealing a complex interplay between the gel energy landscape and embedded activity.

Zhou, Tingtao (ORCID:000000021766719X)↗

NASA NTRS1972

The effect of alloy composition on the mechanism of stress-corrosion cracking of titanium alloys in aqueous environments

A detailed study was made of the relation between the size distribution of Ti3Al particles in a Ti-8Al alloy and the tensile properties measured in air and in saltwater. The size distribution of Ti3Al was varied by isothermal aging for various times at temperatures in the range 770 to 970 K (930 to 1290 F). The aging kinetics were found to be relatively slow. Quantitative measurements of the particle coarsening rate at 920 K (1200 F) showed good agreement with the predicted behavior for coarsening controlled by matrix diffusion, and suggested that the specific free energy of the Ti3Al alpha interface in negligible small. In all cases, the Ti3Al particles were sheared by the glide dislocations. It was concluded that there is a definite correlation between the presence of deformable Ti3Al particles and an alloy's susceptibility to aqueous stress corrosion cracking. Furthermore, the appearance of the surface slip lines and the dislocation substructure in deformed specimens suggest that the specific effect of the Ti3Al particles is to cause a nonhomogeneous planar slip character and an enhanced chemical potential of the slip bands.

Wood, R. A.↗

NASA NTRS2005

Numerical Study of High-Temperature Jet Flow Using RANS/LES and PANS Formulations

Two multi-scale-type turbulence models are implemented in the PAB3D solver. The models are based on modifying the Reynolds Averaged Navier-Stokes (RANS) equations. The first scheme is a hybrid RANS/LES model utilizing the two-equation (k(epsilon)) model with a RANS/LES transition function dependent on grid spacing and the computed turbulence length scale. The second scheme is a modified version of the Partially Averaged Navier-Stokes (PANS) model, where the unresolved kinetic energy parameter (f(sub k)) is allowed to vary as a function of grid spacing and the turbulence length scale. This parameter is estimated based on a novel two-stage procedure to efficiently estimate the level of scale resolution possible for a given flow on a given grid for Partial Averaged Navier-Stokes (PANS). It has been found that the prescribed scale resolution can play a major role in obtaining accurate flow solutions. The parameter f(sub k) varies between zero and one and equal to one in the viscous sub layer, and when the RANS turbulent viscosity becomes smaller than the LES viscosity. The formulation, usage methodology, and validation examples are presented to demonstrate the enhancement of PAB3D's time-accurate and turbulence modeling capabilities. The accurate simulations of flow and turbulent quantities will provide valuable tool for accurate jet noise predictions. Solutions from these models are compared to RANS results and experimental data for high-temperature jet flows. The current results show promise for the capability of hybrid RANS/LES and PANS in simulating such flow phenomena.

Abdol-Hamid, Khaled S.↗

DOE OSTI2025

Structure–Function Relationships in Sequence-Controlled Copolymers for Rare Earth Element Chelation

The ability to tune material function through primary sequence is a defining feature of biological macromolecules, allowing precise control over structure and target interactions in complex aqueous environments. However, translating sequence–structure–function relationships to synthetic macromolecules is challenging due to their dispersity in sequence, conformation, and composition. Here, we report systematic studies of amphiphilic polymer chelators designed to probe how composition and patterning influence binding affinity and selectivity for rare earth elements (REEs), a series of technologically relevant metals with challenging separation profiles. A library of copolymers varying hydrophobic monomer composition and patterning was synthesized via reversible addition–fragmentation chain transfer (RAFT) polymerization, spanning statistical, gradient, and block architectures. REE binding was quantified using a high-throughput colorimetric assay, and reconstruction of polymer ensembles using kinetic stochastic simulations enabled quantitative comparisons of sequence heterogeneity, linking local monomer colocalization to emergent REE binding. Further, we investigated the role of different hydrophobic comonomers in tuning metal coordination, with binding trends linked to structural features that influence binding site desolvation. Complementary dynamic light scattering (DLS) and small-angle X-ray scattering (SAXS) measurements showed that both polymer and monomer architecture modulate metal-induced conformational changes, and that multichain assembly behavior emerges beyond critical hydrophobic thresholds. Sequence control also altered REE selectivity, with nonmonotonic differences observed across compositionally identical polymers with different sequence architectures. Together, these findings establish design principles that connect polymer sequence and structure to binding performance, guiding the design of macromolecular chelators with enhanced affinity and selectivity for applications in separations, sensing, and catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Reactive Transport Modeling of Hydrogen Production from Serpentinization of Olivine

Hydrogen production from serpentinization of ultramafic rocks represents a promising natural pathway for generating carbon-free energy, yet its kinetics and controlling factors remain incompletely understood. A key challenge in advancing serpentinization research lies in the heterogeneity of porosity and permeability in rocks, which leads to nonuniform fluid velocity fields, as well as uncertainties in estimating reactive surface area and identifying appropriate mineral reaction equilibria. Additional complexities arise from the role of dissolved SiO 2 , Fe 2+ /Fe 3+ partitioning, and the limited effect of pH variations within the strongly alkaline regime on hydrogen yields. These challenges hinder straightforward extrapolation from laboratory tests to practical applications of hydrogen production from natural rocks. Here, in this work, we address these questions using a simulation-based reactive transport modeling framework calibrated against controlled laboratory experiments reported elsewhere. The model couples geochemical kinetics, multiphase flow, and mineralogical feedbacks, enabling systematic evaluation of how surface area, dissolved silica concentration, Fe redox state, temperature, and pressure govern serpentinization and H2 generation. We find that surface area exerts the strongest control on reaction rates and hydrogen yields, while Fe 2+ /Fe 3+ ratios act as secondary modulators. Elevated dissolved silica concentrations suppress hydrogen production but accelerate serpentine precipitation, whereas increasing pH beyond 12 within the strongly alkaline regime produces only marginal gains. Finally, we demonstrate that integrating targeted experiments with calibrated simulations offers a powerful and efficient approach for predicting hydrogen yields and assessing parameter trade-offs in industrial-scale applications. This integration can substantially reduce the experimental burden while improving predictive capability, thereby enhancing both the mechanistic understanding and the practical feasibility of hydrogen production from serpentinization.

08 HYDROGEN↗

DOE OSTI2026

Effects of Dissolution Regimes on Flow Channelization and Solute Transport in 3D Fracture Networks: Insights From Graph‐Based Reactive Transport Modeling

We investigate how mineral dissolution reshapes flow pathways and solute transport in three‐dimensional discrete fracture networks using a computationally efficient graph‐based reactive transport model. The DFNs are inspired by field‐site observations of fractured carbonate and represent realistic connectivity and structural heterogeneity. Flow is simulated with the Reynolds equation, and dissolution follows first‐order kinetics with diffusive limitations captured through an effective mass‐transfer coefficient. By systematically varying two key dimensionless parameters, the effective Damköhler number (Da), governing reaction versus advection rates, and a transport parameter (Da), analogous to the Thiele modulus, distinct flow channelization regimes emerge: mildly channelized at low G, highly channelized at intermediate Da, and extreme wormhole formation at high Da and low G. Eulerian and Lagrangian analyses, including breakthrough curves, particle tortuosity, dispersivity, and flow channeling indicators quantitatively characterize the progression of dissolution‐driven network restructuring. Across all regimes, initial fracture heterogeneity persists. The results underscore how the interplay between this initial structure, advection, reaction, and diffusion critically shapes subsurface flow pathways, with implications for applications ranging from groundwater remediation to enhanced geothermal systems.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Finite-element-based simulations of electrodes for CO 2 cascade reduction reactions

The multielectron reduction of CO 2 to liquid fuels could be a path to scalable energy storage, but reaching this goal requires major advances in catalysis and systems engineering. Cascade catalysis, which couples sequential reactions without isolating intermediates, has emerged as a promising route to enhance selectivity and efficiency in CO 2 reduction (CO 2 R). In this review, we examine how finite-element-based simulations of continuum model [finite element method (FEM)] approaches are being used to analyze and guide CO 2 R cascade systems. We first outline the fundamentals of cascade catalysis and recent advances in catalytic materials (metallic, molecular, and hybrid architectures). We then focus on FEM developments at the electrode and device scales, emphasizing how these models capture transport phenomena, local microenvironments, and geometry-dependent effects. To clarify design principles, we present case studies of cascade electrodes organized in systems without and with integrated semiconductors. We further emphasize the integration of FEM with multiscale frameworks (density functional theory, molecular dynamics, kinetic Monte Carlo) and its role in bridging atomic-level insights with device-level performance. Finally, we identify current limitations and future prospects, including improved boundary conditions, coupling with operando experiments, and machine learning-accelerated model development. Together, these insights provide design principles for next-generation CO 2 R cascade systems for efficient solar fuel production.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

On the enhanced sulfur and coking tolerance of Ni-Co-rare earth oxide catalysts for the dry reforming of methane

In this study, sulfur and coking tolerance of Ni-based dry reforming catalysts were examined. Catalysts utilizing both Ce/Zr and Ce/La oxide supports, some with additional Co, were tested. Long-term reaction runs were conducted with and without sulfur in the feed. Catalysts were also characterized by STEM, XPS, XAFS and XANES and CO chemisorption. Only catalysts where Co was also present, and supported on the Ce-Zr oxide, were capable of extended sulfur tolerance at >20 ppm sulfur. This tolerance, along with a greatly reduced coking rate, is linked to Co in intimate contact with Ni, the mixture existing as clusters anchored and influenced electronically by the oxide support. The activation of methane takes place on these sites. Larger metal aggregates formed by ripening during reaction appear to be spectators. The measured activation energies for dry reforming suggest that CO 2 activation takes place at the oxide interface, and is a kinetically significant step.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1995

Computational study of inlet injection for a Pre-Mixed, Shock-Induced Combustion (PM/SIC) engine

A computational simulation of reacting 2-D and 3-D flowfields in a model inlet section of a Pre-Mixed, Shock-Induced Combustion (PM/SIC) engine concept was performed. LARCK, a multi-dimensional Navier-Stokes code with finite-rate kinetics chemistry developed at NASA LaRC by J.A. White, was adapted for this simulation. The flow conditions in the simulation match those envisioned for the PM/SIC engine experiments currently planned at LaRC. The reacting flowfields were Mach 6.3 freestream air and Mach 2 hydrogen at various pressure and temperature conditions injected through a slot injector at the base of the inlet section. In the PM/SIC engine, fuel is injected at the inlet section upstream of the combustor, and reaction is initiated by the shock wave at the inlet which increases the gas temperature and pressure beyond the kinetic limits for reaction. Many challenges exist prior to establishing shock-controlled combustion as a practical engine concept. These challenges include fuel injection schemes that can provide proper fuel-air mixing without creating large losses in the inlet section, and control of the combustion process so that early ignition or combustion propagation through the inlet boundary layer does not occur. For this project, a parametrics study was carried out to model the fuel injection of hydrogen at different flow conditions. It was found that, as the fuel temperature and pressure were increased, the potential for pre-ignition was high at a short distance downstream of the slot injector. The next stage of this work will investigate injection techniques for enhancing mixing of fuel and air in a manner that prevents or reduces the potential for premature ignition observed numerically.

Gonzalez, D. E.↗

DOE OSTI2024

Measurement of the double-differential cross section of muon-neutrino charged-current interactions with low hadronic energy in the NOvA Near Detector

The NOvA collaboration reports cross-section measurements for $\nu_{\mu}$ charged-current interactions with low hadronic energy (maximum kinetic energy of 250 MeV for protons and 175 MeV for pions) in the NOvA Near Detector. The results are presented as a double-differential cross section as a function of the direct observables of the final-state muon kinematics. Results are also presented as a single-differential cross section as a function of the derived square of the four-momentum transfer, $Q^{2}$, and as a function of the derived neutrino energy. The data correspond to an accumulated 8.09$\times10^{20}$ protons-on-target (POT) in the neutrino mode of the NuMI beam, with a narrow band of neutrino energies peaked at 1.8 GeV. The analysis provides a sample of neutrino-nucleus interactions with an enhanced fraction of quasi-elastic and two-particle-two-hole (2p2h) interactions. This enhancement allows quantitative comparisons with various nuclear models. We find strong disagreement between data and theory-based models in various regions of the muon kinematic phase space, especially in the forward muon direction.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2026

Stabilizing Ru-Formyl CO 2 Reduction Intermediates by Formation of Lewis Acid–Base Adducts: A Combined 2DIR and NMR Study

This work demonstrates how the stability of a Ru(II) formyl complex (cis-[Ru(bpy) 2 (CO)(CHO)][PF 6 ]), which is known to decarbonylate to form the metal hydride complex, is greatly enhanced by the addition of a Lewis acidic cation (Li + or Mg 2+ ). Multinuclear NMR is used to measure equilibrium constants of adduct formation, and then ultrafast 2DIR spectroscopy is applied to determine how the electronic structure and local molecular dynamics change upon complexation of the formyl ligand with the Lewis acid. We attribute the stabilization of the complex to the formation of a carbene-like structure where the bond strength between the Ru and the formyl group increases. The change in electronic structure is evidenced by steady-state NMR spectroscopy and DFT calculations. 2DIR measurements further confirm both the complexation of the formyl ligand with the Lewis acid and the change in electronic structure, where changes in vibrational frequency, relaxation kinetics, and a slowdown in the time scale of rotation about the metal-formyl bond are observed in the presence of Lewis acids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Ideal conductor/dielectric model (ICDM): A generalized technique to correct for finite-size effects in molecular simulations of hindered ion transport

Molecular simulations serve as indispensable tools for investigating the kinetics and elucidating the mechanism of hindered ion transport across nanoporous membranes. In particular, recent advancements in advanced sampling techniques have made it possible to access translocation timescales spanning several orders of magnitude. In our prior study [Shoemaker et al., J. Chem. Theory Comput. 18, 7142 (2022)], we identified significant finite size artifacts in simulations of pressure-driven hindered ion transport through nanoporous graphitic membranes. We introduced the ideal conductor model, which effectively corrects for such artifacts by assuming the feed to be an ideal conductor. In the present work, we introduce the ideal conductor dielectric model (Icdm), a generalization of our earlier model, which accounts for the dielectric properties of both the membrane and the filtrate. Using the Icdm model substantially enhances the agreement among corrected free energy profiles obtained from systems of varying sizes, with notable improvements observed in regions proximate to the pore exit. Moreover, the model has the capability to consider secondary ion passage events, including the transport of a co-ion subsequent to the traversal of a counter-ion, a feature that is absent in our original model. We also investigate the sensitivity of the new model to various implementation details. The Icdm model offers a universally applicable framework for addressing finite size artifacts in molecular simulations of ion transport. It stands as a significant advancement in our quest to use molecular simulations to comprehensively understand and manipulate ion transport processes through nanoporous membranes.

Chemistry↗