Conservative Finite Difference Scheme Research Articles

Scalable implicit solvers with dynamic mesh adaptation for a relativistic drift-kinetic Fokker–Planck–Boltzmann model

In this work we consider a relativistic drift-kinetic model for runaway electrons along with a Fokker–Planck operator for small-angle Coulomb collisions, a radiation damping operator, and a secondary knock-on (Boltzmann) collision source. We develop a new scalable fully implicit solver utilizing finite volume and conservative finite difference schemes and dynamic mesh adaptivity. A new data management framework in the PETSc library based on the p4est library is developed to enable simulations with dynamic adaptive mesh refinement (AMR), distributed memory parallelization, and dynamic load balancing of computational work. This framework and the runaway electron solver building on the framework are able to dynamically capture both bulk Maxwellian at the low-energy region and a runaway tail at the high-energy region. To effectively capture features via the AMR algorithm, a new AMR indicator prediction strategy is proposed that is performed alongside the implicit time evolution of the solution. This strategy is complemented by the introduction of computationally cheap feature-based AMR indicators that are analyzed theoretically. Numerical results quantify the advantages of the prediction strategy in better capturing features compared with nonpredictive strategies; and we demonstrate trade-offs regarding computational costs. The robustness with respect to model parameters, algorithmic scalability, and parallel scalability are demonstrated through several benchmark problems including manufactured solutions and solutions of different physics models. We focus on demonstrating the advantages of using implicit time stepping and AMR for runaway electron simulations.

Invariant conservative finite-difference schemes for the one-dimensional shallow water magnetohydrodynamics equations in Lagrangian coordinates

Invariant finite-difference schemes for the one-dimensional shallow water equations in the presence of a magnetic field for various bottom topographies are constructed. Based on the results of the group classification recently carried out by the authors, finite-difference analogues of the conservation laws of the original differential model are obtained. Some typical problems are considered numerically, for which a comparison is made between the cases of a magnetic field presence and when it is absent (the standard shallow water model). The invariance of difference schemes in Lagrangian coordinates and the energy preservation on the obtained numerical solutions are also discussed.

Conservative higher-order finite difference scheme for the coupled nonlinear Schrödinger equations

This paper introduces a conservative higher-order finite difference scheme for solving the coupled nonlinear Schrödinger equations. The Crank–Nicolson method is employed to discretize time derivatives and the sixth-order difference operator is used to discretize space derivatives, correspondingly, the resulting difference scheme has second-order accuracy in time and sixth-order accuracy in space. By utilizing the discrete energy method, the conservation of discrete mass and energy, the boundedness, existence and uniqueness of solution, unconditional stability and the convergence of the new scheme are proved. Then making use of Richardson extrapolation, the time accuracy is increased to the fourth order. Finally, the numerical experiments are conducted to validate the theoretical results presented in the paper.

THE FORMATION AND THE EVOLUTION OF LARGE-SCALE VORTEX STRUCTURES IN STELLAR ACCRETION DISKS

Explaining the causes of angular momentum transfer in accretion stellar disks is an important astrophysical problem, since it is the process that determines the rate of accretion of matter onto the central gravitating body. Previously, within the framework of a two-dimensional approach, it was shown that the introduction of small perturbations into the flow of disk matter leads to the appearance of shear instability. This process is accompanied by the development of large-scale vortex structures. Their movement and evolution lead to a redistribution of angular momentum in the accretion disk. The action of the described mechanism was previously studied numerically only within a two-dimensional approximation, so the goal of the current work is to carry out full-scale three-dimensional modeling. The processes under study are described within the framework of the system of equations of ideal gas dynamics. The article briefly describes the method for their numerical integration, which is based on a conservative finite-difference scheme and the solution of the Riemann problem. The initial data is a stationary toroidal gas state surrounded by a matter with low density and pressure. At the next step, small perturbations of one of the gas-dynamic variables are introduced. The modeling and analysis of the results of numerical calculations show the emergence of vortex structures in the shear flow of a three-dimensional accretion disk. Their movement is accompanied by a redistribution of matter and angular momentum in the volume of the disk, leading to accretion of matter onto the central body.

Hybrid conservative central/WENO finite difference scheme for two-dimensional detonation problems

In the present work, a hybrid finite difference scheme is proposed and tested for the study of two-dimensional detonation problems. The hybrid scheme consists of a fifth-order weighted essentially nonoscillatory (WENO) method for discontinuous regions and a fourth or sixth-order robust conservative central finite difference scheme for the remaining domain. An arbitrary shock sensor is used to identify the discontinuities. The proposed shock sensor is based on the absolute value of the density gradient with an arbitrary threshold computed at each time step using a fast image segmentation technique. The sensor is complemented with a Ducros sensor to avoid the selection of vortices as discontinuities. The hybrid scheme is tested for several benchmark examples, including a two-dimensional transverse detonation wave, detonation wave diffraction, and a multiple obstacles test case. The obtained results show that the proposed sensor detects discontinuities adequately and the hybrid schemes give the same results as the fifth-order WENO in faster computational time.

Analysis of two conservative fourth-order compact finite difference schemes for the Klein-Gordon-Zakharov system in the subsonic limit regime

We propose two conservative fourth-order compact finite difference (CFD4C) schemes and give the rigorous error analysis for the Klein-Gordon-Zakharov system (KGZS) with ε∈(0,1] being a small parameter. In the case 0<ε≪1, i.e., the subsonic limit regime, the solution of KGZS propagates waves with wavelength O(ε) in time and O(1) in space, respectively, with amplitude at O(εα†) with α†=min⁡{α,β+1,2}, where α and β describe the incompatibility between the initial data of the KGZS and the limiting Klein-Gordon equation (KGE) as ε→0+ and satisfy α≥0,β+1≥0. Different from the standard KGZS, the oscillation in time becomes the main difficulty in constructing numerical schemes and giving the error analysis for KGZS in this regime. We prove that the schemes conserve the discrete energies. By defining a discrete ‘energy’, applying cut-off technique and considering the convergence rate of the KGZS to the limiting KGE, we obtain two independent error estimates with the bounds at O(h4/ε1−α⁎+τ2/ε3−α†) and O(h4+τ2+εα⁎) where α⁎=min⁡{1,α,1+β}, τ is time step and h is mesh size. Hence, we get uniformly accurate error estimates with the bounds at O(h4+τ2/(4−α†)) when α≥1 and β≥0, and O(h4α⁎+τ2α⁎/3) when 0<α<1 or −1<β<0. While for α=0 or β=−1, the meshing requirements of CFD4Cs are h=O(ε1/4) and τ=O(ε3/2) for 0<ε≪1. The compact schemes provide much better spatial resolution compared with the usual second order finite difference schemes. Our numerical results support energy-conserving properties and the results of the convergence analysis.

Reformulated Dissipation for the Free-Stream Preserving of the Conservative Finite Difference Schemes on Curvilinear Grids

Reformulated Dissipation for the Free-Stream Preserving of the Conservative Finite Difference Schemes on Curvilinear Grids

A nonlinear heat equation arising from automated-vehicle traffic flow models

In this paper, a new nonlinear heat equation is studied that arises as a model of the collective behaviour of automated vehicles. The properties of the solutions of this equation are studied by introducing the appropriate notion of a weak solution that requires certain entropy-like conditions. To obtain an approximation of the solution of the nonlinear heat equation, a new conservative first-order finite difference scheme is proposed that respects the corresponding entropy conditions, and certain links between the weak solution and the numerical scheme are provided. Finally, a traffic simulation scenario and a comparison with the Lighthill–Whitham–Richards (LWR) model are provided, illustrating the benefits of the use of automated vehicles.

Conservative second-order accurate finite-difference scheme for the coupled Maxwell-Dirac equations

The recent development of nanoelectronic devices that incorporate Dirac materials has highly increased the need for adequate simulation and modelling tools. This paper introduces an accurate, multiphysics finite-difference time-domain method to solve the pertinent Maxwell-Dirac equations. The stability criterion for the Dirac equation with electromagnetic fields is derived, which reduces to the Courant-Friedrichs-Lewy condition in the absence of electromagnetic fields. Validation examples show the second-order accuracy of the novel fully coupled Maxwell-Dirac scheme and illustrate that total norm and energy are caonserved within a relative error of order 10−4. The method is applied to a ZrTe5 waveguide and it is found that even at low field strengths, the charge carriers can be accelerated to 80% of the Fermi velocity. Furthermore, the flexibility of the advocated method allows for the seamless integration into existing computational electromagnetics frameworks and the possible extension to higher-order schemes.

An efficient fourth-order accurate conservative scheme for Riesz space fractional Schrödinger equation with wave operator

In this study, we investigate a high-order accurate conservative finite difference scheme by utilizing a fourth-order fractional central finite difference method for the two-dimensional Riesz space-fractional nonlinear Schrödinger wave equation. The conservation laws of the discrete difference scheme are shown. Meanwhile, the exactness, uniqueness, and prior estimate of the numerical solution are rigorously established. Then, it is proved that the proposed scheme is unconditionally convergent in the discrete L2 and Hγ/2 norm, where γ is a fractional order. Furthermore, we demonstrate that when the fractional order γ and the spatial grid number J increase, the block-Toeplitz coefficient matrix generated by the spatial discretization becomes ill-conditioned. As a result, we adopt an effective linearized iteration method for the nonlinear system, allowing it to be solved efficiently by the Krylov subspace solver with an appropriate circulant preconditioner, in which the fast Fourier transform is applied to speed up the computational cost at each iterative step. Finally, numerical experiments are presented to validate the theoretical findings and the efficiency of the fast algorithm.

Lie group symmetry analysis and invariant difference schemes of the two-dimensional shallow water equations in Lagrangian coordinates

The two-dimensional shallow water equations in Lagrangian coordinates are considered. Lie group classification for the class of the elliptic paraboloid bottom topography is performed.The transformations mapping the two-dimensional shallow water equations with a plane or rotation symmetric bottom into the gas dynamics equations of a polytropic gas with polytropic exponent γ=2 are represented. The group foliation of the two-dimensional shallow water equations in Lagrangian coordinates is discussed.New invariant conservative finite-difference schemes for the equations and their one-dimensional reductions are constructed. The schemes are derived either by extending the known one-dimensional schemes or by direct algebraic construction based on some assumptions on the form of the energy conservation law. Among the proposed schemes there are schemes possessing conservation laws of mass and energy.

Optimal point-wise error estimates of two conservative finite difference schemes for the coupled Gross–Pitaevskii equations with angular momentum rotation terms

This paper is concerned with numerical analysis of two finite difference schemes for solving the coupled Gross–Pitaevskii equations with angular momentum rotation terms which model the dynamics of rotating two-component Bose–Einstein condensates. Both schemes are proved to be uniquely solvable and inherit the conservation laws of mass and energy in the discrete sense. Besides the energy method, the ‘cut-off’ function technique and inductive argument are introduced to establish the optimal error estimates of the numerical methods without any requirement on the grid-ratio or initial values. The convergence order of the numerical solutions to the exact solutions is proved to be of O(τ2+h2) with time-step τ and mesh size h. Several numerical results are reported to validate the error estimates and conservation laws.

A general conservative eighth-order compact finite difference scheme for the coupled Schrödinger-KdV equations

&lt;abstract&gt;&lt;p&gt;In this paper, we present a general conservative eighth-order compact finite difference scheme for solving the coupled Schrödinger-KdV equations numerically. The proposed scheme is decoupled and preserves several physical invariants in discrete sense. The matrices obtained in the eighth-order compact scheme are all circulant symmetric positive definite so that it can be used to solve other similar equations. Numerical experiments on model problems show the better performance of the scheme compared with other numerical schemes.&lt;/p&gt;&lt;/abstract&gt;

Conservative finite difference schemes for dynamical systems

The article presents the implementation of one of the approaches to the integration of dynamical systems, which preserves algebraic integrals in the original fdm for Sage system. This approach, which goes back to the paper by del Buono and Mastroserio, makes it possible, based on any two explicit difference schemes, including any two explicit Runge-Kutta schemes, to construct a new numerical algorithm for integrating a dynamical system that preserves the given integral. This approach has been implemented and tested in the original fdm for Sage system. Details and implementation difficulties are discussed. For testing, two Runge-Kutta schemes were taken having the same order, but different Butcher tables, which does not complicate the method due to paralleling. Two examples are considered - a linear oscillator and a Jacobi oscillator with two quadratic integrals. The second example shows that the preservation of one integral of motion does not lead to the conservation of the other. Moreover, this method allows us to propose a practical application of the well-known ambiguity in the definition of Butcher tables.

Conservative finite difference methods for the Boussinesq paradigm equation

In this paper, we concentrate on developing and analyzing two types of finite difference schemes with energy conservation properties for solving the Boussinesq paradigm equation. Firstly, by introducing the auxiliary potential function ∂u∂t=Δv and lim|x|→∞v=0, the BPE is reformulated as an equivalent system of coupled equations. Then a class of efficient difference schemes are proposed for solving the resulting system, where one proposed scheme is a two-level nonlinear difference scheme and the other is a three-level linearized difference scheme using the invariant energy quadratization technique. Subsequently, we present the theoretical analysis of the proposed energy conservative finite difference schemes, which involves the discrete energy conservation properties, unique solvability and optimal error estimates. By using the discrete energy method, it is proven that the proposed schemes can achieve the optimal convergence rates of O(Δt2+hx2+hy2) in discrete L2-, H1- and L∞-norms. At last, numerical experiments illustrate the physical behaviors and efficiency of the proposed schemes.

Numerical solutions of generalized Rosenau–KDV–RLW equation by using Haar wavelet collocation approach coupled with nonstandard finite difference scheme and quasilinearization

AbstractIn this article, we analyze and propose to compute the numerical solutions of a generalized Rosenau–KDV–RLW (Rosenau‐Korteweg De Vries‐Regularized Long Wave) equation based on the Haar wavelet (HW) collocation approach coupled with nonstandard finite difference (NSFD) scheme and quasilinearization. In the process of the numerical solution, the NSFD scheme is applied to discretize the first‐order time derivative, Haar wavelets are applied on spatial derivatives and the non‐linear term is taken care by quasilinearization technique. To discuss the efficiency of the method we compute error and error. We also use discrete mass and energy conservation to check the accuracy of the proposed methodology. The computed results have been compared with the existing methods, for example, three‐level average implicit finite difference technique, B‐spline collocation, three‐level linear conservative implicit finite difference scheme and conservative fourth‐order stable finite difference scheme.

ON SOME NUMERICAL MODEL TO SOLVING DYNAMICAL EQUATIONS OF NONRELATIVISTIC AND RELATIVISTIC BACKWARD-WAVE TUBE

It is developed an effective computational approach to solution the master corresponding system of differential equations, which describe the nonlinear stationary and non-stationary electromagnetic processes in the nonrelativistic and relativistic backward-wave tubes (carcinotrons) with maximal accounting for the different physical factors such as the relativistic effects, effects of dissipation, the presence of a space charge, wave reflections at the ends of the slowing system, stochastic factors by means including the special elements in a whole system etc as well as the detailed investigation of characteristics (dynamical and topological invariants) of dynamics of a carcinotron in automodulation and chaotic regimes with construction the corresponding bi-furcation diagrams. Below in order to further solve the master system of dynamical equations for carcinotron it is presented in brief the realizing numerical scheme, based on the use of the conservative finite-difference schemes of the "predictor-corrector" type and the sweep algorithm.

Design of a Herringbone-Grooved Bearing for Application in an Electrically Driven Air Compressor

A turbo compressor was investigated to ensure the operational reliability of the charging of fuel cell systems. This study investigated air-lubricated herringbone bearings to support the high-speed rotating shaft. For reliable operation of the rotor bearing system, stable operation in the whole speed range (up to 120 krpm), as well as low lift-off speed, is an important issue. Some publications containing guidelines for an optimized design in terms of stability and lift-off behavior date back to the 1970s, with some simplifying assumptions (such as narrow groove theory and small eccentricity analysis). Many publications have addressed the calculations, as well as the optimization of herringbone-grooved bearings; however, general design guidelines are still missing in the view of the authors. Although the investigations related to bearings for the support of a lightweight rotor for a special compressor of a fuel cell unit, this study could also indicate favorable bearing designs for other high-speed applications. Here, the compressible Reynolds equation was solved in the whole solution domain using a conservative finite difference scheme, and the corresponding bearing characteristics were determined. In a perturbation analysis, the linearized dynamic coefficients of the herringbone bearing are calculated. To compare the suitability and performance of the various herringbone-grooved bearing designs, especially at high speed, the simple model of a Jeffcott rotor airborne with two identical herringbone-grooved journal bearings (HGJBs) was used. The geometrical parameters of the HGJBs were varied, and their effects on bearing characteristics and stability were evaluated. Recommendations concerning favorable geometrical bearing parameters for a sufficiently high stability threshold speed and reasonable low lift-off speed were the result of the parameter study.

Математическое моделирование динамики тепловых ударных волн в нелинейных локально-неравновесных средах

We performed a mathematical simulation of heat transfer in a local non-equilibrium medium whose transfer characteristics are functions of the temperature distribution. A homogeneous polynomial of arbitrary degree represents nonlinearities in thermal conductivity and thermal diffusivity. The mathematical model consists of a hyperbolic nonlinear heat transfer wave equation, initial conditions and nonlinear boundary conditions of the second and first kind. To solve this problem, we used a conservative homogeneous finite-difference scheme along the upper time grid line (implicitly). We then used the tridiagonal matrix algorithm of the second order in the spatial variable and of the first order in time to solve the resulting system of linearised algebraic equations. A periodic series of rectangular temperature or heat flux pulses form the boundary conditions of the first and second kind. Computation results reveal ultimate propagation rates of temperature and heat fronts featuring pronounced first-kind discontinuities with attenuating magnitudes. As the process unfolds, the initial pulses heat the region between the boundary and the heat wave front, while the subsequent pulses traverse this region at a higher velocity due to thermal diffusivity being a function of temperature, their fronts "catching up" with the previous fronts, increasing the discontinuity magnitude at the initial pulse front, that is, forming a thermal shock wave front similar to that of a shock wave in gas dynamics. We obtained such thermal shock waves for boundary conditions of both the first and the second kind. We also analysed kinematic and dynamic characteristics of thermal waves

Advancing Three-Dimensional Coupled Water Quality Model of Marine Ranches: Model Development, Global Sensitivity Analysis, and Optimization Based on Observation System

Marine ranching is a stock enhancement project that has been an important part of aquaculture in China. Due to the lack of scientific management, disasters have occurred, resulting in millions of economic losses. Based on the observation system of marine ranches, a three-dimensional water quality model of marine ranches was developed to study the spatio-temporal variation of key ecological indicators, including the concentrations of chlorophyll-a, zooplankton, detritus, nutrients, and dissolved oxygen (DO). The model is coupled offline with the residual current, temperature, and salinity simulated by a regional oceanic modeling system (ROMS). The conservative characteristic finite difference (C-CFD) scheme is introduced to solve the equations, which guarantees model stability and mass conservation and allows for a larger time step compared to traditional difference schemes. In state-of-the-art water quality models, the biogeochemical processes are parameterized. Due to the complexity of the water quality model, a combination of global sensitivity analysis (GSA) and the adjoint method is introduced as the methodology to optimize the model parameters. Morris’ sampling method is implemented as the GSA method to find out the key factors of the water quality model. The optimization of sensitive parameters with the adjoint method significantly improves the model precision, while the other parameters can be set as empirical values. The results indicate that the combination of GSA and the adjoint method is efficient in the parameter optimization of the water quality model. The model is applicable in marine ranches.