Two-dimensional Transport Model Research Articles

Modeling electrolytic transport for systems with concentration gradients, Ohmic resistance and electrochemical reactions

ABSTRACT This paper describes a two-dimensional multi-component electrolytic transport model that calculates the electric field by applying electroneutrality as an upper bound. This approach avoids directly enforcing electroneutrality in mass transport calculations or using Poisson’s equation. The two coupled equations of this model were numerically solved for cases with no convection. The transport equation was solved using a modified Control Volume method and a Peclet number. The electric field equation was discretized using the finite difference method and solved using the Alternating Direction Implicit method. The model’s results were compared with free-diffusion liquid junction data. Comparisons were also made with one-dimensional transport models. The model was then used to simulate two-dimensional scenarios without prescribed current distributions. The simulations agreed with the comparison data. Hence, the model shows promise in its ability to simulate two-dimensional multi-component electrolytic transport with concentration gradients, Ohmic resistance, and electrochemical reactions.

Sustainable Management of Pollutant Transport in Defective Composite Liners of Landfills: A Semi-Analytical Model

This study presents an analytical model for two-dimensional pollutant transport within a three-layer composite liner system, which comprises a geomembrane (GM), a geosynthetic clay liner (GCL), and a soil liner (SL), with particular attention to defects in the geomembrane. The model integrates key processes such as convection, diffusion, adsorption, and degradation, offering a more accurate prediction of pollutant behavior. Through Laplace and Fourier transforms, pollutant concentration distributions are derived, providing a comprehensive view of pollutant migration in landfill settings. Verification against COMSOL 6.0 simulations underscores the model’s robustness. Results show that there is an optimal thickness for the SL that balances the effectiveness of pollutant containment and material usage, while higher diffusion coefficients and advection velocity accelerate migration. The degradation of organic pollutants reduces concentrations over time, especially with shorter half-lives. These findings not only improve the design of landfill liners but also support more sustainable waste management practices by reducing the risk of environmental contamination. This research contributes to the development of more effective, long-lasting landfill containment systems, enhancing sustainability in waste management infrastructure.

Analytical model for two-dimensional contaminant transport in a cut-off wall and aquitard system

The cut-off wall is an effective countermeasure to prevent contaminant transport from a contaminated site. However, the aquitard beneath the cut-off wall may serve as a preferential path, leading to an earlier breakthrough of contaminants in the cut-off wall system. In this study, we present a 2-D analytical model for characterizing the contaminant transport in the cut-off wall and aquitard system. The Green function method and discretization method are used to investigate contaminant migration in a semi-infinite domain. We show that the one-dimensional cut-off wall model may be insufficient to predict contaminant transport when the aquitard hydraulic conductivity (k2) exceeds 5 × 10-10 m/s. The contaminant mass outflow from the aquitard for the case with k2 = 2 × 10-9 m/s may account for 53 % of the total contaminant mass outflow. Increasing k2 from 5 × 10-10 m/s to 5 × 10-9 m/s leads to a 2.9-fold increase of cumulative mass outflow. An approximate fitting formula is proposed to calculate the required minimum thickness of cut-off wall under different type of aquitards. The thicker cut-off wall is required to ensure the cumulative mass outflow will not exceed the allowable value when k2 is relatively high (e.g., 1 × 10-9 m/s).

Research on the Corrosion Inhibition Effect of Xanthium sibiricum on Reinforced Steel and the Prediction of Reinforced Concrete Performance under a Stray Current and Chloride Environment

This study examined a newly developed environmentally friendly plant-based corrosion inhibitor (Xanthium sibiricum). The natural potential method, linear polarization method, steel weight loss method, and corrosion area method were employed to verify the inhibitor’s effectiveness in chloride-containing concrete. The results indicated that Xanthium sibiricum elevated the natural potential of reinforcing steel in concrete, increased its self-corrosion potential, and reduced the self-corrosion current. After three months of curing, the corrosion rate of steel without an inhibitor was approximately 47.5% faster than the experimental group, with the steel loss rate about 40% more severe. The effectiveness of the inhibitor was influenced by increased chloride content in concrete. A two-dimensional multiphase ion transport model of reinforced concrete with realistic aggregate distribution was established using the finite element method (FEM). This model simulated chloride ion transport under typical civil engineering service environments—the coupled effects of a stray current and chloride environment. A predictive formula for the residual compressive strength of reinforced concrete was derived after corrosion under various voltages and chloride ion concentrations for a specific duration. In conjunction with a pump station project operating in a similar environment, the optimal dosage of the Xanthium sibiricum inhibitor for practical engineering was determined to be 2 g/L. At this dosage, the strength of reinforced concrete specimens increased by approximately 31.1%. Finally, a predictive formula for the residual compressive strength of reinforced concrete with an added inhibitor was obtained after corrosion under various voltages and chloride ion concentrations for a specific duration. The conclusions can enhance the durability and safety of concrete structures in similar projects, showing promising application prospects.

Coupling Microkinetics with Continuum Transport Models to Understand CO2R on Planar and Optimized Patterned Electrodes

Electrochemical CO2 reduction has emerged as an attractive technique to sustainably produce valuable fuels and chemicals, while reducing CO2 emissions. In addition to the specific electrode material, the microenvironment local to the electrode interface – e.g., the electrolyte pH, CO2 concentration, buffer concentration – plays a significant role in determining the selectivity and activity of the electrolyzer. In this work, we demonstrate a multi-scale approach, where microkinetic simulations of Au are coupled to a two-dimensional continuum transport model in a flow reactor configuration. Specifically, concentrations and potentials solved for in the continuum model are fed into microkinetic models along the planar electrode surface, which return concentration fluxes that are used to update the continuum model. We find that local quantities at the electrode interface, such as the CO2 concentration and pH, vary not only with the applied potential and flow rate, but also with position on the electrode. We investigate two ways to increase CO2 availability at the interface: (1) increasing the applied flow rate, thus reducing the boundary layer thickness, and (2) introducing an inert defect in the middle of the electrode, allowing CO2 to partially replenish within the depletion boundary layer. Inspired by the effect of a single inert defect, we explore electrodes patterned with inert defects and multiple catalyst materials, i.e. a generalized tandem system. We show how the specific pattern of the electrode can be optimized in order to target a specific desired product. The different parts of the electrode have different goals that work together: converting CO2 to CO, converting CO to C2+ products, and allowing concentrations to replenish through the boundary layer. Overall, our aim is to demonstrate the need for multi-dimensional, multi-scale approaches to simulating CO2 reduction to elucidate the reaction environment near the electrode.This work was performed under the auspices of the U.S. Department of Energy by Lawrence Livermore National Laboratory under Contract DE-AC52-07NA27344. LLNL release number: LLNL-ABS-857537.

Charge interaction behaviors at interfacial domains in DC GIL insulators

Long-term operation of high voltage direct current at elevated temperatures can result in the accumulation of surface charges in DC gas-insulated transmission line (GIL) insulators. Such a phenomenon leads to localized electric field distortion, increasing the risk of surface discharge. The analysis of interaction behaviors between surface charge and space charge at interfacial domains of GIL insulators is a complex task, which requires a comprehensive understanding of physical mechanisms of the gas–solid interface charging. In this work, a two-dimensional bipolar charge transport and interaction (2D BCTI) model is established, with the consideration of both surface and space charge dynamics. Pulsed electroacoustic tests and surface potential measurements are conducted on DC GIL insulator materials under different electrical-thermal coupling conditions. Experimental results exhibit great consistency with the predictions from the 2D BCTI model. The local accumulation of space charge near interfaces has certain effects on surface potential distribution, which in turn influences charge injection behavior from electrodes. In comparison to traditional surface charge simulation models, the consideration of space charge–surface charge interaction behaviors proves to be essential for estimating the polarity and amplitude of surface potential distribution. This model holds promise for assessing charge characteristics in electrical equipment where direct measurement is challenging.

A Two-Dimensional Transport Model of Bipolar Membrane Electrodialysis for the Electro-Regeneration of Carbon Capture Solvents

Bipolar membrane electrodialysis (BPMED) is a promising technology that uses an electric field to split water into H+ and OH- ions. A repeating stack of cation selective, anion selective, and bipolar membranes are placed between two electrodes, creating a repeating pattern of acid, base, and diluting flow chambers. A major application of BPMED is the electrochemical solvent regeneration of capture solutions used in CO2 point source (PS) or direct air capture (DAC). Protons generated using BPMED can be used to shift the acid-base equilibria back to CO2, providing a more sustainable alternative to traditional thermal regeneration methods when integrated with renewable electricity sources. This method combines electrosynthesis and separation, producing a pure stream of CO2 for sequestration and a concentrated base stream.Despite its potential, current limitations to BPMED implementation for this application stem from a large power consumption and high membrane costs. In addition, gas bubbles in solution can cause high electrical resistances, and experimental studies of this phenomenon are challenging due to the small scale on which it occurs. To address these challenges, we present a two-dimensional model of BPMED for CO2 capture solvent recovery. This model depicts a single repeating cell triplet and was constructed using COMSOL Multiphysics version 6.1. The Navier-Stokes equations were implemented to compute a convection field, and the ionic species concentration was found by solving the Nernst-Planck equation. A source term in the concentration partial differential equations was implemented to account for ionic speciation and water splitting reactions, captured through kinetic rate expressions.Experimental validation was conducted on a PC Cell BED 1-4 recirculating batch system with a stack of seven cell pairs of PC acid 60 (anion selective), PC SK (cation selective), and PC Bip (bipolar) membranes. Good agreement between the model and experimental results across a range of conditions and variables demonstrates accuracy and robustness. The calculated pH and concentration profiles revealed that the majority of the speciation reaction takes place inside, or directly adjacent to, the bipolar membrane, including bubble nucleation. Bubbles form inside membrane pores and displace the electrolyte, greatly increasing the electrical resistance of the membrane, particularly when membrane porosity is low. However, the reaction plane location where bubbles are generated can be manipulated by increasing the channel width or reducing the applied voltage.Future research will focus on expanding to multiphase CFD simulations to investigate the effect on flow. Our two-dimensional model provides valuable insights into the complex and multifaceted nature of BPMED, offering a powerful tool for investigating small-scale processes and advancing the use of BPMED for CO2 capture solvent recovery. Figure 1

Two-Dimensional Heavy Metal Migration in Soil with Adsorption and Instantaneous Injection

This paper studies a two-dimensional transport model to investigate the behavior of heavy metal migration in porous media, specifically considering their transport in soil. The model takes into consideration the combination of adsorption term as well as instantaneous injection at the boundary. Also, the model accounts for both longitudinal and transverse movements, providing a comprehensive understanding of heavy metal transport phenomena. In order to obtain the analytical solution Laplace transform has been implemented. It is found that the peak concentration of heavy metals is found to be greatly affected by the changes in the instantaneous injection value. Additionally, a correlation is observed between retardation factors and heavy metal concentrations, with a decrease in retardation factors resulting in an increase in heavy metal concentration.

Solution of Two-Dimensional Solute Transport Model for Heterogeneous Porous Medium Using Fractional Reduced Differential Transform Method

This study contains a two-dimensional mathematical model of solute transport in a river with temporally and spatially dependent flow, explicitly focusing on pulse-type input point sources with a fractional approach. This model is analyzed by assuming an initial concentration function as a declining exponential function in both the longitudinal and transverse directions. The governing equation is a time-fractional two-dimensional advection–dispersion equation with a variable form of dispersion coefficients, velocities, decay constant of the first order, production rate coefficient for the solute at the zero-order level, and retardation factor. The solution of the present problem is obtained by the fractional reduced differential transform method (FRDTM). The analysis of the initial retardation factor has been carried out via plots. Also, the influence of initial longitudinal and transverse dispersion coefficients and velocities has been examined by graphical analysis. The impact of fractional parameters on pollution levels is also analyzed numerically and graphically. The study of convergence for the FRDTM technique has been conducted to assess its efficacy and accuracy.

Development of Simple Formula for Transverse Dispersion Coefficient in Meandering Rivers

This study aims to develop a straightforward and practical formula to estimate transverse dispersion coefficients in meandering natural rivers, a critical factor for predicting solute transport. We present a novel expression for the transverse dispersion coefficient based on dispersion and hydraulic data sets obtained from tracer experiments conducted in natural rivers. A distinctive feature of the formula is its reliance on one dimensionless hydraulic parameter, u¯u*hRc. To assess the effectiveness and accuracy of our proposed formula, we compare it with previously established equations commonly employed in the field. Furthermore, we apply the formula to natural river bends situated in the Nakdong River of Korea. This equation serves to estimate the initial value of the dispersion coefficient in two-dimensional solute transport modeling. As a result, the calibrated value of the dimensionless transverse dispersion coefficient is 0.97, which is only a 16% difference between the initial value of 0.81 as obtained from the formula. The formula presented in this study, simplifying and utilizing the dimensionless hydraulic parameter, offers a promising approach to estimating transverse dispersion in natural meandering rivers in cases where tracer and secondary flow data are unavailable. Additionally, the formula can be refined with more recent dispersion data, leading to a clearer, more straightforward, and validated formulation that captures the intricate interplay between topography and transverse dispersion.

Two-dimensional modeling of CO2 mineral trapping through the oxalate‑carbonate pathway: Influence of the root system model

Two-dimensional reactive transport models, one with a simplified root system and the other accounting for dynamically evolving root architecture, were constructed to examine the influence of model complexity on capturing the effect of soil-root dynamics relating to the Oxalate Carbonate Pathway (OCP) of the Iroko tree over 170 years. Oxidation of oxalate from fallen tree tissue by soil bacteria enables local soil pH increase, leading to the sequestration of atmospheric carbon in carbonate minerals (calcite) in the shallow soil surrounding the tree. Simulations of both root models corroborate previous one-dimensional models of the OCP focused on Ca and C mass balance, where high weathering rates of Ca-containing silicate minerals in bedrock, along with contributions from groundwater, provided sufficient Ca for precipitation of observed quantities of calcite. Both simulations demonstrate the development of a distinct high pH zone where oxalate is oxidized, Ca accumulates, and calcite precipitates (OCP zone); and a low pH zone where roots collect Ca, later returned to the top soil as calcium oxalate (Total Root Extent/TRE zone) via litterfall. While the extent of OCP zone development near the ground surface was very similar between simulations, differences in localized root water uptake between the two approaches resulted in variation in water and solute transport and influenced the geometry of the OCP zone at depth, with implications for calcite precipitation in the soil. Trends in CO2 and O2 partial pressures in the OCP zone were mirrored in the TRE zone, suggesting linkage between the two zones with regard to gas transport. Near the end of the tree's lifespan, results indicate that soil permeability decreases due to calcite precipitation may limit O2 ingress and availability in the shallow soil, while trapping CO2 released from the oxidation of organics in the shallow soil, with implications for the long-term sustainability of the OCP itself.

Finite-difference model for two-dimensional transport of degradable organic contaminants in a buffer layer-cutoff wall-aquifer system accounting for the consolidation behavior

Consolidation changes the geotechnical properties of a cutoff wall with depth. However, current research in transport theory fails to account for this factor. In this paper, a model for two-dimensional transport of degradable organic contaminants in a buffer layer-cutoff wall-aquifer system is developed, which takes the impacts of consolidation behavior into account for the first time. The finite-difference approach is employed to address this model, and its reasonableness is exhibited by the comparative analyses of three-different scenarios. Following this, a case-study on the barrier performance shows that accounting for the cutoff wall’s consolidation behavior prolongs the defined breakthrough time tb compared to when it is neglected. When a change in the cutoff wall consolidation parameter increases the horizontal and vertical effective stresses, its anti-fouling effectiveness is enhanced. Moreover, reducing the depth of the maximum inflow boundary concentration and increasing the concentration distribution range result in a decrease in tb. When the degradation half-life ≤ 400 years, the degradation effect should be considered to provide an accurate assessment of the cutoff wall’s longevity. For the “enclosure” cutoff walls, the appropriate thicknesses of the buffer layer and the cutoff wall can be determined based on the site environment and the construction cost.

Numerical Simulation of Sediment Transport in Unsteady Open Channel Flow

This paper presented a two-dimensional, well-balanced hydrodynamic and sediment transport model based on the solutions of variable-density shallow-water equations (VDSWEs) and the Exner equation for bed change for simulating sediment transport in unsteady flows. Those equations are solved in a coupled way by the first-order Godunov-type finite volume method. The Harten–Lax–van Leer–Contact (HLLC) Riemann solver is extended to find the local Riemann fluxes to maintain the exact balance between the momentum term and the bed slope term. A well-balanced scheme is superior to an unbalanced scheme to minimize numerical dispersion as demonstrated by the synthetic standing contact-discontinuity test case. Following this, the model is employed to simulate two laboratory experiments and a field case, the 1996 Lake Ha! Ha! flood event in Canada. The results of water surface elevations and bed surface profiles agree well with the measurements. The accuracy and robustness of the numerical schemes make the model a good candidate for practical engineering applications.

Two-dimensional reactive transport model as a new approach for identifying the origins and contribution of arsenic in a soil and water system

Distribution of arsenic (As) in the environment can occur through various paths, depending on the interrelationship between soil and surface water, and groundwater. This study aimed to identify the origins and pathways of As among various potential origins in soil and water systems using two-dimensional reactive transport models (2D RTMs) and to evaluate the contribution of each origin to the geochemical distribution of As in the study area, which included a reservoir, two streams, and two abandoned mines. The 2D RTMs were prepared by combining hydraulic and geochemical reactive models. The MODFLOW-2005 package was used to simulate the water flow in the study area. The geochemical reaction and distribution of As were simulated using Geochemist's Workbench. The concentration of As in the reservoir was ∼50 μg/L, while that in the soils collected around the mine X and Y was 15.0–1853.2 mg/kg (median: 126.1 mg/kg) and 3.6–1629.2 mg/kg (median: 60.9 mg/kg), respectively. Hydraulic and geochemical input data were derived from field surveys and laboratory experiments. Based on the geochemical and hydrogeological input data, 2D RTMs determined that 32.0 kg of As accumulated into the reservoir over 18 years through surface water, which is from mines X (2.7 kg) and Y (29.3 kg). Hence, mine Y would be the primary origin of As to the reservoir. Based on 2D RTM, 2.4 kg of As entered the reservoir via groundwater (DO <0.1 mg/L). Also, precipitation of sulfide minerals (e.g., orpiment; As2S3) influences As dispersion and contribution under subsurface conditions. These modeling results suggest that 2D RTMs could be a new approach for identifying and assessing the contribution of potential As origins in soil and water system.

Effects of macro-scale heterogeneity on the kinetic interface-sensitive tracer test for measuring the fluid–fluid interfacial area in dynamic two-phase flow in porous media

A novel reactive smart tracer method, termed the kinetic interface-sensitive (KIS) tracer test, has been demonstrated in laboratory column experiments to enable measurement of the specific capillary-associated fluid–fluid interfacial area in dynamic two-phase flow displacement processes in porous media. Development of the tracer method towards effective application in real field conditions requires investigation of the influence of the porous media heterogeneity on the front size and the specific interfacial area, and, consequently, in how far a kinetic interface-sensitive tracer experiment, and the corresponding breakthrough curves, are affected. This study employs a two-dimensional Darcy-scale two-phase flow reactive transport model to investigate numerically the KIS tracer transport in heterogeneous porous media. Simulations were carried out for the primary drainage process in a domain formed of fine and coarse porous media. Various heterogeneity patterns, having different numbers of inclusions and different geometrical distributions, were studied. It is shown that the shape of the breakthrough curves can be used as an indicator for quantifying the displacement front roughness, the specific interfacial area in the domain, and the domain heterogeneity, e.g., the existence of preferential flow pathways inside the porous media. The results indicate that when the displacement front roughness is small, the concentration breakthrough curves exhibit a linear increase. The slopes of the breakthrough curves linearly depend on the fraction of the bulk volume occupied by the low-permeability sand inclusions. The volume-averaged specific interfacial area and the size of the transition zone can be determined from the slopes of the breakthrough curves.

Development of a Two-Dimensional Hybrid Sediment-Transport Model

This paper presents the development of a two-dimensional hydrodynamic sediment transport model using the finite volume method based on a collocated unstructured hybrid-mesh system consisting of triangular and quadrilateral cells. The model is a single-phase nonequilibrium sediment-transport model for nonuniform and noncohesive sediments in unsteady turbulent flows that considers multiple sediment-transport processes such as deposition, erosion, transport, and bed sorting. This model features a hybrid unstructured mesh system for easy mesh generation in complex domains. To avoid interpolation from vertices in conventional unstructured models, this model adopted a second-order accurate edge-gradient evaluation method to consider the mesh irregularities based on Taylor’s series expansion. In addition, the multipoint momentum interpolation corrections were integrated to avoid possible nonphysical oscillations during the wetting-and-drying process, common in unsteady sediment transport problems, to ensure both numerical stability and numerical accuracy. The developed sediment transport model was validated by a benchmark degradation case for the erosion process with armoring effects, a benchmark aggradation case for the deposition process, and a naturally meandering river for long-term unsteady sediment-transport processes. Finally, the model was successfully applied to simulate sediment transport in a reservoir that was significantly affected by typhoon events.

Development of a multiphase perfusion model for biomimetic reduced-order dense tumors

Dense fibrous extracellular constitution of solid tumors exerts high resistance to diffusive transport into it; additionally, the scarcity of blood and lymphatic flows hinders convection. The complexity of fluidic transport mechanisms in such tumor environments still presents open questions with translational end goals. For example, clinical diagnosis and targeted drug delivery platforms for such dense tumors can ideally benefit from a quantitative framework on plasma uptake into the tumor. In this study, we present a computational model for physical parameters that may influence blood percolation and penetration into simple biomimetic solid tumor geometry. The model implements three-phase viscous-laminar transient simulation to mimic the transport physics inside a tumor-adhering blood vessel and measures the constituent volume fractions of the three considered phases, viz. plasma, RBCs (red blood cells, also known as “erythrocytes”), and WBCs (white blood cells, also known as “leukocytes”) at three different flow times, while simultaneously recording the plasma pressure and velocity at the entry point to the tumor’s extracellular space. Subsequently, to quantify plasma perfusion within the tumor zone, we proposed a reduced-order two-dimensional transport model for the tumor entry zone and its extracellular space for three different fenestra diameters: 0.1, 0.3, and 0.5 µm; the simulations were two-phase viscous-laminar transient. The findings support the hypothesis that plasma percolation into the tumor is proportional to the leakiness modulated by the size of fenestra openings, and the rate of percolation decays with the diffusion distance.

Theoretical Analysis of a Nonequilibrium Transport Model of Two-Dimensional Nonisothermal Reactive Chromatography Accounting for Bi-Langmuir Isotherm.

The current study investigates a nonequilibrium and nonlinear two-dimensional lumped kinetic transport model of nonisothermal reactive liquid chromatography, considering the Bi-Langmuir adsorption isotherm, heterogeneous reaction rates, radial and axial concentration variations, and the adsorption and reaction enthalpies. The mathematical models of packed bed chromatographic processes are expressed by a highly nonlinear system of coupled partial differential algebraic equations connecting the phenomena of convection, diffusion, and reaction, for mass and energy balance, the differential algebraic equations for mass balance in the solid phase, and the algebraical expressions for the adsorption isotherms and for the reaction rates. The nonlinearity of the reaction term and the adsorption isotherm preclude the derivation of an analytical solution for the model equations. For this reason, a semidiscrete, high-resolution, finite-volume technique is extended and employed in this study to obtain the numerical solution. Several consistency checks are performed to evaluate the model predictions and analyze the precision of the proposed numerical scheme. A number of heterogeneously catalyzed stoichiometric reactions are numerically simulated to examine reactor performance under the influence of temperature and Bi-Langmuir adsorption dynamics, the level of coupling between mass and energy fronts, and to study the effects of various critical parameters. The numerical results obtained are beneficial for optimal predictive control and process optimization during production and the development of methods for systematic design and fault detection of nonisothermal liquid chromatographic reactors, and hence constitute the first step to provide deeper insight into the overall evaluation of integrated reaction and separation processes.

Numerical Study of Aluminum Pitting Corrosion Using a Two-Dimensional Multi-Species Transport Model

Numerical Study of Aluminum Pitting Corrosion Using a Two-Dimensional Multi-Species Transport Model

Space charge characteristics in epoxy/nano-MgO composites: Experiment and two-dimensional model simulation

Space charge accumulation in polymer dielectrics may lead to serious electric field distortion and even insulation failure during long-term operations of power equipment and electronic devices, especially under conditions of high temperature and direct current electric stress. The addition of nanoparticles into polymer matrices has been found effective in suppressing space charge accumulation and alleviating electric field distortion issues. Yet, the underlying mechanisms of nanoparticle doping remain a challenge to explore, especially from multi-dimensional composite insights. Here, a two-dimensional bipolar charge transport model with consideration of interface zones between organic/inorganic phases is proposed for the investigation into space charge behaviors of nanodielectrics. To validate the effectiveness and feasibility of the model, pulsed electroacoustic experiments are performed on epoxy/nano-MgO composites with different doping ratios of nanoparticles. Experimental observations match well with simulation anticipations, i.e., higher doping ratios of nanoparticles below the percolation threshold exhibit better capabilities to inhibit space charge accumulation. The deep traps (∼1.50 eV) generated in the interface zones are demonstrated to capture free charges, forming a reverse electric field in the region adjacent to electrodes and impeding the space charge migration toward the interior of the composite. This model is anticipated to provide theoretical insight for understanding space charge characteristics in polymer nanodielectrics and computing charge dynamics in extreme conditions where experiments are challenging to perform.