scholarly journals On the Numerical Solution of One-Dimensional Nonlinear Nonhomogeneous Burgers’ Equation

2014 ◽  
Vol 2014 ◽  
pp. 1-15 ◽  
Author(s):  
Maryam Sarboland ◽  
Azim Aminataei

The nonlinear Burgers’ equation is a simple form of Navier-Stocks equation. The nonlinear nature of Burgers’ equation has been exploited as a useful prototype differential equation for modeling many phenomena. This paper proposes two meshfree methods for solving the one-dimensional nonlinear nonhomogeneous Burgers’ equation. These methods are based on the multiquadric (MQ) quasi-interpolation operatorℒ𝒲2and direct and indirect radial basis function networks (RBFNs) schemes. In the present schemes, the Taylors series expansion is used to discretize the temporal derivative and the quasi-interpolation is used to approximate the solution function and its spatial derivatives. In order to show the efficiency of the present methods, several experiments are considered. Our numerical solutions are compared with the analytical solutions as well as the results of other numerical schemes. Furthermore, the stability analysis of the methods is surveyed. It can be easily seen that the proposed methods are efficient, robust, and reliable for solving Burgers’ equation.

2002 ◽  
Vol 13 (3) ◽  
pp. 283-320 ◽  
Author(s):  
M. J. WARD ◽  
J. WEI

Equilibrium solutions to the one-dimensional Gierer–Meinhardt model in the form of sequences of spikes of different heights are constructed asymptotically in the limit of small activator diffusivity ε. For a pattern with k spikes, the construction yields k1 spikes that have a common small amplitude and k2 = k− k1 spikes that have a common large amplitude. A k- spike asymmetric equilibrium solution is obtained from an arbitrary ordering of the small and large spikes on the domain. It is shown that such solutions exist when the inhibitor diffusivity D is less than some critical value Dm that depends upon k1, on k2, and on other parameters associated with the Gierer–Meinhardt model. It is also shown that these asymmetric k-spike solutions bifurcate from the symmetric solution branch sk, for which k spikes have equal height. These asymmetric solutions provide connections between the branch sk and the other symmetric branches sj , for j = 1,…, k− 1. The stability of the asymmetric k-spike patterns with respect to the large O(1) eigenvalues and the small O(ε2) eigenvalues is also analyzed. It is found that the asymmetric patterns are stable with respect to the large O(1) eigenvalues when D > De, where De depends on k1 and k2, on certain parameters in the model, and on the specific ordering of the small and large spikes within a given k-spike sequence. Numerical values for De are obtained from numerical solutions of a matrix eigenvalue problem. Another matrix eigenvalue problem that determines the small eigenvalues is derived. For the examples considered, it is shown that the bifurcating asymmetric branches are all unstable with respect to these small eigenvalues.


2013 ◽  
Vol 18 (4) ◽  
pp. 463-479
Author(s):  
Natalia Boal ◽  
Francisco J. Gaspar ◽  
Francisco J. Lisbona ◽  
Petr N. Vabishchevich

This paper deals with the stabilization of the poroelasticity system, in the incompressible fully dynamic case. The stabilization term is a perturbation of the equilibrium equation that allows us to use central difference schemes to approximate the first order spatial derivatives, yielding numerical solutions without oscillations independently of the chosen discretization parameters. The perturbation term is a discrete Laplacian of the forward time difference, affected by a stabilization parameter depending on the mesh size and the properties of the porous medium. In the one dimensional case, this parameter is shown to be optimal. Some numerical experiments are presented to show the efficiency of the proposed stabilization technique.


Mathematics ◽  
2019 ◽  
Vol 7 (2) ◽  
pp. 113 ◽  
Author(s):  
Muaz Seydaoğlu

An efficient technique is proposed to solve the one-dimensional Burgers’ equation based on multiquadric radial basis function (MQ-RBF) for space approximation and a Lie-Group scheme for time integration. The comparisons of the numerical results obtained for different values of kinematic viscosity are made with the exact solutions and the reported results to demonstrate the efficiency and accuracy of the algorithm. It is shown that the numerical solutions concur with existing results and the proposed algorithm is efficient and can be easily implemented.


2021 ◽  
Vol 83 (3) ◽  
Author(s):  
Ginger Egberts ◽  
Fred Vermolen ◽  
Paul van Zuijlen

AbstractTo deal with permanent deformations and residual stresses, we consider a morphoelastic model for the scar formation as the result of wound healing after a skin trauma. Next to the mechanical components such as strain and displacements, the model accounts for biological constituents such as the concentration of signaling molecules, the cellular densities of fibroblasts and myofibroblasts, and the density of collagen. Here we present stability constraints for the one-dimensional counterpart of this morphoelastic model, for both the continuous and (semi-) discrete problem. We show that the truncation error between these eigenvalues associated with the continuous and semi-discrete problem is of order $${{\mathcal {O}}}(h^2)$$ O ( h 2 ) . Next we perform numerical validation to these constraints and provide a biological interpretation of the (in)stability. For the mechanical part of the model, the results show the components reach equilibria in a (non) monotonic way, depending on the value of the viscosity. The results show that the parameters of the chemical part of the model need to meet the stability constraint, depending on the decay rate of the signaling molecules, to avoid unrealistic results.


2021 ◽  
Vol 22 (4) ◽  
pp. 2030
Author(s):  
Hela Ferjani ◽  
Hammouda Chebbi ◽  
Mohammed Fettouhi

The new organic–inorganic compound (C6H9N2)2BiCl5 (I) has been grown by the solvent evaporation method. The one-dimensional (1D) structure of the allylimidazolium chlorobismuthate (I) has been determined by single crystal X-ray diffraction. It crystallizes in the centrosymmetric space group C2/c and consists of 1-allylimidazolium cations and (1D) chains of the anion BiCl52−, built up of corner-sharing [BiCl63−] octahedra which are interconnected by means of hydrogen bonding contacts N/C–H⋯Cl. The intermolecular interactions were quantified using Hirshfeld surface analysis and the enrichment ratio established that the most important role in the stability of the crystal structure was provided by hydrogen bonding and H···H interactions. The highest value of E was calculated for the contact N⋯C (6.87) followed by C⋯C (2.85) and Bi⋯Cl (2.43). These contacts were favored and made the main contribution to the crystal packing. The vibrational modes were identified and assigned by infrared and Raman spectroscopy. The optical band gap (Eg = 3.26 eV) was calculated from the diffuse reflectance spectrum and showed that we can consider the material as a semiconductor. The density functional theory (DFT) has been used to determine the calculated gap, which was about 3.73 eV, and to explain the electronic structure of the title compound, its optical properties, and the stability of the organic part by the calculation of HOMO and LUMO energy and the Fukui indices.


2003 ◽  
Vol 2003 (43) ◽  
pp. 2735-2746 ◽  
Author(s):  
Ekaterina T. Kolkovska

We consider the one-dimensional Burgers equation perturbed by a white noise term with Dirichlet boundary conditions and a non-Lipschitz coefficient. We obtain existence of a weak solution proving tightness for a sequence of polygonal approximations for the equation and solving a martingale problem for the weak limit.


2003 ◽  
Vol 14 (08) ◽  
pp. 1087-1105 ◽  
Author(s):  
ZHONGCHENG WANG ◽  
YONGMING DAI

A new twelfth-order four-step formula containing fourth derivatives for the numerical integration of the one-dimensional Schrödinger equation has been developed. It was found that by adding multi-derivative terms, the stability of a linear multi-step method can be improved and the interval of periodicity of this new method is larger than that of the Numerov's method. The numerical test shows that the new method is superior to the previous lower orders in both accuracy and efficiency and it is specially applied to the problem when an increasing accuracy is requested.


Sign in / Sign up

Export Citation Format

Share Document