Weak Approximations of the Wright–Fisher Process

In this paper, we construct first- and second-order weak split-step approximations for the solutions of the Wright–Fisher equation. The discretization schemes use the generation of, respectively, two- and three-valued random variables at each discretization step. The accuracy of constructed approximations is illustrated by several simulation examples.

An EHL Extension of the Unsteady FBNS Algorithm

Many engineering applications rely on lubricated gaps where the hydrodynamic pressure distribution is influenced by cavitation phenomena and elastic deformations. To obtain details about the conditions within the lubricated gap, solvers are required that can model cavitation and elastic deformation effects efficiently when a large amount of discretization cells is employed. The presented unsteady EHL-FBNS solver can compute the solution of such large problems that require the consideration of both mass-conserving cavitation and elastic deformation. The execution time of the presented algorithm scales almost with N log(N) where N is the number of computational grid points. A detailed description of the algorithm and the discretized equations is presented. The MATLAB© code is provided in the supplements along with a maintained version on GitHub to encourage its usage and further development. The output of the solver is compared to and validated with simulated and experimental results from the literature to provide a detailed comparison of different discretization schemes of the Couette term in presence of gap height discontinuities. As a final result, the most favourable scheme is identified for the unsteady study of surface textures in ball-on-disc tribometers under severe EHL conditions.

Consistency and convergence for a family of finite volume discretizations of the Fokker–Planck operator

We introduce a family of various finite volume discretization schemes for the Fokker–Planck operator, which are characterized by different Stolarsky weight functions on the edges. This family particularly includes the well-established Scharfetter–Gummel discretization as well as the recently developed square-root approximation (SQRA) scheme. We motivate this family of discretizations both from the numerical and the modeling point of view and provide a uniform consistency and error analysis. Our main results state that the convergence order primarily depends on the quality of the mesh and in second place on the choice of the Stolarsky weights. We show that the Scharfetter–Gummel scheme has the analytically best convergence properties but also that there exists a whole branch of Stolarsky means with the same convergence quality. We show by numerical experiments that for small convection the choice of the optimal representative of the discretization family is highly non-trivial while for large gradients the Scharfetter–Gummel scheme stands out compared to the others.

A monolithic approach toward coupled electrodynamic–thermomechanical problems with regard to weak formulations

Weak formulations of boundary value problems are the basis of various numerical discretization schemes. They are classically derived applying the method of weighted residuals or a variational principle. For electrodynamical and caloric problems, variational approaches are not straightforwardly obtained from physical principles like in mechanics. Weak formulations of Maxwell’s equations and of energy or charge balances thus are frequently derived from the method of weighted residuals or tailored variational approaches. Related formulations of multiphysical problems, combining mechanical balance equations and the axioms of electrodynamics with those of heat conduction, however, raise the additional issue of lacking consistency of physical units, since fluxes of charge and heat intrinsically involve time rates and temperature is only included in the heat balance. In this paper, an energy-based approach toward combined electrodynamic–thermomechanical problems is presented within a classical framework, merging Hamilton’s and Jourdain’s variational principles, originally established in analytical mechanics, to obtain an appropriate basis for a multiphysical formulation. Complementing the Lagrange function by additional potentials of heat flux and electric current and appropriately defining generalized virtual powers of external fields including dissipative processes, a consistent formulation is obtained for the four-field problem and compared to a weighted residuals approach.

Systematic Validation Study of an Unsteady Cavitating Flow over a Hydrofoil Using Conditional Averaging: LES and PIV

We present results of Large-eddy simulations (LES) modeling of steady sheet and unsteady cloud cavitation on a two-dimensional hydrofoil which are validated against Particle image velocimetry (PIV) data. The study is performed for the angle of attack of 9∘ and high Reynolds numbers ReC of the order of 106 providing a strong adverse pressure gradient along the surface. We employ the Schnerr–Sauer and Kunz cavitation models together with the adaptive mesh refinement in critical flow regions where intensive phase transitions occur. Comparison of the LES and visualization results confirms that the flow dynamics is adequately reproduced in the calculations. To correctly match averaged velocity distributions, we propose a new methodology based on conditional averaging of instantaneous velocity fields measured by PIV which only provides information on the liquid phase. This approach leads to an excellent overall agreement between the conditionally averaged fields of the mean velocity and turbulence intensity obtained experimentally and numerically. The benefits of second-order discretization schemes are highlighted as opposed to the lower-order TVD scheme.

Comparison of Various Discretization Schemes for Simulation of Large Field Case Reservoirs Using Unstructured Grids

Solving the equations governing multiphase flow in geological formations involves the generation of a mesh that faithfully represents the structure of the porous medium. This challenging mesh generation task can be greatly simplified by the use of unstructured (tetrahedral) grids that conform to the complex geometric features present in the subsurface. However, running a million-cell simulation problem using an unstructured grid on a real, faulted field case remains a challenge for two main reasons. First, the workflow typically used to construct and run the simulation problems has been developed for structured grids and needs to be adapted to the unstructured case. Second, the use of unstructured grids that do not satisfy the K-orthogonality property may require advanced numerical schemes that preserve the accuracy of the results and reduce potential grid orientation effects. These two challenges are at the center of the present paper. We describe in detail the steps of our workflow to prepare and run a large-scale unstructured simulation of a real field case with faults. We perform the simulation using four different discretization schemes, including the cell-centered Two-Point and Multi-Point Flux Approximation (respectively, TPFA and MPFA) schemes, the cell- and vertex-centered Vertex Approximate Gradient (VAG) scheme, and the cell- and face-centered hybrid Mimetic Finite Difference (MFD) scheme. We compare the results in terms of accuracy, robustness, and computational cost to determine which scheme offers the best compromise for the test case considered here.

AN EFFICIENT APPROACH FOR SOLVING FRACTIONAL VARIABLE ORDER REACTION SUB-DIFFUSION BASED ON HERMITE FORMULA

Anomalous Reaction-Sub-diffusion equations play an important role transferred in a lot of our daily applications in our life, especially in applied chemistry. In the presented work, a modified type of these models is considered which is the Reaction-Sub-diffusion equation of variable order, the linear and nonlinear models and we will refer to it by VORSDE. An accurate technique depends on a mix of the finite difference methods (FDM) together with Hermite formula is introduced to study these important types of anomalous equations. Regarding the analysis of the stability for the mentioned, it is done using the variable Von-Neumann technique; also the convergent analysis is introduced. As a result of the previous steps, we derived a stability condition which will be held for many discretization schemes of the variable order derivative and some other parameters and we checked it numerically.