Perturbative corrections for the scaling of heat transport in a Hele-Shaw geometry and its application to geological vertical fractures

2019 ◽  
Vol 864 ◽  
pp. 746-767 ◽  
Author(s):  
Juvenal A. Letelier ◽  
Nicolás Mujica ◽  
Jaime H. Ortega
Keyword(s):  
Heat Transport ◽  
Dissipation Rate ◽  
Periodic Boundary ◽  
Periodic Boundary Conditions ◽  
Thermal Dissipation ◽  
Anisotropy Ratio ◽  
Theoretical Relation ◽  
Fracture Scale ◽  
Perturbative Corrections ◽  
Cell Anisotropy

In this work, we investigate numerically the perturbative effects of cell aperture in heat transport and thermal dissipation rate for a vertical Hele-Shaw geometry, which is used as an analogue representation of a planar vertical fracture at the laboratory scale. To model the problem, we derive a two-dimensional set of equations valid for this geometry. For Hele-Shaw cells heated from below and above, with periodic boundary conditions in the horizontal direction, the model gives new nonlinear scalings for both the time-averaged Nusselt number $\langle Nu\rangle _{\unicode[STIX]{x1D70F}}$ and dimensionless mean thermal dissipation rate $\langle \unicode[STIX]{x1D717}\rangle _{\unicode[STIX]{x1D70F}}$ in the high-Rayleigh regime. We demonstrate that $\langle Nu\rangle _{\unicode[STIX]{x1D70F}}$ and $\langle \unicode[STIX]{x1D717}\rangle _{\unicode[STIX]{x1D70F}}$ depend upon the cell anisotropy ratio $\unicode[STIX]{x1D716}$, which measures the ratio between the cell gap and height. We show that $\langle Nu\rangle _{\unicode[STIX]{x1D70F}}$ values in the high-Rayleigh regime decrease when $\unicode[STIX]{x1D716}$ grows, supporting the field observations at the fracture scale. When $\unicode[STIX]{x1D716}\ll 1$, our results are in agreement with the scalings found using the Darcy model. The numerical results satisfy the theoretical relation $\langle Nu\rangle _{\unicode[STIX]{x1D70F}}=Ra\langle \unicode[STIX]{x1D717}\rangle _{\unicode[STIX]{x1D70F}}$, which is obtained from the model. This latter relation is valid for all values of Rayleigh number considered. The perturbative effects of cell aperture are observed only in the exponents of the scalings $\langle Nu\rangle _{\unicode[STIX]{x1D70F}}\sim Ra^{\unicode[STIX]{x1D6FE}(\unicode[STIX]{x1D716})}$ and $\langle \unicode[STIX]{x1D717}\rangle _{\unicode[STIX]{x1D70F}}\sim Ra^{\unicode[STIX]{x1D6FE}(\unicode[STIX]{x1D716})-1}$.

An algorithm for solving the Navier–Stokes equations with shear-periodic boundary conditions and its application to homogeneously sheared turbulence

2017 ◽  
Vol 833 ◽  
pp. 687-716 ◽  
Author(s):  
M. Houssem Kasbaoui ◽  
Ravi G. Patel ◽  
Donald L. Koch ◽  
Olivier Desjardins
Keyword(s):  
Reynolds Number ◽  
Boundary Conditions ◽  
Dissipation Rate ◽  
Periodic Boundary ◽  
Stokes Equations ◽  
Transverse Velocity ◽  
Periodic Boundary Conditions ◽  
Integration Time ◽  
Asymptotic State ◽  
Self Similar

Simulations of homogeneously sheared turbulence (HST) are conducted until a universal self-similar state is established at the long non-dimensional time $\unicode[STIX]{x1D6E4}t=20$, where $\unicode[STIX]{x1D6E4}$ is the shear rate. The simulations are enabled by a new robust and discretely conservative algorithm. The method solves the governing equations in physical space using the so-called shear-periodic boundary conditions. Convection by the mean homogeneous shear flow is treated implicitly in a split step approach. An iterative Crank–Nicolson time integrator is chosen for robustness and stability. The numerical strategy captures without distortion the Kelvin modes, rotating waves that are fundamental to homogeneously sheared flows and are at the core of rapid distortion theory. Three direct numerical simulations of HST with the initial Taylor scale Reynolds number $Re_{\unicode[STIX]{x1D706}0}=29$ and shear numbers of $S_{0}^{\ast }=\unicode[STIX]{x1D6E4}q^{2}/\unicode[STIX]{x1D716}=3$, 15 and 27 are performed on a $2048\times 1024\times 1024$ grid. Here, $\unicode[STIX]{x1D716}$ is the dissipation rate and $1/2q^{2}$ is the turbulent kinetic energy. The long integration time considered allows the establishment of a self-similar state observed in experiments but often absent from simulations conducted over shorter times. The asymptotic state appears to be universal with a long time production to dissipation rate ${\mathcal{P}}/\unicode[STIX]{x1D716}\sim 1.5$ and shear number $S^{\ast }\sim 10$ in agreement with experiments. While the small scales exhibit strong anisotropy increasing with initial shear number, the skewness of the transverse velocity derivative decreases with increasing Reynolds number.

Molecular Dynamics Using Non-Variational Polarizable Force Fields: Theory, Periodic Boundary Conditions Implementation and Application to the Bond Capacity Model

2019 ◽  
Author(s):  
Pier Paolo Poier ◽  
Louis Lagardere ◽  
Jean-Philip Piquemal ◽  
Frank Jensen
Keyword(s):  
Boundary Conditions ◽  
Periodic Boundary ◽  
Force Fields ◽  
Periodic Boundary Conditions ◽  
Computational Effort ◽  
Polarizable Force Fields ◽  
Energy Models ◽  
Capacity Model ◽  
Energy Gradient ◽  
Polarization Models

<div> <div> <div> <p>We extend the framework for polarizable force fields to include the case where the electrostatic multipoles are not determined by a variational minimization of the electrostatic energy. Such models formally require that the polarization response is calculated for all possible geometrical perturbations in order to obtain the energy gradient required for performing molecular dynamics simulations. </p><div> <div> <div> <p>By making use of a Lagrange formalism, however, this computational demanding task can be re- placed by solving a single equation similar to that for determining the electrostatic variables themselves. Using the recently proposed bond capacity model that describes molecular polarization at the charge-only level, we show that the energy gradient for non-variational energy models with periodic boundary conditions can be calculated with a computational effort similar to that for variational polarization models. The possibility of separating the equation for calculating the electrostatic variables from the energy expression depending on these variables without a large computational penalty provides flexibility in the design of new force fields. </p><div><div><div> </div> </div> </div> <p> </p><div> <div> <div> <p>variables themselves. Using the recently proposed bond capacity model that describes molecular polarization at the charge-only level, we show that the energy gradient for non-variational energy models with periodic boundary conditions can be calculated with a computational effort similar to that for variational polarization models. The possibility of separating the equation for calculating the electrostatic variables from the energy expression depending on these variables without a large computational penalty provides flexibility in the design of new force fields. </p> </div> </div> </div> </div> </div> </div> </div> </div> </div>

scholarly journals Global classical solutions to the elastodynamic equations with damping

2021 ◽  
Vol 2021 (1) ◽  
Author(s):  
Mengmeng Liu ◽  
Xueyun Lin
Keyword(s):  
Initial Data ◽  
Periodic Boundary ◽  
Energy Estimate ◽  
Periodic Boundary Conditions ◽  
Deformation Tensor ◽  
Classical Solutions ◽  
Small Initial Data ◽  
Weighted Energy ◽  
Global Classical Solutions ◽  
Weighted Energy Estimate

AbstractIn this paper, we show the global existence of classical solutions to the incompressible elastodynamics equations with a damping mechanism on the stress tensor in dimension three for sufficiently small initial data on periodic boxes, that is, with periodic boundary conditions. The approach is based on a time-weighted energy estimate, under the assumptions that the initial deformation tensor is a small perturbation around an equilibrium state and the initial data have some symmetry.

scholarly journals Normal periodic solutions for the fractional abstract Cauchy problem

2021 ◽  
Vol 2021 (1) ◽  
Author(s):  
Jennifer Bravo ◽  
Carlos Lizama
Keyword(s):  
Banach Space ◽  
Cauchy Problem ◽  
Periodic Solution ◽  
Periodic Boundary ◽  
Abstract Cauchy Problem ◽  
Periodic Boundary Conditions ◽  
Sectorial Operator ◽  
Closed Linear Operator ◽  
Resolvent Set ◽  
Spectral Angle

AbstractWe show that if A is a closed linear operator defined in a Banach space X and there exist $t_{0} \geq 0$ t 0 ≥ 0 and $M>0$ M > 0 such that $\{(im)^{\alpha }\}_{|m|> t_{0}} \subset \rho (A)$ { ( i m ) α } | m | > t 0 ⊂ ρ ( A ) , the resolvent set of A, and $$ \bigl\Vert (im)^{\alpha }\bigl(A+(im)^{\alpha }I \bigr)^{-1} \bigr\Vert \leq M \quad \text{ for all } \vert m \vert > t_{0}, m \in \mathbb{Z}, $$ ∥ ( i m ) α ( A + ( i m ) α I ) − 1 ∥ ≤ M  for all  | m | > t 0 , m ∈ Z , then, for each $\frac{1}{p}<\alpha \leq \frac{2}{p}$ 1 p < α ≤ 2 p and $1< p < 2$ 1 < p < 2 , the abstract Cauchy problem with periodic boundary conditions $$ \textstyle\begin{cases} _{GL}D^{\alpha }_{t} u(t) + Au(t) = f(t), & t \in (0,2\pi ); \\ u(0)=u(2\pi ), \end{cases} $$ { D t α G L u ( t ) + A u ( t ) = f ( t ) , t ∈ ( 0 , 2 π ) ; u ( 0 ) = u ( 2 π ) , where $_{GL}D^{\alpha }$ D α G L denotes the Grünwald–Letnikov derivative, admits a normal 2π-periodic solution for each $f\in L^{p}_{2\pi }(\mathbb{R}, X)$ f ∈ L 2 π p ( R , X ) that satisfies appropriate conditions. In particular, this happens if A is a sectorial operator with spectral angle $\phi _{A} \in (0, \alpha \pi /2)$ ϕ A ∈ ( 0 , α π / 2 ) and $\int _{0}^{2\pi } f(t)\,dt \in \operatorname{Ran}(A)$ ∫ 0 2 π f ( t ) d t ∈ Ran ( A ) .

scholarly journals Sharp Lyapunov-type inequalities for second-order half-linear difference equations with different kinds of boundary conditions

2021 ◽  
Vol 115 (3) ◽  
Author(s):  
Robert Stegliński
Keyword(s):  
Difference Equation ◽  
Boundary Conditions ◽  
Difference Equations ◽  
Periodic Boundary ◽  
Periodic Boundary Conditions ◽  
Second Order ◽  
Linear Difference Equations ◽  
Order Difference ◽  
Second Order Difference Equation ◽  
Order Difference Equation

AbstractIn this work, we establish optimal Lyapunov-type inequalities for the second-order difference equation with p-Laplacian $$\begin{aligned} \Delta (\left| \Delta u(k-1)\right| ^{p-2}\Delta u(k-1))+a(k)\left| u(k)\right| ^{p-2}u(k)=0 \end{aligned}$$ Δ ( Δ u ( k - 1 ) p - 2 Δ u ( k - 1 ) ) + a ( k ) u ( k ) p - 2 u ( k ) = 0 with Dirichlet, Neumann, mixed, periodic and anti-periodic boundary conditions.

Sign in / Sign up

Export Citation Format

Share Document