Search results

This study aims to investigate the dual solutions for axisymmetric flow and heat transfer due to a permeable radially shrinking disk in copper oxide (CuO) and silver (Ag) hybrid nanofluids with radiation effect.

The partial differential equations that governed the problem will undergo a transformation into a set of similarity equations. Following this transformation, a numerical solution will be obtained using the boundary value problem solver, bvp4c, built in the MATLAB software. Later, analysis and discussion are conducted to specifically examine how various physical parameters affect both the flow characteristics and the thermal properties of the hybrid nanofluid.

Dual solutions are discovered to occur for the case of shrinking disk (λ < 0). Stronger suction triggers the critical values’ expansion and delays the boundary layer separation. Through stability analysis, it is determined that one of the solutions is stable, whereas the other solution exhibits instability, over time. Moreover, volume fraction upsurge enhances skin friction and heat transfer in hybrid nanofluid. The hybrid nanofluid’s heat transfer also heightened with the influence of radiation.

Flow over a shrinking disk has received limited research focus, in contrast to the extensively studied axisymmetric flow problem over a diverse set of geometries such as flat surfaces, curved surfaces and cylinder. Hence, this study highlights the axisymmetric flow due to a shrinking disk under radiation influence, using hybrid nanofluids containing CuO and Ag. Upon additional analysis, it is evidently shows that only one of the solutions exhibits stability, making it a physically dependable choice in practical applications. The authors are very confident that the findings of this study are novel, with several practical uses of hybrid nanofluids in modern industry.

The purpose of this paper is to analyze numerically the blowup in finite time of the solutions to a one-dimensional, bidirectional, nonlinear wave model equation for the propagation of small-amplitude waves in shallow water, as a function of the relaxation time, linear and nonlinear drift, power of the nonlinear advection flux, viscosity coefficient, viscous attenuation, and amplitude, smoothness and width of three types of initial conditions.

An implicit, first-order accurate in time, finite difference method valid for semipositive relaxation times has been used to solve the equation in a truncated domain for three different initial conditions, a first-order time derivative initially equal to zero and several constant wave speeds.

The numerical experiments show a very rapid transient from the initial conditions to the formation of a leading propagating wave, whose duration depends strongly on the shape, amplitude and width of the initial data as well as on the coefficients of the bidirectional equation. The blowup times for the triangular conditions have been found to be larger than those for the Gaussian ones, and the latter are larger than those for rectangular conditions, thus indicating that the blowup time decreases as the smoothness of the initial conditions decreases. The blowup time has also been found to decrease as the relaxation time, degree of nonlinearity, linear drift coefficient and amplitude of the initial conditions are increased, and as the width of the initial condition is decreased, but it increases as the viscosity coefficient is increased. No blowup has been observed for relaxation times smaller than one-hundredth, viscosity coefficients larger than ten-thousandths, quadratic and cubic nonlinearities, and initial Gaussian, triangular and rectangular conditions of unity amplitude.

The blowup of a one-dimensional, bidirectional equation that is a model for the propagation of waves in shallow water, longitudinal displacement in homogeneous viscoelastic bars, nerve conduction, nonlinear acoustics and heat transfer in very small devices and/or at very high transfer rates has been determined numerically as a function of the linear and nonlinear drift coefficients, power of the nonlinear drift, viscosity coefficient, viscous attenuation, and amplitude, smoothness and width of the initial conditions for nonzero relaxation times.

We present an approach for pricing American put options with a regime-switching volatility. Our method reveals that the option price can be expressed as the sum of two components: the price of a European put option and the premium associated with the early exercise privilege. Our analysis demonstrates that, under these conditions, the perpetual put option consistently commands a higher price during periods of high volatility compared to those of low volatility. Moreover, we establish that the optimal exercise boundary is lower in high-volatility regimes than in low-volatility regimes. Additionally, we develop an analytical framework to describe American puts with an Erlang-distributed random-time horizon, which allows us to propose a numerical technique for approximating the value of American puts with finite expiry. We also show that a combined approach involving randomization and Richardson extrapolation can be a robust numerical algorithm for estimating American put prices with finite expiry.

Two-dimensional (2D) problems are governed by unsteady anisotropic modified-Helmholtz equation of time–space dependent coefficients are considered. The problems are transformed into a boundary-only integral equation which can be solved numerically using a standard boundary element method (BEM). Some examples are solved to show the validity of the analysis and examine the accuracy of the numerical method.

The 2D problems which are governed by unsteady anisotropic modified-Helmholtz equation of time–space dependent coefficients are solved using a combined BEM and Laplace transform. The time–space dependent coefficient equation is reduced to a time-dependent coefficient equation using an analytical transformation. Then, the time-dependent coefficient equation is Laplace transformed to get a constant coefficient equation, which can be written as a boundary-only integral equation. By utilizing a BEM, this integral equation is solved to find numerical solutions to the problems in the frame of the Laplace transform. These solutions are then inversely transformed numerically to obtain solutions in the original time–space frame.

The main finding of this research is the derivation of a boundary-only integral equation for the solutions of initial-boundary value problems governed by a modified-Helmholtz equation of time–space dependent coefficients for anisotropic functionally graded materials with time-dependent properties.

The originality of the research lies on the time dependency of properties of the functionally graded material under consideration.

Fused Filament Fabrication (FFF) is an extrusion-based manufacturing process using fused thermoplastics. Despite its low cost, the FFF is not extensively used in high-value industrial sectors mainly due to parts' anisotropy (related to the deposition strategy) and residual stresses (caused by successive heating cycles). Thus, this study aims to investigate the process improvement and the optimization of the printed parts.

In this work, a meshless technique – the Radial Point Interpolation Method (RPIM) – is used to numerically simulate the viscoplastic extrusion process – the initial phase of the FFF. Unlike the FEM, in meshless methods, there is no pre-established relationship between the nodes so the nodal mesh will not face mesh distortions and the discretization can easily be modified by adding or removing nodes from the initial nodal mesh. The accuracy of the obtained results highlights the importance of using meshless techniques in this field.

Meshless methods show particular relevance in this topic since the nodes can be distributed to match the layer-by-layer growing condition of the printing process.

Using the flow formulation combined with the heat transfer formulation presented here for the first time within an in-house RPIM code, an algorithm is proposed, implemented and validated for benchmark examples.

This study aims to derive a novel spatial numerical method based on multidimensional local Taylor series representations for solving high-order advection-diffusion (AD) equations.

The parabolic AD equations are reduced to the nonhomogeneous elliptic system of partial differential equations by utilizing the Chebyshev spectral collocation method (ChSCM) in the temporal variable. The implicit-explicit local differential transform method (IELDTM) is constructed over two- and three-dimensional meshes using continuity equations of the neighbor representations with either explicit or implicit forms in related directions. The IELDTM yields an overdetermined or underdetermined system of algebraic equations solved in the least square sense.

The IELDTM has proven to have excellent convergence properties by experimentally illustrating both h-refinement and p-refinement outcomes. A distinctive feature of the IELDTM over the existing numerical techniques is optimizing the local spatial degrees of freedom. It has been proven that the IELDTM provides more accurate results with far fewer degrees of freedom than the finite difference, finite element and spectral methods.

This study shows the derivation, applicability and performance of the IELDTM for solving 2D and 3D advection-diffusion equations. It has been demonstrated that the IELDTM can be a competitive numerical method for addressing high-space dimensional-parabolic partial differential equations (PDEs) arising in various fields of science and engineering. The novel ChSCM-IELDTM hybridization has been proven to have distinct advantages, such as continuous utilization of time integration and optimized formulation of spatial approximations. Furthermore, the novel ChSCM-IELDTM hybridization can be adapted to address various other types of PDEs by modifying the theoretical derivation accordingly.

This paper aims to focus on the finite element method in the frequency domain (FD-FEM) for the transient electric field in the non-sinusoidal steady state under the non-sinusoidal periodic voltage excitation.

Firstly, the boundary value problem of the transient electric field in the frequency domain is described, and the finite element equation of the FD-FEM is derived by Galerkin’s method. Secondly, the constrained electric field equation on the boundary in the frequency domain (FD-CEFEB) is also derived, which can solve the electric field intensity on the boundary and the dielectric interface with high accuracy. Thirdly, the calculation procedures of the FD-FEM with FD-CEFEB are introduced in detail. Finally, a numerical example of the press-packed insulated gate bipolar transistor under the working condition of the repetitive turn-on and turn-off is given.

The FD-CEFEB improves numerical accuracy of electric field intensity on the boundary and interfacial charge density, which can be achieved by modifying the existing FD-FEMs’ code in appropriate steps. Moreover, the proposed FD-FEM and the FD-CEFEB will only increase calculation costs by a little compared with the traditional FD-FEMs.

The FD-CEFEB can directly solve the electric field intensity on the boundary and the dielectric interface with high accuracy. This paper provides a new FD-FEM for the transient electric field in the non-sinusoidal steady state with high accuracy, which is suitable for combined insulation structure with a long time constant.

This study aims to propose and numerically assess different ways of discretising a very weak formulation of the Poisson problem.

We use integration by parts twice to shift smoothness requirements to the test functions, thereby allowing low-regularity data and solutions.

Various conforming discretisations are presented and tested, with numerical results indicating good accuracy and stability in different types of problems.

This is one of the first articles to propose and test concrete discretisations for very weak variational formulations in primal form. The numerical results, which include a problem based on real MRI data, indicate the potential of very weak finite element methods for tackling problems with low regularity.

This paper aims to present dual solutions for the two-dimension copper oxide with silver (CuO–Ag) and zinc oxide with silver (ZnO–Ag) hybrid nanofluid flow past a permeable shrinking sheet in a dusty fluid with velocity slip.

The governing partial differential equations for the two dust particle phases are reduced to the pertinent ordinary differential equations using a similarity transformation. Closed-form analytical solutions for the reduced skin friction and reduced Nusselt number, as well as for the velocity and temperature profiles, were presented, both graphically and in tables, under specific non-dimensional physical parameters such as the suction parameter, Prandtl number, slip parameter and shrinking parameter, which are also presented in both figures and tables.

The results indicate that for the shrinking flow, the wall skin friction is higher in the dusty fluid when compared with the clear (viscous) fluid. In addition, the effect of the fluid–particle interaction parameter to the fluid phase can be seen more clearly in the shrinking flow. Furthermore, multiple (dual, upper and lower branch solutions) are found for the governing similarity equations and the upper branch solution expanded with higher values of the suction parameter. It can be confirmed that the lower branch solution is unstable.

In practice, the study of the stretching/shrinking flow is crucially important and useful. Both the problems of steady and unsteady flow of a dusty fluid have a wide range of possible applications in practice, such as in the centrifugal separation of particles, sedimentation and underground disposal of radioactive waste materials.

Even though the problem of dusty fluid has been broadly investigated, very limited results can be found for a shrinking sheet. Indeed, this paper has succeeded to obtain analytically dual solutions. The stability analysis can be performed by following many published papers on stretching/shrinking sheets. Finally, the critical values and plotting curves for obtaining single or dual solution are successfully presented.

A novel type of heat transfer fluid known as hybrid nanofluids is used to improve the efficiency of heat exchangers. It is observed from literature evidence that hybrid nanofluids outperform single nanofluids in terms of thermal performance. This study aims to address the stagnation point flow induced by Williamson hybrid nanofluids across a vertical plate. This fluid is drenched under the influence of mixed convection in a Darcy–Forchheimer porous medium with heat source/sink and entropy generation.

By applying the proper similarity transformation, the partial differential equations that represent the leading model of the flow problem are reduced to ordinary differential equations. For the boundary value problem of the fourth-order code (bvp4c), a built-in MATLAB finite difference code is used to tackle the flow problem and carry out the dual numerical solutions.

The shear stress decreases, but the rate of heat transfer increases because of their greater influence on the permeability parameter and Weissenberg number for both solutions. The ability of hybrid nanofluids to strengthen heat transfer with the incorporation of a porous medium is demonstrated in this study.

The findings may be highly beneficial in raising the energy efficiency of thermal systems.

The originality of the research lies in the investigation of the Darcy–Forchheimer stagnation point flow of a Williamson hybrid nanofluid across a vertical plate, considering buoyancy forces, which introduces another layer of complexity to the flow problem. This aspect has not been extensively studied before. The results are verified and offer a very favorable balance with the acknowledged papers.