Pith. sign in

REVIEW 4 major objections 3 minor 1 cited by

An adaptive two-grid preconditioner and linearly implicit scheme for shale gas transport in fractured porous media

T0 review · 4 major / 3 minor · reviewed 2026-08-12 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Shale gas transport in fractured media can be advanced in time with a fixed linear operator and solved with an adaptive spectral two-grid preconditioner whose iteration count is independent of fracture-matrix permeability contrast.

desk verdict Useful integration of fixed-operator ImEx with adaptive spectral two-grid; the stability proof has a fixable hypothesis mismatch and an unverified concentration-range assumption. read the letter →

arxiv 2411.17903 v2 pith:C64XQUB7 submitted 2024-11-26 math.NA cs.NA

Maria Vasilyeva , Ben S. Southworth , Shubin Fu This is my paper · ORCID
classification math.NAcs.NA MSC 65M6065F0865N5576S0565M12
keywords shalegastransportfracturedporousmediamixed-dimensionalmodellinearlyimplicitschemeexplicit-implicit-nullmethodadaptivespectralcoarsespacetwo-gridpreconditionerhigh-contrastpermeability
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

Shale gas moves through a porous matrix and a network of thin fractures with very different permeabilities, which makes simulation stiff and expensive. The paper tries to establish two claims: that the nonlinear transport equations can be advanced in time with a fixed linear operator chosen once, giving an unconditionally stable linearly implicit scheme under a dominance condition; and that the resulting linear system can be solved by an adaptive spectral two-grid preconditioner whose iteration count does not grow with fracture-matrix permeability contrast. If both claims hold, each time step reduces to one preconditioned conjugate gradient solve with a preconditioner built once offline, instead of rebuilding a multiscale solver whenever the concentration field changes. The paper supports the claims with tests on 30 and 160 fracture networks, contrasts from $10^{3}$ to $10^{9}$, and both homogeneous and heterogeneous matrix properties.

What carries the argument

The load-bearing mechanism is the combination of two ingredients. First, the additive split $S(c)=S^{(\mathrm{lin})}+S^{(\mathrm{nl})}(c)$, $D(c)=D^{(\mathrm{lin})}+D^{(\mathrm{nl})}(c)$, with $S^{(\mathrm{lin})}, D^{(\mathrm{lin})}$ chosen as coefficient-wise upper bounds over $[c_{\min},c_{\max}]$, converts the nonlinear evolution into one linear system $A=S^{(\mathrm{lin})}+\tau D^{(\mathrm{lin})}$ per time step; Theorem 1 and the appendix's canonical three-step representation give the unconditional-stability criterion $S^{(\mathrm{lin})}>S^{(\mathrm{nl})}$, $D^{(\mathrm{lin})}>D^{(\mathrm{nl})}$. Second, the two-grid preconditioner uses a symmetric Gauss-Seidel smoother, spectrally equivalent to the diagonal of $A$, and a spectral coarse space generated by the local generalized eigenvalue problem (14), with adaptive thresholding $\lambda_{m_i}<\delta_\lambda$; partition-of-unity multiplication gives continuous basis functions, and the local approximation bound (16) feeds the global estimate $\|v-Pv\|_D^2 \le (H^2/\Lambda^*)\|v\|_A^2$ that yields $K_{TG}=c_m H^2/\Lambda^*$.

What would settle it

At each time step of the reported 30- and 160-fracture runs, compute the smallest eigenvalue of $S^{(\mathrm{lin})}-S^{(\mathrm{nl})}(c^n)$ and of $D^{(\mathrm{lin})}-D^{(\mathrm{nl})}(c^n)$; if any eigenvalue is negative while $c^n$ lies inside $[c_{\min}, c_{\max}]$, the unconditional-stability claim fails. Separately, fix the threshold $\delta_\lambda=10^{-3}$ and increase the fracture-matrix permeability contrast beyond $10^{9}$; if the average number of PCG iterations rises systematically with contrast, the contrast-independence claim fails.

Watch

Extended reading notes

Core claim

On its own terms, the central discovery is that an additive split of the operators, $S(c)=S^{(\mathrm{lin})}+S^{(\mathrm{nl})}(c)$ and $D(c)=D^{(\mathrm{lin})}+D^{(\mathrm{nl})}(c)$, can be tuned so the linear part is fixed for all time steps and yet dominates the nonlinear part. The authors choose $a^*_m(x)=a_m(x,c_{\min})$, $b^*_m(x)=b_{m1}(x,c_{\min})+b_{m2}(x,c_{\max})$, and $b^*_f(x)=b_f(x,c_{\max})$, then prove in Theorem 1 that the linearly implicit scheme (10) is unconditionally stable whenever $S^{(\mathrm{lin})}>S^{(\mathrm{nl})}(c^n)$ and $D^{(\mathrm{lin})}>D^{(\mathrm{nl})}(c^n)$ at every step; the appendix verifies positivity of $R_n - \tfrac{1}{4}D_n$ through a canonical three-step energy estimate. For the linear solve, the coarse space is built from the local generalized eigenvalue problem $A^{\omega_i}\psi^{\omega_i}_k = \lambda^{\omega_i}_k D^{\omega_i}\psi^{\omega_i}_k$, keeping eigenvectors below a threshold $\delta_\lambda$ and multiplying them by partition-of-unity functions to form the interpolation operator $P$; the resulting two-grid method has convergence factor controlled by $K_{TG}=c_m H^2/\Lambda^*$ from estimate (19). The paper reports that with this adaptive coarse space, PCG converges to a $10^{-9}$ residual in O(1) iterations across all tested contrasts, fracture counts, and matrix heterogeneities, while the implicit scheme remains first-order accurate in time.

Load-bearing premise

The argument assumes that the fixed linear operator chosen as an upper bound over the concentration interval $[c_{\min}, c_{\max}]$ indeed dominates the true nonlinear operator at every time step; the paper sets the coefficients this way but never checks the resulting matrix inequalities during the simulations, and the upper-bound construction does not exactly match the positive-definiteness hypotheses stated in the stability theorem.

Editorial extensions

If this is right

  • Each nonlinear time step becomes one linear solve with the same matrix $A=S^{(\mathrm{lin})}+\tau D^{(\mathrm{lin})}$, so spectral multiscale basis functions are computed once offline and reused for the entire simulation.
  • With the adaptive coarse space, PCG reaches a residual tolerance of $10^{-9}$ in O(1) iterations across contrasts $10^{3}$ to $10^{9}$, for 30 and 160 fractures, on both $10\times10$ and $20\times20$ coarse grids.
  • A single threshold $\delta_\lambda=10^{-3}$ automatically places more basis functions where fractures create small local eigenvalues, matching the performance of 4 to 8 fixed basis functions per node while keeping the coarse system smaller.
  • The linearly implicit scheme shows first-order time accuracy: relative $L^2$ errors shrink as the number of time steps increases in the reported tests.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Beyond the paper: the dominance condition could be monitored cheaply at runtime, and if the concentration approaches the assumed bounds, the upper-bound operator and local basis functions could be refreshed only in the affected subdomains.
  • Beyond the paper: the same additive-split-plus-spectral-coarse-space construction should transfer to other degenerate nonlinear parabolic problems with monotone coefficient dependence, such as unsaturated flow or gas-condensate transport, provided the local eigenvalue decay still separates the stiff modes.
  • Beyond the paper: the reported tables suggest that what drives the required number of basis functions is the number of small-eigenvalue modes created by fracture geometry rather than the permeability contrast itself; a quantitative formula linking coarse-space dimension to fracture density would be a natural next step.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 3 minor

Summary. The paper proposes a linearly implicit (ImEx) time integration scheme for a mixed-dimensional shale gas transport model with explicit fracture networks, combined with an adaptive spectral two-grid preconditioner for the resulting fixed-matrix linear systems. The linear operator is built once from pointwise coefficient upper bounds, so multiscale basis functions are computed offline only once, and each time step reduces to solving a linear system with a previously constructed preconditioner. The paper claims unconditional stability of the time scheme under conditions on the operator splitting, and claims iterative convergence independent of the fracture/matrix permeability contrast. Numerical experiments cover 30 and 160 fracture networks with homogeneous and heterogeneous matrices over three permeability contrasts.

Significance. If the stability and convergence claims hold, the fixed-operator, one-time-preconditioner approach would be a practically valuable advance: it avoids rebuilding multiscale basis functions every nonlinear iteration or time step, which is the dominant cost in GMsFEM-type solvers for time-dependent problems. The numerical study is extensive, with two fracture geometries, homogeneous and heterogeneous backgrounds, three permeability contrasts, two coarse grids, and both fixed and adaptive numbers of spectral basis functions; the observed iteration counts are low and largely robust across these variations. The adaptive eigenvalue threshold gives a single-parameter mechanism for choosing local coarse-space dimension, which is an attractive practical feature. However, as detailed below, the stability theorem as stated does not apply to the constructed splitting, and the concentration-range condition on which the argument relies is never verified; the theoretical two-grid bound also does not itself establish contrast independence. These gaps affect the two central claims of the paper, so substantial revision is needed.

major comments (4)
  1. [Section 3.2, Theorem 1, and Appendix] Theorem 1 states hypotheses for symmetric positive definite S(nl)(c^n) and D(nl)(c^n) and gives condition (11) as sufficient for unconditional stability. The constructed splitting, however, uses coefficient-wise upper bounds a*_m = a_m(x,cmin), b*_m = b_m1(x,cmin)+b_m2(x,cmax), and b*_f = b_f(x,cmax), which implies S(c^n) <= S(lin) and D(c^n) <= D(lin) in the Loewner order. Consequently S(nl) = S(c^n) - S(lin) and D(nl) = D(c^n) - D(lin) are negative semidefinite, not positive definite. The theorem as stated is therefore inapplicable to the scheme actually analyzed and tested. The appendix's proof derives the correct sufficient condition, S(lin)-S(nl)(c^n) > 0 and D(lin)-D(nl)(c^n) > 0, which is weaker and is what should appear in the theorem statement. This is a load-bearing mismatch because the abstract and conclusion assert unconditional stability under the stated splitting.
  2. [Section 3.2 and Section 5] The stability condition ultimately rests on the assumption that all nodal concentrations satisfy cmin <= c_m^n, c_f^n <= cmax for every time step n. The paper states this assumption but never proves a discrete maximum principle, nor does it verify the resulting discrete matrix inequalities (11) for the computed solutions in the numerical tests. The fracture source/sink term f_f(c_f) = c_w ZRT κ_w (c_w - c_f)/mu is nonlinear and can act as a source when c_f < c_w, and the explicit treatment of the nonlinear residual provides no mechanism to enforce the bound. Since the invariance of the concentration range underlies both the unconditional stability claim and the offline construction of the fixed preconditioner A = S(lin) + tau D(lin), the authors should either prove that the scheme preserves the interval [cmin, cmax] or add a direct numerical verification for all reported test cases, and if necessary introduce a projection or clipping step that preserves the stability argument.
  3. [Section 4.3, Eq. (19)] The two-grid convergence bound K_TG = c_m H^2 / Lambda* is asserted for the matrix A in (13), which is the fixed time-stepping matrix S(lin)+tau D(lin) of the mixed-dimensional system. The derivation combines the local projection estimate (16) with a partition-of-unity argument, but it does not establish that the local eigenproblems (14) with a 'modified diagonal part' are consistent with the global matrix A, nor does it quantify how Lambda* and the constant c_m depend on the fracture permeability contrast and on the time step tau. As written, inequality (19) does not prove contrast-independent convergence; such independence is only an empirical observation in Tables 3 and 4. The authors should either make the dependence on contrast explicit in the bound or present the theoretical result as conditional on the spectral coarse space capturing the relevant fracture-contrast modes.
  4. [Abstract and Conclusion, Tables 3-4] The claim of 'iterative convergence independent of the contrast of fracture and porous matrix permeability' is stronger than the data and theory support. For instance, Table 4, Test 2a on the 20x20 coarse grid with m=4 basis functions shows average iterations of 12.5, 28.0, and 72.0 for contrasts 10^3, 10^6, and 10^9, respectively, and the adaptive rows also show a noticeable (albeit milder) increase with contrast (e.g., 14.5, 16.7, 21.0 for delta_lambda=10^-4). The authors should qualify the claim, for example as 'robust' or 'weakly dependent on contrast' in the adaptive regime, or state precisely which configurations demonstrate independence.
minor comments (3)
  1. [Section 5] In the parameters paragraph, 'kf = 10^3, 10^6 and 10^6' appears to be a typo; the subsequent text and tables indicate the intended contrasts are 10^3, 10^6, and 10^9.
  2. [Appendix] The symbol S_n is used in the sentence 'for X n = Bn and X n = Sn' but is never defined; it should presumably be S(lin)+S(nl)(c^n) or the corresponding mass-type matrix from (10).
  3. [Throughout] There are several typographical errors, for example 'particulaly' (Section 3.2), 'porious media' (Section 4), and 'we proof Theorem 1' (Appendix); these should be corrected in revision.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity found: the stability condition is an explicit design constraint, the spectral coarse space is a standard solver component, and the numerical claims are validated independently of any fitted prediction.

full rationale

The paper's derivation chain does not reduce to its inputs. The linearly implicit scheme (10) is analyzed in the Appendix via a canonical three-step form; stability is derived from the sufficient conditions S(lin)-S(nl)>0 and D(lin)-D(nl)>0. The linear operators are then constructed from pointwise coefficient upper bounds (a*_m=a_m(x,cmin), b*_m=b_m1(x,cmin)+b_m2(x,cmax), b*_f=b_f(x,cmax)) precisely so that these dominance conditions hold when concentrations remain in [cmin,cmax]. Constructing operators to satisfy a sufficient condition is not circular: the dominance condition is an input, and unconditional stability is a derived consequence. The adaptive coarse space is built from eigenvectors of the same matrix A that the two-grid preconditioner inverts; this is the standard spectral-AMG/GMsFEM design principle, not a renamed prediction, and the convergence bound (19) depends on the chosen spectral gaps in an explicit, standard way. The threshold delta_lambda is acknowledged as a tunable parameter, and robustness is demonstrated across three thresholds and fixed numbers of basis functions. The most serious concern in the paper is a correctness gap, not circularity: under the upper-bound construction S(nl) and D(nl) are nonpositive, so the stated hypotheses of Theorem 1 (symmetric positive definite S(nl),D(nl)) do not match the constructed splitting, and the required discrete matrix inequalities are never checked on the mixed-dimensional system. That affects the validity of the unconditional-stability claim, but it is not an instance of a prediction being equivalent to an input by construction. The numerical experiments provide independent evidence for the solver claims across 30/160 fractures, homogeneous/heterogeneous matrices, and contrasts 1e3-1e9.

Assumptions & free parameters 2 free parameters · 4 assumptions · 0 invented entities

The method rests on standard numerical linear algebra tools (FEM, GMsFEM spectral coarse spaces, two-grid convergence theory) and on the physical shale gas model from prior literature. No new physical entities are introduced. The only tunable method parameters are the eigenvalue threshold delta_lambda and the smoothing count nu; the paper's qualitative conclusions are robust to delta_lambda but not completely parameter-free.

free parameters (2)
  • delta_lambda (eigenvalue threshold) = 10^-3 (primary), 10^-2 and 10^-4 explored
    Adaptive basis selection uses m_i = max{i: lambda_i < delta_lambda}; this threshold controls coarse-space size and iteration counts and is chosen by the user. The paper demonstrates robustness across a range but gives no automatic selection rule.
  • nu (number of Gauss-Seidel smoothing iterations) = 5
    Smoothing count in the two-grid preconditioner is fixed at 5 for all experiments; convergence results depend on this choice.
assumptions (4)
  • standard math Finite element discretization on a fracture-conforming unstructured mesh yields symmetric positive semidefinite stiffness and mass matrices.
    Used throughout Section 3.1 to define the discrete system (9) and the later two-grid analysis.
  • standard math The two-grid convergence framework of Falgout, Vassilevski, and Zikatanov, and the spectral equivalence of symmetric Gauss-Seidel to the matrix diagonal, are valid for this mixed-dimensional system.
    Basis for Section 4.3 convergence estimate K_TG = c_m H^2/Lambda*.
  • domain assumption The mixed-dimensional shale gas model with Langmuir adsorption and Darcy/diffusive transport is the correct physical description.
    Equations (1)-(7) are taken from prior shale gas modeling literature; the practical value of the method depends on this model.
  • domain assumption Concentrations stay in [cmin, cmax], with cmin from boundary or production pressure and cmax from initial concentration.
    Used to fix the upper-bound coefficients a*_m, b*_m, and b*_f that define the stable linear operator in Section 3.2. Violation of this range would break the stated stability condition.

how reviews work

0 comments
Cite this review

Pith. "Pith review of An adaptive two-grid preconditioner and linearly implicit scheme for shale gas transport in fractured porous media." pith.science (2026). https://pith.science/paper/C64XQUB7

@misc{pith2026241117903,
  author       = {Pith},
  title        = {Pith review of: An adaptive two-grid preconditioner and linearly implicit scheme for shale gas transport in fractured porous media},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/C64XQUB7}},
  note         = {Machine review of arXiv:2411.17903}
}
read the original abstract

We consider a nonlinear mixed-dimensional model for simulating gas transport in shale formation. The mathematical model consists of a coupled system of nonlinear equations, where flow within fractures is represented using a lower-dimensional representation. For the numerical solution of the coupled transport problem, we construct an unstructured mesh that resolves lower dimensional fractures on the grid level and use the finite element approximation to build a discrete system. To construct an efficient scheme for the resulting nonlinear problem, we use an explicit-implicit method for time integration, where we carefully choose an additive partition of the nonlinear operators to separate the stiff linear component and integrate it implicitly to ensure the stability of the time integration. Next, we invert the linear partition of the operator by constructing an efficient two-grid preconditioner for shale gas transport in fractured porous media. We use a local pointwise smoother on the fine grid and carefully design an adaptive multiscale space for coarse grid approximation based on local generalized eigenvalue problems. We utilize an adaptive thresholding to automatically identify local dominant modes which correspond to the very small eigenvalues in local domains. We remark that such spatial features are automatically captured through our local spectral problems, and connect these to fracture information in the global formulation of the problem. Approximation properties of the local spectral space with convergence of the proposed two-grid algorithm are given. Numerical results are presented for two fracture distributions with 30 and 160 fractures, demonstrating iterative convergence independent of the contrast of fracture and porous matrix permeability.

Figures

Figures reproduced from arXiv: 2411.17903 by the authors.

Figure 1
Figure 1. Illustration of a fine-scale triangulation conformed with fracture grid. [PITH_FULL_IMAGE:figures/full_fig_p006_1.png] view at source ↗
Figure 2
Figure 2. Illustration of nonlinear coefficients am(c) and bm(c) for shale gas transport. To satisfy the stability estimate, we should use an upper bound for operators S n and Dn in an additive representation. For S n, we have a block diagonal matrix with linear block Sf and linearized block S n m that depends on coefficients af and am, respectively. Similarly, for operator D, we have a block structure with blocks that rely o… view at source ↗
Figure 3
Figure 3. Illustration of a coarse grid T H with local domains (ω35 and ω39) and coarse elements (K1 and K87). where ϕ ωi k = χiψ ωi k . In order to construct the multiscale space V ω off, we solve local spectral problems in each ω. The analysis in [30, 76, 31] motivates the following generalized eigenvalue problem: A ωiψ ωi k = λ ωi k Dωiψ ωi k , (14) where A ωi =  Aωi mm A ωi mf A ωi fm A ωi f f  , denotes a principle sub… view at source ↗
Figures from the paper (7 more)
Figure 4
Figure 4. Figure 4: Eigenvalues and eigenvectors in local domains [PITH_FULL_IMAGE:figures/full_fig_p012_4.png]
Figure 5
Figure 5. Figure 5: Computational domains and heterogeneous coefficients. Fracture lines (white color) and fine-scale [PITH_FULL_IMAGE:figures/full_fig_p015_5.png]
Figure 6
Figure 6. Figure 6: Solution at final time for homogeneous (first and third plots) and heterogeneous porous matrix [PITH_FULL_IMAGE:figures/full_fig_p015_6.png]
Figure 7
Figure 7. Figure 7: Number of nonlinear iterations for implicit scheme with Picard iterations. Contrast 10 [PITH_FULL_IMAGE:figures/full_fig_p016_7.png]
Figure 8
Figure 8. Figure 8: Relative error dynamics for linearly implicit scheme. Contrast 10 [PITH_FULL_IMAGE:figures/full_fig_p016_8.png]
Figure 9
Figure 9. Figure 9: Number of time steps (Nt) vs relative error at final time. Contrast 109 10 × 10 and 20 × 20 coarse grids. The influence of permeability contrast is investigated for 103 , 106 , and 109 contrasts. The study uses 100 in an iterative solver as the maximum number of iterat…
Figure 10
Figure 10. Figure 10: Illustration of adaptivity in spectral space construction for Test 1a and 2a. First column: first [PITH_FULL_IMAGE:figures/full_fig_p020_10.png]

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

  1. Multiscale approximation and two-grid preconditioner for extremely anisotropic heat flow

    math.NA 2024-12 conditional novelty 6.0 of 10

    Spectral multiscale basis functions built from local anisotropic Laplacians give accurate coarse models and O(1) two-grid convergence for heat flow with anisotropy up to 10^12.

Reference graph

Works this paper leans on

89 extracted references · 68 canonical work pages · cited by 1 Pith paper

  1. [1]

    Multiscale domain decomposition methods for elliptic problems with high aspect ratios

    Jørg Aarnes and Thomas Y Hou. Multiscale domain decomposition methods for elliptic problems with high aspect ratios. Acta Mathematicae Applicatae Sinica , 18(1):63–76, 2002

    2002

  2. [2]

    A convergence analysis of generalized multiscale finite element methods

    Eduardo Abreu, Ciro D ´ ıaz, and Juan Galvis. A convergence analysis of generalized multiscale finite element methods. Journal of Computational Physics , 396:303–324, 2019

    2019

  3. [3]

    Unconditionally stable schemes for convection-diffusion problems

    Nadezhda Afanas’eva, Petr Vabishchevich, and Maria Vasilyeva. Unconditionally stable schemes for convection-diffusion problems. Russian Mathematics, 57:1–11, 2013

    2013

  4. [4]

    Multiscale model reduction for shale gas transport in a coupled discrete fracture and dual-continuum porous media

    I Yucel Akkutlu, Yalchin Efendiev, Maria Vasilyeva, and Yuhe Wang. Multiscale model reduction for shale gas transport in a coupled discrete fracture and dual-continuum porous media. Journal of Natural Gas Science and Engineering , 48:65–76, 2017

    2017

  5. [5]

    Multiscale model reduction for shale gas transport in poroelastic fractured media

    I Yucel Akkutlu, Yalchin Efendiev, Maria Vasilyeva, and Yuhe Wang. Multiscale model reduction for shale gas transport in poroelastic fractured media. Journal of Computational Physics , 353:356–376, 2018

    2018

  6. [6]

    Multiscale gas transport in shales with local kerogen heterogeneities

    I Yucel Akkutlu and Ebrahim Fathi. Multiscale gas transport in shales with local kerogen heterogeneities. SPE Journal , 17(04):1–002, 2012

    2012

  7. [7]

    Multiscale model reduction for shale gas transport in fractured media

    Ibrahim Y Akkutlu, Yalchin Efendiev, and Maria Vasilyeva. Multiscale model reduction for shale gas transport in fractured media. Computational Geosciences, 20:953–973, 2016

    2016

  8. [8]

    Multiscale model reduction for shale gas transport in fractured media

    IY Akkutlu, Yalchin Efendiev, and Maria Vasilyeva. Multiscale model reduction for shale gas transport in fractured media. Computational Geosciences, pages 1–21, 2015. 23

    2015

Show all 89 references
  1. [9]

    Efficient algebraic two-level schwarz preconditioner for sparse matrices

    Hussam Al Daas, Pierre Jolivet, and Tyrone Rees. Efficient algebraic two-level schwarz preconditioner for sparse matrices. SIAM Journal on Scientific Computing , 45(3):A1199–A1213, 2023

    2023

  2. [10]

    Derivation of the double porosity model of single phase flow via homogenization theory

    Todd Arbogast, Jim Douglas, Jr, and Ulrich Hornung. Derivation of the double porosity model of single phase flow via homogenization theory. SIAM Journal on Mathematical Analysis , 21(4):823–836, 1990

    1990

  3. [11]

    Implicit-explicit runge-kutta methods for time-dependent partial differential equations

    Uri M Ascher, Steven J Ruuth, and Raymond J Spiteri. Implicit-explicit runge-kutta methods for time-dependent partial differential equations. Applied Numerical Mathematics , 25(2-3):151–167, 1997

    1997

  4. [12]

    Implicit-explicit methods for time-dependent partial differential equations

    Uri M Ascher, Steven J Ruuth, and Brian TR Wetton. Implicit-explicit methods for time-dependent partial differential equations. SIAM Journal on Numerical Analysis , 32(3):797–823, 1995

    1995

  5. [13]

    Basic concepts in the theory of seepage of homogeneous liquids in fissured rocks [strata]

    GI Barenblatt, Iu P Zheltov, and IN Kochina. Basic concepts in the theory of seepage of homogeneous liquids in fissured rocks [strata]. Journal of applied mathematics and mechanics , 24(5):1286–1303, 1960

    1960

  6. [14]

    Spectral upscaling for graph laplacian problems with application to reservoir simulation

    Andrew T Barker, Chak S Lee, and Panayot S Vassilevski. Spectral upscaling for graph laplacian problems with application to reservoir simulation. SIAM Journal on Scientific Computing , 39(5):S323– S346, 2017

    2017

  7. [15]

    Pyamg: Algebraic multigrid solvers in python

    Nathan Bell, Luke N Olson, and Jacob Schroder. Pyamg: Algebraic multigrid solvers in python. Journal of Open Source Software , 7(72):4142, 2022

    2022

  8. [16]

    Pyamg: algebraic multigrid solvers in python

    Nathan Bell, Luke N Olson, Jacob Schroder, and Ben Southworth. Pyamg: algebraic multigrid solvers in python. Journal of Open Source Software , 8(87):5495, 2023

    2023

  9. [17]

    An iterative method with convergence rate chosen a priori

    Marian Brezina, Caroline Heberton, Jan Mandel, and Petr Vanek. An iterative method with convergence rate chosen a priori. UCD/CCM Report, 140:513–521, 1999

    1999

  10. [18]

    Smoothed aggregation spectral element agglomeration amg: Sa-ρ amge

    Marian Brezina and Panayot S Vassilevski. Smoothed aggregation spectral element agglomeration amg: Sa-ρ amge. In International Conference on Large-Scale Scientific Computing , pages 3–15. Springer, 2011

    2011

  11. [19]

    Linearly implicit runge–kutta methods for advection–reaction–diffusion equations

    Mar ´ ıa Paz Calvo, Javier De Frutos, and Julia Novo. Linearly implicit runge–kutta methods for advection–reaction–diffusion equations. Applied Numerical Mathematics , 37(4):535–549, 2001

    2001

  12. [20]

    Spectral amge ( ρ amge)

    Tim Chartier, Robert D Falgout, VE Henson, J Jones, T Manteuffel, S McCormick, J Ruge, and Panayot S Vassilevski. Spectral amge ( ρ amge). SIAM Journal on Scientific Computing , 25(1):1–26, 2003

    2003

  13. [21]

    Spectral element agglomerate amge

    Timothy Chartier, Robert Falgout, Van Emden Henson, Jim E Jones, Tom A Manteuffel, John W Ruge, Steve F McCormick, and Panayot S Vassilevski. Spectral element agglomerate amge. In Domain decomposition methods in science and engineering XVI , pages 513–521. Springer, 2007

    2007

  14. [22]

    Residual-driven online generalized multiscale finite element methods

    Eric T Chung, Yalchin Efendiev, and Wing Tat Leung. Residual-driven online generalized multiscale finite element methods. Journal of Computational Physics , 302:176–190, 2015

    2015

  15. [23]

    Non-local multi-continua upscaling for flows in heterogeneous fractured media

    Eric T Chung, Yalchin Efendiev, Wing Tat Leung, Maria Vasilyeva, and Yating Wang. Non-local multi-continua upscaling for flows in heterogeneous fractured media. Journal of Computational Physics , 372:22–34, 2018

    2018

  16. [24]

    Chung, Y

    E.T. Chung, Y. Efendiev, and G. Li. An adaptive GMsFEM for high-contrast flow problems. Journal of Computational Physics , 273:54–76, 2014

    2014

  17. [25]

    A robust two-level schwarz preconditioner for sparse matrices

    Hussam Al Daas, Pierre Jolivet, Fr´ ed´ eric Nataf, and Pierre-Henri Tournier. A robust two-level schwarz preconditioner for sparse matrices. arXiv preprint arXiv:2401.03915 , 2024. 24

    2024 arXiv

  18. [26]

    On the coupling of 1d and 3d diffusion-reaction equations: applica- tion to tissue perfusion problems

    Carlo D’angelo and Alfio Quarteroni. On the coupling of 1d and 3d diffusion-reaction equations: applica- tion to tissue perfusion problems. Mathematical Models and Methods in Applied Sciences , 18(08):1481– 1504, 2008

    2008

  19. [27]

    A mixed finite element method for darcy flow in fractured porous media with non-matching grids

    Carlo D’Angelo and Anna Scotti. A mixed finite element method for darcy flow in fractured porous media with non-matching grids. ESAIM: Mathematical Modelling and Numerical Analysis , 46(2):465– 489, 2012

    2012

  20. [28]

    The explicit–implicit–null method: Removing the numerical insta- bility of pdes

    Laurent Duchemin and Jens Eggers. The explicit–implicit–null method: Removing the numerical insta- bility of pdes. Journal of Computational Physics , 263:37–52, 2014

    2014

  21. [29]

    Efendiev, J

    Y. Efendiev, J. Galvis, and P.S. Vassilevski. Spectral element agglomerate algebraic multigrid methods for elliptic problems with high-contrast coefficients. In Domain decomposition methods in science and engineering XIX , volume 78 of Lect. Notes Comput. Sci. Eng. , pages 407...

    2011

  22. [30]

    Efendiev, J

    Y. Efendiev, J. Galvis, and X.H. Wu. Multiscale finite element methods for high-contrast problems using local spectral basis functions. Journal of Computational Physics , 230:937–955, 2011

    2011

  23. [31]

    Generalized multiscale finite element methods (gmsfem)

    Yalchin Efendiev, Juan Galvis, and Thomas Y Hou. Generalized multiscale finite element methods (gmsfem). Journal of computational physics , 251:116–135, 2013

    2013

  24. [32]

    Spectral element agglomerate algebraic multi- grid methods for elliptic problems with high-contrast coefficients

    Yalchin Efendiev, Juan Galvis, and Panayot S Vassilevski. Spectral element agglomerate algebraic multi- grid methods for elliptic problems with high-contrast coefficients. In Domain decomposition methods in science and engineering XIX , pages 407–414. Springer, 2010

    2010

  25. [33]

    Multiscale spectral amge solvers for high- contrast flow problems

    Yalchin Efendiev, Juan Galvis, and Panayot S Vassilevski. Multiscale spectral amge solvers for high- contrast flow problems. ISC-Preprint, Texas A&M University , 2012

    2012

  26. [34]

    Multiscale finite element and domain decomposition methods for high-contrast problems using local spectral basis functions

    Yalchin Efendiev, Juan Galvis, and Xiao-Hui Wu. Multiscale finite element and domain decomposition methods for high-contrast problems using local spectral basis functions. J Comput Phys, 230(4):937–955, 2009

    2009

  27. [35]

    Multiscale finite element methods for high-contrast problems using local spectral basis functions

    Yalchin Efendiev, Juan Galvis, and Xiao-Hui Wu. Multiscale finite element methods for high-contrast problems using local spectral basis functions. Journal of Computational Physics , 230(4):937–955, 2011

    2011

  28. [36]

    On generalizing the algebraic multigrid framework

    Robert D Falgout and Panayot S Vassilevski. On generalizing the algebraic multigrid framework. SIAM Journal on Numerical Analysis , 42(4):1669–1693, 2004

    2004

  29. [37]

    On two-grid convergence estimates

    Robert D Falgout, Panayot S Vassilevski, and Ludmil T Zikatanov. On two-grid convergence estimates. Numerical linear algebra with applications , 12(5-6):471–494, 2005

    2005

  30. [38]

    Matrix heterogeneity effects on gas transport and adsorption in coalbed and shale gas reservoirs

    Ebrahim Fathi and I Y¨ ucel Akkutlu. Matrix heterogeneity effects on gas transport and adsorption in coalbed and shale gas reservoirs. Transport in porous media, 80(2):281–304, 2009

    2009

  31. [39]

    Mass transport of adsorbed-phase in stochastic porous medium with fluctuating porosity field and nonlinear gas adsorption kinetics

    Ebrahim Fathi and I Yucel Akkutlu. Mass transport of adsorbed-phase in stochastic porous medium with fluctuating porosity field and nonlinear gas adsorption kinetics. Transport in porous media, 91(1):5–33, 2012

    2012

  32. [40]

    Multi-component gas transport and adsorption effects during co2 injection and enhanced shale gas recovery

    Ebrahim Fathi and I Yucel Akkutlu. Multi-component gas transport and adsorption effects during co2 injection and enhanced shale gas recovery. International Journal of Coal Geology , 123:52–61, 2014

    2014

  33. [41]

    A reduced model for darcy’s problem in networks of fractures

    Luca Formaggia, Alessio Fumagalli, Anna Scotti, and Paolo Ruffo. A reduced model for darcy’s problem in networks of fractures. ESAIM: Mathematical Modelling and Numerical Analysis , 48(4):1089–1116, 2014. 25

    2014

  34. [42]

    A numerical study of microscale flow behavior in tight gas and shale gas reservoir systems

    CM Freeman, GJ Moridis, and TA Blasingame. A numerical study of microscale flow behavior in tight gas and shale gas reservoir systems. Transport in porous media, 90(1):253–268, 2011

    2011

  35. [43]

    An efficient multiscale preconditioner for large-scale highly heterogeneous flow

    Shubin Fu, Eric Chung, and Lina Zhao. An efficient multiscale preconditioner for large-scale highly heterogeneous flow. SIAM Journal on Scientific Computing , pages S352–S377, 2023

    2023

  36. [44]

    An adaptive preconditioner for three-dimensional single-phase compressible flow in highly heterogeneous porous media

    Shubin Fu, Eric Chung, and Lina Zhao. An adaptive preconditioner for three-dimensional single-phase compressible flow in highly heterogeneous porous media. Multiscale Modeling & Simulation , 22(1):155– 177, 2024

    2024

  37. [45]

    Domain decomposition preconditioners for multiscale flows in high- contrast media

    Juan Galvis and Yalchin Efendiev. Domain decomposition preconditioners for multiscale flows in high- contrast media. Multiscale Modeling & Simulation , 8(4):1461–1483, 2010

    2010

  38. [46]

    Domain decomposition preconditioners for multiscale flows in high contrast media: reduced dimension coarse spaces

    Juan Galvis and Yalchin Efendiev. Domain decomposition preconditioners for multiscale flows in high contrast media: reduced dimension coarse spaces. Multiscale Model. Simul. , 8(5):1621–1644, 2010

    2010

  39. [47]

    Hairer and G

    E. Hairer and G. Wanner. Solving Ordinary Differential Equations II: Stiff and Differential-Algebraic Problems. Springer Series in Computational Mathematics. Springer Berlin Heidelberg, 2010

    2010

  40. [48]

    Iterative multiscale finite- volume method

    Hadi Hajibeygi, Giuseppe Bonfigli, Marc Andre Hesse, and Patrick Jenny. Iterative multiscale finite- volume method. Journal of Computational Physics , 227(19):8604–8621, 2008

    2008

  41. [49]

    Adaptive iterative multiscale finite volume method

    Hadi Hajibeygi and Patrick Jenny. Adaptive iterative multiscale finite volume method. Journal of Computational Physics , 230(3):628–643, 2011

    2011

  42. [50]

    A hierarchical fracture model for the iterative multiscale finite volume method

    Hadi Hajibeygi, Dimitris Karvounis, and Patrick Jenny. A hierarchical fracture model for the iterative multiscale finite volume method. Journal of Computational Physics , 230(24):8729–8743, 2011

    2011

  43. [51]

    Removing the stiffness from interfacial flows with surface tension

    Thomas Y Hou, John S Lowengrub, and Michael J Shelley. Removing the stiffness from interfacial flows with surface tension. Journal of Computational Physics , 114(2):312–338, 1994

    1994

  44. [52]

    Nanoscale gas flow in shale gas sediments

    Farzam Javadpour, D Fisher, M Unsworth, et al. Nanoscale gas flow in shale gas sediments. Journal of Canadian Petroleum Technology, 46(10), 2007

    2007

  45. [53]

    Additive runge–kutta schemes for convection–diffusion– reaction equations

    Christopher A Kennedy and Mark H Carpenter. Additive runge–kutta schemes for convection–diffusion– reaction equations. Applied numerical mathematics , 44(1-2):139–181, 2003

    2003

  46. [54]

    Higher-order additive runge–kutta schemes for ordinary differential equations

    Christopher A Kennedy and Mark H Carpenter. Higher-order additive runge–kutta schemes for ordinary differential equations. Applied numerical mathematics , 136:183–205, 2019

    2019

  47. [55]

    Pressure and fluid-driven fracture propagation in porous media using an adaptive finite element phase field model

    Sanghyun Lee, Mary F Wheeler, and Thomas Wick. Pressure and fluid-driven fracture propagation in porous media using an adaptive finite element phase field model. Computer Methods in Applied Mechanics and Engineering , 305:111–132, 2016

    2016

  48. [56]

    Morphology, genesis, and distribution of nanometer-scale pores in siliceous mudstones of the mississippian barnett shale

    Robert G Loucks, Robert M Reed, Stephen C Ruppel, and Daniel M Jarvie. Morphology, genesis, and distribution of nanometer-scale pores in siliceous mudstones of the mississippian barnett shale. Journal of Sedimentary Research, 79(12):848–861, 2009

    2009

  49. [57]

    Modeling fractures and barriers as interfaces for flow in porous media

    Vincent Martin, J´ erˆ ome Jaffr´ e, and Jean E Roberts. Modeling fractures and barriers as interfaces for flow in porous media. SIAM Journal on Scientific Computing , 26(5):1667–1691, 2005

    2005

  50. [58]

    The partition of unity finite element method: basic theory and applications

    Jens M Melenk and Ivo Babuˇ ska. The partition of unity finite element method: basic theory and applications. Computer methods in applied mechanics and engineering , 139(1-4):289–314, 1996

    1996

  51. [59]

    Convergence analysis of perturbed two-grid and multigrid methods

    Yvan Notay. Convergence analysis of perturbed two-grid and multigrid methods. SIAM journal on numerical analysis, 45(3):1035–1044, 2007. 26

    2007

  52. [60]

    Implicit–explicit runge–kutta schemes and applications to hy- perbolic systems with relaxation

    Lorenzo Pareschi and Giovanni Russo. Implicit–explicit runge–kutta schemes and applications to hy- perbolic systems with relaxation. Journal of Scientific computing , 25:129–155, 2005

    2005

  53. [61]

    Algebraic multigrid

    John W Ruge and Klaus St¨ uben. Algebraic multigrid. InMultigrid methods, pages 73–130. SIAM, 1987

    1987

  54. [62]

    Additive difference schemes and iteration methods for problems of mathematical physics

    AA Samarskii and PN Vabishchevich. Additive difference schemes and iteration methods for problems of mathematical physics. Journal of Mathematical Sciences , 104(6):1657–1660, 2001

    2001

  55. [63]

    The theory of difference schemes , volume 240

    Alexander A Samarskii. The theory of difference schemes , volume 240. CRC Press, 2001

    2001

  56. [64]

    Dimensionally re- duced flow models in fractured porous media: crossings and boundaries

    Nicolas Schwenck, Bernd Flemisch, Rainer Helmig, and Barbara I Wohlmuth. Dimensionally re- duced flow models in fractured porous media: crossings and boundaries. Computational Geosciences, 19(6):1219–1230, 2015

    2015

  57. [65]

    On numerical homogenization of shale gas transport

    Alexey Talonov and Maria Vasilyeva. On numerical homogenization of shale gas transport. Journal of Computational and Applied Mathematics , 301:44–52, 2016

    2016

  58. [66]

    The high-order variable-coefficient explicit-implicit-null method for diffusion and dispersion equations

    Meiqi Tan, Juan Cheng, and Chi-Wang Shu. The high-order variable-coefficient explicit-implicit-null method for diffusion and dispersion equations. Communications on Applied Mathematics and Compu- tation, pages 1–36, 2024

    2024

  59. [67]

    Additive operator-difference schemes

    Petr N Vabishchevich. Additive operator-difference schemes. In Additive Operator-Difference Schemes. de Gruyter, 2013

    2013

  60. [68]

    Explicit-implicit schemes for convection-diffusion-reaction problems

    Petr Nikolaevich Vabishchevich and Mariya Vasil’evna Vasil’eva. Explicit-implicit schemes for convection-diffusion-reaction problems. Numerical Analysis and Applications , 5:297–306, 2012

    2012

  61. [69]

    Explicit–implicit schemes for first-order evolution equations

    PN Vabishchevich. Explicit–implicit schemes for first-order evolution equations. Differential Equations, 56:882–889, 2020

    2020

  62. [70]

    Efficient decoupling schemes for multiscale multicontinuum problems in fractured porous media

    Maria Vasilyeva. Efficient decoupling schemes for multiscale multicontinuum problems in fractured porous media. Journal of Computational Physics , 487:112134, 2023

    2023

  63. [71]

    Implicit-explicit schemes for decoupling multicontinuum problems in porous media

    Maria Vasilyeva. Implicit-explicit schemes for decoupling multicontinuum problems in porous media. arXiv preprint arXiv:2404.16576 , 2024

    2024 arXiv

  64. [72]

    Generalized multiscale finite element method for discrete network (graph) models

    Maria Vasilyeva. Generalized multiscale finite element method for discrete network (graph) models. Journal of Computational and Applied Mathematics , 457:116275, 2025

    2025

  65. [73]

    Multiscale modeling of heat and mass transfer in fractured media for enhanced geothermal systems applications

    Maria Vasilyeva, Masoud Babaei, Eric T Chung, and Denis Spiridonov. Multiscale modeling of heat and mass transfer in fractured media for enhanced geothermal systems applications. Applied Mathematical Modelling, 67:159–178, 2019

    2019

  66. [74]

    Upscaled model for mixed dimensional coupled flow problem in fractured porous media using non-local multicontinuum (nlmc) method

    Maria Vasilyeva, Eric T Chung, Yalchin Efendiev, Wing Tat Leung, and Yating Wang. Upscaled model for mixed dimensional coupled flow problem in fractured porous media using non-local multicontinuum (nlmc) method. In Finite Difference Methods. Theory and Applications: 7th Intern...

    2018

  67. [75]

    Learning macro- scopic parameters in nonlinear multiscale simulations using nonlocal multicontinua upscaling techniques

    Maria Vasilyeva, Wing T Leung, Eric T Chung, Yalchin Efendiev, and Mary Wheeler. Learning macro- scopic parameters in nonlinear multiscale simulations using nonlocal multicontinua upscaling techniques. Journal of Computational Physics , 412:109323, 2020

    2020

  68. [76]

    Multilevel block factorization preconditioners: Matrix-based analysis and algo- rithms for solving finite element equations

    Panayot S Vassilevski. Multilevel block factorization preconditioners: Matrix-based analysis and algo- rithms for solving finite element equations . Springer Science & Business Media, 2008. 27

    2008

  69. [77]

    Coarse spaces by algebraic multigrid: Multigrid convergence and upscaling error estimates

    Panayot S Vassilevski. Coarse spaces by algebraic multigrid: Multigrid convergence and upscaling error estimates. Advances in Adaptive Data Analysis , 3(01n02):229–249, 2011

    2011

  70. [78]

    Partially explicit splitting scheme with explicit-implicit-null method for nonlinear multiscale flow problems

    Yating Wang and Wing Tat Leung. Partially explicit splitting scheme with explicit-implicit-null method for nonlinear multiscale flow problems. Communications in Nonlinear Science and Numerical Simula- tion, page 108094, 2024

    2024

  71. [79]

    The behavior of naturally fractured reservoirs

    JE Warren, P Jj Root, et al. The behavior of naturally fractured reservoirs. Society of Petroleum Engineers Journal, 3(03):245–255, 1963

    1963

  72. [80]

    Permeability of organic-rich shale

    Asana Wasaki, I Yucel Akkutlu, et al. Permeability of organic-rich shale. In SPE Annual Technical Conference and Exhibition . Society of Petroleum Engineers, 2014

    2014

  73. [81]

    A multiple-continuum model for simulating single-phase and multiphase flow in naturally fractured vuggy reservoirs

    Yu-Shu Wu, Yuan Di, Zhijiang Kang, and Perapon Fakcharoenphol. A multiple-continuum model for simulating single-phase and multiphase flow in naturally fractured vuggy reservoirs. Journal of Petroleum Science and Engineering , 78(1):13–22, 2011

    2011

  74. [82]

    Algebraic multigrid methods

    Jinchao Xu and Ludmil Zikatanov. Algebraic multigrid methods. Acta Numerica, 26:591–721, 2017

    2017

  75. [83]

    Convergence analysis of inexact two-grid methods: A theoretical framework

    Xuefeng Xu and Chen-Song Zhang. Convergence analysis of inexact two-grid methods: A theoretical framework. SIAM Journal on Numerical Analysis , 60(1):133–156, 2022

    2022

  76. [84]

    Galvis Y

    J. Galvis Y. Efendiev and P. Vassilevski. Spectral element agglomerate algebraic multigrid methods for elliptic problems with high-contrast coefficients. In Y. Huang, R. Kornhuber, O. Widlund, and J. Xu, editors, Domain Decomposition Methods in Science and Engineering XIX , vo...

    2011

  77. [85]

    A two-grid preconditioner with an adaptive coarse space for flow simulations in highly heterogeneous media

    Yanfang Yang, Shubin Fu, and Eric T Chung. A two-grid preconditioner with an adaptive coarse space for flow simulations in highly heterogeneous media. Journal of Computational Physics , 391:1–13, 2019

    2019

  78. [86]

    An adaptive generalized multiscale finite element method based two-grid preconditioner for large scale high-contrast linear elasticity problems

    Yanfang Yang, Shubin Fu, and Eric T Chung. An adaptive generalized multiscale finite element method based two-grid preconditioner for large scale high-contrast linear elasticity problems. Journal of Scien- tific Computing , 92(1):21, 2022

    2022

  79. [87]

    Discrete fracture-vug network model for modeling fluid flow in fractured vuggy porous media

    Jun Yao, Zhaoqin Huang, Yajun Li, Chenchen Wang, Xinrui Lv, et al. Discrete fracture-vug network model for modeling fluid flow in fractured vuggy porous media. In International oil and gas conference and exhibition in China . Society of Petroleum Engineers, 2010

    2010

  80. [88]

    Numerical simulation of gas transport mechanisms in tight shale gas reservoirs

    Jun Yao, Hai Sun, Dong-yan Fan, Chen-chen Wang, and Zhi-xue Sun. Numerical simulation of gas transport mechanisms in tight shale gas reservoirs. Petroleum Science, 10(4):528–537, 2013

    2013

  81. [89]

    Triple-continuum modeling of shale gas reservoirs considering the effect of kerogen

    Min Zhang, Jun Yao, Hai Sun, Jian-lin Zhao, Dong-yan Fan, Zhao-qin Huang, and Yue-ying Wang. Triple-continuum modeling of shale gas reservoirs considering the effect of kerogen. Journal of Natural Gas Science and Engineering , 24:252–263, 2015. 28

    2015

Pith tools

Reviewed August 12, 2026 · model on record in the stance chip above.