scholarly journals A modified perturbation method for mathematical models with randomness: An analysis through the steady-state solution to Burgers' partial differential equation

Author(s):  
Julia Calatayud ◽  
Juan Carlos Cortés ◽  
Marc Jornet
Entropy ◽  
2020 ◽  
Vol 22 (1) ◽  
pp. 84
Author(s):  
Matthew Woolway ◽  
Byron A. Jacobs ◽  
Ebrahim Momoniat ◽  
Charis Harley ◽  
Dieter Britz

This work investigates the convergence dynamics of a numerical scheme employed for the approximation and solution of the Frank–Kamenetskii partial differential equation. A framework for computing the critical Frank–Kamenetskii parameter to arbitrary accuracy is presented and used in the subsequent numerical simulations. The numerical method employed is a Crank–Nicolson type implicit scheme coupled with a fourth order spatial discretisation as well as a Newton–Raphson update step which allows for the nonlinear source term to be treated implicitly. This numerical implementation allows for the analysis of the convergence of the transient solution toward the steady-state solution. The choice of termination criteria, numerically dictating this convergence, is interrogated and it is found that the traditional choice for termination is insufficient in the case of the Frank–Kamenetskii partial differential equation which exhibits slow transience as the solution approaches the steady-state. Four measures of convergence are proposed, compared and discussed herein.


Author(s):  
A. Yusnaeni ◽  
Kasbawati Kasbawati ◽  
Toaha Syamsuddin

AbstractIn this paper, we study a mathematical model of an immune response system consisting of a number of immune cells that work together to protect the human body from invading tumor cells. The delay differential equation is used to model the immune system caused by a natural delay in the activation process of immune cells. Analytical studies are focused on finding conditions in which the system undergoes changes in stability near a tumor-free steady-state solution. We found that the existence of a tumor-free steady-state solution was warranted when the number of activated effector cells was sufficiently high. By considering the lag of stimulation of helper cell production as the bifurcation parameter, a critical lag is obtained that determines the threshold of the stability change of the tumor-free steady state. It is also leading the system undergoes a Hopf bifurcation to periodic solutions at the tumor-free steady-state solution.Keywords: tumor–immune system; delay differential equation; transcendental function; Hopf bifurcation. AbstrakDalam makalah ini, dikaji model matematika dari sistem respon imun yang terdiri dari sejumlah sel imun yang bekerja sama untuk melindungi tubuh manusia dari invasi sel tumor. Persamaan diferensial tunda digunakan untuk memodelkan sistem kekebalan yang disebabkan oleh keterlambatan alami dalam proses aktivasi sel-sel imun. Studi analitik difokuskan untuk menemukan kondisi di mana sistem mengalami perubahan stabilitas di sekitar solusi kesetimbangan bebas tumor. Diperoleh bahwa solusi kesetimbangan bebas tumor dijamin ada ketika jumlah sel efektor yang diaktifkan cukup tinggi. Dengan mempertimbangkan tundaan stimulasi produksi sel helper sebagai parameter bifurkasi, didapatkan lag kritis yang menentukan ambang batas perubahan stabilitas dari solusi kesetimbangan bebas tumor. Parameter tersebut juga mengakibatkan sistem mengalami percabangan Hopf untuk solusi periodik pada solusi kesetimbangan bebas tumor.Kata kunci: sistem tumor–imun; persamaan differensial tundaan; fungsi transedental; bifurkasi Hopf.


2019 ◽  
Vol 6 (4) ◽  
pp. 647-656 ◽  
Author(s):  
Takayuki Yamada

Abstract A unified method for extracting geometric shape features from binary image data using a steady-state partial differential equation (PDE) system as a boundary value problem is presented in this paper. The PDE and functions are formulated to extract the thickness, orientation, and skeleton simultaneously. The main advantage of the proposed method is that the orientation is defined without derivatives and thickness computation is not imposed a topological constraint on the target shape. A one-dimensional analytical solution is provided to validate the proposed method. In addition, two-dimensional numerical examples are presented to confirm the usefulness of the proposed method. Highlights A steady state partial differential equation for extraction of geometrical shape features is formulated. The functions for geometrical shape features are formulated by the solution of the proposed PDE. Analytical solution is provided in one-dimension. Numerical examples are provided in two-dimension.


2011 ◽  
Vol 18 (1-2) ◽  
pp. 281-287 ◽  
Author(s):  
Hu Ding ◽  
Li-Qun Chen

Nonlinear models of transverse vibration of axially moving viscoelastic beams subjected external transverse loads via steady-state periodical response are numerically investigated. An integro-partial-differential equation and a partial-differential equation of transverse motion can be derived respectively from a model of the coupled planar vibration for an axially moving beam. The finite difference scheme is developed to calculate steady-state response for the model of coupled planar and the two models of transverse motion under the simple support boundary. Numerical results indicate that the amplitude of the steady-state response for the model of coupled vibration and two models of transverse vibration predict qualitatively the same tendencies with the changing parameters and the integro-partial-differential equation gives results more closely to the coupled planar vibration.


Author(s):  
S. S. Okoya

This paper is devoted to closed-form solutions of the partial differential equation:θxx+θyy+δexp(θ)=0, which arises in the steady state thermal explosion theory. We find simple exact solutions of the formθ(x,y)=Φ(F(x)+G(y)), andθ(x,y)=Φ(f(x+y)+g(x-y)). Also, we study the corresponding nonlinear wave equation.


2021 ◽  
Vol 2021 ◽  
pp. 1-11
Author(s):  
Muhammad Sinan ◽  
Kamal Shah ◽  
Zareen A. Khan ◽  
Qasem Al-Mdallal ◽  
Fathalla Rihan

In this study, we investigate the semianalytic solution of the fifth-order Kawahara partial differential equation (KPDE) with the approach of fractional-order derivative. We use Caputo-type derivative to investigate the said problem by using the homotopy perturbation method (HPM) for the required solution. We obtain the solution in the form of infinite series. We next triggered different parametric effects (such as x, t, and so on) on the structure of the solitary wave propagation, demonstrating that the breadth and amplitude of the solitary wave potential may alter when these parameters are changed. We have demonstrated that He’s approach is highly effective and powerful for the solution of such a higher-order nonlinear partial differential equation through our calculations and simulations. We may apply our method to an additional complicated problem, particularly on the applied side, such as astrophysics, plasma physics, and quantum mechanics, to perform complex theoretical computation. Graphical presentation of few terms approximate solutions are given at different fractional orders.


Author(s):  
Anna-Carina Kurth ◽  
Kevin Schmidt ◽  
Oliver Sawodny

Abstract Through chemostat reactors, organisms can be observed under laboratory conditions. Hereby, the reactor contains the biomass, whose growth can be controlled via the dilution rate respectively the speed of a pump. Due to physical limitations, input constraints need to be considered. The population density in the reactor can be described by a hyperbolic nonlinear integro partial differential equation of first order. The steady-states and generalized eigenvalues and -modes of these integro partial differential equation are determined. In order to track a desired reference trajectory an optimal and an inversion-based feedforward control are designed. For the optimal feedforward control, the singular arc of the control is calculated and a switching strategy is stated, which explicitly considers the input constraints. For the inversion-based feedforward control, the integro partial differential equation is first linearized around the steady-state. To comply with the input constraints a control system simulator is designed. For the simulation model, the integro partial differential equation is approximated using Galerkin's method. Simulations show the functionality of the designed controls and provide the basis for comparison. The inversion-based feedforward control operates well near the steady-state, whereas the performance of the optimal feedforward control is not bounded to the proximity to the steady-state.


2010 ◽  
Vol 133 (2) ◽  
Author(s):  
Hossein Rastgoftar ◽  
Mohammad Eghtesad ◽  
Alireza Khayatian

In this paper, an analytical method and a partial differential equation based solution to control temperature distribution for functionally graded (FG) plates is introduced. For the rectangular FG plate under consideration, it is assumed that the material properties such as thermal conductivity, density, and specific heat capacity vary in the width direction, and the governing heat conduction equation of the plate is a second-order partial differential equation. Using Lyapunov’s theorem, it is shown that by applying controlled heat flux through the boundary of the domain, the temperature distribution of the plate will approach a desired steady-state distribution. Numerical simulation is provided to verify the effectiveness of the proposed method such that by applying the boundary transient heat flux, in-domain distributed temperature converges to its desired steady-state temperature.


Sign in / Sign up

Export Citation Format

Share Document