Finite Difference Method Research Articles

Mixed, quantum-classical description of electron density transfer in the collision process

In this work, we investigate an ion-atom model describing the time-dependent evolution of electron density during the collision. For a S 3 + − H system, numerical simulations are based on classical trajectory calculations, and the electron density behaviour is described with the time-dependent Schrödinger equation. We apply the finite difference method to obtain quantitative insights into the charge transfer dynamics, providing detailed information about the spatial and temporal evolution of the collision process. The results are given for representative examples of the collision, from eV to keV range of energies, in head-on collision as well as for different values of impact parameter. A validity and precision of the proposed model and interpretation of the particle collision in terms of eigenstates are also discussed.

Non-central m-point formula in method of lines for solving the Korteweg-de Vries (KdV) equation

AbstractThe present study is committed to devising efficient spatial discretization with two non-central difference formulae incorporated in the method of lines (MOL). The method is then implemented numerically on the renowned dispersive evolution equation, the Korteweg-de Vries (KdV) model while infusing Euler and fourth-order Rung-Kutta (RK4) methods, respectively. The resulting schemes are proven to be numerically stable using Fourier’s stability approach, with the MOL matrix admitting negative real eigenvalues. Moreover, this proposal has been assessed on certain initial-boundary value problems of the KdV model by examining so many factors, like the percentage errors, absolute error differences, the error norms, and the invariants $$I_{k},$$ I k , for $$k=1,2,3,$$ k = 1 , 2 , 3 , among others. In this work, we solve KdV by applying a new approach of MOL to improve the results presented in some published articles. Lastly, the assessment revealed that our proposal was found to be better than those of the exponential finite-difference method (Exponential) and the Heat Balance Integral (HBI) methods, respectively, which serve as competing approaches for validation, in terms of both the accuracy and efficiency. The numerical examples are obtained by using MATLAB.

Numerical modeling of mixed convective nanofluid flow with fractal stochastic heat and mass transfer using finite differences

This study presents the first comprehensive numerical simulation of heat and mass transfer in fractal-like mixed convective nanofluid flows. The flow of non-Newtonian nanofluids over flat and oscillating sheets is modelled mathematically, and a finite difference scheme is used to solve this model. The two-stage scheme can tackle fractal and fractal stochastic mathematical models of partial differential equations. The consistency in the mean square is proved, and Fourier series stability analysis is adopted to find stability conditions for fractal stochastic partial differential equation. The scheme is applied to solve the unsteady Casson nanofluid flow over the flat and oscillatory sheet, which affects thermal radiation, heat source, and chemical reaction. The existence of the solution is also provided for the Navier-Stokes equation of the considered flow model using fractal time derivative. The graph illustrates that the proposed fractal technique achieves faster convergence than the Crank-Nicolson approach. Applications in energy systems, materials science, and environmental engineering are just a few of the domains that could benefit from a better understanding of mixed convective nanofluid flows with fractal features, and that is what this research study hopes to accomplish. Scientists and engineers may better develop efficient and environmentally friendly systems by simulating and analyzing these complicated processes with the suggested finite difference technique.

A new efficient two-step iterative method for solving absolute value equations

PurposeThis study presents a two-step numerical iteration method specifically designed to solve absolute value equations. The proposed method is valuable and efficient for solving absolute value equations. Several numerical examples were taken to demonstrate the accuracy and efficiency of the proposed method.Design/methodology/approachWe present a two-step numerical iteration method for solving absolute value equations. Our two-step method consists of a predictor-corrector technique. The new method uses the generalized Newton method as the predictor step. The four-point open Newton-Cotes formula is considered the corrector step. The convergence of the proposed method is discussed in detail. This new method is highly effective for solving large systems due to its simplicity and effectiveness. We consider the beam equation, using the finite difference method to transform it into a system of absolute value equations, and then solve it using the proposed method.FindingsThe paper provides empirical insights into how to solve a system of absolute value equations.Originality/valueThis paper fulfills an identified need to study absolute value equations.

Analysis of a Narrow-Stencil Finite Difference Method for Approximating Viscosity Solutions of Fully Nonlinear Second Order Parabolic PDEs

Analysis of a Narrow-Stencil Finite Difference Method for Approximating Viscosity Solutions of Fully Nonlinear Second Order Parabolic PDEs

Numerical study of shock wave boundary layer interaction flows using a novel turbulence model

This study proposes a novel partial average fluctuation velocity (PAFV) turbulence based on the PAFV and the average strain rate tensor, which is used to calculate the shock wave boundary layer interaction flows, including cascade flow and transonic compressor flow. In order to illustrate the adequacy, the calculation results of several other turbulence models were also compared with PAFV, including k-ω, shear stress transport, SA-neg (negative Spalart-Allmaras), and stress-BSL (ω-based Reynolds stress model). The suitability of the PAFV turbulence model was initially tested by calculating the transonic flow in the L030-4 cascade. Subsequently, numerical computations were performed for the 3D (three-dimensional) National Aeronautics and Space Administration Rotor 67 transonic compressor rotor. In the numerical simulation, the flow control equations were transformed in a general curvilinear coordinate system for precision enhancement, and the spatial discretization of control equations was executed using the finite difference method. The convection term was processed using the Steger–Warming flux vector splitting method, the diffusion term was discretized using a second-order central difference scheme, and the time derivative term was discretized using the third-order Runge–Kutta method. A parallel computing strategy with multi-block partitioning was employed to speed up the calculations and facilitate faster convergence alongside a local time step method. All numerical computation procedures were written in Python, revealing the following: (1) The PAFV turbulence model aptly simulates the transonic flow within the cascade channel. The shock wave pattern in the cascade channel aligns well with the experimental schlieren images. (2) During the Rotor 67 transonic compressor rotor calculation, the flow–pressure ratio, flow efficiency, and inlet and outlet total temperatures and pressures align well with the experimental data. Furthermore, the resolution of the flow field structure related to the interaction between the tip leakage vortex and channel shock wave is high. (3) The turbulence model uses PAFV as the turbulence velocity scale, thereby suppressing the fluctuation velocity and turbulence viscosity near shock wave areas, preventing excessive turbulence viscosity. Moreover, PAFV inherently exhibits transport properties and anisotropic characteristics, making it well-suited for handling transonic compressor flows.

Bioconvective flow of thermally radiative Casson nanofluid with aerobic microorganisms over a stretching surface

A mathematical model is proposed to analyze the flow characteristics of a thermally radiative, chemically reacting Casson nanofluid incorporating aerobic microorganisms over a stretching sheet. Nonlinear dimensionless ordinary differential equations are formulated through suitable similarity transformations. The resulting coupled system of nonlinear ordinary differential equations is numerically solved using the finite difference method with the bvp4c routine in MATLAB. The verification of the numerical solution’s accuracy is conducted by comparing it with previously established results under specific conditions pertinent to the present investigation. Skin friction coefficients, local Nusselt number, and Sherwood number are computed for various parameters and their implications for engineering applications are discussed. The investigation also explores the impact of different dimensionless parameters on velocity, temperature, mass concentration, motile microorganism density, and oxygen concentration profiles for both Casson and Newtonian nanofluids, offering a comprehensive analysis. The study highlights that suspending non-Newtonian nanofluids with fast-moving aerobic microorganisms and nanoparticles results in enhanced heat transfer rates, thus presenting potential benefits for practical applications such as the manufacturing of biofilms.

Multiple hydraulic fracture propagation simulation in deep shale gas reservoir considering thermal effects

Deep and ultra-deep shale gas will gradually become the focus of unconventional oil and gas resources in the petroleum industry. In the development of shallow shale reservoirs, the multi-cluster hydraulic fracturing of horizontal wells is the commonly applied technique to create a network of fractures and achieve reservoir stimulation. However, for deep shale reservoirs, the behavior of multiple hydraulic fracture propagation is affected by stress interference between fractures. On the other hand, it is also influenced by the temperature disturbances caused by the injection of low-temperature fracturing fluids into high-temperature formations and the resulting thermal effects. Currently, most research on multiple hydraulic fracture propagation in shale only considers the interaction between fluids and solids while neglecting thermal effects. Therefore, based on the displacement discontinuity method, finite volume method, and finite difference method, this study established a numerical model for the non-planar propagation of multi-cluster hydraulic fractures in deep shale, considering the thermal effects. And the model was validated through numerical solutions of the temperature field and double-fracture propagation. Moreover, the number of fracture clusters, pumping rate, and formation temperature were selected as influencing factors, and the differences in multiple hydraulic fracture propagation with and without considering thermal effects were analyzed. The results show that as the formation temperature increases, the thermal effects become more significant, leading to greater differences in the morphology of multi-cluster fractures. Moreover, the thermal effects are more significant when the formation temperature exceeds 90 °C. Additionally, appropriately reducing the number of fracture clusters or increasing the pumping rate can to some extent mitigate stress interference between fractures and promote uniform fracture propagation. The findings of this study contribute to understanding the influence of thermal effects on multiple hydraulic fracture propagation and provide theoretical basis for the efficient development of deep shale gas.

Nonlinear dynamics of cantilevered hyperelastic pipes conveying fluid: comparative study of linearelasticity and hyperelasticity

In light of the softness and lightweight characteristics, hyperelastic pipes are prone to large deformation. Intending to provide accurate predictions of their large deformation vibrations, the present study aims to establish a rotation-angle-based geometrically exact model for cantilevered pipes conveying fluid composed of hyperelastic materials. Additionally, given the presence of nonlinear geometric relations, it is imperative to assess the role of hyperelasticity in the nonlinear dynamics of fluid-conveying pipes in comparison to linearelastic pipes. The modeling process integrates the Mooney-Rivlin hyperelastic model with a geometrically exact description of the pipe. And the final governing equation is developed through energy analysis and the extended Hamilton's principle. Subsequently, Galerkin's method and the finite difference method are employed to address the governing equation. Following the verification of the proposed model, parametric studies are conducted in the form of planar oscillation diagrams, time histories, phase trajectories, and bifurcation diagrams. The findings of the present study highlight three critical insights. Firstly, the hyperelastic model introduces numerous complicated nonlinear terms as function of the rotation angle with high-order nonlinearity. In the present study, at least fifth-order nonlinear terms should be preserved to adequately predict large deformation vibrations. Secondly, the Mooney-Rivlin hyperelastic model imparts a hardening effect in comparison to the linearelastic model, significantly reducing vibration amplitudes. Lastly, the third-order approximate model is demonstrated to be inadequate in predicting large deformation vibrations at high flow velocity, as it contains the approximation on not only the geometric description but the hyperelasticity-related terms.

Thermal-magnetic coupling inhomogeneous temperature field analysis of cylindrical eddy current brakes

With the increasing requirements for safety and reliability, it is now necessary to analyze the thermal characteristics to the same extent as the electromagnetic design for eddy current brakes(ECB), because the working characteristic will be greatly influenced by temperature. In this paper, the transient thermal-magnetic coupling analysis of permanent magnet cylindrical eddy current brake, which is used as the recoil mechanism in artillery, is presented. First, the structure and working principle of this ECB are described in detail. Second, the analytical electromagnetic model is established by multilayer theory to derive the heat source of this ECB. Thirdly, the thermal-magnetic coupling model is established to calculate the dynamic inhomogeneous temperature field of the ECB. The heat conduction equation is solved in different regions with inhomogeneous internal heat source by Finite Difference method(FDM) in two-dimensional cylindrical coordinate. Finally, the validity of the proposed model is verified by commercial Finite Element software. The thermal characteristic and its effects on braking force is analyzed during working. The maximum value of temperature approximately rises by 16℃ in each firing. The braking force will decline by 8% considering temperature.

Theoretical study of electron mobility in AlGaN back-barrier Al2O3/InAlN/GaN heterostructures under optical phonon scattering

Considering the strong built-in electric field caused by polarization in wurtzite nitride heterostructures, the energy levels and wave functions of the two-dimensional electrons in AlyGa1−yN back-barrier Al2O3/InxAl1−xN/GaN heterostructures are calculated using the finite element difference method. The dispersion relations and electrostatic potentials of optical phonons are obtained by the transfer matrix method. The electron mobility under optical phonon scattering is studied based on the theory of Lei–Ting force balance equation. The influences of AlyGa1−yN back barrier, ternary mixed crystal effect, and size effect are also analyzed by comparing with the Al2O3/InxAl1−xN/GaN heterostructure without back barrier. It is found that the introduction of a back barrier can attenuate the effect of gate-insulating layer, but enhance the effect of barrier and channel layers on electronic states. Then, the interaction between electrons and optical phonons is weakened, so the electron mobility in AlyGa1−yN back-barrier Al2O3/InxAl1−xN/GaN heterostructure is higher than that in Al2O3/InxAl1−xN/GaN structure under the same conditions. It is also found that using a thinner AlyGa1−yN film with y = 0.25 as the back-barrier layer of Al2O3/InxAl1−xN/GaN heterostructure is more conducive to improving 2DEG mobility. These conclusions can provide references for the preparation of InxAl1−xN/GaN heterojunction devices.

Convergence analysis of a second order numerical scheme for the Flory–Huggins-Cahn-Hilliard-Navier–Stokes system

We present an optimal rate convergence analysis for a second order accurate in time, fully discrete finite difference scheme for the Cahn-Hilliard-Navier–Stokes (CHNS) system, combined with logarithmic Flory–Huggins energy potential. The numerical scheme has been recently proposed, and the positivity-preserving property of the logarithmic arguments, as well as the total energy stability, have been theoretically justified. In this paper, we rigorously prove second order convergence of the proposed numerical scheme, in both time and space. Since the CHNS is a coupled system, the standard ℓ∞(0,T;ℓ2)∩ℓ2(0,T;Hh2) error estimate could not be easily derived, due to the lack of regularity to control the numerical error associated with the coupled terms. Instead, the ℓ∞(0,T;Hh1)∩ℓ2(0,T;Hh3) error analysis for the phase variable and the ℓ∞(0,T;ℓ2) analysis for the velocity vector, which shares the same regularity as the energy estimate, is more suitable to pass through the nonlinear analysis for the error terms associated with the coupled physical process. Furthermore, the highly nonlinear and singular nature of the logarithmic error terms makes the convergence analysis even more challenging, since a uniform distance between the numerical solution and the singular limit values of is needed for the associated error estimate. Many highly non-standard estimates, such as a higher order asymptotic expansion of the numerical solution (up to the third order accuracy in time and fourth order in space), combined with a rough error estimate (to establish the maximum norm bound for the phase variable), as well as a refined error estimate, have to be carried out to conclude the desired convergence result. To our knowledge, it will be the first work to establish an optimal rate convergence estimate for the Cahn-Hilliard-Navier–Stokes system with a singular energy potential.

Compact difference scheme for the two-dimensional semilinear wave equation

In this article, we investigate the use of a high-order compact finite difference method for solving a two-dimensional damped semilinear wave equation. We use a new iterative methods that employs compact finite difference operators to approximate the second-order spatial derivatives and achieve fourth-order convergence. We establish the stability and convergence of the compact finite difference scheme in the sense of the discrete H1-norm. Finally, we present numerical results to demonstrate the accuracy and efficiency of the proposed scheme.

Numerical method for second order singularly perturbed delay differential equations with fractional order in time via fitted computational method

The numerical solution of time fractional parabolic differential equations with singular perturbations and delay is the subject of this article. An arbitrarily small perturbation parameter ɛ(0<ɛ<<1) is multiplied the highest order derivative term. The problem’s solution shows an exponential boundary layer on the right side of the spatial domain for small values of ɛ. The solution and its derivatives are discussed together with their properties and bounds. The Crank–Nicolson approach for time discretization and the fitted operator non-standard finite difference method for space discretization are used to solve the time fractional singularly perturbed time delay problem under consideration. The comparison principle and solution bound are used to examine and analyze the scheme’s stability. The scheme’s uniform convergence is examined and demonstrated. The proposed scheme has a uniform convergence order of Mx−1+(Δt)2−γ. Two numerical illustrations for various values of ɛ are taken into consideration to validate the theoretical analysis of the scheme.The constructed scheme provides a better order of convergence and more accurate results than methods found in the literature.

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.

Accurate Description of Coulombic Interactions in Organic Field‐Effect Transistors Enabled by Efficient 3D Poisson's Equation Solver with Mixed Boundary Conditions

Accurate Description of Coulombic Interactions in Organic Field‐Effect Transistors Enabled by Efficient 3D Poisson's Equation Solver with Mixed Boundary Conditions

Simulating time-harmonic acoustic wave effects induced by periodic holes/inclusions on surfaces

This paper introduces a localized meshless method to analyze time-harmonic acoustic wave propagation on curved surfaces with periodic holes/inclusions. In particular, the generalized finite difference method is used as a localized meshless technique to discretize the surface gradient and Laplace-Beltrami operators defined extrinsically in the governing equations. An absorbing boundary condition is introduced to reduce reflections from boundaries and accurately simulate wave propagation on unclosed surfaces with periodic inclusions. Several benchmark examples demonstrate the efficiency and accuracy of the proposed method in simulating acoustic wave propagation on surfaces with diverse geometries, including complex shapes and periodic holes or inclusions.

Electro-fluid-dynamics (EFD) of soft-bodied organisms swimming through mucus having dilatant, viscous, and pseudo-plastic properties

The sperm propelling mechanism has been proposed as a possible resource for soft micro-robots in confined spaces, with potential applications in biomedical engineering. Human sperm cells essentially swim through the non-Newtonian liquid (cervical mucus) to reach their target. Thus, sperm cells swimming through non-Newtonian fluids is not vital only for physiology, but also for the fabrication of swimming micro-robots. Inspired by these remarkable applications, we examine the basic mechanics of spermatozoa motility using an undulating sheet model. This undulating sheet is bounded between two rigid walls which is self-propeling in the negative axial direction. The liquid around the spermatozoa is taken as Carreau fluid with electro-osmotic properties. The application of the lubrication approximation results in the reduction of momentum equations. The resulting ODE is solved numerically via the finite difference method and MATLAB’s built-in routine bvp5c. The unknowns that are present in the boundary conditions are refined by the root-finding algorithm. Power losses, cell speed, flow rate, velocity of the fluid, and streamline pattern are visualized by graphs. The findings of this study have important implications for the designing and optimization of electrically controlled microswimmers.

Fully conservative difference schemes for the rotation-two-component Camassa–Holm system with smooth/nonsmooth initial data

This paper derives a semi-discrete conservative difference scheme for the rotation-two-component Camassa–Holm system based on its Hamiltonian invariants. Mass, momentum and energy are preserved for the semi-discrete scheme. Furthermore, a fully discrete finite difference scheme is proposed without destroying any one of the conservative laws. Combining a nonlinear iteration with a threshold strategy, the accuracy of the scheme is guaranteed. Meanwhile, this scheme captures the formation and propagation of solitary wave solutions in long time behavior under smooth/nonsmooth initial data. Remarkably, a new type of asymmetric wave breaking phenomenon is revealed in the case of the nonzero rotational parameter.

Analysis of the Sapphire Fiber Fabry-Perot Interferometer Fringe Visibility under Different Excitation Conditions

In this paper, a theoretical analysis of how the excitation conditions affect the sapphire fiber Fabry-Perot interferometer (SFPI) visibility was performed. The conditions were considered, in which an SFPI was excited by a single-mode fiber (SMF), a multimode fiber (MMF), and a fiber collimator. The finite difference method (FDM) was used to realize the numerical solution of the modal electric fields, and then, the modal excited distributions in the sapphire fiber and the SFPI visibility were calculated. The results showed that different numbers of modes were excited in sapphire fibers under different excitation conditions and finally affected the fringe visibility of the SFPI. The fiber collimator excited the fewest modes and the visibility remained at the highest level. Finally, an experiment was performed, and the experimental results agreed well with the theoretical results.