Linear Functional Equations and Twisted Polynomials

A certain variety of non-switched polynomials provides a uni-figure representation for a wide range of linear functional equations. This is properly adapted for the calculations. We reinterpret from this point of view a number of algorithms.

Deep brain implantable microelectrode arrays for detection and functional localization of the subthalamic nucleus in rats with Parkinson’s disease

The subthalamic nucleus(STN)is considered the best target for deep brain stimulation treatments of Parkinson’s disease(PD).It is difficult to localize the STN due to its small size and deep location.Multichannel microelectrode arrays(MEAs)can rapidly and precisely locate the STN,which is important for precise stimulation.In this paper,16-channel MEAs modified with multiwalled carbon nanotube/poly(3,4-ethylenedioxythiophene):poly(styrene sulfonate)(MWCNT/PEDOT:PSS)nanocomposites were designed and fabricated,and the accurate and rapid identification of the STN in PD rats was performed using detection sites distributed at different brain depths.These results showed that nuclei in 6-hydroxydopamine hydrobromide(6-OHDA)-lesioned brains discharged more intensely than those in unlesioned brains.In addition,the MEA simultaneously acquired neural signals from both the STN and the upper or lower boundary nuclei of the STN.Moreover,higher values of spike firing rate,spike amplitude,local field potential(LFP)power,and beta oscillations were detected in the STN of the 6-OHDA-lesioned brain,and may therefore be biomarkers of STN localization.Compared with the STNs of unlesioned brains,the power spectral density of spikes and LFPs synchronously decreased in the delta band and increased in the beta band of 6-OHDA-lesioned brains.This may be a cause of sleep and motor disorders associated with PD.Overall,this work describes a new cellular-level localization and detection method and provides a tool for future studies of deep brain nuclei.

THE EXACT MEROMORPHIC SOLUTIONS OF SOME NONLINEAR DIFFERENTIAL EQUATIONS

We find the exact forms of meromorphic solutions of the nonlinear differential equations■,n≥3,k≥1,where q,Q are nonzero polynomials,Q■Const.,and p_(1),p_(2),α_(1),α_(2)are nonzero constants withα_(1)≠α_(2).Compared with previous results on the equation p(z)f^(3)+q(z)f"=-sinα(z)with polynomial coefficients,our results show that the coefficient of the term f^((k))perturbed by multiplying an exponential function will affect the structure of its solutions.

A deep learning method based on prior knowledge with dual training for solving FPK equation

The evolution of the probability density function of a stochastic dynamical system over time can be described by a Fokker–Planck–Kolmogorov(FPK) equation, the solution of which determines the distribution of macroscopic variables in the stochastic dynamic system. Traditional methods for solving these equations often struggle with computational efficiency and scalability, particularly in high-dimensional contexts. To address these challenges, this paper proposes a novel deep learning method based on prior knowledge with dual training to solve the stationary FPK equations. Initially, the neural network is pre-trained through the prior knowledge obtained by Monte Carlo simulation(MCS). Subsequently, the second training phase incorporates the FPK differential operator into the loss function, while a supervisory term consisting of local maximum points is specifically included to mitigate the generation of zero solutions. This dual-training strategy not only expedites convergence but also enhances computational efficiency, making the method well-suited for high-dimensional systems. Numerical examples, including two different two-dimensional(2D), six-dimensional(6D), and eight-dimensional(8D) systems, are conducted to assess the efficacy of the proposed method. The results demonstrate robust performance in terms of both computational speed and accuracy for solving FPK equations in the first three systems. While the method is also applicable to high-dimensional systems, such as 8D, it should be noted that computational efficiency may be marginally compromised due to data volume constraints.

Two-Stream Approximation to the Radiative Transfer Equation:A New Improvement and Comparative Accuracy with Existing Methods

Mathematical modeling of the interaction between solar radiation and the Earth's atmosphere is formalized by the radiative transfer equation(RTE), whose resolution calls for two-stream approximations among other methods. This paper proposes a new two-stream approximation of the RTE with the development of the phase function and the intensity into a third-order series of Legendre polynomials. This new approach, which adds one more term in the expression of the intensity and the phase function, allows in the conditions of a plane parallel atmosphere a new mathematical formulation of γparameters. It is then compared to the Eddington, Hemispheric Constant, Quadrature, Combined Delta Function and Modified Eddington, and second-order approximation methods with reference to the Discrete Ordinate(Disort) method(δ –128 streams), considered as the most precise. This work also determines the conversion function of the proposed New Method using the fundamental definition of two-stream approximation(F-TSA) developed in a previous work. Notably,New Method has generally better precision compared to the second-order approximation and Hemispheric Constant methods. Compared to the Quadrature and Eddington methods, New Method shows very good precision for wide domains of the zenith angle μ 0, but tends to deviate from the Disort method with the zenith angle, especially for high values of optical thickness. In spite of this divergence in reflectance for high values of optical thickness, very strong correlation with the Disort method(R ≈ 1) was obtained for most cases of optical thickness in this study. An analysis of the Legendre polynomial series for simple functions shows that the high precision is due to the fact that the approximated functions ameliorate the accuracy when the order of approximation increases, although it has been proven that there is a limit order depending on the function from which the precision is lost. This observation indicates that increasing the order of approximation of the phase function of the RTE leads to a better precision in flux calculations. However, this approach may be limited to a certain order that has not been studied in this paper.

Implementation of the Integrated Green’s Function Method for 3D Poisson’s Equation in a Large Aspect Ratio Computational Domain

The solution of Poisson’s Equation plays an important role in many areas, including modeling high-intensity and high-brightness beams in particle accelerators. For the computational domain with a large aspect ratio, the integrated Green’s function method has been adopted to solve the 3D Poisson equation subject to open boundary conditions. In this paper, we report on the efficient implementation of this method, which can save more than a factor of 50 computing time compared with the direct brute force implementation and its improvement under certain extreme conditions.

Euler’s First-Order Explicit Method–Peridynamic Differential Operator for Solving Population Balance Equations of the Crystallization Process

Using Euler’s first-order explicit(EE)method and the peridynamic differential operator(PDDO)to discretize the time and internal crystal-size derivatives,respectively,the Euler’s first-order explicit method–peridynamic differential operator(EE–PDDO)was obtained for solving the one-dimensional population balance equation in crystallization.Four different conditions during crystallization were studied:size-independent growth,sizedependent growth in a batch process,nucleation and size-independent growth,and nucleation and size-dependent growth in a continuous process.The high accuracy of the EE–PDDO method was confirmed by comparing it with the numerical results obtained using the second-order upwind and HR-van methods.The method is characterized by non-oscillation and high accuracy,especially in the discontinuous and sharp crystal size distribution.The stability of the EE–PDDO method,choice of weight function in the PDDO method,and optimal time step are also discussed.

Energy Stable Nodal DG Methods for Maxwell’s Equations of Mixed-Order Form in Nonlinear Optical Media

In this work,we develop energy stable numerical methods to simulate electromagnetic waves propagating in optical media where the media responses include the linear Lorentz dispersion,the instantaneous nonlinear cubic Kerr response,and the nonlinear delayed Raman molecular vibrational response.Unlike the first-order PDE-ODE governing equations considered previously in Bokil et al.(J Comput Phys 350:420–452,2017)and Lyu et al.(J Sci Comput 89:1–42,2021),a model of mixed-order form is adopted here that consists of the first-order PDE part for Maxwell’s equations coupled with the second-order ODE part(i.e.,the auxiliary differential equations)modeling the linear and nonlinear dispersion in the material.The main contribution is a new numerical strategy to treat the Kerr and Raman nonlinearities to achieve provable energy stability property within a second-order temporal discretization.A nodal discontinuous Galerkin(DG)method is further applied in space for efficiently handling nonlinear terms at the algebraic level,while preserving the energy stability and achieving high-order accuracy.Indeed with d_(E)as the number of the components of the electric field,only a d_(E)×d_(E)nonlinear algebraic system needs to be solved at each interpolation node,and more importantly,all these small nonlinear systems are completely decoupled over one time step,rendering very high parallel efficiency.We evaluate the proposed schemes by comparing them with the methods in Bokil et al.(2017)and Lyu et al.(2021)(implemented in nodal form)regarding the accuracy,computational efficiency,and energy stability,by a parallel scalability study,and also through the simulations of the soliton-like wave propagation in one dimension,as well as the spatial-soliton propagation and two-beam interactions modeled by the two-dimensional transverse electric(TE)mode of the equations.

Pseudo S-Asymptotically(ω,c)-Periodic Solutions to Fractional Differential Equations of Sobolev Type

In this paper,we firstly recall some basic results on pseudo S-asymptotically(ω,c)-periodic functions and Sobolev type fractional differential equation.We secondly investigate some existence of pseudo S-asymptotically(ω,c)-periodic solutions for a semilinear fractional differential equations of Sobolev type.We finally present a simple example.

Normalized Ground State Solutions of Nonlinear Choquard Equations with Nonconstant Potential

In this paper,we mainly focus on the following Choquard equation-{△u-V(x)(I_(a*)|u|^(p))|u|^(p-2)u=λu,x∈R^(N),u∈H^(1)(R^(N))where N≥1,λ∈R will arise as a Lagrange multiplier,0<a<N and N+a/N<p<N+a+2/N Under appropriate hypotheses on V(x),we prove that the above Choquard equation has a normalized ground state solution by utilizing variational methods.

From Hölder Continuous Solutions of 3D Incompressible Navier-Stokes Equations to No-Finite Time Blowup on [ 0,∞ ]

This article gives a general model using specific periodic special functions, that is, degenerate elliptic Weierstrass P functions composed with the LambertW function, whose presence in the governing equations through the forcing terms simplify the periodic Navier Stokes equations (PNS) at the centers of arbitrary r balls of the 3-Torus. The continuity equation is satisfied together with spatially periodic boundary conditions. The yicomponent forcing terms consist of a function F as part of its expression that is arbitrarily small in an r ball where it is associated with a singular forcing expression both for inviscid and viscous cases. As a result, a significant simplification occurs with a v3(vifor all velocity components) only governing PDE resulting. The extension of three restricted subspaces in each of the principal directions in the Cartesian plane is shown as the Cartesian product ℋ=Jx,t×Jy,t×Jz,t. On each of these subspaces vi,i=1,2,3is continuous and there exists a linear independent subspace associated with the argument of the W function. Here the 3-Torus is built up from each compact segment of length 2R on each of the axes on the 3 principal directions x, y, and z. The form of the scaled velocities for non zero scaled δis related to the definition of the W function such that e−W(ξ)=W(ξ)ξwhere ξdepends on t and proportional to δ→0for infinite time t. The ratio Wξis equal to 1, making the limit δ→0finite and well defined. Considering r balls where the function F=(x−ai)2+(y−bi)2+(z−ci)2−ηset equal to −1e+rwhere r>0. is such that the forcing is singular at every distance r of centres of cubes each containing an r-ball. At the centre of the balls, the forcing is infinite. The main idea is that a system of singular initial value problems with infinite forcing is to be solved for where the velocities are shown to be locally Hölder continuous. It is proven that the limit of these singular problems shifts the finite time blowup time ti∗for first and higher derivatives to t=∞thereby indicating that there is no finite time blowup. Results in the literature can provide a systematic approach to study both large space and time behaviour for singular solutions to the Navier Stokes equations. Among the references, it has been shown that mathematical tools can be applied to study the asymptotic properties of solutions.

The Jaffa Transform for Hessian Matrix Systems and the Laplace Equation

Hessian matrices are square matrices consisting of all possible combinations of second partial derivatives of a scalar-valued initial function. As such, Hessian matrices may be treated as elementary matrix systems of linear second-order partial differential equations. This paper discusses the Hessian and its applications in optimization, and then proceeds to introduce and derive the notion of the Jaffa Transform, a new linear operator that directly maps a Hessian square matrix space to the initial corresponding scalar field in nth dimensional Euclidean space. The Jaffa Transform is examined, including the properties of the operator, the transform of notable matrices, and the existence of an inverse Jaffa Transform, which is, by definition, the Hessian matrix operator. The Laplace equation is then noted and investigated, particularly, the relation of the Laplace equation to Poisson’s equation, and the theoretical applications and correlations of harmonic functions to Hessian matrices. The paper concludes by introducing and explicating the Jaffa Theorem, a principle that declares the existence of harmonic Jaffa Transforms, which are, essentially, Jaffa Transform solutions to the Laplace partial differential equation.

Existence of a Sigh-Changing Solution Result for Logarithmic Schrödinger Equations with Weight Function

This paper is devoted to studying the existence of solutions for the following logarithmic Schrödinger problem: −div(a(x)∇u)+V(x)u=ulogu2+k(x)| u |q1−2u+h(x)| u |q2−2u, x∈ℝN.(1)We first prove that the corresponding functional I belongs to C1(HV1(ℝN),ℝ). Furthermore, by using the variational method, we prove the existence of a sigh-changing solution to problem (1).

Derivation of a Revised Tsiolkovsky Rocket Equation That Predicts Combustion Oscillations

Our study identifies a subtle deviation from Newton’s third law in the derivation of the ideal rocket equation, also known as the Tsiolkovsky Rocket Equation (TRE). TRE can be derived using a 1D elastic collision model of the momentum exchange between the differential propellant mass element (dm) and the rocket final mass (m1), in which dm initially travels forward to collide with m1 and rebounds to exit through the exhaust nozzle with a velocity that is known as the effective exhaust velocity ve. We observe that such a model does not explain how dm was able to acquire its initial forward velocity without the support of a reactive mass traveling in the opposite direction. We show instead that the initial kinetic energy of dm is generated from dm itself by a process of self-combustion and expansion. In our ideal rocket with a single particle dm confined inside a hollow tube with one closed end, we show that the process of self-combustion and expansion of dm will result in a pair of differential particles each with a mass dm/2, and each traveling away from one another along the tube axis, from the center of combustion. These two identical particles represent the active and reactive sub-components of dm, co-generated in compliance with Newton’s third law of equal action and reaction. Building on this model, we derive a linear momentum ODE of the system, the solution of which yields what we call the Revised Tsiolkovsky Rocket Equation (RTRE). We show that RTRE has a mathematical form that is similar to TRE, with the exception of the effective exhaust velocity (ve) term. The ve term in TRE is replaced in RTRE by the average of two distinct exhaust velocities that we refer to as fast-jet, vx1, and slow-jet, vx2. These two velocities correspond, respectively, to the velocities of the detonation pressure wave that is vectored directly towards the exhaust nozzle, and the retonation wave that is initially vectored in the direction of rocket propagation, but subsequently becomes reflected from the thrust surface of the combustion chamber to exit through the exhaust nozzle with a time lag behind the detonation wave. The detonation-retonation phenomenon is supported by experimental evidence in the published literature. Finally, we use a convolution model to simulate the composite exhaust pressure wave, highlighting the frequency spectrum of the pressure perturbations that are generated by the mutual interference between the fast-jet and slow-jet components. Our analysis offers insights into the origin of combustion oscillations in rocket engines, with possible extensions beyond rocket engineering into other fields of combustion engineering.

Enhancing Proficiency in Quadratic Equations and Functions Through the MILAPlus Strategy

This study investigates the efficacy of the Mathematics Independent Learning Activity Practice and Play Unite Scheme(MILAPlus)as an instructional strategy to improve the proficiency levels of Grade 9 students in quadratic equations and functions through a study carried out at Quezon National High School.The research involved 116 Grade 9 students and utilized a quantitative approach,incorporating both pre-assessment and post-assessment measures.The research utilizes a quasi-experimental design,examining the academic performance of students before and after the introduction of MILAPlus.The pre-assessment establishes a baseline,and the subsequent post-assessment measures the impact of the instructional strategy.Statistical analyses,including t-tests,assess the significance of differences in mean scores and mean percentage scores,providing quantitative insights into the effectiveness of MILAPlus.Findings from the study revealed a statistically significant improvement in both mean scores and mean percentage scores after the utilization of MILAPlus,indicating enhanced proficiency in quadratic equations and functions.The Mean Proficiency Scores(MPS)also showed a substantial increase,demonstrating a marked improvement in overall proficiency levels among Grade 9 students.In light of the results,recommendations were given including the continued utilization of MILAPlus as an instructional strategy and aligning its development with prescribed learning competencies.Emphasizing the consistent adherence to policies and guidelines for MILAPlus implementation is suggested for sustaining positive effects on students’long-term performance in mathematics.This research contributes valuable insights into the practical application and effectiveness of MILAPlus within the context of Grade 9 mathematics education at Quezon National High School.

High-Order Spatial FDTD Solver of Maxwell’s Equations for Terahertz Radiation Production

We applied a spatial high-order finite-difference-time-domain (HO-FDTD) scheme to solve 2D Maxwell’s equations in order to develop a fluid model employed to study the production of terahertz radiation by the filamentation of two femtosecond lasers in air plasma. We examined the performance of the applied scheme, in this context, we implemented the developed model to study selected phenomena in terahertz radiation production, such as the excitation energy and conversion efficiency of the produced THz radiation, in addition to the influence of the pulse chirping on properties of the produced radiation. The obtained numerical results have clarified that the applied HO-FDTD scheme is precisely accurate to solve Maxwell’s equations and sufficiently valid to study the production of terahertz radiation by the filamentation of two femtosecond lasers in air plasma.

Numerical Treatments of Functional Fredholm Integral Equation in 2D with Discontinuous Kernels

This work proposes a new definition of the functional Fredholm integral equation in 2D of the second kind with discontinuous kernels (FT-DFIE). Furthermore, the work is concerned to study this new equation numerically. The existence of a unique solution of the equation is proved. In addition, the approximate solutions are obtained by two powerful methods Toeplitz Matrix Method (TMM) and Product Nystr?m Methods (PNM). The given numerical examples showed the efficiency and accuracy of the introduced methods.

New Numerical Integration Formulations for Ordinary Differential Equations

An entirely new framework is established for developing various single- and multi-step formulations for the numerical integration of ordinary differential equations. Besides polynomials, unconventional base-functions with trigonometric and exponential terms satisfying different conditions are employed to generate a number of formulations. Performances of the new schemes are tested against well-known numerical integrators for selected test cases with quite satisfactory results. Convergence and stability issues of the new formulations are not addressed as the treatment of these aspects requires a separate work. The general approach introduced herein opens a wide vista for producing virtually unlimited number of formulations.

Exact Solutions of Forced Schrödinger Equation and How to Choose the External Forces

Schrödinger equations are very common equations in physics and mathematics for nonlinear physics to model the dynamics of wave propagation in waveguides such as power lines, atomic chains, optical fibers, and even in quantum mechanics. But all these equations are most often studied without worrying about what would happen if this equation were maintained, that is to say, had a second member synonymous with an external force. It is true that on a physical level, such equations can be considered as describing the generation of waves on a waveguide using an external force. However, the in-depth analysis of this aspect is not at the center of our reflection in this article, but for us, it is a question of proposing exact solutions to this type of equation and above all proposing the general form of the external force so that the obtaining exact solutions is possible.

高压管汇材料疲劳性能测试及P-S-N模型曲线的拟合

高压管汇作为压裂设备中的主要易损件之一,其失效危害较大。它的失效原因主要是疲劳、冲蚀、腐蚀或者材料缺陷引起的刺漏和爆裂,其中尤以疲劳失效最不可预估。目前,对于高压管汇材料的疲劳性能研究不够深入,为解决高压管汇材料疲劳寿命的准确描述问题,以某国产高压管汇材料为例,进行了一系列疲劳试验,并基于试验数据,采用多种分布模型和不同S-N模型进行拟合分析,得出综合评价拟合能力最强的P-S-N模型。结果表明,该材料在中长疲劳寿命区,Weibull三参数模型在7级应力水平下综合评价能力最好;在存活率分别为50%、90%、99%、99.9%时,指数S-N模型的拟合系数均大于0.98,拟合能力最好。得出的P-S-N模型曲线可以为高压管汇的疲劳寿命以及安全设计提供依据。