Cauchy problems with periodic controls.
Central schemes and contact discontinuities
Central schemes and contact discontinuities
We introduce a family of new second-order Godunov-type central schemes for one-dimensional systems of conservation laws. They are a less dissipative generalization of the central-upwind schemes, proposed in [A. Kurganov et al., submitted to SIAM J. Sci. Comput.], whose construction is based on the maximal one-sided local speeds of propagation. We also present a recipe, which helps to improve the resolution of contact waves. This is achieved by using the partial characteristic decomposition,...
Central WENO schemes for hyperbolic systems of conservation laws
Central WENO schemes for hyperbolic systems of conservation laws
We present a family of high-order, essentially non-oscillatory, central schemes for approximating solutions of hyperbolic systems of conservation laws. These schemes are based on a new centered version of the Weighed Essentially Non-Oscillatory (WENO) reconstruction of point-values from cell-averages, which is then followed by an accurate approximation of the fluxes via a natural continuous extension of Runge-Kutta solvers. We explicitly construct the third and fourth-order scheme and demonstrate...
Central-upwind schemes for the Saint-Venant system
We present one- and two-dimensional central-upwind schemes for approximating solutions of the Saint-Venant system with source terms due to bottom topography. The Saint-Venant system has steady-state solutions in which nonzero flux gradients are exactly balanced by the source terms. It is a challenging problem to preserve this delicate balance with numerical schemes. Small perturbations of these states are also very difficult to compute. Our approach is based on extending semi-discrete central schemes...
Central-Upwind Schemes for the Saint-Venant System
We present one- and two-dimensional central-upwind schemes for approximating solutions of the Saint-Venant system with source terms due to bottom topography. The Saint-Venant system has steady-state solutions in which nonzero flux gradients are exactly balanced by the source terms. It is a challenging problem to preserve this delicate balance with numerical schemes. Small perturbations of these states are also very difficult to compute. Our approach is based on extending semi-discrete central...
Comparative Study of a Solid Film Dewetting in an Attractive Substrate Potentials with the Exponential and the Algebraic Decay
We compare dewetting characteristics of a thin nonwetting solid film in the absence of stress, for two models of a wetting potential: the exponential and the algebraic. The exponential model is a one-parameter (r) model, and the algebraic model is a two-parameter (r, m) model, where r is the ratio of the characteristic wetting length to the height of the unperturbed film, and m is the exponent of h (film height) in a smooth function that interpolates the system's surface energy above and below...
Comparison of explicit and implicit difference methods for quasilinear functional differential equations
We give a theorem on error estimates of approximate solutions for explicit and implicit difference functional equations with unknown functions of several variables. We apply this general result to investigate the stability of difference methods for quasilinear functional differential equations with initial boundary condition of Dirichlet type. We consider first order partial functional differential equations and parabolic functional differential problems. We compare the properties of explicit...
Comparison of explicit and implicit difference schemes for parabolic functional differential equations
Initial-boundary value problems of Dirichlet type for parabolic functional differential equations are considered. Explicit difference schemes of Euler type and implicit difference methods are investigated. The following theoretical aspects of the methods are presented. Sufficient conditions for the convergence of approximate solutions are given and comparisons of the methods are presented. It is proved that the assumptions on the regularity of the given functions are the same for both methods. It...
Compressible two-phase flows by central and upwind schemes
This paper concerns numerical methods for two-phase flows. The governing equations are the compressible 2-velocity, 2-pressure flow model. Pressure and velocity relaxation are included as source terms. Results obtained by a Godunov-type central scheme and a Roe-type upwind scheme are presented. Issues of preservation of pressure equilibrium, and positivity of the partial densities are addressed.
Compressible two-phase flows by central and upwind schemes
This paper concerns numerical methods for two-phase flows. The governing equations are the compressible 2-velocity, 2-pressure flow model. Pressure and velocity relaxation are included as source terms. Results obtained by a Godunov-type central scheme and a Roe-type upwind scheme are presented. Issues of preservation of pressure equilibrium, and positivity of the partial densities are addressed.
Computation of bifurcated branches in a free boundary problem arising in combustion theory
Computation of bifurcated branches in a free boundary problem arising in combustion theory
We consider a parabolic 2D Free Boundary Problem, with jump conditions at the interface. Its planar travelling-wave solutions are orbitally stable provided the bifurcation parameter does not exceed a critical value . The latter is the limit of a decreasing sequence of bifurcation points. The paper deals with the study of the 2D bifurcated branches from the planar branch, for small k. Our technique is based on the elimination of the unknown front, turning the problem into a fully nonlinear...
Conditions aux limites pour un système strictement hyperbolique fournies, par le schéma de Godunov
Conservative numerical methods for a two-temperature resistive MHD model with self-generated magnetic field term
We propose numerical methods on Cartesian meshes for solving the 2-D axisymmetric two-temperature resistivive magnetohydrodynamics equations with self-generated magnetic field and Braginskii’s [1] closures. These rely on a splitting of the complete system in several subsystems according to the nature of the underlying mathematical operator. The hyperbolic part is solved using conservative high-order dimensionally split Lagrange-remap schemes whereas...
Constraint preserving schemes using potential-based fluxes. III. Genuinely multi-dimensional schemes for MHD equations∗∗∗
We design efficient numerical schemes for approximating the MHD equations in multi-dimensions. Numerical approximations must be able to deal with the complex wave structure of the MHD equations and the divergence constraint. We propose schemes based on the genuinely multi-dimensional (GMD) framework of [S. Mishra and E. Tadmor, Commun. Comput. Phys. 9 (2010) 688–710; S. Mishra and E. Tadmor, SIAM J. Numer. Anal. 49 (2011) 1023–1045]. The schemes are formulated in terms of vertex-centered potentials....
Constraint preserving schemes using potential-based fluxes. III. Genuinely multi-dimensional schemes for MHD equations∗∗∗
We design efficient numerical schemes for approximating the MHD equations in multi-dimensions. Numerical approximations must be able to deal with the complex wave structure of the MHD equations and the divergence constraint. We propose schemes based on the genuinely multi-dimensional (GMD) framework of [S. Mishra and E. Tadmor, Commun. Comput. Phys. 9 (2010) 688–710; S. Mishra and E. Tadmor, SIAM J. Numer. Anal. 49 (2011) 1023–1045]. The schemes are formulated in terms of vertex-centered potentials....
Constructing Orthogonal Curvilinear Meshes by Solving Initial Value Problems.
Contaminant transport with adsorption in dual-well flow
Numerical approximation schemes are discussed for the solution of contaminant transport with adsorption in dual-well flow. The method is based on time stepping and operator splitting for the transport with adsorption and diffusion. The nonlinear transport is solved by Godunov’s method. The nonlinear diffusion is solved by a finite volume method and by Newton’s type of linearization. The efficiency of the method is discussed.