HIV testing patterns for United States Air Force personnel, 2008-2012.

R. Kugblenu,P. Paulin,K. Tastad,J. Okulicz

DOI: https://doi.org/10.1016/j.puhe.2015.11.019

IF: 4.984

2016-04-01

Public Health

Abstract:

What problem does this paper attempt to address?

Time discretization of a nonlinear initial value problem

M.S Mock

DOI: https://doi.org/10.1016/0021-9991(76)90018-8

IF: 4.645

1976-05-01

Journal of Computational Physics

Abstract:Time-dependent models of semiconductor devices are readily formulated as a system of three partial-differential equations in three variables. Each of the equations is linear and well posed for one of the variables; however, if the equations are solved sequentially at each step, a nonlinear instability arises from the coupling terms, unless prohibitively small time steps are used. Three methods are proposed for avoiding such an instability. Convergence and stability results are obtained for each method, and the relative merits of each method are discussed from a practical viewpoint. An example of a computation using one of these methods is included.

physics, mathematical,computer science, interdisciplinary applications
High accuracy least-squares solutions of nonlinear differential equations

Daniele Mortari,Hunter Johnston,Lidia Smith

DOI: https://doi.org/10.1016/j.cam.2018.12.007

2019-05-15

Abstract:This study shows how to obtain least-squares solutions to initial and boundary value problems of ordinary nonlinear differential equations. The proposed method begins using an approximate solution obtained by any existing integrator. Then, a least-squares fitting of this approximate solution is obtained using a constrained expression, derived from Theory of Connections. In this expression, the differential equation constraints are embedded and are always satisfied. The resulting constrained expression is then used as an initial guess in a Newton iterative process that increases the solution accuracy to machine error level in no more than two iterations for most of the problems considered. An analysis of speed and accuracy has been conducted for this method using two nonlinear differential equations. For non-smooth solutions or for long integration times, a piecewise approach is proposed. The highly accurate value estimated at the final time is then used as the new initial guess for the next time range, and this process is repeated for subsequent time ranges. This approach has been applied and validated solving the Duffing oscillator obtaining a final solution error on the order of 10-12. To complete the study, a final numerical test is provided for a boundary value problem with a known solution.
Solving higher-order nonlinear Volterra integro-differential equations using two discretization methods

Boutheina Tair,Walid Slimani

DOI: https://doi.org/10.1007/s12190-024-02075-7

2024-04-13

Journal of Applied Mathematics and Computing

Abstract:The primary objective of this article is to devis an effective method for solving a high-degree nonlinear Volterra integro-differential equation. First, we conduct a comprehensive literature review to underscore the significance of such mathematical equations in scientific research. Subsequently, we establish sufficient conditions to demonstrate the existence and uniqueness of the exact solution. Then, we employ two different discretization methods, namely Nyström and collocation. Then, we obtain two approximate systems, each of which yields a solution converging to the exact solution. The convergence of these systems is rigorously demonstrated through two theorems. Finally, we present numerical tests to evaluate the effectiveness and simplicity of the proposed numerical processes.

mathematics, applied
Transforming Stiffness and Chaos

Jan Scheffel

2024-02-27

Abstract:Stiff and chaotic differential equations are challenging for time-stepping numerical methods. For explicit methods, the required time step resolution significantly exceeds the resolution associated with the smoothness of the exact solution for specified accuracy. In order to improve efficiency, the question arises whether transformation to asymptotically stable solutions can be performed, for which neighbouring solutions converge towards each other at a controlled rate. Employing the concept of local Lyapunov exponents, it is demonstrated that chaotic differential equations can be successfully transformed to obtain high accuracy, whereas stiff equations cannot. For instance, the accuracy of explicit fourth order Runge-Kutta solution of the Lorenz chaotic equations can be increased by two orders of magnitude. Alternatively, the time step can be significantly extended with retained accuracy.

Numerical Analysis,Mathematical Physics,Chaotic Dynamics
Adaptive timestepping strategies for nonlinear stochastic systems

Cónall Kelly,Gabriel J. Lord

DOI: https://doi.org/10.48550/arXiv.1610.04003

2016-10-13

Abstract:We introduce a class of adaptive timestepping strategies for stochastic differential equations with non-Lipschitz drift coefficients. These strategies work by controlling potential unbounded growth in solutions of a numerical scheme due to the drift. We prove that the Euler-Maruyama scheme with an adaptive timestepping strategy in this class is strongly convergent. Specific strategies falling into this class are presented and demonstrated on a selection of numerical test problems. We observe that this approach is broadly applicable, can provide more dynamically accurate solutions than a drift-tamed scheme with fixed stepsize, and can improve MLMC simulations.

Numerical Analysis,Probability
Variational Quantum Crank-Nicolson and Method of Lines for the Solution of Initial Value Problems

Francisco Guzman-Cajica,Francisco S. Guzman

2024-10-04

Abstract:In this paper we use a Variational Quantum Algorithm to solve Initial Value Problems with the Implicit Crank-Nicolson and the Method of Lines (MoL) evolution schemes. The unknown functions use a spectral decomposition with the Fourier basis. The examples developed to illustrate the implementation are the Advection equation, the wave equation written as a system of first order coupled equations and the viscous Burgers equation as a non-linear case. The problems are solved using: i) standard Finite Differences as the solution to compare with, ii) the State Vector Formalism (SVF), and iii) the Sampling Error Formalism (SEF). The contributions of this paper include: 1) cost functions for generic first order in time PDEs using the implicit Crank-Nicholson and the MoL, 2) detailed convergence or self-convergence tests are presented for all the equations solved, 3) a system of three coupled PDEs is solved, 4) solutions using sampling error are presented, which highlights the importance of simulating the sampling process and 5) a fast version of the SVF and SEF was developed which can be used to test different optimizers faster.

Quantum Physics,Computational Physics
Parallel-in-Time Probabilistic Numerical ODE Solvers

Nathanael Bosch,Adrien Corenflos,Fatemeh Yaghoobi,Filip Tronarp,Philipp Hennig,Simo Särkkä

2024-09-11

Abstract:Probabilistic numerical solvers for ordinary differential equations (ODEs) treat the numerical simulation of dynamical systems as problems of Bayesian state estimation. Aside from producing posterior distributions over ODE solutions and thereby quantifying the numerical approximation error of the method itself, one less-often noted advantage of this formalism is the algorithmic flexibility gained by formulating numerical simulation in the framework of Bayesian filtering and smoothing. In this paper, we leverage this flexibility and build on the time-parallel formulation of iterated extended Kalman smoothers to formulate a parallel-in-time probabilistic numerical ODE solver. Instead of simulating the dynamical system sequentially in time, as done by current probabilistic solvers, the proposed method processes all time steps in parallel and thereby reduces the span cost from linear to logarithmic in the number of time steps. We demonstrate the effectiveness of our approach on a variety of ODEs and compare it to a range of both classic and probabilistic numerical ODE solvers.

Numerical Analysis,Machine Learning
Efficient Quantum Algorithms for Nonlinear Stochastic Dynamical Systems

Abeynaya Gnanasekaran,Amit Surana,Tuhin Sahai

2023-07-30

Abstract:In this paper, we propose efficient quantum algorithms for solving nonlinear stochastic differential equations (SDE) via the associated Fokker-Planck equation (FPE). We discretize the FPE in space and time using two well-known numerical schemes, namely Chang-Cooper and implicit finite difference. We then compute the solution of the resulting system of linear equations using the quantum linear systems algorithm. We present detailed error and complexity analyses for both these schemes and demonstrate that our proposed algorithms, under certain conditions, provably compute the solution to the FPE within prescribed $\epsilon$ error bounds with polynomial dependence on state dimension $d$. Classical numerical methods scale exponentially with dimension, thus, our approach, under the aforementioned conditions, provides an \emph{exponential speed-up} over traditional approaches.

Dynamical Systems,Data Structures and Algorithms,Numerical Analysis,Quantum Physics
Nonlinear Stability and Convergence of ERKN Integrators for Solving Nonlinear Multi-Frequency Highly Oscillatory Second-Order ODEs with Applications to Semi-Linear Wave Equations

Changying Liu,Xinyuan Wu

DOI: https://doi.org/10.1016/j.apnum.2020.02.020

IF: 2.994

2020-01-01

Applied Numerical Mathematics

Abstract:In this paper, we first present the nonlinear stability and convergence of extended Runge-Kutta-Nyström (ERKN) integrators for nonlinear multi-frequency highly oscillatory second-order ODEs with a takanami number and dominant frequency matrix. We then rigorously analyse the global errors of the blend of the ERKN time integrators and the Fourier pseudospectral spatial discretisation (ERKN-FP) when applied to semi-linear wave equations. The theoretical results based on the energy techniques show that the nonlinear stability and the global error bounds are entirely independent of the takanami number, the dominant frequency matrix, and the spatial mesh size. The analysis also provides a new perspective on the ERKN time integrators. That is, the ERKN-FP methods are free from the restriction on the CFL condition when applied to semi-linear wave equations. Numerical experiments verify our theoretical analysis results and demonstrate the remarkable superiority of the ERKN time integrators in comparison with traditional time-integration methods in the literature.
A reliable technique for nonlinear evolution equations with applications in mathematical physics

Mahmoud AE Abdelrahman,Mohammed B Almatrafi

DOI: https://doi.org/10.1177/14613484241301478

2024-11-27

Journal of Low Frequency Noise Vibration and Active Control

Abstract:Journal of Low Frequency Noise, Vibration and Active Control, Ahead of Print. This study presents a reliable and precise solver for solving several kinds of nonlinear evolution equations. The subequation approach is the basic ingredient of this solver. In this regard, a solver is created to precisely resolve and depict the nonlinear evolution equations' whole wave structures. This analysis offers the whole information obtained in the representation of solutions to these equations. The robust solver offers hyperbolic, trigonometric, and rational forms of solutions. A set of test problems drawn from the literature are provided in order to validate the accuracy of the proposed solver. The solutions that were discovered could be useful for a vital recent applied science observations. The method that was demonstrated is a tool that mathematics, engineers, and physicists can use as a box solver. Theoretical analysis and reported solutions show that the suggested solver is suitable and efficacious.

acoustics
Stable numerical methods for a stochastic nonlinear Schrödinger equation with linear multiplicative noise

Xiaobing Feng,Shu Ma

DOI: https://doi.org/10.3934/dcdss.2021071

2022-01-01

Abstract:This paper is concerned with fully discrete finite element approximations of a stochastic nonlinear Schrödinger (sNLS) equation with linear multiplicative noise of the Stratonovich type. The goal of studying the sNLS equation is to understand the role played by the noises for a possible delay or prevention of the collapsing and/or blow-up of the solution to the sNLS equation. In the paper we first carry out a detailed analysis of the properties of the solution which lays down a theoretical foundation and guidance for numerical analysis, we then present a family of three-parameters fully discrete finite element methods which differ mainly in their time discretizations and contains many well-known schemes (such as the explicit and implicit Euler schemes and the Crank-Nicolson scheme) with different combinations of time discetization strategies. The prototypical \begin{document}$ \theta $\end{document} -schemes are analyzed in detail and various stability properties are established for its numerical solution. An extensive numerical study and performance comparison are also presented for the proposed fully discrete finite element schemes.
Algebraically Stable SDIRK Methods with Controllable Numerical Dissipation for First/second-Order Time-Dependent Problems

Yazhou Wang,Xiaodai Xue,Kumar K. Tamma,Nikolaus A. Adams

DOI: https://doi.org/10.1016/j.jcp.2024.113032

IF: 4.645

2024-01-01

Journal of Computational Physics

Abstract:In this paper, a family of four-stage singly diagonally implicit Runge-Kutta methods are proposed to solve first-/second-order time-dependent problems, exhibiting the following numerical properties: fourth-order accuracy in time, unconditional stability, controllable numerical dissipation, and adaptive time step selection. The BN-stability condition is employed as a constraint to optimize parameters in the Butcher table, having significant benefits, and hence is recommended for nonlinear dynamics problems in contrast to existing methods. Numerical examples involving both first- and second-order linear/nonlinear dynamics problems validate the proposed method, and numerical results reveal that the proposed methods are free from the order reduction phenomenon when applied to nonlinear dynamics problems. The performance of adaptive time-stepping using the embedded scheme is further illustrated by the phase-field modeling problem. Additionally, the advantages and disadvantages of three-stage third-order accurate algebraically stable methods are discussed. The proposed high-order time integration can be readily integrated into high-order spatial discretization methods, such as the high-order spectral element method employed in this paper, to obtain high-order discretization in space and time dimensions.
Design and dynamical behavior of a fourth order family of iterative methods for solving nonlinear equations

Alicia Cordero,Arleen Ledesma,Javier G. Maimó,Juan R. Torregrosa

DOI: https://doi.org/10.3934/math.2024415

2024-01-01

AIMS Mathematics

Abstract:In this paper, a new fourth-order family of iterative schemes for solving nonlinear equations has been proposed. This class is parameter-dependent and its numerical performance depends on the value of this free parameter. For studying the stability of this class, the rational function resulting from applying the iterative expression to a low degree polynomial was analyzed. The dynamics of this rational function allowed us to better understand the performance of the iterative methods of the class. In addition, the critical points have been calculated and the parameter spaces and dynamical planes have been presented, in order to determine the regions with stable and unstable behavior. Finally, some parameter values within and outside the stability region were chosen. The performance of these methods in the numerical section have confirmed not only the theoretical order of convergence, but also their stability. Therefore, the robustness and wideness of the attraction basins have been deduced from these numerical tests, as well as comparisons with other existing methods of the same order of convergence.

mathematics, applied
On the stability analysis of numerical schemes for solving non-linear polynomials arises in engineering problems

Mudassir Shams,Nasreen Kausar,Serkan Araci,Liang Kong

DOI: https://doi.org/10.3934/math.2024433

2024-01-01

AIMS Mathematics

Abstract:This study shows the link between computer science and applied mathematics. It conducts a dynamics investigation of new root solvers using computer tools and develops a new family of single-step simple root-finding methods. The convergence order of the proposed family of iterative methods is two, according to the convergence analysis carried out using symbolic computation in the computer algebra system CAS-Maple 18. Without further evaluations of a given nonlinear function and its derivatives, a very rapid convergence rate is achieved, demonstrating the remarkable computing efficiency of the novel technique. To determine the simple roots of nonlinear equations, this paper discusses the dynamic analysis of one-parameter families using symbolic computation, computer animation, and multi-precision arithmetic. To choose the best parametric value used in iterative schemes, it implements the parametric and dynamical plane technique using CAS-MATLAB$ ^{@}R2011b. $ The dynamic evaluation of the methods is also presented utilizing basins of attraction to analyze their convergence behavior. Aside from visualizing iterative processes, this method illustrates not only iterative processes but also gives useful information regarding the convergence of the numerical scheme based on initial guessed values. Some nonlinear problems that arise in science and engineering are used to demonstrate the performance and efficiency of the newly developed method compared to the existing method in the literature.

mathematics, applied
A high-order implicit time integration method for linear and nonlinear dynamics with efficient computation of accelerations

Daniel O'Shea,Xiaoran Zhang,Shayan Mohammadian,Chongmin Song

2024-09-20

Abstract:An algorithm for a family of self-starting high-order implicit time integration schemes with controllable numerical dissipation is proposed for both linear and nonlinear transient problems. This work builds on the previous works of the authors on elastodynamics by presenting a new algorithm that eliminates the need for factorization of the mass matrix providing benefit for the solution of nonlinear problems. The improved algorithm directly obtains the acceleration at the same order of accuracy of the displacement and velocity using vector operations (without additional equation solutions). The nonlinearity is handled by numerical integration within a time step to achieve the desired order of accuracy. The new algorithm fully retains the desirable features of the previous works: 1. The order of accuracy is not affected by the presence of external forces and physical damping; 2. numerical dissipation in the algorithm is controlled by a user-specified parameter, leading to schemes ranging from perfectly nondissipative A-stable to L-stable; 3. The effective stiffness matrix is a linear combination of the mass, damping, and stiffness matrices as in the trapezoidal rule. The proposed algorithm is shown to replicate the numerical results demonstrated on linear problems in previous works. Additional numerical examples of linear and nonlinear vibration and wave propagation are presented herein. Notably, the proposed algorithms show the same convergence rates for nonlinear problems as linear problems, and very high accuracy. Second-order time integration methods commonly used in commercial software produce significantly polluted acceleration responses for a common class of wave propagation problems. The high-order time integration schemes presented here perform noticably better at suppressing spurious high-frequency oscillations and producing reliable and useable acceleration responses.

Numerical Analysis
A novel optimal fourth-order iteration scheme for solving nonlinear problems in applied sciences

Sunil Kumar,R. Ishwariya,Moin-ud-Din Junjua,Saima Akram

DOI: https://doi.org/10.1007/s12190-024-02259-1

2024-10-06

Journal of Applied Mathematics and Computing

Abstract:Numerous real life models in applied sciences and engineering commonly deal with nonlinear equations that have to be solved using reliable numerical methods. The computational science research is continuously growing, which is a marked by either introduction of new iteration algorithms or improvement of the existing ones. Nevertheless, these numerical methods may have high computational cost, but they do have a faster rate of convergence. Hence, the aim of this paper is create new two-step fourth order iteration algorithms for finding simple roots of nonlinear equations. To achieve this goal we have taken Steffensen's method in the first sub-step and a general quadratic interpolation polynomial in the second sub-step of a two-step iterative scheme. The main theorem demonstrates the fourth-order convergence of the proposed scheme. The proposed scheme is derivative free and achieves optimal fourth order convergence by using only three function evaluations per iteration satisfying the conjecture of Kung and Traub. To validate the theoretical results and to demonstrate the effectiveness of the presented methods, we explore some nonlinear models in medical sciences such as the law of blood flow, blood rheology, fluid permeability in biogels, and thermal regulation of the human body. The numerical results of the proposed scheme are presented in terms of number of iterations, errors in consecutive iterations, computational convergence order (CCO) and CPU time (sec). Several researchers have investigated basins of attraction on simple polynomials of the form in the complex plan. However, we have studied the fractal behavior of different methods on real world problems stated above in terms of basins of attraction for comparison of their convergence regions. The proposed scheme generate the basins of attraction in less time with wider regions of convergence in comparison with the existing methods of similar kind. Dynamical analysis of the proposed derivative free scheme illustrate its superiority in comparison with some existing optimal fourth order methods involving derivatives.

mathematics, applied
Explicit numerical methods for solving singular initial value problems for systems of second-order nonlinear ODEs

B. Y. Datsko,M. V. Kutniv,Kutniv, M. V.

DOI: https://doi.org/10.1007/s11075-024-01820-0

IF: 2.37

2024-03-28

Numerical Algorithms

Abstract:A new approach for numerical solving initial value problems for systems of second-order nonlinear ordinary differential equations with a singularity of the first kind at the start point is proposed. By substitution of the independent variable , we reduce the original initial value problem on the interval [0, a ] to the equivalent one on the interval . For solving this initial value problem at the grid node of finite grid , new fourth-order explicit Runge-Kutta-type methods have been constructed. For finding the solution in other nodes of the grid, we can apply any of the standard Runge-Kutta methods or linear multistep ones, using the solution at the point , calculated by the constructed in this article methods, as an initial condition. For the proposed approach, a new effective numerical algorithm with a given tolerance has been developed.

mathematics, applied
An iterative Monte Carlo method to solve nonlinear second-order differential equations

Martín Chávez-Páez,Enrique González-Tovar,Guillermo Iván Guerrero-García

2024-08-02

Abstract:The Monte Carlo method is a thriving and mathematically beautiful numerical technique used extensively, nowadays, to deal with many demanding problems in diverse fields. Here, we present an iterative Monte Carlo algorithm to work out very general nonlinear second-order differential equations, with Dirichlet boundary conditions. An example of its usage is, also, reported.

Computational Physics
Asymptotic Behaviour of Time Stepping Methods for Phase Field Models

Xinyu Cheng,Dong Li,Keith Promislow,Brian Wetton

DOI: https://doi.org/10.1007/s10915-020-01391-x

2021-01-16

Journal of Scientific Computing

Abstract:Abstract Adaptive time stepping methods for metastable dynamics of the Allen–Cahn and Cahn–Hilliard equations are investigated in the spatially continuous, semi-discrete setting. We analyse the performance of a number of first and second order methods, formally predicting step sizes required to satisfy specified local truncation error $$\sigma $$ σ in the limit of small length scale parameter $$\epsilon \rightarrow 0$$ ε → 0 during meta-stable dynamics. The formal predictions are made under stability assumptions that include the preservation of the asymptotic structure of the diffuse interface, a concept we call profile fidelity. In this setting, definite statements about the relative behaviour of time stepping methods can be made. Some methods, including all so-called energy stable methods but also some fully implicit methods, require asymptotically more time steps than others. The formal analysis is confirmed in computational studies. We observe that some provably energy stable methods popular in the literature perform worse than some more standard schemes. We show further that when Backward Euler is applied to meta-stable Allen–Cahn dynamics, the energy decay and profile fidelity properties for these discretizations are preserved for much larger time steps than previous analysis would suggest. The results are established asymptotically for general interfaces, with a rigorous proof for radial interfaces. It is shown analytically and computationally that for most reaction terms, Eyre type time stepping performs asymptotically worse due to loss of profile fidelity.

mathematics, applied
A multistep algorithm for ODEs

Marius-F. Danca

DOI: https://doi.org/10.48550/arXiv.1103.1227

2011-03-07

Abstract:The objective of this paper is to prove the convergence of a linear implicit multi-step numerical method for ordinary differential equations. The algorithm is obtained via Taylor approximations. The convergence is proved following the Dahlquist theory. As an additional topic, the time stability is established too. Comparative tests between some of the most known numerical methods and this method are presented.

Chaotic Dynamics,Numerical Analysis

HIV testing patterns for United States Air Force personnel, 2008-2012.

Time discretization of a nonlinear initial value problem

High accuracy least-squares solutions of nonlinear differential equations

Solving higher-order nonlinear Volterra integro-differential equations using two discretization methods

Transforming Stiffness and Chaos

Adaptive timestepping strategies for nonlinear stochastic systems

Variational Quantum Crank-Nicolson and Method of Lines for the Solution of Initial Value Problems

Parallel-in-Time Probabilistic Numerical ODE Solvers

Efficient Quantum Algorithms for Nonlinear Stochastic Dynamical Systems

Nonlinear Stability and Convergence of ERKN Integrators for Solving Nonlinear Multi-Frequency Highly Oscillatory Second-Order ODEs with Applications to Semi-Linear Wave Equations

A reliable technique for nonlinear evolution equations with applications in mathematical physics

Stable numerical methods for a stochastic nonlinear Schrödinger equation with linear multiplicative noise

Algebraically Stable SDIRK Methods with Controllable Numerical Dissipation for First/second-Order Time-Dependent Problems

Design and dynamical behavior of a fourth order family of iterative methods for solving nonlinear equations

On the stability analysis of numerical schemes for solving non-linear polynomials arises in engineering problems

A high-order implicit time integration method for linear and nonlinear dynamics with efficient computation of accelerations

A novel optimal fourth-order iteration scheme for solving nonlinear problems in applied sciences

Explicit numerical methods for solving singular initial value problems for systems of second-order nonlinear ODEs

An iterative Monte Carlo method to solve nonlinear second-order differential equations

Asymptotic Behaviour of Time Stepping Methods for Phase Field Models

A multistep algorithm for ODEs