Numerically Research Articles

Dispersivity calculation in digital twins of multiscale porous materials using the micro-continuum approach

The micro-continuum method is a novel approach to simulate flow and transport in multiscale porous materials. For such materials, the domain can be divided into three sub-domains depending on the local porosity ε: fully resolved solid phase, for which ε = 0, fully resolved pores, for which ε = 1.0, and unresolved pores, for which 0 < ε < 1.0. For such domains, the flow can be solved using the Darcy-Brinkman-Stokes (DBS) equation, which offers a seamless transition between unresolved pores, where flow is described by Darcy’s law, and resolved pores, where flow is described by the Stokes equations. Species transport can then be modelled using a volume-averaged equation. In this work, we present a derivation of the closure problem for the micro-continuum approach. Effective dispersivity tensors can then be calculated through a multi-stage process. First, high resolution images are chosen for characterizing the structure of the unresolved pores. Porosity, permeability and effective dispersivity for the unresolved parts are calculated by solving a closure problem based on Direct Numerical Simulation (DNS) in the high-resolution images. The effective dispersivity is then expressed as a function of the Péclet number, which describes the ratio of advective to diffusive transport. This relationship, along with porosity and permeability, is then integrated into the multiscale domain and the effective dispersivity tensor for the full image is calculated. Our novel method is validated by comparison with the numerical solution obtained for a fully-resolved simulation in a multiscale 2D micromodel. It is then applied to obtain an effective dispersivity model in digital twins for two multiscale materials: hierarchical ceramic foams and microporous carbonate rocks.

Type-2 fuzzy initial value problems under granular differentiability

This article investigates type-2 fuzzy initial value problems and introduces a novel strategy that capitalises on granular differentiability. Incorporating type-2 fuzzy numbers to depict the problem’s uncertainty may be advantageous from a practical standpoint. This work employs triangularly perfect quasi type-2 fuzzy numbers (TPQT2FNs) and defines the granular differentiability of TPQT2FN-valued functions. In addition, the solution approach for initial value problems with type-2 fuzzy initial conditions is discussed in the context of granular differentiability by transforming the type-2 fuzzy problem into a type-1 fuzzy problem using the lower membership function (LMF) and upper membership function (UMF) concepts. A couple of numerical examples are then examined to determine the applicability of the proposed method, and comparisons are made with existing type-2 fuzzy results and, in a special case, type-1 fuzzy results. In order to aid readers’ comprehension and study the behaviour of the numerical solution, three-dimensional graphical results are also shown.

On the influence of heat-induced evaporation over interfacial atomization driven by surface acoustic waves

In this work, the influence of evaporative flux over the atomization driven by surface acoustic waves (SAWs) of a Newtonian water drop is studied. The drop is placed on a heated substrate at constant temperature, higher than the saturation temperature at a given vapor pressure. In this manner, an interfacial temperature distribution arises along the drop free surface in terms evaporative mass flux and vapor recoil, which repercussion over aerosol size is studied by determining the asymptotic evolution equation governing the acoustically driven free surface. At such scenario, the connection between surface tension and temperature is also considered; thus, thermocapillary flow is incorporated into our drop model, described in terms of fundamental parameters, like the evaporation number, Marangoni number, and acoustic capillary number. Numerical solution of the evolution equation led us to obtain a simplified representation of the drop interfacial deformation mechanism, capable of predicting atomization and portraying the influence of evaporation over atomization. Subsequent analysis shows that the incorporation of evaporation at SAW atomization traduces in normal stresses counteracting the acoustic and thermocapillary effect, leading to the development of smaller drop aspect ratios with respect to the no-evaporative case. Being aware that the aerosol size is deeply related to the aspect ratio, we propose an analytical expression to estimate aerosol diameter under evaporative conditions. The results show that aspect ratio reduction leads to a decrement on aerosol size, up to two orders of magnitude, with respect to the no-evaporative case. Our study is a first approach providing insight about the importance of evaporation on aerosol regulation at SAW atomization.

Breakdown of laminar regime in flat plate boundary layer seen in numerical solutions of the Navier–Stokes equations obtained with 16th-order scheme

Numerical solutions of the Navier–Stokes equations describing the onset and the development of the boundary layer instability of a subsonic flow about a finite width flat plate are described. They were obtained via exciters-free calculations with the 16th-order multioperators-based scheme optimized to resolve very small solutions scales. Some relevant details of the scheme are briefly outlined. The occurrence of the Tollmien-Schlichting waves near the leading edge, their downstream propagation and breakdowns some distance away are displayed. The role of the oblique waves seen in the solutions is emphasized. The nonlinear stage in the form of structured vortices followed by the solutions randomization are described.

Unconditional superconvergence analysis of a new energy stable nonconforming BDF2 mixed finite element method for BBM-Burgers equation

The focus of this paper is to establish a new energy stable 2-step backword differentiation formula (BDF2) fully-discrete mixed finite element method (MFEM) for the Benjamin–Bona–Mahony-Burgers (BBM-Burgers) equation with the nonconforming rectangular EQ1rot element and zero-order Raviart–Thomas (R-T) element (EQ1rot/Q10×Q01). Based on the energy stable property, the existence and uniqueness of the numerical solution are proved by the Brouwer fixed point theorem. Subsequently, with the assistance of some typical properties of this element pair and the interpolation post-processing approach, the unconditional superclose and superconvergence results with O(h2+τ2) are obtained directly without any constraints between the spatial partition parameter h and the time step τ. It is worthy to mention that the method and analysis presented herein are very different from the time-space splitting approach utilized in the previous studies and simplify the implement. Finally, two numerical experiments are executed to validate the theoretical analysis.

Application of the Cubic Exponential B-spline Collocation Operator Splitting Method for Numerical solutions of Burgers Equation

In this study, the cubic exponential B-spline collocation method has been proposed for the numerical solutions of the Burgers equation with the operator splitting. To apply the operator splitting method, the Burgers' equation has decomposed into two sub-equations based on the time term: the linear part (diffusion) and the nonlinear part (convection). Subsequently, for each sub-equation, Crank-Nicolson finite difference schemes in the temporal direction and cubic exponential B-spline functions and their derivatives, have applied at the x_m nodal points in the spatial direction. The algebraic equation systems obtained have been solved numerically using the Lie-Trotter and Strang splitting schemes to get the solutions of the main equation. Some advantages of the splitting methods include preserving the physical characteristics of the solution, yielding more convergent results over long time intervals, enabling simpler algorithms, and facilitating the storage of solution vectors on computer. To assess the accuracy of the computed numerical results the L_2 and L_∞ error norms have been used. Additionally, the obtained results have been compared with some studies in the literature. The stability analysis of the applied method has been investigated using the von Neumann Fourier series method.

Sinc method in spectrum completion and inverse Sturm–Liouville problems

Cardinal series representations for solutions of the Sturm–Liouville equation with a complex‐valued potential are obtained, by using the corresponding transmutation operator. Consequently, partial sums of the series approximate the solutions uniformly with respect to in any strip of the complex plane. This property of the obtained series representations leads to their applications in a variety of spectral problems. In particular, we show their applicability to the spectrum completion problem, consisting in computing large sets of the eigenvalues from a reduced finite set of known eigenvalues, without any information on the potential as well as on the constants from boundary conditions. Among other applications this leads to an efficient numerical method for computing a Weyl function from two finite sets of the eigenvalues. This possibility is explored in the present work and illustrated by numerical tests. Finally, based on the cardinal series representations obtained, we develop a method for the numerical solution of the inverse two‐spectra Sturm–Liouville problem and show its numerical efficiency.

An easy-to-use analytical model for standing column wells operating with bleed

Accurate assessment of subsurface heat transport is vital for the design of standing column well systems. When bleed is utilized, the mathematical problem is governed by coupled heat transfer and groundwater flow within a radially convergent flow field, making the development of models challenging. While numerical models exist to simulate these transient processes, the absence of simple and accessible analytical solutions limits their broader application. This study addresses this gap by developing a novel easy-to-use analytical model to accurately represent the heat advection-diffusion problem in standing column wells operating with bleed. The proposed model combines the well-known infinite line source model with an innovative scaling function, inspired from the field of solute transport, through a simple convolution product. Notably, the developed model depends on only three dimensionless parameters: dimensionless advection time, Péclet number, and bleed ratio. Rigorous validation against two distinct sets of reference numerical solutions demonstrated the model’s efficiency, accuracy, and reliability across a broad spectrum of nine physical parameters. Key results include relative root mean square errors on the order of 10−3 across 200 reference solutions, confirming the model’s robustness. These findings highlight the model’s potential to significantly advance both research and practical applications in the design and optimization of standing column wells.

Influence of space conjugate temperature varying non-uniform heat sink/source on hydromagnetic slip water-EG (50:50) nanofluid

Significant cooling is an essential requirement in the field of aeronautical and bio-medical engineering, nuclear reactor system, solar collectors, and in development of electronic chips etc. Involving rotating conical geometries. In view of this, flow and heat transfer over conical geometries subject to different constraints of motion are highly needed. Implementation of magnetic field subject to flow pattern controls the fluid motion thereby imparting better cooling. Consideration of nanofluid instead of regular fluid yields prominent cooling of the associated surface. Such relevance has motivated the authors to work on magnetohydrodynamics and heat transfer investigation of water-EG (50:50) mixture based Al2O3 and Fe3O4 past heated and rotating down-pointing upright cone subject to impact of space and temperature varying non-uniform heat source or sink. Numerical solution of dimensionless governing equations is accomplished by implementing Runge-Kutta method. The findings indicate that swirl and axial velocities peter out with rise in magnetic parameter while both exhibit opposite impact in response to slip parameter. Temperature profiles upgrade due to amplification of space and temperature dependent parameters. Skin friction and heat transportation upsurge with growth of solid volume fraction of nanoparticle.

A Second-Order Numerical Method for a Class of Optimal Control Problems

The numerical solution of optimal control problems through second-order methods is examined in this paper. Controlled processes are described by a system of nonlinear ordinary differential equations. There are two specific characteristics of the class of control actions used. The first one is that controls are searched for in a given class of functions, which depend on unknown parameters to be found by minimizing an objective functional. The parameter values, in general, may be different at different time intervals. The second feature of the considered problem is that the boundaries of time intervals are also optimized with fixed values of the parameters of the control actions in each of the intervals. The special cases of the problem under study are relay control problems with optimized switching moments. In this work, formulas for the gradient and the Hessian matrix of the objective functional with respect to the optimized parameters are obtained. For this, the technique of fast differentiation is used. A comparison of numerical experiment results obtained with the use of first- and second-order optimization methods is presented.

Exploring the effects of finite size and indenter shape on the contact behavior of functionally graded thermoelectric materials

The performance enhancement of functionally graded thermoelectric (FGTE) devices is significantly influenced by contact studies of the FGTE materials. It is unclear how the finite thickness and the punch geometry influence the FGTE materials’ contact behaviors. This paper investigates the frictionless contact problem between three types of rigid punches (flat, triangular, and cylindrical) and the FGTE strip with finite thickness. The electric-thermo-elastic parameters of the FGTE strip vary in the thickness direction according to an exponential function. Based on the Fourier integral transform and the transformation matrix method, the problem is transformed into the numerical solution of three sets of singular integral equations. The presence of singular features on either side of the punch demands the adoption of specific collocation strategies. The distribution of the normal current density, the normal energy flux, and the normal contact stress is obtained by adjusting multiple electric-thermo-elastic parameters. The contact stresses in the case of punches with varying shapes can be effectively controlled by modulating the coefficient of thermal expansion and the strip thickness, whereas the effect of the electrical conductivity, the shear modulus, and the thermoelectric load on these stresses depends on whether they are increased or decreased

Lattice Boltzmann simulations of unsteady Bingham fluid flows

Transient flows of viscoplastic fluids have very peculiar characteristics. The startup and cessation flows of viscoplastic materials have been subject to many theoretical and numerical investigations. The most challenging aspect of numerical solutions of viscoplastic fluids is the viscosity singularity during the transition from yielded to unyielded material. Hence, the proper representation of yield surfaces is the most critical aspect of numerical methods in viscoplastic fluid flow. In the present work, we use a lattice Boltzmann scheme to solve an ideal Bingham fluid’s startup and cessation flows. This numerical scheme advantage is that can represent infinite viscosity without noticeable numerical instabilities, producing yield surfaces with more accuracy and quality. Theoretical solutions for the startup flow are available in the literature. However, it is unclear which is more accurate and what their validity ranges are. Nonetheless, these solutions served as a reference for the present simulations. The overall aspect of the numerical solutions agreed with the theoretical models. The cessation flow of the Bingham fluid was also simulated. Unlike a Newtonian fluid, this type of flow is known to have a finite period until cessation. The simulations correctly reproduced this behavior. The transient yield surfaces matched very well with augmented Lagrangian solutions.

Variable density and heat generation impact on chemically reactive Carreau nanofluid heat-mass transfer over stretching sheet with convective heat condition

The present study focuses on the physical significance of heat generation and chemical reaction on Carreau nanofluid with convective heat conditions. Heat transfer is characterized using convective boundary conditions. The governing partial differential equations (PDEs) are transformed into ordinary differential equations (ODEs) by using well define stream functions and similarity transformations. Using a shooting methodology, the Keller-box method with Newton Raphson scheme is used to elaborate the numerical solutions of physical phenomena. Utilizing a similar technique to find the impact of physical parameter such as the production of heat , the rate of reaction , Biot numbers , Brownian motion variable , the thermophoresis parameters , the Weissenberg quantity , Prandtl number , and Lewis number on velocity profile, temperature profile and mass transmission profile are determined graphically. The skin-friction coefficient , local Nusselt , and Sherwood numbers are analyzed numerically. Increment in fluid velocity and slip temperature are depicted with high Biot number. Maximum magnitude of fluid temperature and fluid concentration function are depicted at high value of temperature dependent density. The magnitude of heat and mass transportation enhanced with maximum choice of Brownian motion.

A novel binary data classification algorithm based on the modified reaction-diffusion predator-prey system with Holling-II function.

In this study, we propose a modified reaction-diffusion prey-predator model with a Holling-II function for binary data classification. In the model, we use u and v to represent the densities of prey and predators, respectively. We modify the original equation by substituting the term v with f-v to obtain a stable and clear nonlinear decision surface. By employing a finite difference method for numerical solution of the original model, we conduct various experiments in two-dimensional and three-dimensional spaces to validate the feasibility of the classifier. Additionally, with consideration for wide real applications, we perform classification experiments on electroencephalogram signals, demonstrating the effectiveness and robustness of the classifier in binary data classification.

A Variant of Accelerated Ramadan Group Adomian Decomposition Method for Numerical Solution of Fractional Riccati Differential Equations

A Variant of Accelerated Ramadan Group Adomian Decomposition Method for Numerical Solution of Fractional Riccati Differential Equations

Utilizing numerical techniques for modeling radiating Casson nanofluid flow with thermophoretic phenomenon on a heated stretching surface

This work aims to use mathematical modeling to investigate the mechanics of heat and mass transport in dissipative Casson nanofluid flows over a linear rough sheet. This study considers various elements such as thermal radiation, magnetic fields, heat generation, the varied properties of porous media, and the thermophoretic impact. Besides that, it looks into the variations in viscosity, diffusivity, and thermal conductivity, as well as the energy dissipation from viscous internal friction and fluid temperature modifications. The method involves coming up with and changing the boundary layer equations into a group of linked nonlinear ordinary differential equations that use variables that do not have any dimensions. The foundation of the solution strategy for these equations is the Hermite collocation method (HCM), which is renowned for its precision and adaptability. It offers an organized approach to solving the complex differential equations, allowing for accurate numerical solutions. The use of graphical representations ensures thorough data analysis and clarity, while also providing insightful information about the computed outcomes. The code validation method uses numerical comparisons with recent research to confirm the algorithm’s accuracy and dependability, as well as its resilience in comparison to the existing literature. Important conclusions from the study show that thermal boundary layers and nanofluid velocity decrease with increases in the porosity parameter, slip parameter, and magnetic field parameter.

Modeling Brittle-to-Ductile Transitions in Rock Masses: Integrating the Geological Strength Index with the Hoek–Brown Criterion

Many studies focus on brittle–ductile transition stress in intact rocks; however, in real life, we deal with rock mass which contains many discontinuities. To fill this gap, this research focuses on the brittle–ductile transition stress of rock mass by considering the influence of different Geological Strength Index (GSI) values on the brittle–ductile transition stress of rock mass. In other words, the Hoek–Brown failure criteria for rock mass were reformulated mathematically including the ductility parameter (d), which is defined as the ratio of differential stress to minor stress. Then, the results were analyzed and plotted between σ3*σc and GSI, considering different (d) and Hoek–Brown material constant (mi) values. The brittle–ductile transition stress, σ3*, was determined by intersecting the Hoek–Brown failure envelope with Mogi’s line, with ductility parameters d ranging from 3.4 (silicate rocks) to 5.0 (carbonate rocks). Numerical solutions were derived for σ3*σc as a function of GSI using Matlab, and the results were fitted with an exponential model. The analysis revealed an exponential relationship between σ3*σc and GSI for values above 32, with accuracy better than 3%. Increased ductility reduces rock mass strength, with higher d values leading to lower σ3*σc. The diminishing returns in confinement strength at higher GSI values suggest that rock masses with higher GSI can sustain more confinement but with reduced effectiveness as GSI increases. These findings provide a framework for predicting brittle–ductile transitions in rock engineering.

Hybrid finite element and laplace transform method for efficient numerical solutions of fractional PDEs on graphics processing units

Fractional Partial Differential equations (FPDEs) are essential for modeling complex systems across various scientific and engineering areas. However, efficiently solving FPDEs presents significant computational challenges due to their inherent memory effects, often leading to increased execution times for numerical solutions. This study proposes a highly parallelizable hybrid computational approach that combines the Finite Element Method (FEM) for spatial discretization with Numerical Inversion of the Laplace Transform (NILT) for time-domain solutions, optimized for execution on Graphics Processing Units (GPUs). The NILT method’s high parallelizability, stemming from the independence of its series terms, combined with the robust spatial discretization provided by FEM, enables the efficient and accurate solution of FPDEs on GPUs, demonstrating substantial performance improvements over traditional CPU-based implementations. We observe a generalized pattern in execution time behavior that accounts for both the number of nodes and the number of NILT terms. Specifically, execution time scales quadratically with the number of nodes, while showing only a logarithmic marginal increase with the number of NILT terms These behaviors not only enables consistent performance assessment but also highlights potential areas for algorithm optimization. Validation against exact solutions of fractional diffusion and wave equations, employing Caputo, modified Caputo-Fabrizio, and modified Atangana-Baleanu derivatives, demonstrates the accuracy and convergence of the hybrid FEM-NILT method. Notably, the exact solutions of wave equation are novel in literature. The results highlight the method’s potential for enabling high-precision, large-scale simulations in fractional calculus applications, thereby advancing computational capabilities and efficiency in the field.

Numerical Study of the Reaction Diffusion Prey–Predator Model Having Holling II Increasing Function in the Predator Under Noisy Environment

The current study deals with the stochastic diffusive prey–predator model with Holling II increasing function in the predator. As it is a population dynamics its population can be affected by various factors such as disease in prey–predator species, environmental fluctuation, etc. So, the underlying model is considered under the impact of a temporal multiplicative noise. The numerical solution of the governing model is obtained by finite difference schemes. To check the accuracy of the numerical solution, the stability and consistency of the schemes in the mean square sense are derived. All the schemes are consistent with the model and von Neumann stable. Several simulations are drawn to check the efficacy of the scheme. A test problem is investigated with numerical schemes and for the different values of the parameters, 3D and 2D plots are drawn. A graphical representation of the given system is drawn which shows the behavior of both species and the efficiency of the novel schemes. Hope so, these results will motivate the researcher to consider the stochastic prey–predator model for a better understanding of the ecosystem.

Fluctuation-driven dynamics of liquid nano-threads with external hydrodynamic perturbations

Instability and rupture dynamics of a liquid nano-thread, subjected to external hydrodynamic perturbations, are captured by a stochastic lubrication equation (SLE) incorporating thermal fluctuations via Gaussian white noise. Linear instability analysis of the SLE is conducted to derive the spectra and distribution functions of thermal capillary waves influenced by external perturbations and thermal fluctuations. The SLE is also solved numerically using a second-order finite difference method with a correlated noise model. Both theoretical and numerical solutions, validated through molecular dynamics, indicate that surface tension forces due to specific external perturbations overcome the random effects of thermal fluctuations, determining both the thermal capillary waves and the evolution of perturbation growth. The results also show two distinct regimes: (i) the hydrodynamic regime, where external perturbations dominate, leading to uniform ruptures, and (ii) the thermal-fluctuation regime, where external perturbations are surpassed by thermal fluctuations, resulting in non-uniform ruptures. The transition between these regimes, modelled by a criterion developed from linear instability theory, exhibits a strong dependence on the amplitudes and wavenumbers of the external perturbations.