Some Results on Number Theory and Differential Equations

A multimode perturbation method for frequency response analysis of nonlinearly vibrational beams with periodic distribution parameters is proposed. The partial differential equation with spatial varying parameters for nonlinear vibration of beams with periodic parameters under harmonic excitations is derived. The procedure of the multimode perturbation method includes three main steps: first, the nonlinear partial differential equation is transformed into linear partial differential equations with varying parameters by applying perturbation method; second, the linear partial differential equations are transformed into ordinary differential equations with multimode coupling by applying Galerkin method, where multiple vibration modes of the beams are used and the equations are suitable to nonlinear vibration of periodic structures with high parameter-varying wave in wide frequency band; third, the ordinary differential equations are solved by applying harmonic balance method to obtain vibration response of the nonlinear periodic beam, which is used for characteristics analysis of frequency response and spatial mode. Furthermore, the stability problem of nonlinear harmonic vibration as multidegree-of-freedom system with periodic time-varying parameters is solved by applying direct eigenvalue analysis approach. The proposed method can incorporate multiple vibration modes into response analysis of nonlinear periodic structures and consider mode-coupling effects due to structural nonlinearity and parametric periodicity. Finally, a nonlinear beam with periodic supports under harmonic excitations is studied. Numerical results on frequency response of the beam are given to illustrate an application of the proposed method, new frequency response characteristics, and influences of periodic parameters on structural response. The results have potential application to nonlinear structural vibration control and support damage detection of nonlinear structures with periodic supports.

In this work, a novel quantum Fourier ordinary differential equation (ODE) solver is proposed to solve both linear and nonlinear partial differential equations (PDEs). Traditional quantum ODE solvers transform a PDE into an ODE system via spatial discretization and then integrate it, thereby converting the task of solving the PDE into computing the integral for the driving function f(x). These solvers rely on the quantum amplitude estimation algorithm, which requires the driving function f(x) to be within the range of [0, 1] and necessitates the construction of a quantum circuit for the oracle R that encodes f(x). This construction can be highly complex, even for simple functions like f(x) = x. An important exception arises for the specific case of f(x) = sin2(mx + c), which can be encoded more efficiently using a set of Ry rotation gates. To address these challenges, we expand the driving function f(x) as a Fourier series and propose the quantum Fourier ODE Solver. This approach not only simplifies the construction of the oracle R but also removes the restriction that f(x) must lie within [0,1]. The proposed method was evaluated by solving several representative linear and nonlinear PDEs, including the Navier–Stokes (N–S) equations. The results show that the quantum Fourier ODE solver produces results that closely match both analytical and reference solutions.

The method of Lyapunov functions is used to investigate the stability of systems with distributed and lumped parameters, described by linear partial and ordinary differential equations. The original high-order partial differential equations are represented by a first-order system of partial differential evolution equations and constraint equations; to do that, we introduce additional variables. Passing to the first-order system of partial differential equations and representing ordinary differential equations in the normal Cauchy form, we obtain a possibility to construct the Lyapunov function as a sum of integral and classical quadratic forms and develop general methods of the investigation of the stability of a broad class of systems with distributed and lumped parameters. For example, we consider the stability of the work of the wind-driven lift pump and take into account the elasticity of the shaft transmitting the torque from the wind engine to the pump.

On the smoothing property of linear delay partial differential equations

We derive a method of obtaining approximate numerical solution of linear variable-coefficient partial differential equations (PDEs) with two independent variables from Kida's optimum approximation theory. It is shown that a certain generalized filter bank implements linear PDEs. By applying generalized discrete orthogonality of Kida's optimum approximation to this filter bank, we prove that our approximate solution satisfies a given linear PDEs and the corresponding initial or boundary conditions at all given sample points, simultaneously.

This paper presents a semi-analytical solution for one-dimensional advection-diffusion equation coupled with a linear partial differential equation with constant coefficients. The mathematical model describes a grain-fumigant-air system during fumigation processes, where fumigant gas transports through a storage silo. The coupled system considers both diffusion and advection mechanisms with constant velocity and diffusivity parameters. The solution methodology employs the Laplace transformation technique to convert the partial differential equations into ordinary differential equations in the Laplace domain. The Stehfest numerical algorithm is subsequently applied to invert the Laplace transforms and obtain the time-domain solution. Numerical computations are performed using MATLAB software to simulate the fumigant concentration distributions. Graphical results illustrate the fumigant gas concentration in air versus vertical height within the silo for different time intervals. Additional plots demonstrate the fumigant concentration absorbed by grain particles over time. The analysis examines effects of varying initial gas concentration and flow velocity on the transport process. Results indicate that higher initial concentrations and increased velocities accelerate the fumigation process, requiring less time to fill the silo completely. The proposed solution provides a mathematical framework for optimizing fumigation parameters in agricultural storage applications.

Homotopy perturbation method (HPM) is a semi-analytical technique for solving linear as well as nonlinear ordinary/partial differential equations. The method may also be used to solve a system of coupled linear and nonlinear differential equations. This chapter first illustrates the basic idea of the HPM. It uses the present method to solve a nonlinear ordinary differential equation, a nonlinear partial differential equation, and a coupled nonlinear system of partial differential equations. It is worth mentioning that the HPM can also be used for handling linear ordinary and partial differential equations. The test problems demonstrated in the chapter confirm that the HPM is an efficient method for solving linear/nonlinear ordinary, partial, and coupled differential equations. The series usually converges with an increase in the number of terms, but one may not expect the compact form solution always.

17 - Engineering mathematics

In this work, we develop a method for dynamic output feedback covariance control of the state covariance of linear dissipative stochastic partial differential equations (PDEs) using spatially distributed control actuation and sensing with noise. Such stochastic PDEs arise naturally in the modeling of surface height profile evolution in thin film growth and sputtering processes. We begin with the formulation of the stochastic PDE into a system of infinite stochastic ordinary differential equations (ODEs) by using modal decomposition. A finite-dimensional approximation is then obtained to capture the dominant mode contribution to the surface roughness profile (i.e., the covariance of the surface height profile). Subsequently, a state feedback controller and a Kalman-Bucy filter are designed on the basis of the finite-dimensional approximation. The dynamic output feedback covariance controller is subsequently obtained by combining the state feedback controller and the state estimator. The steady-state expected surface covariance under the dynamic output feedback controller is then estimated on the basis of the closed-loop finite-dimensional system. An analysis is performed to obtain a theoretical estimate of the expected surface covariance of the closed-loop infinite-dimensional system. Applications of the linear dynamic output feedback controller to the linearized stochastic Kuramoto-Sivashinsky equation are presented.

The aim of this paper is to derive the dependence of the nature of a solution of a class of differential equations of n-th order with polynomial coefficients on the solutions of the corresponding characteristic algebraic equation of n-th degree.

The aim of this research is to study the Z transformation to solve linear partial difference equations. This is due to the fact that linear partial difference equations cannot be solved directly using the Z transformation, and the research concluded that direct partial difference conditions cannot be treated directly by the Z transformation except after converting the partial difference equations to linear difference equations using one of the transformation methods. The research recommended using Binomial transformation to solve partial difference equations.

Differential equations (DEs) are commonly used to describe dynamic systems evolving in one (ordinary differential equations or ODEs) or in more than one dimensions (partial differential equations or PDEs). In real data applications, the parameters involved in the DE models are usually unknown and need to be estimated from the available measurements together with the state function. In this paper, we present frequentist and Bayesian approaches for the joint estimation of the parameters and of the state functions involved in linear PDEs. We also propose two strategies to include state (initial and/or boundary) conditions in the estimation procedure. We evaluate the performances of the proposed strategy through simulated examples and a real data analysis involving (known and necessary) state conditions.

The Fourier spectral method which uses the Fourier transform and the Fourier inverse transform are applied to the nonlinear and linear partial differential equations. The partial differential equation becomes the ordinary differential equation after the Fouier transform used and it is solved by the Runge-Kutta method. The solutions of linear Burgers' equation are carefully checked by being compared with the analytical solutions and this method was applied to nonlinear Burgers' equation. The initial conditions are cos x, sin x and Gaussian distribution.

The notion of bisimulation is an important tool for analysis and synthesis of complex dynamical systems. While the notion of bisimulation for systems described by ordinary differential and difference equations is well studied in the literature, a notion of bisimulation for systems of partial differential and difference equations is practically absent. In this paper, we consider the problem of formulating a notion of bisimulation for a class of dynamical systems described by linear partial difference equations with real constant coefficients. We define two notions of bisimulation, namely, S-bisimulation and Π-induced bisimulation, and provide algebraic characterizations for the same.

Recently a new method for solving linear and integrable nonlinear partial differential equations (PDEs) in two dimensions has been introduced. For linear PDEs, this method involves the construction of appropriate integral representations for both the solution and the associated spectral function. Here we present an alternative approach for constructing these integral representations. This approach is based on the introduction of what we call a fundamental differential form, and on a reformulation of Green's formula. The new approach is illustrated for the Laplace and the modified Helmholtz equations in a convex polygon, for a dispersive evolution equation with spatial derivatives of arbitrary order on the half-line, and for the heat equation in a domain involving a moving boundary. In addition we show that the new approach can be used to construct appropriate integral representations for multidimensional linear PDEs.