Using extensive numerical simulation of the Navier–Stokes equations, we study the transition from the Darcy’s law for slow flow of fluids through a disordered porous medium to the nonlinear flow regime in which the effect of inertia cannot be neglected. The porous medium is represented by two-dimensional slices of a three-dimensional image of a sandstone. We study the problem over wide ranges of porosity and the Reynolds number, as well as two types of boundary conditions, and compute essential features of fluid flow, namely, the strength of the vorticity, the effective permeability of the pore space, the frictional drag, and the relationship between the macroscopic pressure gradient
For polyhedral constrained optimization problems and a feasible point
- Award ID(s):
- 2031213
- NSF-PAR ID:
- 10415087
- Publisher / Repository:
- Springer Science + Business Media
- Date Published:
- Journal Name:
- Mathematical Programming
- Volume:
- 205
- Issue:
- 1-2
- ISSN:
- 0025-5610
- Format(s):
- Medium: X Size: p. 107-134
- Size(s):
- p. 107-134
- Sponsoring Org:
- National Science Foundation
More Like this
-
Abstract and the fluid velocity$${\varvec{\nabla }}P$$ v . The results indicate that when the Reynolds number Re is low enough that the Darcy’s law holds, the magnitude of the vorticity is nearly zero. As Re increases, however, so also does$$\omega _z$$ , and its rise from nearly zero begins at the same Re at which the Darcy’s law breaks down. We also show that a nonlinear relation between the macroscopic pressure gradient and the fluid velocity$$\omega _z$$ v , given by, , provides accurate representation of the numerical data, where$$-{\varvec{\nabla }}P=(\mu /K_e)\textbf{v}+\beta _n\rho |\textbf{v}|^2\textbf{v}$$ and$$\mu$$ are the fluid’s viscosity and density,$$\rho$$ is the effective Darcy permeability in the linear regime, and$$K_e$$ is a generalized nonlinear resistance. Theoretical justification for the relation is presented, and its predictions are also compared with those of the Forchheimer’s equation.$$\beta _n$$ -
Abstract It has been recently established in David and Mayboroda (Approximation of green functions and domains with uniformly rectifiable boundaries of all dimensions.
arXiv:2010.09793 ) that on uniformly rectifiable sets the Green function is almost affine in the weak sense, and moreover, in some scenarios such Green function estimates are equivalent to the uniform rectifiability of a set. The present paper tackles a strong analogue of these results, starting with the “flagship degenerate operators on sets with lower dimensional boundaries. We consider the elliptic operators associated to a domain$$L_{\beta ,\gamma } =- {\text {div}}D^{d+1+\gamma -n} \nabla $$ with a uniformly rectifiable boundary$$\Omega \subset {\mathbb {R}}^n$$ of dimension$$\Gamma $$ , the now usual distance to the boundary$$d < n-1$$ given by$$D = D_\beta $$ for$$D_\beta (X)^{-\beta } = \int _{\Gamma } |X-y|^{-d-\beta } d\sigma (y)$$ , where$$X \in \Omega $$ and$$\beta >0$$ . In this paper we show that the Green function$$\gamma \in (-1,1)$$ G for , with pole at infinity, is well approximated by multiples of$$L_{\beta ,\gamma }$$ , in the sense that the function$$D^{1-\gamma }$$ satisfies a Carleson measure estimate on$$\big | D\nabla \big (\ln \big ( \frac{G}{D^{1-\gamma }} \big )\big )\big |^2$$ . We underline that the strong and the weak results are different in nature and, of course, at the level of the proofs: the latter extensively used compactness arguments, while the present paper relies on some intricate integration by parts and the properties of the “magical distance function from David et al. (Duke Math J, to appear).$$\Omega $$ -
Abstract For a given material,
controllable deformations are those deformations that can be maintained in the absence of body forces and by applying only boundary tractions. For a given class of materials,universal deformations are those deformations that are controllable for any material within the class. In this paper, we characterize the universal deformations in compressible isotropic implicit elasticity defined by solids whose constitutive equations, in terms of the Cauchy stress and the left Cauchy-Green strain$$\varvec{\sigma }$$ , have the implicit form$$\textbf{b}$$ . We prove that universal deformations are homogeneous. However, an important observation is that, unlike Cauchy (and Green) elasticity, not every homogeneous deformation is constitutively admissible for a given implicit-elastic solid. In other words, the set of universal deformations is material-dependent, yet it remains a subset of homogeneous deformations.$$\varvec{\textsf{f}}(\varvec{\sigma },\textbf{b})=\textbf{0}$$ -
Abstract We introduce a distributional Jacobian determinant
in dimension two for the nonlinear complex gradient\det DV_{\beta}(Dv) for anyV_{\beta}(Dv)=\lvert Dv\rvert^{\beta}(v_{x_{1}},-v_{x_{2}}) , whenever\beta>-1 andv\in W^{1\smash{,}2}_{\mathrm{loc}} .This is new when\beta\lvert Dv\rvert^{1+\beta}\in W^{1\smash{,}2}_{\mathrm{loc}} .Given any planar ∞-harmonic function 𝑢, we show that such distributional Jacobian determinant\beta\neq 0 is a nonnegative Radon measure with some quantitative local lower and upper bounds.We also give the following two applications.\det DV_{\beta}(Du) Applying this result with
, we develop an approach to build up a Liouville theorem, which improves that of Savin.Precisely, if 𝑢 is an ∞-harmonic function in the whole\beta=0 with\mathbb{R}^{2} then\liminf_{R\to\infty}\inf_{c\in\mathbb{R}}\frac{1}{R}\barint_{B(0,R)}\lvert u(x)-c\rvert\,dx<\infty, for someu=b+a\cdot x andb\in\mathbb{R} .a\in\mathbb{R}^{2} Denoting by
the 𝑝-harmonic function having the same nonconstant boundary condition as 𝑢, we show thatu_{p} as\det DV_{\beta}(Du_{p})\to\det DV_{\beta}(Du) in the weak-⋆ sense in the space of Radon measure.Recall thatp\to\infty is always quasiregular mappings, butV_{\beta}(Du_{p}) is not in general.V_{\beta}(Du) -
Abstract A classical parking function of length
n is a list of positive integers whose nondecreasing rearrangement$$(a_1, a_2, \ldots , a_n)$$ satisfies$$b_1 \le b_2 \le \cdots \le b_n$$ . The convex hull of all parking functions of length$$b_i \le i$$ n is ann -dimensional polytope in , which we refer to as the classical parking function polytope. Its geometric properties have been explored in Amanbayeva and Wang (Enumer Combin Appl 2(2):Paper No. S2R10, 10, 2022) in response to a question posed by Stanley (Amer Math Mon 127(6):563–571, 2020). We generalize this family of polytopes by studying the geometric properties of the convex hull of$${\mathbb {R}}^n$$ -parking functions for$${\textbf{x}}$$ , which we refer to as$${\textbf{x}}=(a,b,\dots ,b)$$ -parking function polytopes. We explore connections between these$${\textbf{x}}$$ -parking function polytopes, the Pitman–Stanley polytope, and the partial permutahedra of Heuer and Striker (SIAM J Discrete Math 36(4):2863–2888, 2022). In particular, we establish a closed-form expression for the volume of$${\textbf{x}}$$ -parking function polytopes. This allows us to answer a conjecture of Behrend et al. (2022) and also obtain a new closed-form expression for the volume of the convex hull of classical parking functions as a corollary.$${\textbf{x}}$$