基礎講座 電気泳動(Electrophoresis)

A criterion is given for the convergence of numerical solutions of the Navier-Stokes equations in two dimensions under steady conditions. The criterion applies to all cases, of steady viscous flow in two dimensions and shows that if the local ' mesh Reynolds number ', based on the size of the mesh used in the solution, exceeds a certain fixed value, the numerical solution will not converge.

On the Navier-Stokes equations

The potential for sea ice-albedo feedback to give rise to nonlinear climate change in the Arctic Ocean – defined as a nonlinear relationship between polar and global temperature change or, equivalently, a time-varying polar amplification – is explored in IPCC AR4 climate models. Five models supplying SRES A1B ensembles for the 21 st century are examined and very linear relationships are found between polar and global temperatures (indicating linear Arctic Ocean climate change), and between polar temperature and albedo (the potential source of nonlinearity). Two of the climate models have Arctic Ocean simulations that become annually sea ice-free under the stronger CO 2 increase to quadrupling forcing. Both of these runs show increases in polar amplification at polar temperatures above-5 o C and one exhibits heat budget changes that are consistent with the small ice cap instability of simple energy balance models. Both models show linear warming up to a polar temperature of-5 o C, well above the disappearance of their September ice covers at about-9 o C. Below-5 o C, surface albedo decreases smoothly as reductions move, progressively, to earlier parts of the sunlit period. Atmospheric heat transport exerts a strong cooling effect during the transition to annually ice-free conditions. Specialized experiments with atmosphere and coupled models show that the main damping mechanism for sea ice region surface temperature is reduced upward heat flux through the adjacent ice-free oceans resulting in reduced atmospheric heat transport into the region.

Geophysical fluid dynamics.

Finite element methods for Maxwell's equations.

The classical Poisson-Boltzmann (PB) theory of electrolytes assumes a dilute solution of point charges with mean-field electrostatic forces. Even for very dilute solutions, however, it predicts absurdly large ion concentrations (exceeding close packing) for surface potentials of only a few tenths of a volt, which are often exceeded, e.g. in microfluidic pumps and electrochemical sensors. Since the 1950s, several modifications of the PB equation have been proposed to account for the finite size of ions in equilibrium, but in this two-part series, we consider steric effects on diffuse charge dynamics (in the absence of electro-osmotic flow). In this first part, we review the literature and analyze two simple models for the charging of a thin double layer, which must form a condensed layer of close-packed ions near the surface at high voltage. A surprising prediction is that the differential capacitance typically varies non-monotonically with the applied voltage, and thus so does the response time of an electrolytic system. In PB theory, the capacitance blows up exponentially with voltage, but steric effects actually cause it to decrease above a threshold voltage where ions become crowded near the surface. Other nonlinear effects in PB theory are also strongly suppressed by steric effects: The net salt adsorption by the double layers in response to the applied voltage is greatly reduced, and so is the tangential "surface conduction" in the diffuse layer, to the point that it can often be neglected compared to bulk conduction (small Dukhin number).

Steric effects in the dynamics of electrolytes at large applied voltages: I. Double-layer charging

The divergence constraint of the incompressible Navier--Stokes equations is revisited in the mixed finite element framework. While many stable and convergent mixed elements have been developed throughout the past four decades, most classical methods relax the divergence constraint and only enforce the condition discretely. As a result, these methods introduce a pressure-dependent consistency error which can potentially pollute the computed velocity. These methods are not robust in the sense that a contribution from the right-hand side, which influences only the pressure in the continuous equations, impacts both velocity and pressure in the discrete equations. This article reviews the theory and practical implications of relaxing the divergence constraint. Several approaches for improving the discrete mass balance or even for computing divergence-free solutions will be discussed: grad-div stabilization, higher order mixed methods derived on the basis of an exact de Rham complex, $H(div)$-conforming finite ...

On the Divergence Constraint in Mixed Finite Element Methods for Incompressible Flows

Pressure-robustness and discrete Helmholtz projectors in mixed finite element methods for the incompressible Navier–Stokes equations

Classical inf-sup stable mixed finite elements for the incompressible (Navier--)Stokes equations are not pressure-robust, i.e., their velocity errors depend on the continuous pressure. However, a modification only in the right-hand side of a Stokes discretization is able to reestablish pressure-robustness, as shown recently for several inf-sup stable Stokes elements with discontinuous discrete pressures. In this contribution, this idea is extended to low and high order Taylor--Hood and mini elements, which have continuous discrete pressures. For the modification of the right-hand side a velocity reconstruction operator is constructed that maps discretely divergence-free test functions to exactly divergence-free ones. The reconstruction is based on local $H(\mathrm{div})$-conforming flux equilibration on vertex patches, and fulfills certain orthogonality properties to provide consistency and optimal a priori error estimates. Numerical examples for the incompressible Stokes and Navier--Stokes equations conf...

Divergence-free Reconstruction Operators for Pressure-Robust Stokes Discretizations with Continuous Pressure Finite Elements

On velocity errors due to irrotational forces in the Navier-Stokes momentum balance

Recently, it was understood how to repair a certain L 2 -orthogonality of discretely divergence-free vector fields and gradient fields such that the velocity error of inf–sup stable discretizations for the incompressible Stokes equations becomes pressure-independent. These new ‘pressure-robust’ Stokes discretizations deliver a small velocity error, whenever the continuous velocity field can be well approximated on a given grid. On the contrary, classical inf–sup stable Stokes discretizations can guarantee a small velocity error only when both the velocity and the pressure field can be approximated well, simultaneously. In this contribution, ‘pressure-robustness’ is extended to the time-dependent Navier–Stokes equations. In particular, steady and time-dependent potential flows are shown to build an entire class of benchmarks, where pressure-robust discretizations can outperform classical approaches significantly. Speedups will be explained by a new theoretical concept, the ‘discrete Helmholtz projector’ of an inf–sup stable discretization. Moreover, different discrete nonlinear convection terms are discussed, and skew-symmetric pressure-robust discretizations are proposed.

Christian Merdon

Papers

On the Divergence Constraint in Mixed Finite Element Methods for Incompressible Flows

Pressure-robustness and discrete Helmholtz projectors in mixed finite element methods for the incompressible Navier–Stokes equations

Divergence-free Reconstruction Operators for Pressure-Robust Stokes Discretizations with Continuous Pressure Finite Elements

On velocity errors due to irrotational forces in the Navier-Stokes momentum balance

Pressure-robustness and discrete Helmholtz projectors in mixed finite element methods for the incompressible Navier--Stokes equations