Linear B-spline finite element method for the improved Boussinesq equation

In this paper, we develop and validate a numerical procedure for solving a class of initial boundary value problems for the improved Boussinesq equation. The finite element method with linear B-spline basis functions is used to discretize the nonlinear partial differential equation in space and derive a second order system involving only ordinary derivatives. It is shown that the coefficient matrix for the second order term in this system is invertible. Consequently, for the first time, the initial boundary value problem can be reduced to an explicit initial value problem to which many accurate numerical methods are readily applicable. Various examples are presented to validate this technique and demonstrate its capacity to simulate wave splitting, wave interaction and blow-up behavior

The well-posedness and solutions of Boussinesq-type equations

We develop well-posedness theory and analytical and numerical solution techniques for Boussinesq-type equations. Firstly, we consider the Cauchy problem for a generalized Boussinesq equation. We show that under suitable conditions, a global solution for this problem exists. In addition, we derive sufficient conditions for solution blow-up in finite time.Secondly, a generalized Jacobi/exponential expansion method for finding exact solutions of non-linear partial differential equations is discussed. We use the proposed expansion method to construct many new, previously undiscovered exact solutions for the Boussinesq and modified Korteweg-de Vries equations. We also apply it to the shallow water long wave approximate equations. New solutions are deduced for this system of partial differential equations.Finally, we develop and validate a numerical procedure for solving a class of initial boundary value problems for the improved Boussinesq equation. The finite element method with linear B-spline basis functions is used to discretize the equation in space and derive a second order system involving only ordinary derivatives. It is shown that the coefficient matrix for the second order term in this system is invertible. Consequently, for the first time, the initial boundary value problem can be reduced to an explicit initial value problem, which can be solved using many accurate numerical methods. Various examples are presented to validate this technique and demonstrate its capacity to simulate wave splitting, wave interaction and blow-up behavior

Numerical approximations of second-order matrix differential equations using higher-degree splines

Many studies of mechanical systems in engineering are based on second-order matrix models. This work discusses the second-order generalization of previous research on matrix differential equations dealing with the construction of approximate solutions for non-stiff initial problems Y 00(x) = f(x, Y (x), Y 0 (x)) using higher-degree matrix splines without any dimensional increase. An estimation of the approximation error for some illustrative examples are presented by using Mathematica. Several MatLab functions have also been developed, comparing, under equal conditions, accuracy and execution times with built-in MatLab functions. Experimental results show the advantages of solving the above initial problem by using the implemented MatLab functions.The authors wish to thank for financial support by the Universidad Politecnica de Valencia [grant number PAID-06-11-2020].Defez Candel, E.; Tung ., MM.; Solis Lozano, FJ.; Ibáñez González, JJ. (2015). Numerical approximations of second-order matrix differential equations using higher-degree splines. Linear and Multilinear Algebra. 63(3):472-489. https://doi.org/10.1080/03081087.2013.873427S472489633Loscalzo, F. R., & Talbot, T. D. (1967). Spline Function Approximations for Solutions of Ordinary Differential Equations. SIAM Journal on Numerical Analysis, 4(3), 433-445. doi:10.1137/0704038Al-Said, E. A. (2001). The use of cubic splines in the numerical solution of a system of second-order boundary value problems. Computers & Mathematics with Applications, 42(6-7), 861-869. doi:10.1016/s0898-1221(01)00204-8Al-Said, E. A., & Noor, M. A. (2003). Cubic splines method for a system of third-order boundary value problems. Applied Mathematics and Computation, 142(2-3), 195-204. doi:10.1016/s0096-3003(02)00294-1Kadalbajoo, M. K., & Patidar, K. C. (2002). Numerical solution of singularly perturbed two-point boundary value problems by spline in tension. Applied Mathematics and Computation, 131(2-3), 299-320. doi:10.1016/s0096-3003(01)00146-1Micula, G., & Revnic, A. (2000). An implicit numerical spline method for systems for ODEs. Applied Mathematics and Computation, 111(1), 121-132. doi:10.1016/s0096-3003(98)10111-xDefez, E., Soler, L., Hervás, A., & Santamaría, C. (2005). Numerical solution ofmatrix differential models using cubic matrix splines. Computers & Mathematics with Applications, 50(5-6), 693-699. doi:10.1016/j.camwa.2005.04.012Defez, E., Hervás, A., Soler, L., & Tung, M. M. (2007). Numerical solutions of matrix differential models using cubic matrix splines II. Mathematical and Computer Modelling, 46(5-6), 657-669. doi:10.1016/j.mcm.2006.11.027Ascher, U., Pruess, S., & Russell, R. D. (1983). On Spline Basis Selection for Solving Differential Equations. SIAM Journal on Numerical Analysis, 20(1), 121-142. doi:10.1137/0720009Brunner, H. (2004). On the Divergence of Collocation Solutions in Smooth Piecewise Polynomial Spaces for Volterra Integral Equations. BIT Numerical Mathematics, 44(4), 631-650. doi:10.1007/s10543-004-3828-5Tung, M. M., Defez, E., & Sastre, J. (2008). Numerical solutions of second-order matrix models using cubic-matrix splines. Computers & Mathematics with Applications, 56(10), 2561-2571. doi:10.1016/j.camwa.2008.05.022Defez, E., Tung, M. M., Ibáñez, J. J., & Sastre, J. (2012). Approximating and computing nonlinear matrix differential models. Mathematical and Computer Modelling, 55(7-8), 2012-2022. doi:10.1016/j.mcm.2011.11.060Claeyssen, J. R., Canahualpa, G., & Jung, C. (1999). A direct approach to second-order matrix non-classical vibrating equations. Applied Numerical Mathematics, 30(1), 65-78. doi:10.1016/s0168-9274(98)00085-3Froese, C. (1963). NUMERICAL SOLUTION OF THE HARTREE–FOCK EQUATIONS. Canadian Journal of Physics, 41(11), 1895-1910. doi:10.1139/p63-189Marzulli, P. (1991). Global error estimates for the standard parallel shooting method. Journal of Computational and Applied Mathematics, 34(2), 233-241. doi:10.1016/0377-0427(91)90045-lShore, B. W. (1973). Comparison of matrix methods applied to the radial Schrödinger eigenvalue equation: The Morse potential. The Journal of Chemical Physics, 59(12), 6450-6463. doi:10.1063/1.1680025ZHANG, J. F. (2002). OPTIMAL CONTROL FOR MECHANICAL VIBRATION SYSTEMS BASED ON SECOND-ORDER MATRIX EQUATIONS. Mechanical Systems and Signal Processing, 16(1), 61-67. doi:10.1006/mssp.2001.1441Flett, T. M. (1980). Differential Analysis. doi:10.1017/cbo978051189719

A class of cubic and quintic spline modified collocation methods for the solution of two-point boundary value problems

This paper is concerned with the study of a class of methods for solving second and fourth-order two-point boundary-value problems. The methods under consideration are modifications of the standard cubic and quintic spline collocation techniques, and are derived by making use of recent results con- cerning the a posteriori correction of cubic and quintic interpolating spline

A discrete least squares collocation method for two-dimensional nonlinear time-dependent partial differential equations

In this paper, we develop regularized discrete least squares collocation and finite volume methods for solving two-dimensional nonlinear time-dependent partial differential equations on irregular domains. The solution is approximated using tensor product cubic spline basis functions defined on a background rectangular (interpolation) mesh, which leads to high spatial accuracy and straightforward implementation, and establishes a solid base for extending the computational framework to three-dimensional problems. A semi-implicit time-stepping method is employed to transform the nonlinear partial differential equation into a linear boundary value problem. A key finding of our study is that the newly proposed mesh-free finite volume method based on circular control volumes reduces to the collocation method as the radius limits to zero. Both methods produce a large constrained least-squares problem that must be solved at each time step in the advancement of the solution. We have found that regularization yields a relatively well-conditioned system that can be solved accurately using QR factorization. An extensive numerical investigation is performed to illustrate the effectiveness of the present methods, including the application of the new method to a coupled system of time-fractional partial differential equations having different fractional indices in different (irregularly shaped) regions of the solution domain

Un método Wavelet-Galerkin para ecuaciones diferenciales parciales parabólicas

In this paper an Adaptive Wavelet-Galerkin method for the solution ofparabolic partial differential equations modeling physical problems withdifferent spatial and temporal scales is developed. A semi-implicit timedifference scheme is applied andB-spline multiresolution structure on theinterval is used. As in many cases these solutions are known to presentlocalized sharp gradients, local error estimators are designed and an ef-ficient adaptive strategy to choose the appropriate scale for each time isdeveloped. Finally, experiments were performed to illustrate the applica-bility and efficiency of the proposed method.En este trabajo se desarrolla un método Wavelet-Galerkin Adaptativopara la resolución de ecuaciones diferenciales parabólicas que modelanproblemas físicos, con diferentes escalas en el espacio y en el tiempo. Seutiliza un esquema semi-implícito en diferencias temporales y la estructuramultirresolución de las B-splines sobre intervalo.Como es sabido que enmuchos casos las soluciones presentan gradientes localmente altos, se handiseñado estimadores locales de error y una estrategia adaptativa eficientepara elegir la escala apropiada en cada tiempo. Finalmente, se realizaronexperimentos que ilustran la aplicabilidad y la eficiencia del método pro-puestoFil: Vampa, Victoria Cristina. Universidad Nacional de La Plata. Facultad de Ingeniería; ArgentinaFil: Martín, María Teresa. Consejo Nacional de Investigaciones Científicas y Técnicas; Argentina. Universidad Nacional de La Plata. Facultad de Ingeniería; Argentin