Finite Difference Method Research Articles

Analytical Properties and Numerical Preservation of an Age-Group Susceptible-Infected-Recovered Model: Application to the Diffusion of Information

Abstract This paper analyzes an age-group susceptible-infected-recovered (SIR) model. Theoretical results concerning the conservation of the total population, the positivity of the analytical solution, and the final size of the epidemic are derived. Since the model is a nonlinear system of ordinary differential equations (ODEs), a numerical approximation is considered, based on Standard and non-Standard Finite Difference methods, and on a Modified Patankar-Runge–Kutta (MPRK) method. The numerical preservation of the qualitative properties of the analytical solution is studied. The obtained results are applied to the diffusion of information in social networks, and the effectiveness of the different numerical approaches is shown through several numerical tests on real data.

DeHNSSo: The Delft Harmonic Navier-Stokes Solver for nonlinear stability problems with complex geometric features

A nonlinear Harmonic Navier-Stokes (HNS) framework is introduced for simulating instabilities in laminar spanwise-invariant shear layers, featuring sharp and smooth wall surface protuberances. While such cases play a critical role in the process of laminar-to-turbulent transition, classical stability theory analyses such as parabolized or local stability methods fail to provide (accurate) results, due to their underlying assumptions. The generalized incompressible Navier-Stokes (NS) equations are expanded in perturbed form, using a spanwise and temporal Fourier ansatz for flow perturbations. The resulting equations are discretized using spectral collocation in the wall-normal direction and finite-difference methods in the streamwise direction. The equations are then solved using a direct sparse-matrix solver. The nonlinear mode interaction terms are converged iteratively. The solution implementation makes use of a generalized domain transformation to account for geometrical smooth surface features, such as humps. No-slip conditions can be embedded in the interior domain to account for the presence of sharp surface features such as forward- or backward-facing steps. Common difficulties with Navier-Stokes solvers, such as the treatment of the outflow boundary and convergence of nonlinear terms, are considered in detail. The performance of the developed solver is evaluated against several cases of representative boundary layer instability growth, including linear and nonlinear growth of Tollmien-Schlichting waves in a Blasius boundary layer and stationary crossflow instabilities in a swept flat-plate boundary layer. The latter problem is also treated in the presence of a geometrical smooth hump and a sharp forward-facing step at the wall. HNS simulation results, such as perturbation amplitudes, growth rates, and shape functions, are compared to benchmark flow stability analysis methods such as Parabolized Stability Equations (PSE), Adaptive Harmonic Linearized Navier-Stokes (AHLNS), or Direct Numerical Simulations (DNS). Good agreement is observed in all cases. The HNS solver is subjected to a grid convergence study and a simple performance benchmark, namely memory usage and computational cost. The computational cost is found to be considerably lower than high-fidelity DNS at comparable grid resolutions. Program summaryProgram Title: DeHNSSoCPC Library link to program files:https://doi.org/10.17632/9bnms99kk2.1Developer's repository link:https://github.com/SvenWesterbeek/DeHNSSoLicensing provisions: GPLv3Programming language: MatlabSupplementary material: The supplementary material contains the code as well as a user manual.Nature of problem: Fluid flows are subject to laminar-to-turbulent transition following the growth of instabilities. To avoid computationally demanding Direct Numerical Simulations (DNS), perturbation theory is often applied to their analysis. However, classical stability methods based on the Orr-Sommerfeld equation or the Parabolized Stability Equations neglect the influence of streamwise gradients in varying degrees. The validity of these assumptions is difficult to estimate a priori.Solution method: The Delft Harmonic Navier-Stokes Solver (DeHNSSo) solves the harmonic Navier-Stokes equations nonlinearly on domains featuring sharp and smooth spanwise invariant surface features using a generalized grid approach in combination with an embedded boundary method. This allows the user to include the effects of streamwise gradients on flow at a fraction of the cost of DNS.Additional comments including restrictions and unusual features: In DeHNSSo, the equations are solved using direct matrix solvers. As such, memory is treated as a scarce resource. The problem is formulated in a mode-independent manner such that left-hand side matrices need only be computed and stored once despite incorporating nonlinear terms. Additionally, DeHNSSo offers the user the possibility to prescribe inhomogeneous boundary conditions to introduce and solve the receptivity problem or manipulate instabilities. Due to the double Fourier expansion, the solver is restricted to spanwise and temporally periodic problems.

A future scenario earthquake and ground motion hazards for Kathmandu, Nepal

Since Nepal is an earthquake-prone country due to the collision of the Indian and Eurasian plates, it is crucial for its capital of Kathmandu to evaluate ground motion hazards from a future large earthquake. For this purpose, we constructed a realistic scenario earthquake with realistic rupture parameters in a likely location. To obtain the location, a new distribution of the rupture zones of large historical earthquakes along the Main Himalayan Thrust (MHT) was obtained by integrating the distribution from a previous study and the results of trench surveys. Global Navigation Satellite System observations indicated that the plate boundary is strongly coupled from the southern boundary of the MHT to a depth of approximately 10 km and there is almost no lateral change in the coupling. This implies that all regions along the MHT have similar rates of strain increase. Therefore, it is most probable that the rupture zone of the oldest previous event will rupture as a scenario earthquake. In the new distribution, the 1255 earthquake is the oldest. However, large earthquakes occurred in 1934 and 2015 within its rupture zone. Thus, we adopted the area obtained by removing the 1934 and 2015 rupture zones from the western part of the 1255 rupture zone. The relationship between the rupture area size and seismic moment of the 2015 earthquake lies between the scaling formulas for crustal earthquakes and plate-boundary earthquakes, but is closer to the former. Therefore, using this and the scheme for characterized source models, we determined realistic rupture parameters. We then simulated broadband ground motions in Kathmandu using these rupture parameters, our 3-D velocity structure models, and a hybrid method combining the finite-difference method and the stochastic Green’s function method. We obtained the peak ground accelerations (PGAs) of simulated ground motions, and calculated the seismic intensities in the Modified Mercalli intensity scale from the PGAs as indexes of hazards for Kathmandu. Intensities IX coincide with the center of the Kathmandu Valley, and intensities VIII and VII are found in the area surrounded by the sedimentary boundary and the southernmost part of the valley.Graphical abstract

A quantitative method of calculating transient nonlinear heat partition coefficient between clutch friction discs with deep learning

In the current clutch temperature field study, the generally used constant heat partition coefficient tends to overestimate the separator disc temperature and underestimate the friction disc temperature. Although some researchers have found the characteristics of the time-varying heat partition coefficient, a suitable method is still needed to apply it to temperature calculations. This study provides a quantitative method for the application of the transient nonlinear heat partition coefficient to temperature calculations. The finite difference method is adopted to figure out the time-varying curve of the heat partition coefficient by coupling the contact temperature of the friction components. The numerical results show that the heat partition coefficient is independent of rotation speed with three stages: initial value, rapid time-varying, and steady-state. Different from the analytical method, we apply a deep learning method to train the quantisation function to characterise these three stages, avoiding complex formula derivation. As a result, the quantitative function can characterise the time-varying heat partition coefficient accurately, with an average error of 0.19%, 3.05% and 0.62% for the inert, time-varying, and steady-state stages, respectively. In addition, the accuracy of applying the quantisation function in temperature simulation is verified by friction experiments, and the error is less than 8%. This is superior to the results of solving the temperature field by a constant heat partition coefficient.

Topological optimization in 3D-magnetostatics: development of adjoint methods using the equations of magnetic moments

PurposeThe purpose of this paper is mainly to develop the adjoint method within the method of magnetic moment (MMM) and thus, to provide an efficient new way to solve topology optimization problems in magnetostatic to design 3D-magnetic circuits.Design/methodology/approachFirst, the MMM is recalled and the optimization design problem is reformulated as a partial derivative equation-constrained optimization problem where the constraint is the Maxwell equation in magnetostatic. From the Karush–Khun–Tucker optimality conditions, a new problem is derived which depends on a Lagrangian parameter. This problem is called the adjoint problem and the Lagrangian parameter is called the adjoint parameter. Thus, solving the direct and the adjoint problems, the values of the objective function as well as its gradient can be efficiently obtained. To obtain a topology optimization code, a semi isotropic material with penalization (SIMP) relaxed-penalization approach associated with an optimization based on gradient descent steps has been developed and used.FindingsIn this paper, the authors provide theoretical results which make it possible to compute the gradient via the continuous adjoint of the MMMs. A code was developed and it was validated by comparing it with a finite difference method. Thus, a topology optimization code associating this adjoint based gradient computations and SIMP penalization technique was developed and its efficiency was shown by solving a 3D design problem in magnetostatic.Research limitations/implicationsThis research is limited to the design of systems in magnetostatic using the linearity of the materials. The simple examples, the authors provided, are just done to validate our theoretical results and some extensions of our topology optimization code have to be done to solve more interesting design cases.Originality/valueThe problem of design is a 3D magnetic circuit. The 2D optimization problems are well known and several methods of resolution have been introduced, but rare are the problems using the adjoint method in 3D. Moreover, the association with the MMMs has never been treated yet. The authors show in this paper that this association could provide gains in CPU time.

Compact Low Loss Ribbed Asymmetric Multimode Interference Power Splitter

Optical power splitters (OPSs) are utilized extensively in integrated photonic circuits, drawing significant interest in research on power splitters with adjustable splitting ratios. This paper introduces a compact, low-loss 1 × 2 asymmetric multimode interferometric (MMI) optical power splitter on a silicon-on-insulator (SOI) platform. The device is simulated using the finite difference method (FDM) and eigenmode expansion solver (EME). It is possible to attain various output power splitting ratios by making the geometry of the MMI central section asymmetric relative to the propagation axis. Six distinct optical power splitters are designed with unconventional splitting ratios in this paper, which substantiates that the device can achieve any power splitter ratios (PSRs) in the range of 95:5 to 50:50. The dimensions of the multimode section were established at 2.9 × (9.5–10.9) μm. Simulation results show a range of unique advantages of the device, including a low extra loss of less than 0.4 dB, good fabrication tolerance, and power splitting ratio fluctuation below 3% across the 1500 nm to 1600 nm wavelength span.

Finite element method for fractional order parabolic obstacle problem with nonlinear source term

This paper considers the finite element method for a fractional order parabolic obstacle problem in time with a nonlinear source term. An error estimation of the spatially semi-discrete problem is discussed. A fully discrete problem is proposed by introducing the finite difference approximation for the time variable. The L2-unconditional stability and convergence of the approximate solution have been proven.

High-Q Multiband Narrowband Absorbers Based on Two-Dimensional Graphene Metamaterials

In this paper, an absorber with multi-band, tunable, high Q, and high sensitivity, based on terahertz periodic two-dimensional patterned graphene surface plasmon resonance (SPR), is proposed. The absorber consists of a bottom metal film separated by a periodically patterned graphene metamaterial structure and a SiO2 dielectric layer, where the patterned graphene layer is etched by “+” and “L” shapes and circles. It has simple structural features that can greatly simplify the fabrication process. We have analyzed the optical properties of a graphene surface plasmon perfect metamaterial absorber based on graphene in the terahertz region using the finite-difference method in time domain (FDTD). The results show that the absorber device exhibits three perfect absorption peaks in the terahertz bands of f1 = 1.55 THz, f2 = 4.19 THz, and f3 = 6.92 THz, with absorption rates as high as 98.70%, 99.63%, and 99.42%, respectively. By discussing the effects of parameters such as the geometrical dimensions of patterned graphene metamaterial structure “+” width W1, “L” width W2, circular width R, and the thickness of the dielectric layer on the absorption performance of absorber, as well as investigating the chemical potential and relaxation time of patterned-layer graphene material, it was found that the amplitude of the absorption peaks and the frequency of resonance of absorber devices can be dynamically adjusted. Finally, we simulated the spectra as the surrounding refractive index n varied to better evaluate the sensing performance of the structure, yielding structural sensitivities up to 382 GHz/RIU. Based on this study, we find that the results of our research will open new doors for the use of multi-band, tunable, polarization-independent metamaterial absorbers that are insensitive to large-angle oblique incidence.

Numerical Solution of Nonlinear Reaction-Advection-Diffusion Equation in Sense of Caputo-Fabrizio Derivative with Finite Difference and Collocation Method

Numerical Solution of Nonlinear Reaction-Advection-Diffusion Equation in Sense of Caputo-Fabrizio Derivative with Finite Difference and Collocation Method

Ultra-high Q-factor quasi-BIC BaTiO3 metasurface for electro-optic modulation

Metasurfaces play a crucial role in trapping electromagnetic waves with specific wavelengths, serving as a significant platform for enhancing light–matter interactions. In all kinds of dynamic modulation metasurfaces, electro-optic modulation metasurfaces have attracted much attention due to its advantages of fast, stable and high efficiency. In order to respond to the extremely weak refractive index change of the electro-optical effect of the materials, the metasurfaces are required to support optical signals with high Q values. The quasi-bound state in the continuum (Q-BIC) is often used to enhance the light-field modulation capability of metasurfaces and to improve the modulation sensitivity of electro-optic modulators due to its ability to generate high Q-factor resonances. However, the design of an electro-optic modulation metasurface that facilitates the application of voltage and achieves modulation efficiency of nearly 100% is still in urgent need of development. In this study, single-crystal BTO metasurfaces are modeled using finite-difference time-domain method, and the structural symmetry is broken to introduce a Q-BIC resonance to generate a high Q-factor optical signal of 2.45 × 104 for high-depth electro-optic modulation. By simulating an applied electric field of 143 V/mm on the metasurface, a slight refractive index change of BTO of 8 × 10−4 was produced, leading to an electro-optical intensity modulation depth of 100%. Furthermore, the nanostructure of the metasurface was carefully designed to facilitate nano-fabrication and voltage application, and it is ideal for the development of low-power, CMOS-compatible, and miniaturized electro-optic modulation devices. Although the results of this study are based on simulations, they provide a crucial theoretical basis and guidance for the realization of efficient and realistic design of dynamic metasurfaces.

An efficient computational scheme for solving coupled time-fractional Schrödinger equation via cubic B-spline functions.

The time fractional Schrödinger equation contributes to our understanding of complex quantum systems, anomalous diffusion processes, and the application of fractional calculus in physics and cubic B-spline is a versatile tool in numerical analysis and computer graphics. This paper introduces a numerical method for solving the time fractional Schrödinger equation using B-spline functions and the Atangana-Baleanu fractional derivative. The proposed method employs a finite difference scheme to discretize the fractional derivative in time, while a θ-weighted scheme is used to discretize the space directions. The efficiency of the method is demonstrated through numerical results, and error norms are examined at various values of the non-integer parameter, temporal directions, and spatial directions.

Finite Difference Method for Infection Model of HPV with Cervical Cancer under Caputo Operator

In this paper, a fractional model in the Caputo sense is used to characterize the dynamics of HPV with cervical cancer. Generalized mean value theorem has been used to examine whether the infection model has a unique positive solution. The model has two equilibrium points: the disease-free point and the endemic point. The examination of the system’s local and global stability is provided in terms of the basic reproductive number Rp°. The global stability analysis has been carried out using an appropriate Lyapunov function and the LaSalle invariant principle. The results demonstrate that in the infection model, if Rp°<1, then the solution converges to the disease-free equilibrium, which is both locally and globally asymptotically stable. Whilst Rp°>1, the endemic equilibrium is considered to exist. Simulations are implemented via a finite difference method with Grünwald-Letnikov discretization approach for Caputo derivative operator to define how changes in parameters impact the dynamic behavior of the system using Matlab.

Nonlinear Dynamic Properties of Rigid–Elastic–Liquid Coupled Ball Bearings Considering External Load Excitation

In this study, a novel dynamic model of a ball bearing based on rigid–elastic–liquid (REL) coupling structure considering external load excitation is established, and the dynamic response of its outer ring under external load excitation is examined. First, the differential equation of the journal motion (i.e. the outer ring) is derived based on Newton’s law, and the oil film force induced by the REL coupling structure is determined by using the central finite difference method, planar bending theory of thin ring, and Simpson’s rule. Further, the fourth and fifth-order variable step length Runge–Kutta method is used to solve the kinematic differential equation. Second, the dynamic characteristics of the outer ring of the bearing under different speeds and elastic support structural parameters are examined. Finally, the feasibility of the theoretical model is verified through experiments. The results show that compared with the traditional squeeze film damper (SFD) system, the ball bearing with REL coupling structure is more stable, and the vibration damping effect is better under higher bearing speed, lower squirrel cage stiffness, and smaller elastic ring flexibility coefficient. The above findings provide useful insights into the dynamic characteristics of REL coupled ball bearings under unbalanced excitation, which can serve as a reference for the optimization design of such ball bearings.

Threshold current density for diffusion-controlled stability of electrolytic surface nanobubbles

Understanding the stability mechanism of surface micro/nanobubbles adhered to gas-evolving electrodes is essential for improving the efficiency of water electrolysis, which is known to be hindered by the bubble coverage on electrodes. Using molecular simulations, the diffusion-controlled evolution of single electrolytic nanobubbles on wettability-patterned nanoelectrodes is investigated. These nanoelectrodes feature hydrophobic islands as preferential nucleation sites and allow the growth of nanobubbles in the pinning mode. In these simulations, a threshold current density distinguishing stable nanobubbles from unstable nanobubbles is found. When the current density remains below the threshold value, nucleated nanobubbles grow to their equilibrium states, maintaining their nanoscopic size. However, for the current density above the threshold value, nanobubbles undergo unlimited growth and can eventually detach due to buoyancy. Increasing the pinning length of nanobubbles increases the degree of nanobubble instability. By connecting the current density with the local gas oversaturation, an extension of the stability theory for surface nanobubbles [Lohse and Zhang, Phys. Rev. E 91, 031003(R) (2015)] accurately predicts the nanobubble behavior found in molecular simulations, including equilibrium contact angles and the threshold current density. For larger systems that are not accessible to molecular simulations, continuum numerical simulations with the finite difference method combined with the immersed boundary method are performed, again demonstrating good agreement between numerics and theories.

A two-grid characteristic block-centered finite difference method for Darcy-Forchheimer slightly compressible miscible displacement problem

In this paper, a two-grid characteristic block-centered finite difference method is developed for solving a nonlinear Darcy-Forchheimer slightly compressible miscible displacement problem in porous media. The method is proposed for solving the nonlinear system in two steps. In the first step, the nonlinear equations are approximated on a coarse grid using the characteristic block-centered finite difference method. In the second step, the nonlinear system is linearized using Newton's method with a small positive parameter to ensure the differentiability of the nonlinear term |u|u in the Darcy-Forchheimer equation. The proposed two-grid method is rigorously analyzed, in which a priori error estimates of the velocity, pressure, concentration and its flux are provided, and the rates of convergence are O(ε+Δt+H3+h2). Finally, numerical experiments are conducted to demonstrate the effectiveness of the proposed methods by comparing the efficiency, especially in terms of CPU time.

Finite element modeling of thermal‐hydro‐mechanical coupled processes in unsaturated freezing soils considering air‐water capillary pressure and cryosuction

AbstractThis paper presents a comprehensive computational model for analyzing thermo‐hydro‐mechanical coupled processes in unsaturated porous media under frost actions. The model employs the finite element method to simulate multiphase fluid flows, heat transfer, phase change, and deformation behaviors. A new soil freezing characteristic curve model is proposed to consider the suctions from air‐water capillary pressure and water‐ice cryosuction. A total pore pressure with components from liquid water pressure, air pressure, and ice pressure is used in the effective stress law. Vapor and dry air are considered miscible gases, utilizing the ideal gas law and Dalton's law. The governing equations encompass the linear momentum balance equation, the energy balance equation, and mass conservation equations for water species (ice, liquid, and vapor) and dry air. Weak forms are formulated based on primary variables of displacement, water pressure, air pressure, and temperature. The spatial discretization is achieved through the finite element method, while temporal discretization employs the fully implicit finite difference method, resulting in a system of fully coupled nonlinear equations. To verify the proposed computational model, a numerical implementation is developed and validated against a set of experimental data from the literature. The successful verification demonstrates the robustness of the model. A detailed discussion of the contributions from phase change strain and different sources of pore pressure is also addressed.

Calculation of Plates of Different Shapes Using Generalized Equations of the Finite Difference Method

Statement of the problem. The article reviewed nonaxisymmetric tasks for calculating slabs in shape of a circle sector, wedge and a slab formed by a connection of rectangular and circular areas. The solution is based on generalized equations of the finite differences method. The algorithm allows us to take into account finite discontinuities of the desired function, the right side of the original differential equations, as well as discontinuities of derivatives of these functions without accumulation of the grid and using points outside the contour. It resolves oneself into a system of differential algebraic equations generated for field and contour points. In this paper the initial differential equations and their numerical approximation are provided. Results. Calculations of a semicircle and a wedge under different boundary conditions are presented, as well as an uncut slab in the form of a semicircle with an annular support and a composite slab are consi-dered. The solution of the above problems has been carried out on the minimum possible computational grid, comparison with the available solutions has been provided, and static checks of the calculation have been carried out. Conclusions. Along with the most common numerical method for calculating buildings and structures, it is possible and quite effective to use these equations. It is proved by means of some examples that such an approximation has a number of advantages: the results allow us to talk about the sufficient accuracy of the solution with a small number of partitions; fundamental simplicity and universality of the method.

Solving Time-fractional Order Radon Diffusion Equation in Water by Finite Difference Method

Solving Time-fractional Order Radon Diffusion Equation in Water by Finite Difference Method

Analysis of Rectangular Plate with Opening by Finite Difference Method

Steel plates are commonly used to support lateral or vertical loads. Before the design of such a plate, analysis is performed to check the stability of plate for the proposed load. There are several methods for this analysis. In this research, a finite difference method (FDM) is proposed to analyze a rectangular steel plate. The plate is considered to be subjected to an arbitrary transverse uniformly distributed loading and is considered to be clamped at the two opposite edges and free at the other two edges. At first the plate is analyzed considering it to be completely solid and then it was re-analyzed considering a rectangular opening at its centre. The ordinary finite difference method is used to solve the governing differential equation of the plate deflection. The proposed method can be easily programmed to readily apply on a plate problem. The work covers the determination of displacement components at different points of the plate and checking the result by software (STAAD.Pro) analysis.

Numerical simulation of unsteady blood based Carreau nanofluid flow in an inclined porous saturated artery having overlapping stenosis with electro-osmosis

The present paper focuses on the computational modeling of the unsteady non-Newtonian electro-magnetohydrodynamic flow of nano-blood through the inclined and tapered saturated porous artery with overlapping stenosis. The Carreau fluid model is adopted to simulate the non-Newtonian behavior of the blood. This study also uses a modified form of Darcy’s law applied to the Carreau model. Blood is doped with a suspension of homogeneous biocompatible nanoparticles of different shapes. In addition, a combination of a magnetic field and an electric field has also been taken into account. The significance of this work rests in the ability to enhance comprehension of the complicated behavior of hemodynamic flow when exposed to external fields and porous media with the addition of nanoparticles, which has substantial implications for the regulation and treatment of cardiovascular diseases. The finite-difference method is used to obtain the solution of the resulting equations. The simulations are applicable to targeted magnetic therapy for stenotic artery disease and the delivery of nanomedicines.