Search

The numerical shock instability problem significantly limits the application of the finite volume method in hypersonic flow simulations, especially for high-order schemes, which will encounter more severe shock instability problems. In the previous work, the MSAT (matrix stability analysis tool) was introduced to quantitatively analyze the numerical shock instability problem. However, the utilization of the three-point reconstruction stencil in the original version of MSAT presents challenges in conducting investigations into high-order schemes. The current work provides an extension of MSAT by extending the reconstruction stencil to five points, a challenging but necessary endeavor for high-order schemes. This extension enables MSAT to provide analyses for schemes with high-order reconstruction methods. Furthermore, the updated MSAT can be utilized to investigate the impact of characteristic variable reconstruction on the shock instability problem. Results demonstrate the reliability of the updated MSAT and reveal novel characteristics in the shock instability of high-order schemes. Consequently, the updated MSAT will contribute to the investigation of shock instabilities for high-order schemes, ultimately aiding in mitigating this challenging problem.

The simulation of supersonic or hypersonic flows often suffers from numerical shock instabilities if the flow field contains strong shocks, limiting the further application of shock-capturing schemes. In this paper, we develop the unified matrix stability analysis method for schemes with three-point stencils and present MSAT (matrix stability analysis tool), an open-source tool to quantitatively analyze the shock instability problem. Based on the finite-volume approach on the structured grid, MSAT can be employed to investigate the mechanism of the shock instability problem, evaluate the robustness of numerical schemes, and then help to develop robust schemes. Also, MSAT has the ability to analyze the practical simulation of supersonic or hypersonic flows, evaluate whether it will suffer from shock instabilities, and then assist in selecting appropriate numerical schemes accordingly. As a result, MSAT is a helpful tool that can investigate the shock instability problem and help to cure it.

In present paper, a new technique is developed named Uniform Algebraic Trigonometric (UAT) tension B-spline based DQM, for the better numerical approximation of coupled 1D and coupled 2D Burgers’ equation. The formula of UAT tension B-spline of order 4 is developed in this paper, by using above-mentioned basis function. Strong stability preserving SSP-RK 43 method is implemented to solve obtained ordinary differential equations. To check the effectiveness and accuracy of the proposed scheme, L2 and L∞error norms are used. Four numerical experiments are discussed to put light upon applicability of scheme. By means of tables and figures, results obtained are given. Obtained results are in a good agreement with existing ones. Stability of proposed scheme is also checked by implementing the Matrix stability analysis method. This novel scheme will surely open some new dimensions in the field of numerical approximation, for other researchers in their future work.

This paper proposes a new scheme termed as modified extended cubic B-Spline differential quadrature (mECDQ) method for time dependent partial differential equations. Specially, the numerical computation of the Burgers’ equation is obtained using mECDQ method. First of all the modified extended cubic B-splines are used as a set of basis function in DQ to evaluate the weighting coefficients. The mECDQ method converts the initial boundary value system of Burgers’ equation into a initial value system of ordinary differential equations (ODEs). The resulting system is solved by using an optimal five stage four order strong stability preserving Runge–Kutta method (SSP-RK54). The accuracy and efficiency of the method is illustrated by considering five test problems. The proposed results are compared with the exact solutions in terms of L2 and L∞ error norms and the existing results. The mECDQ scheme produces better results than the results due to almost all the existing schemes. The stability analysis of the scheme is also carried out using the matrix stability analysis method for various grid points. This shows that mECDQ scheme is conditionally stable.

The positive-definite stability analysis (PDSA) is presented as a technique complementary to the companion-matrix stability analysis (CMSA). The PDSA is used to analyze the stability of marching-on-in-time (MOT) schemes. The heart of the PDSA is formed by the analysis on particular linear combinations of interaction matrices from an MOT scheme, which are assumed to be real-valued. If these are all positive definite, then the PDSA guarantees the stability of the scheme. The PDSA can be of a lower complexity than the full CMSA. The construction of the PDSA is shown and applied to two numerical examples. [ABSTRACT FROM AUTHOR]

The aim of this study is to develop compact difference method to approximate parabolic PDEs of fourth order equipped with Dirichlet and Neumann boundary conditions involving half-step points. The proposed method converges quaternary and quadratically in space and time, respectively. The imbedding technique has been applied to approximate derivative terms of lower order by means of the governing differential equation to deduce the high-order method. The primary utility of this new discretization is that it can be straightaway applied to problems with singularities without necessitating fictitious nodes or special approach which has consequently lowered computational complicacy. We have examined linear stability of the proposed three-level implicit difference scheme using matrix stability analysis. In addition, we also obtained the solution of the first-order spatial derivative which is of significance in several physical problems. The efficacy of the proposed approximation is confirmed through numerical tests performed on a collection of physically relevant problems comprising the Euler Bernoulli beam equation and the highly nonlinear good Boussinesq equation. Numerical experiments evidently exhibit that the method provides more accurate results in contrast with the existing numerical techniques. The present method is able to simulate well the complex and intriguing long time dynamics of the good Boussinesq equation. [ABSTRACT FROM AUTHOR]

The contrast current density volume integral equation, discretized with piecewise constant spatial basis and test functions and Dirac-delta temporal test functions and the piecewise polynomial temporal basis functions, results in a causal implicit marching-on-intime scheme that we refer to as the marching-on-in-time contrast current density volume integral equation (MOT-JVIE). The companion matrix stability analysis of the MOT-JVIE solver shows that for a fixed spatial and temporal step size, the stability is independent of the scatterer's dielectric contrast for quadratic spline temporal basis functions. Whereas, Lagrange and cubic spline exhibit instabilities at higher contrast. We relate this stability performance to the expansion and testing procedure in time. We further illustrate the capabilities of the MOT-JVIE based on quadratic spline temporal basis functions by: comparing the MOT-JVIE solution to time-domain results from literature and frequency-domain results from a commercial combined field integral equation solver. Finally, we present a long time sequence for a high-contrast scatterer discretized with 24,000 spatial unknowns. [ABSTRACT FROM AUTHOR]

Modern shock-capturing schemes often suffer from numerical shock instabilities when simulating strong shocks, limiting their application in supersonic or hypersonic flow simulations. In the current study, we devote our efforts to investigating the shock instability problem for second-order schemes, which has not gotten enough attention in previous research but is crucial to address. To this end, we develop the matrix stability analysis method for the finite-volume Monotone Upstream-centered Schemes for Conservation Laws (MUSCL) approach, taking into account the influence of reconstruction. With the help of this newly developed method, the shock instability of second-order schemes is investigated quantitatively and efficiently. The results demonstrate that when second-order schemes are employed, whether shock instabilities will occur is closely related to the property of Riemann solvers, just like the first-order case. However, enhancing spatial accuracy still impacts the shock instability problem, and the impact can be categorized into two types depending on the dissipation of Riemann solvers. Furthermore, the research emphasizes the impact of the numerical shock structure, highlighting both its role as the source of instability and the influence of its state on the occurrence of instability. One of the most significant contributions of this study is the confirmation of the multidimensional coupled nature of shock instability. Both one-dimensional and multidimensional instabilities are proven to influence the instability problem, and they have different properties. Moreover, this paper reveals that increasing the aspect ratio and distortion angle of the computational grid can help mitigate shock instabilities. The current work provides an effective tool for quantitatively investigating the shock instabilities for second-order schemes, revealing the inherent mechanism and thereby contributing to the elimination of shock instability. [ABSTRACT FROM AUTHOR]

This study presents a numerical solution of Fisher's equation. For time integration, Crank–Nicolson and fourth-order one-step implicit schemes are used and for space discretization, quintic B-spline Collocation and quintic B-spline Galerkin methods are employed. The truncation error is analyzed and the stability of the suggested methods is discussed matrix stability analysis. Three examples are studied to compare the present results with existing numerical results by computing error norm L ∞ and the order of convergence. The obtained results show that the proposed methods are satisfactorily efficient in terms of accuracy. [ABSTRACT FROM AUTHOR]

In the present paper, a Barycentric Lagrange interpolation polynomial-based Differential quadrature method is applied to get the numerical solution of non-linear coupled 2D But-gers' equation, For this process, first of all, the coupled Burgers' equation is dealt with generalized formula, By using the obtained generalized formula, required weighting coefficients in the differential quadrature method are easily calculated, Afterward, the reduced set of ordinary differential equations is solved by using the SSP-RK43 scheme, To test the efficiency of the newly applied method, the L., error norm has been computed, To test the stability of the proposed method, the Matrix stability analysis method is implemented, which shows that the developed method is unconditionally stable. [ABSTRACT FROM AUTHOR]

In December 2019, the newest infection, named novel coronavirus COVID-19, appeared, and it became a challenge to the world. In this paper, we investigated a mathematical model that explains the transmission of COVID- 19 with the help of ordinary differential equations. Here, we divided the entire model into eight different compartments, namely SEI1I2I3QHR model based on three different age groups. Equilibrium analysis with positivity and boundedness are also investigated. The reproduction number is analyzed using Vandedressche's next-generation matrix. Stability analysis for disease-free equilibrium (DFE) is studied. We have taken some basic parameters with their primary conditions from the Ethiopian Ministry of Health Institute from 29 February 2021 to 07 June 2021. [ABSTRACT FROM AUTHOR]

Present paper deals with the numerical solution of coupled 1D Burgers' equation by implementing the Non-Uniform Algebraic Hyperbolic (NUAH) B-spline Differential Quadrature Method. In the present paper, the spatial variable discretization is done using NUAH B-spline, and the obtained system of ODE is dealt with using the SSP-RK43 scheme. To get the improvised results, the concept of modified cubic NUAH B-spline is incorporated. To test the effectiveness and accuracy of the proposed scheme, numerical examples are discussed. Stability analysis of the proposed scheme is investigated by the matrix stability analysis method. The present regime is worthwhile to deal with some complex natured PDEs, where finding the analytical solution is cumbersome. [ABSTRACT FROM AUTHOR]

We propose a differential quadrature method (DQM) based on cubic hyperbolic B-spline basis functions for computing 3D wave equations. This method converts the problem into a system of ODEs. We use an optimum five-stage and order four SSP Runge-Kutta (SSPRK-(5,4)) scheme to solve the obtained system of ODEs. The matrix stability analysis is also investigated. The accuracy and efficiency of the proposed method are demonstrated via three numerical examples. It has been found that the proposed method gives more accurate results than the existing methods. The main purpose of this work is to present an accurate, economically easy-to-implement, and stable technique for solving hyperbolic partial differential equations. [ABSTRACT FROM AUTHOR]

A simple HLLE-type scheme is proposed for all Mach number flows. In the proposed scheme, no extra wave structure is added in the HLLE scheme to resolve the shear wave while the contact wave is resolved by adding a wave structure similar to the HLLEM scheme. The resolution of the shear layers and the flow features at low Mach numbers is achieved by a velocity reconstruction method based on the face normal Mach number. Robustness against the numerical instabilities is achieved by scaling the velocity reconstruction function in the vicinity of shock with a multi-dimensional pressure sensor. The ability of the proposed scheme to resolve low Mach flow features is demonstrated through asymptotic analysis while the stability of the proposed scheme for strong shock is demonstrated through linear perturbation and matrix stability analyses. A set of numerical test cases are solved to show that the proposed scheme is free from numerical shock instability problems at high speeds and is capable of resolving the flow features at very low Mach numbers. • A new velocity reconstruction method for an all Mach number HLLE scheme. • Velocity reconstruction based on face normal Mach number and pressure function. • Velocity reconstruction resolves shear layers and low Mach flow features while presrving shock stability. • HLLEM type method for resolving contact wave while preserving shock stability. • Analytical and numerical demonstration of all Mach capability of the proposed scheme. [ABSTRACT FROM AUTHOR]

In this work, a numerical regime based upon the Non-Uniform Algebraic Hyperbolic Bspline is developed with the aid of the differential quadrature method named NUAH B-spline DQM to fetch the numerical approximation of non-linear coupled ZD Burgers equation. In the first step, the spatial discretization is mentioned using the NUAH B-spline of order 4: for improvised results, modified formulae are implemented, In the second step, the reduced system of ordinary differential equations is dealt using SSP-RK43 regime. The proposed regime is analyzed for stability also. Finally, the effectiveness and adaptability of the proposed regime are confirmed with the aid of numerical experiments, and fetched results are compared with the existing outcomes in the literature. It is observed that the present scheme is unconditionally stable by implementing the concept of the matrix stability analysis method. It is found that the present regime is easy to implement and produces accurate results. [ABSTRACT FROM AUTHOR]

The aim of present study is to develop a numerical scheme by using the notion of Barycentric Lagrange interpolation based Differential quadrature method to solve the coupled 1D Burgers' equation. This method reduced the mentioned partial differential equation into the set of ordinary differential equations, which can be dealt by the SSP-RK43 scheme. The proposed method has been implemented upon the different numerical examples in order to test the accuracy and effectiveness of the proposed method. It is observed that obtained results are in good compatibility with the exact solution and are better than the previous results. The stability of the proposed method is also discussed by using the matrix stability analysis method, which represents that the proposed method is unconditionally stable. [ABSTRACT FROM AUTHOR]

In this, research work. we present a novel technique named Hyperbolic B-spline based Differential quadrature method to get the numerical approximations of the 1 D and 2D coupled nonlinear Schrodinger equations. The coupled non-linear Schrodinger equation represents the range of physical phenomena with the aid of a model equation for the optical fiber with the linear birefringence. In present paper. Differential quadrature method is applied over the mentioned PDEs and reduced set of ordinary differential equations is solved by SSP-RK43 scheme. Accuracy and effectiveness of the proposed scheme is checked by L∞ error norm and conservation of mass properly is also tested. In present paper. elastic property is also checked. whether the wave amplitudes have restored their shapes or not for the different used parameters and stability of the proposed scheme is checked via matrix stability analysis method. which indicates that the present method has the property of unconditionally stable. [ABSTRACT FROM AUTHOR]

Present paper deals with the numerical solution of coupled 2D non-linear Burgers equation by implementing the notion of modified cubic Hyperbolic B-spline based Differential quadrature method. Feasibility and accuracy of proposed scheme is checked by using L2 and L∞ error norms. Comparison of errors is shown for the mentioned examples. Reduced set of ordinary differential equations is tackled by using SSP-RK43 scheme. Stability of the proposed scheme is checked by using the Matrix stability analysis method. It is noticed that proposed numerical scheme produced the results with accuracy and it is quite easy to implement. [ABSTRACT FROM AUTHOR]

The paper devises a new robust carbuncle-free Roe Riemann solver for strong shock, different from hybrid method, entropy fix, Liou’s conjecture (J Comput Phys 160:623-648, 2000) and artificial viscosity (Rodionov in J Comput Phys 345:308-329, 2017). Roe scheme encounters carbuncle phenomenon, violates entropy condition and lacks positivity property. The remedy integrates shear viscosity into momentum flux to damp undesirable perturbation and prevents unstable vorticity mode from triggering shock instability. Compared to HLL scheme, shear viscosity is properly established through dimensional analysis and analogy method. The non-linear wave speeds are slightly modified by incorporating the neighboring cell information for positive conservation. The pressure-based sensing function is applied to preserve shear layer while keeping shock robustness. The resulting scheme is very easily implemented by converting the existing Roe code. The matrix stability analysis confirms that this approach is shock-stable and more robust than entropy fix. A series of numerical results demonstrate its potential features: positivity-preserving property, entropy-satisfying property, accurate boundary-layer resolution and high robustness against shock instability. In addition, it signifies that the reason to cause shock instability for Roe scheme may be not the pressure difference term in mass flux, but the inadequate shear viscosity. [ABSTRACT FROM AUTHOR]

This article concerns with the numerical investigations of the reaction–diffusion systems (RDSs) arising in the study of pattern formation in biological and chemical systems with the employment of the quartic-trigonometric B-spline functions. The computationally numerical scheme uses collocation method which is established by a relatively new B-splines for the spatial discretizations and, for time integration Crank–Nicolson technique is adapted. Therefore, solutions of the RDSs are assembled by the wholly discretized space-time scheme. A matrix stability analysis is performed for the numerical scheme after linearization process. Experimental cases include Brusselator model, Gray–Scott model, Schnakenberg model as well as a linear problem in one-dimensional domain. Numerical solutions are compared to the existing studies. Spatial pattern formation is demonstrated by present computational algorithm. [ABSTRACT FROM AUTHOR]

A multiple locally 1-D (MLOD) finite-difference time-domain (FDTD) method for inhomogeneous coupled transmission lines and stability analyses are presented. The method is aptly called the MLOD coupled line (CL)-FDTD method. Various split matrices are proposed, and the corresponding update equations are formulated and discussed. All the proposed split matrices yield implicit electric field update equations with tridiagonal or block tridiagonal matrices on the left-hand sides. For more efficiency, the block tridiagonal matrices for implicit electric field may be reformulated and replaced with tridiagonal matrices for implicit magnetic field. The stability analysis is first performed using the von Neumann method in the Fourier domain. It is shown that the von Neumann method alone may not be sufficient to ascertain stability for inhomogeneous media. To include media inhomogeneity, the two-media reduced-matrix stability analysis is proposed. It allows us to efficiently analyze the key stability characteristics in inhomogeneous media and is useful for quick detection of any potential instability that is not apparent via the von Neumann method. The stability characteristics with variation of media parameters are also investigated and discussed. The numerical results are provided to validate the stability and accuracy of the proposed MLOD CL-FDTD method with various split matrices. [ABSTRACT FROM AUTHOR]

This chapter is concerned with the problem of the functional observerbased sliding mode control (SMC) design for parametric uncertain discrete-time delayed stochastic systems including mismatched parameter uncertainty in the state matrix and in the delayed state matrix. Stability analysis of the sliding function is presented in a time delayed stochastic system with the linear matrix inequality (LMI) approach. Moreover, it is shown that the state trajectories can be driven onto the specified sliding surface despite the presence of state delay, unmatched parameter uncertainty, and stochastic noise in the system. The research is motivated by the fact that system states are not always accessible for the state feedback. Therefore, SMC is estimated using the functional observer technique. To mitigate the side effect of parametric uncertainty on the estimation error, a sufficient condition of stability is proposed based on Gershgorin's circle theorem. The claims made are validated through numerical simulations. [ABSTRACT FROM AUTHOR]

In this paper, an attempt is made to understand the transmission dynamics of Ebola virus disease (EVD) incorporating fuzziness in all biological parameters due to its natural variability. To characterize the transmission trajectories of Ebola outbreak, we propose and analyze two SEIR and SEIRHD type transmission models. Using triangular fuzzy numbers for the imprecise parameters, we first study the existence of the equilibria and their stability. Both of the model have two equilibria, namely the disease-free and endemic. Stability of the disease-free and endemic equilibria is related with basic reproduction number that has been calculated from next generation matrix. Stability analysis of the system shows that the disease free equilibrium is locally as well as globally asymptotically stable when the basic reproduction number is less than unity. Under some additional conditions, the model system becomes locally asymptotically stable at unique endemic equilibrium when basic reproduction number is greater than unity. Finally, we perform some numerical experiments to justify the theoretical estimate. [ABSTRACT FROM AUTHOR]

This study is concerned with the problem of the functional observer‐based sliding mode control (SMC) design for parametric uncertain discrete‐time delayed stochastic systems includes mismatched parameter uncertainty in the state matrix and in the delayed state matrix. Stability analysis of sliding function is presented in the time delayed stochastic system with a linear matrix inequality approach. Moreover, it is shown that the state trajectories can be driven onto the specified sliding surface despite the presence of state delay, unmatched parameter uncertainty and stochastic noise in the system. The research is motivated by the fact that the system states are not always accessible for the state feedback. Therefore, SMC is estimated using the functional observer technique. To mitigate the side effect of the parameter uncertainty on the estimation error, a sufficient condition of stability is proposed based on Gershgorin disc theorem. The claims made are validated through numerical simulations. [ABSTRACT FROM AUTHOR]

In this article, we discuss modified three level implicit difference methods of order two in time and four in space for the numerical solution of two‐ and three‐dimensional telegraphic equation with Robin boundary conditions. Ghost points are introduced to obtain fourth‐order approximations for boundary conditions. Matrix stability analysis is carried out to prove stability of the method for telegraphic equations in two and three dimensions with Neumann boundary conditions. Numerical experiments are carried out and the results are found to be better when compared with the results obtained by other existing methods. [ABSTRACT FROM AUTHOR]

Purpose In this study, a second-order standard wave equation extended to a two-dimensional viscous wave equation with timely differentiated advection-diffusion terms has been solved by differential quadrature methods (DQM) using a modification of cubic B-spline functions. Two numerical schemes are proposed and compared to achieve numerical approximations for the solutions of nonlinear viscous wave equations.Design/methodology/approach Two schemes are adopted to reduce the given system into two systems of nonlinear first-order partial differential equations (PDE). For each scheme, modified cubic B-spline (MCB)-DQM is used for calculating the spatial variables and their derivatives that reduces the system of PDEs into a system of nonlinear ODEs. The solutions of these systems of ODEs are determined by SSP-RK43 scheme. The CPU time is also calculated and compared. Matrix stability analysis has been performed for each scheme and both are found to be unconditionally stable. The results of both schemes have been extensively discussed and compared. The accuracy and reliability of the methods have been successfully tested on several examples.Findings A comparative study has been carried out for two different schemes. Results from both schemes are also compared with analytical solutions and the results available in literature. Experiments show that MCB-DQM with Scheme II yield more accurate and reliable results in solving viscous wave equations. But Scheme I is comparatively less expensive in terms of CPU time. For MCB-DQM, less depository requirements lead to less aggregation of approximation errors which in turn enhances the correctness and readiness of the numerical techniques. Approximate solutions to the two-dimensional nonlinear viscous wave equation have been found without linearizing the equation. Ease of implementation and low computation cost are the strengths of the method.Originality/value For the first time, a comparative study has been carried out for the solution of nonlinear viscous wave equation. Comparisons are done in terms of accuracy and CPU time. It is concluded that Scheme II is more suitable. [ABSTRACT FROM AUTHOR]

The Central Inspection for Ecological and Environmental Protection (CIEEP) and Environmental Vertical Management Reform (EVMR) policies have reshaped the relationships among stakeholders and the power structure in local environmental governance systems. To reveal the current local environmental collaborative governance mechanism, this study placed CIEEP and EVMR policies within the same research framework and considered multiple governance entities, including the central government, higher-level environmental protection departments, local environmental protection departments, enterprises, and the public, employing evolutionary game theory and numerical simulation techniques to explore the boundary conditions, key influencing factors, influence mechanisms, and realization paths of local environmental collaborative governance from a dynamic process perspective. The results show that: 1) The cost of enterprise environmental behavior, intensity of reward and punishment for enterprise environmental behavior, government subsidies, supervision cost of local environmental protection departments, and intensity of reward and punishment for local environmental protection departments all have important effects on local environmental governance; 2) Compared with the cost of enterprise environmental behavior, enterprise environmental behavior is more sensitive to changes in supervision cost of local environmental protection departments, and the intensity of reward and punishment for enterprise environmental behavior has the same major impact on enterprise environmental behavior as that of local environmental protection departments; 3) The influence of CIEEP policy on enterprise environmental behavior is unstable, showing signs of strategic response. Therefore, local environmental governance systems should be further optimized by encouraging local environmental protection departments to actively supervise and designing differentiated policy tools. This study provides a reference for solving local environmental governance problems under current environmental regulation policies. [ABSTRACT FROM AUTHOR]

A local differential quadrature method based on Fourier series expansion numerically solves the generalized nonlinear Schrodinger equation. For time integration, a Runge-Kutta fourth-order method is used. Matrix stability analysis is used to examine the method's stability. Three test problems involving the motion of a single solitary wave, the interaction of two solitary waves, and a solution that blows up in finite time, respectively, demonstrate the accuracy and efficiency of the provided method. Finally, the numerical results obtained from the presented method are compared with the exact solution and those obtained in earlier works available in the literature.

In present paper, a new technique is developed named Uniform Algebraic Trigonometric (UAT) tension B-spline based DQM, for the better numerical approximation of coupled 1D and coupled 2D Burgers' equation. The formula of UAT tension B-spline of order 4 is developed in this paper, by using above-mentioned basis function. Strong stability preserving SSP-RK 43 method is implemented to solve obtained ordinary differential equations. To check the effectiveness and accuracy of the proposed scheme, L 2 and L ∞ error norms are used. Four numerical experiments are discussed to put light upon applicability of scheme. By means of tables and figures, results obtained are given. Obtained results are in a good agreement with existing ones. Stability of proposed scheme is also checked by implementing the Matrix stability analysis method. This novel scheme will surely open some new dimensions in the field of numerical approximation, for other researchers in their future work. [ABSTRACT FROM AUTHOR]

In this paper, a new numerical method named Barycentric Lagrange interpolation-based differential quadrature method is implemented to get numerical solution of 1D and 2D coupled nonlinear Schrödinger equations. In the present study, spatial discretization is done with the aid of Barycentric Lagrange interpolation basis function. After that, a reduced system of ordinary differential equations is solved using strong stability, preserving the Runge-Kutta 43 method. In order to check the accuracy of the proposed scheme, we have used the formula of L ∞ error norm. The matrix stability analysis method is implemented to test the proposed method's stability, which confirms that the proposed scheme is unconditionally stable. The present scheme produces better results, and it is easy to implement to obtain numerical solutions of a class of partial differential equations. [ABSTRACT FROM AUTHOR]

In this paper, two-dimensional (2D) heat equations on disjoint rectangles are considered. The solutions are connected by interface Robin's-type internal conditions. The problem has external Dirichlet boundary conditions that, in the forward (direct) formulation, are given functions. In the inverse problem formulation, the Dirichlet conditions are unknown functions, and the aim is to be reconstructed upon integral observations. Well-posedness both for direct and inverse problems is established. Using the given 2D integrals of the unknown solution on each of the domains and the specific interface boundary conditions, we reduce the 2D inverse problem to a forward heat 1D one. The resulting 1D problem is solved using the explicit Saul'yev finite difference method. Numerical test examples are discussed to illustrate the efficiency of the approach. [ABSTRACT FROM AUTHOR]

An implicit causal space-time Galerkin scheme applied to the contrast current density volume integral equation gives rise to amarching-on-in-time scheme known as MOT-JVIE, which is accelerated and stabilized via a fully embedded FIR filter to compute the electromagnetic scattering from high permittivity dielectric objects discretized with over a million voxels. A review of two different acceleration approaches, previously developed for two-dimensional time-domain surface integral equations based on fast Fourier transforms (FFTs), leads to an understanding why these schemes obtain the same order of acceleration and the extension of this FFTacceleration to a three-dimensional MOT-JVIE. The positive definite stability analysis (PDSA) for the MOT-JVIE shows that the number of voxels for a stable MOT-JVIE discretization is restricted by the finite precision of the matrix elements. The application of the PDSA provides the insight that stability can be enforced through regularization, at the cost of accuracy. To minimize the impact in accuracy, FIR-regularization is introduced, which is based on low group-delay linear-phase high-pass FIR-filters. We demonstrate the capabilities of the FFT-accelerated FIR-regularized MOT-JVIE for a number of numerical experiments with high permittivity dielectric scatterers. [ABSTRACT FROM AUTHOR]

This paper presents a comprehensive study on the spectral properties of mimetic finite-difference operators and their application in the robust fluid–structure interaction (FSI) analysis of aircraft wings under uncertain operating conditions. By delving into the eigenvalue behavior of mimetic Laplacian operators and extending the analysis to stochastic settings, we develop a novel stochastic mimetic framework tailored for addressing uncertainties inherent in the fluid dynamics and structural mechanics of aircraft wings. The framework integrates random matrix theory with mimetic discretization methods, enabling the incorporation of uncertainties in fluid properties, structural parameters, and coupling conditions at the fluid–structure interface. Through spectral and localization analysis of the coupled stochastic mimetic operator, we assess the system's stability, sensitivity to perturbations, and computational efficiency. Our results highlight the potential of the stochastic mimetic approach for enhancing reliability and robustness in the design of aircraft wings, paving the way for optimization algorithms that integrate uncertainties directly into the design process. Our findings reveal a significant impact of stochastic perturbations on the spectral radius and eigenfunction localization, indicating heightened system sensitivity. The introduction of randomized singular value decomposition (RSVD) within our framework not only enhances computational efficiency but also preserves accuracy in low-rank approximations, which is critical for handling large-scale systems. Moreover, Monte Carlo simulations validate the robustness of our stochastic mimetic framework, showcasing its efficacy in capturing the nuanced dynamics of FSI under uncertainty. This study contributes to the fields of numerical methods and aerospace engineering by offering a rigorous and scalable approach for conducting uncertainty-aware FSI analysis, which is crucial for the development of safer and more efficient aircraft. [ABSTRACT FROM AUTHOR]

The crucial role of urban industrial land redevelopment in sustainable urban renewal has garnered widespread attention. While some scholars have explored the interest game among stakeholders in industrial land redevelopment, they primarily focus on the government-led model. Moreover, there remains a research gap concerning the impact of government intervention on the redevelopment of industrial land. This article utilizes evolutionary game theory to investigate the interest game between local governments and original land users in the model of urban industrial land redevelopment dominated by original land users. We establish evolutionary game models considering incentives and the combination of incentives and regulations, explore the interest balance strategy, and examine the impact of positive incentives and mandatory regulations on industrial land redevelopment. Furthermore, we employ a numerical simulation to unveil the impact of initial strategies and parameter adjustments on game strategy. The research results are as follows: (1) Under the original land user-led redevelopment model, only two evolutionary stability strategies exist: either the original land users implement industrial land redevelopment with positive responses from local governments, or neither party advances the process. (2) Government intervention is pivotal in facilitating the redevelopment of inefficient industrial land as economic subsidies and punitive measures motivate more participants to adopt proactive strategies. (3) The increase in government support positively correlates with the likelihood of industrial land redevelopment implementation by original land users. (4) The interests and costs of original land users emerge as crucial parameters influencing strategic decisions. This study enriches the understanding of the interests of core participants in industrial land redevelopment and provides valuable insights for sustainable urban renewal. [ABSTRACT FROM AUTHOR]

We study parameterizations of Lyapunov--Krasovskii functionals (LKFs) to analyze the stability of linear time-delay systems. We discuss the solution to the delay Lyapunov matrix, which constructs an LKF associated with a prescribed time derivative, and relate it to the approaches commonly used in the numerical computation of LKFs. We then compare two approaches for the stability analysis of time-delay systems based on semidefinite programming, namely the method based on integral inequalities and the method based on sum-of-squares programming, which have recently emerged as optimization-based methods to compute LKFs. We discuss their main assumptions and establish connections between both methods. Finally, we formulate a projection-based method allowing us to use general sets of functions to parameterize LKFs, thus encompassing the sets of polynomial functions in the literature. The solutions of the proposed stability conditions and the construction of the corresponding LKFs as stability certificates are illustrated with numerical examples. [ABSTRACT FROM AUTHOR]

The Scheduled Relaxation Jacobi (SRJ) method is a promising technique for the solution of the Poisson equation since it is possible to improve the spectral radius of the method through a proper choice of the relaxation parameters and it is also embarrassingly parallel. This is accomplished in the present study by first developing an expression for the spectral radius of the method using matrix stability analysis. The analysis is shown to provide interesting insights into the method such as the narrowing of the spectrum of the eigenvalues of the iteration matrix as well as bounds on the relaxation factors. Next, a fast and robust procedure is proposed for optimizing the values for the relaxation factors and the corresponding number of iterations. The predicted optimal values as well as those given in an earlier work are used to solve a model Poisson problem with Neumann boundary conditions and a comparison of the performance is presented. The method has been extended to cases where a non-uniform mesh is used. Optimal values for the iteration parameters are obtained for a sample non-uniform mesh. Rapid convergence is demonstrated when the SRJ method is used with these values to solve the transformed Poisson equation for the model problem. Extension of the method to levels beyond 5 is shown to be straightforward and optimal values for this case are reported. These show that further reductions in the value of the spectral radius are possible. [ABSTRACT FROM AUTHOR]

In this paper, we propose a modified cubic B-spline differential quadrature method (MCB-DQM) to solve three-dimensional (3D) coupled viscous Burger equation with appropriate initial and boundary conditions. In this method, modified cubic B-spline is treated as a basis function in the differential quadrature method (DQM) to compute the weighting coefficients. In this way, the Burger equation is reduced into a system of ordinary differential equations. An optimal strong stability-preserving Runge-Kutta (SSP-RK) method is employed to solve the resulting system of ordinary differential equations. In order to illustrate the accuracy and efficiency of the proposed method, a numerical problem is considered. From the numerical experiment, it is found that the computed result is in good agreement with the exact solution. Stability analysis of the method is also carried out using the matrix stability analysis method and found to be unconditionally stable. [ABSTRACT FROM AUTHOR]