Search

Chemical reactions occur everywhere in both natural and artificial systems. Some of the reactions occur during the flow of a fluid (such a process is referred to as a reactive flow). Given the hazardous nature of some reactive flows, computer simulations (rather than physical experiments) are necessary for ascertaining or enhancing our understanding of such systems. The process of simulation involves mathematical and numerical modeling of the reactive flows. Mathematical models for reactive flow problems are complicated partial differential equations that often lack exact solutions, thus, numerical solutions are employed. Numerical methods must preserve almost all the relevant properties of the problem for accuracy reasons. Dispersion relations are important properties of wave propagation problems and numerical methods that satisfy them are called dispersion preserving methods. Furthermore, stiff transport models are wave propagation problems that cannot be solved efficiently with explicit methods. However, fully implicit methods are computationally expensive. A combination of implicit and explicit methods called implicit–explicit methods is usually employed to efficiently resolve stiffness. An example of problems of interest in this regard are the advection–diffusion–reaction (ADR) models. In this discussion, spectral analysis is performed on two implicit–explicit methods to ascertain their dispersion preserving abilities in order to determine their suitability for simulating general stiff reactive flow problems. The analysis shows that both implicit–explicit methods are dispersion preserving, however, one particular method is more suitable for general wave propagation problems. [ABSTRACT FROM AUTHOR]

In this paper, we consider simultaneous reconstruction of diffusion coefficient and initial state for a one-dimensional heat equation through boundary control and measurement. The boundary measurement is proposed to make the system approximately observable, and both the coefficient and initial state are shown to be identifiable by this measurement under a boundary switch on/off control. By a Dirichlet series representation for the observation, we can transform the problem into an inverse process of reconstruction of the spectrum and coefficients for the Dirichlet series in terms of observation. This happens to be the reconstruction of spectral data for an exponential sequence with measurement error, and it enables us to develop an algorithm based on the matrix pencil method in signal analysis. A theoretical error analysis for the algorithm concerning the coefficient reconstruction is carried out for the proposed method. The numerical simulations are presented to verify the proposed algorithm. [ABSTRACT FROM AUTHOR]

In this paper, numerical boundary stabilisation of a non-uniform hyperbolic system of balance laws is studied. For the numerical discretisation of the balance laws, a first order explicit upwind scheme is used for the spatial discretisation; and for the temporal discretisation a splitting technique is employed. A discrete L 2 -Lyapunov function is employed to investigate conditions for the stability of the system. After constructing discrete numerical Lyapunov functionals, we prove an asymptotic exponential stability result for a class of non-uniform linear hyperbolic systems of balance laws. Convergence of the solution to its equilibrium is proved. Further application of the approach to practical problems through concrete examples is presented together with suggestions for numerical implementation. The numerical computations also demonstrate the stability of the numerical scheme with parameters chosen to satisfy the stability requirements. [ABSTRACT FROM AUTHOR]

In this paper, simultaneous identification of damping or anti-damping coefficient and initial value for some Riesz spectral systems is considered. An identification algorithm is proposed based on the fact that the output of the systemhappens to be decomposed into a product of an exponential function and a periodic function. The former contains information of the damping coefficient, while the latter does not. The convergence and error analysis are also developed. Three examples, namely an anti-stablewave equation with boundary anti-damping, the Schrödinger equation with internal antidamping and two connected strings with middle joint anti-damping, are investigated and demonstrated by numerical simulations to show the effectiveness of the proposed algorithm. [ABSTRACT FROM AUTHOR]

This paper deals with an inverse problem of determining the diffusion coefficient, spacewise dependent source term, and the initial value simultaneously for a one-dimensional heat equation based on the boundary control, boundary measurement, and temperature distribution at a given single instant in time. By a Dirichlet series representation for the boundary observation, the identification of the diffusion coefficient and initial value can be transformed into a spectral estimation problem of an exponential series with measurement error, which is solved by the matrix pencil method. For the identification of the source term, a finite difference approximation method in conjunction with the truncated singular value decomposition is adopted, where the regularization parameter is determined by the generalized cross-validation criterion. Numerical simulations are performed to verify the result of the proposed algorithm. Copyright © 2016 John Wiley & Sons, Ltd. [ABSTRACT FROM AUTHOR]

In this paper higher order implicit Runge-Kutta schemes are applied to fluid-structure interaction (FSI) simulations. A staggered approach with a structural predictor is applied to an FSI problem. The equations governing the dynamics of the structure are integrated in time by the Explicit Single Diagonal Implicit Runge-Kutta (ESDIRK) schemes and the arbitrary high order finite volume scheme is taken as the fluid solver. The performance of the ESDIRK scheme of order of convergence three to five is tested. Comparative studies with other time integration schemes which have been successfully applied to FSI problems are undertaken. Comparisons to test the performance of the scheme are carried out. The staggered approach is applied to couple the structure and the compressible fluid, therefore the added mass effect is not considered. However the influence of the structural predictors is analyzed through energy conservation. [ABSTRACT FROM AUTHOR]

Suitable numerical discretizations for boundary control problems of systems of nonlinear hyperbolic partial differential equations are presented. Using a discrete Lyapunov function, exponential decay of the discrete solutions of a system of hyperbolic equations for a family of first-order finite volume schemes is proved. The decay rates are explicitly stated. The theoretical results are accompanied by computational results. [ABSTRACT FROM AUTHOR]

Dynamics of multiphase flows through networks are considered. The dynamics of flow through the connected arcs are governed by an isothermal no-slip drift-flux model. Such problems arise in the context of multicomponent flows or in gas transport in pipe networks in which a phase change takes place due to geometrical or physical forces. Coupling conditions for the vertices (joints) in a network have been proposed. We present conditions at and introduce a mathematical representation of the vertex flow for the no-slip drift-flux case of multiphase flows. Mathematical analysis of coupling conditions at the vertices as well as numerical simulations and comparative studies with theoretical predictions are undertaken. [ABSTRACT FROM AUTHOR]