Flow Problems Research Articles

The influence of heterogeneous catalytic processes on the heat flux to the surface and the chemical composition of the shock layer at high-speed flow around blunt bodies

Numerical modeling of the high-speed flow of a multicomponent dissociated gas around a blunt body was carried out for the flow regime corresponding to the thermally stressed point of the gliding trajectory of entry into the Earth's atmosphere in the surface temperature range of 800–2000 K. As boundary conditions, a stage-by-stage model of heterogeneous catalysis of the interaction of dissociated air with the thermal protection material SiO2 was used, previously developed by the authors on the basis of transition state theory and quantum mechanical calculations. The dependence of the modes of heterogeneous chemical reactions on the surface temperature and the density of surface adsorption centers and their influence on the chemical composition of the shock layer and the convective heat flux to the streamlined surface are analyzed. The contributions of individual mechanisms of heterogeneous recombination to the total rate of formation of molecules on the surface are estimated. The possibility of using the density of adsorption sites as an effective parameter when setting boundary conditions on a real thermal protection surface in flow problems is shown.

A hybrid algorithm for quadratically constrained quadratic optimization problems

Quadratically Constrained Quadratic Programs (QCQPs) are an important class of optimization problems with diverse real-world applications. In this work, we propose a variational quantum algorithm for general QCQPs. By encoding the variables in the amplitude of a quantum state, the requirement for the qubit number scales logarithmically with the dimension of the variables, which makes our algorithm suitable for current quantum devices. Using the primal-dual interior-point method in classical optimization, we can deal with general quadratic constraints. Our numerical experiments on typical QCQP problems, including Max-Cut and optimal power flow problems, demonstrate better performance of our hybrid algorithm over classical counterparts.

Optimal power dispatch in hybrid power system for medium- and large-scale practical power systems using self-adaptive bonobo optimizer algorithm

Incorporating renewable energy sources (RESs) introduces a notable amount of uncertainty in the optimal planning and operation of electrical power grids. Under these circumstances, this paper proposes the application of a recently introduced metaheuristic optimization technique to solve the stochastic optimal power flow (OPF) problem involving wind and solar power sources. The self-adaptive bonobo optimizer (SaBO) is used to minimize three distinct objective functions: (i) Total generation cost (TGC) minimization, including both thermal and wind/solar generation costs, (ii) Power loss minimization, (iii) Combined generation cost and emissions effect minimization. The costs associated with the stochastic generation of wind and solar power included direct costs, reserves and penalty costs from the overestimation and underestimation of available wind and solar power, respectively. The performance of the proposed algorithm is evaluated on two power systems: the modified IEEE 30-bus and the Algerian DZA 114-bus test systems. To demonstrate the efficacy of the SaBO, the obtained results have been compared with those obtained from the Kepler optimization algorithm (KOA) and other recently published optimizers under the same case studies and constraints. The comparative results clearly show the superiority of the SaBO algorithm over all other well-known optimization algorithms provided in the literature for solving the OPF problem. This is evidenced by minimizing total generation costs of 781.2363 $/h for the modified IEEE 30-bus and 16,706.1630 $/h for the Algerian DZA-114-bus system. Furthermore, the integration of RES led to a notable 2.33% and 11.67% reduction in total generation cost for the IEEE 30-bus and Algerian DZA 114-bus systems, respectively, compared to their initial configurations without RESs. The promising findings highlight the powerful of the optimizer to solve non-linear and complex optimization problems in power systems.

Well-Balanced conservative central upwind scheme for solving the dam-break flow problem over erodible bed

This work deals with the numerical solution of dam-break flow over an erodible bed. The mathematical model is a combination of the shallow water, the transport diffusion and the bed morphology change equations. The system is solved by a well-Balanced central upwind scheme with conservative property. Several tests are illustrated in order to validate the accuracy and the performance of the model. A comparison of central upwind scheme and Roe scheme is presented.

Impact of electro-osmotic, activation energy and chemical reaction on Sisko fluid over Darcy–Forchheimer porous stretching cylinder

This study provides numerical solution for two-dimensional electro-osmotically motivated electro-magneto-hydrodynamic Sisko fluid through an elongating porous cylinder. The electro-osmotic, activation energy with chemical reaction are also interpreted for this flow model. The numerical equations that govern this flow problem are renewed into dimensionless form by applying appropriate transformations exploiting Runge–Kutta–Fehlberg fifth order through shooting technique method. The after-effects are illustrated as graphs for the concerned physical parameters and their impact on convection and conduction is also studied. The velocity rises for electro-osmotic parameter whereas drops for magnetic and porous parameter. The temperature shows an increasing profile for the curvature, electric and thermal conductivity parameter while decreases for Prandtl number. The concentration of nanoparticles in the fluid boosts for activation energy but deflates for curvature parameter. The present findings appear to be in good accord when compared to earlier published studies. Numerous opportunities and applications are presented by applying electro-osmotic forces to non-Newtonian fluid flow, especially in the fields of nanotechnology, electro-kinetics and micro-fluids. The current study can be applied to design effective electro-magnetic devices, particularly in certain thermal transport characteristic regime. The main findings demonstrate the great utility of electro-osmosis in micro-fluidic devices, chemical analysis, soil analysis and cement slurries for managing flow and heat transmission.

Chaotic-quasi-oppositional-phasor based multi populations gorilla troop optimizer for optimal power flow solution

This paper proposes an efficient modified Gorilla Troops Optimizer (MGTO) to address the problem of prone to local optima and stagnation possibility in the standard Gorilla Troops Optimizer (GTO). The proposed MGTO is based on multiple strategies including chaotic initialization, quasi-oppositional based learning, phasor motion, and multi-populations strategy. The proposed MGTO is applied to resolve the non-convex single objectives optimal power flow (OPF) problem to reduce the fuel cost without and with vple, the emissions, the fuel cost considering prohibted zones (POZs) and piecewise functions, while for multiple objective cases the fuel cost is optimized with the power losses, the voltage deviation and the emissions. To validate the effectiveness of the proposed optimizer, it is tested on IEEE 30-bus and 57-bus systems, as well as 23 standard benchmark functions and the obtained outcomes are compared with the other seven well-known optimization techniques such as SCSO, GWO, WOA, SCA, PF, BWO and GTO. The simulation results verified that the proposed MGTO provides a remarkable and robust solution for the OPF problem solution compared to the other state-of-the-arts algorithms.

Experimental and numerical investigation of shelf flow crossing over a strait

AbstractMotivated by the phenomenon of Scotian Shelf Crossover events, the problem of a shelf flow that is interrupted by a strait is considered. Laboratory experiments in a rotating tank with barotropic and baroclinic flow over flat and sloping shelves confirm that the flow is steered by the bathymetric contours and mainly circumnavigates the gulf. In order to jump across the strait, as suggested by earlier theories, the flow must have unrealistically high Rossby numbers. However, the near bottom friction relaxes the bathymetric constraint and causes the formation of a peculiar jet crossing the strait diagonally. For the dissipation values such that a half of the transport goes around the gulf and half crosses the strait diagonally, the diagonal crossover jet becomes most evident. Numerical solutions for realistic values of the frictional parameter reproduce the results of the laboratory experiments and consideration of the actual Gulf of Maine bathymetry reproduces patterns similar to those observed by drift trajectories and in the satellite derived sea surface temperature fields.

The Maximum Attainable Flow and Minimal Cost Problem in a Network

The Maximum Attainable Flow and Minimal Cost Problem in a Network

Surface tension and wetting at free surfaces in smoothed particle hydrodynamics

Surface tension and wetting are dominating physical effects in microscale and nanoscale flows. We present an efficient and reliable model of surface tension and equilibrium contact angles in smoothed particle hydrodynamics for free-surface problems. We demonstrate its robustness and accuracy by simulating several three-dimensional free-surface flow problems driven by interfacial tension.

A systematic analysis of three-dimensional Riemann problems for verification of compressible-flow solvers

Numerical simulation is a well-established way of analyzing compressible flows. Due to high computational demands of solvers for such flow problems, their verification is typically limited to two-dimensional (2D) cases. However, 2D simulations suppress fundamental three-dimensional (3D) aspects of the flow evolution in (compressible) turbulent flows or shock-bubble and shock-drop interactions. With increase of computational power, 3D simulations become more feasible for routine analyses. The verification of 3D simulation frameworks is often limited to transformations of lower dimensional test cases in the 3D space. There is a lack of strictly 3D reference test cases for gas dynamics. In this work, we present a set of genuine 3D Riemann problems in order to validate and verify numerical solvers for compressible flows. The problems are designed such that each octant’s constant initial data connects two neighboring states by an elementary wave only. The problem design is inspired by well-established 2D Riemann problems most prominently posted by Lax and Liu (1998). In contrast to the twenty published 2D cases, more than 300 distinct combinations can be found in 3D. We provide example solutions for the particularly interesting ones of these case combinations and show how the cases help to expose shortcomings of numerical solvers. We provide reference data from computations with an open-source compressible multiresolution flow solver. For the reference solutions, we employ the Harten-Lax-van Leer contact (HLLC) Riemann solver and a weighted essentially non-oscillatory (WENO) reconstruction stencil of fifth order. The reference solutions use an effective resolution of one billion cells. We additionally make the full compute pipeline of this work publicly available, so interested researchers may reproduce and extend the current work.

Implementation and Validation of an Original OpenFOAM Code for Fluid–Structure Interaction Problems in Compressible Flow

Implementation and Validation of an Original OpenFOAM Code for Fluid–Structure Interaction Problems in Compressible Flow

Robust reverse 1-center problems on trees with interval costs

This paper concerns the problem of reducing the edge lengths of a tree in a minimum cost way such that the 1-center function at a predetermined vertex is at most a given threshold value. Moreover, it is assumed that the costs of reducing edge lengths are not exactly determined, but they are estimated within intervals and their total deviation is limited within a given bound. This problem is called the robust reverse 1-center problem on trees with interval costs. We first reformulate the problem as a Stackelberg game and translate it into a univariate optimization problem. Based on the special properties of the feasible set, we can state the corresponding problem as a parameterized minimum cost flow problem. Then we develop an algorithm that solves the problem in O ( n 2 log ⁡ n ) time due to the convexity of the univariate objective function. We finally discuss the problem of improving the international transportation network in Iran as a real-life application.

Optimal convergence analysis of two RPC-SAV schemes for the unsteady incompressible magnetohydrodynamics equations

Abstract In this paper, we present and analyze two linear and fully decoupled schemes for solving the unsteady incompressible magnetohydrodynamics equations. The rotational pressure-correction (RPC) approach is adopted to decouple the system, and the recently developed scalar auxiliary variable (SAV) method is used to treat the nonlinear terms explicitly and keep energy stability. One is the first-order RPC-SAV-Euler and the other one is generalized Crank–Nicolson-type scheme: GRPC-SAV-CN. For the RPC-SAV-Euler scheme, both unconditionally energy stability and optimal convergence are derived. The new GRPC-SAV-CN is constructed and can be regarded as a parameterized scheme, which includes PC-SAV-CN when the parameter $\beta =0$ and RPC-SAV-CN when $\beta \in (0,\frac {1}{2}]$; see Algorithm 3.2. However, Jiang and Yang (Jiang, N. & Yang, H. (2021) SIAM J. Sci. Comput., 43, A2869–A2896) point out that the SAV method has low accuracy by several commonly tested benchmark flow problem when solving Navier–Stokes equations. To improve the accuracy, we added two stabilization $-\alpha _{1}\varDelta t\nu \varDelta (\widetilde {\textbf {u}}^{n+1}-{\textbf {u}}^{n})$ and $\alpha _{2}\varDelta t\sigma ^{-1}\mbox {curl}\mbox {curl} (\textbf {H}^{n+1}-\textbf {H}^{n})$ in the GRPC-SAV-CN scheme, which play decisive roles in giving optimal error estimates. The unconditionally energy stability of the proposed scheme is given. We prove that the PC-SAV-CN scheme has second-order convergence speed, and the RPC-SAV-CN one has 1.5-order convergence rate. Finally, some numerical examples are presented to verify the validity and convergence of the numerical schemes.

A comparison study of spatial and temporal schemes for flow and transport problems in fractured media with large parameter contrasts on small length scales

AbstractIn this work, various high-accuracy numerical schemes for transport problems in fractured media are further developed and compared. Specifically, to capture sharp gradients and abrupt changes in time, schemes with low order of accuracy are not always sufficient. To this end, discontinuous Galerkin up to order two, Streamline Upwind Petrov-Galerkin, and finite differences, are formulated. The resulting schemes are solved with sparse direct numerical solvers. Moreover, time discontinuous Galerkin methods of order one and two are solved monolithically and in a decoupled fashion, respectively, employing finite elements in space on locally refined meshes. Our algorithmic developments are substantiated with one regular fracture network and several further configurations in fractured media with large parameter contrasts on small length scales. Therein, the evaluation of the numerical schemes and implementations focuses on three key aspects, namely accuracy, monotonicity, and computational costs.

Graph‐based solution for smart grid real‐time operation and control

AbstractUnder the envisioned smart grid paradigm, there is an increasing demand for a fast, accurate, and efficient power flow solution for distribution system operation and control. Various solution techniques have been proposed, each with its own unique formulation, solution methodology, advantages, and drawbacks. Motivated by challenges associated with the integration of renewable distributed energy resources and electric vehicles into distribution systems and further by the speed and convergence limitations of existing tools, this paper presents a novel graph‐based power flow solution for smart grid's real‐time operation and control, named Flow‐AugmentationPF algorithm. The proposed method formulates a power flow problem as a network‐flow problem and solves it by using a maximum‐flow algorithm, inspired by the push‐relabel max‐flow technique. The performance of the proposed algorithm is tested and validated using several benchmark networks of different sizes, topologies, and parameters and compared against the most commonly used solution techniques and commercial software packages, namely PSS/E and PSCAD. The proposed formulation is simple, accurate, fast, yet computationally efficient, as it is based on matrix‐vector multiplication, and is also scalable, considering the formulation works as a graph‐based method, which, inherently, allows for parallel computation for added computational speed.

Numerical approximations of groundwater flow problem using fractional variational iteration method with fractional derivative of singular and nonsingular kernels

This study aims to present the numerical approximations of groundwater flow problem, namely the time-fractional Boussinesq equation using the fractional variational iteration method. Fractional derivatives in time variable are taken in Caputo and Caputo–Fabrizio sense. This work compares the solutions of the time-fractional Boussinesq equation for the Caputo and Caputo–Fabrizio fractional operators. Graphical representation of the water table head’s behavior for different time values is shown.

Analytical and numerical study of water-based silver nanofluid (Ag) across a Riga plate with nonlinear radiation and viscous dissipation: A three-dimensional study

The primary endeavor of this investigation is to determine the Darcy–Forchheimer flow of water-based silver nanofluid (Ag) across a bidirectional Riga plate with viscous dissipation and nonlinear radiation. By manipulating the pertinent variables, the governing equations undergo a transformation into nonlinear ordinary differential equations. The remodeled flow problems are investigated numerically and analytically using the MATLAB (bvp4c built-in function) algorithm and the homotopy analysis method scheme. The repercussions of the appropriate factors on the x− and y− direction velocity profiles, temperature profile, skin friction coefficient of the x− and y− directions, and local Nusselt number are examined via tables and graphs. It is found that both directional fluid velocities decline in the presence of the suction/injection parameter. Increasing the radiation parameter and the Eckert number exalts the fluid temperature. The surface shear stress decays as the suction/injection parameter and the Forchheimer number grow larger. In addition, a greater heat transfer rate appears on the Riga plate compared to the stationary plate.

Second-order, positive, and unconditional energy dissipative scheme for modified Poisson–Nernst–Planck equations

First-order energy dissipative schemes in time are available in literature for the Poisson–Nernst–Planck (PNP) equations, but second-order ones are still in lack. This work proposes novel second-order discretization in time and finite volume discretization in space for modified PNP equations that incorporate effects arising from ionic steric interactions and dielectric inhomogeneity. A multislope method on unstructured meshes is proposed to reconstruct positive, accurate approximations of mobilities on faces of control volumes. Numerical analysis proves that the proposed numerical schemes are able to unconditionally ensure the existence of positive numerical solutions, original energy dissipation, mass conservation, and preservation of steady states at discrete level. Extensive numerical simulations are conducted to demonstrate numerical accuracy and performance in preserving properties of physical significance. Applications in ion permeation through a 3D nanopore show that the modified PNP model, equipped with the proposed schemes, has promising applications in the investigation of ion selectivity and rectification. The proposed second-order discretization can be extended to design temporal second-order schemes with original energy dissipation for a type of gradient flow problems with entropy.

Simulation of hybrid boiling nano fluid flow with convective boundary conditions through a porous stretching sheet through Levenberg Marquardt artificial neural networks approach

In the domain of artificial neural networks, the Levenberg-Marquardt technique stands out for its innovative approach and convergent stability. It generates a numerical approach for the progression of a hybrid Cu–Al2CO3/water nanofluid on a porous stretched sheet (HNF-PSS), employing regression plots, state transition measures, histogram representations, and mean squared errors. This work explores the analysis of the fluid flow problem associated with HNF-PSS, introducing a novel application of an intelligent computing system that utilizes the effectiveness of neural networks trained by the Levenberg-Marquardt algorithm (NN-TLMA). The initial mathematical expression, expressed in terms of partial differential equations (PDEs) for HNF-PSS, undergoes transformation into dimensionless nonlinear ordinary differential equations (ODEs). The collection of data for the envisioned NN-TLMA involves parameters influencing the velocity of the HNF-PSS system model, generated using the Lobatto IIIa formula. The efficacy of the suggested NN-TLMA for solving the HNF-PSS is robustly verified by examinations of state transition dynamics, regression analyses, mean square error, and error histogram studies. The consistent alignment of the obtained results with corresponding solutions attests to the legitimacy of the structure, achieving a high level of accuracy within the range of 10–6 to 10–8. Insightful findings demonstrate that the thermal Biot number contributes significantly to an upsurge in fluid temperature. Furthermore, an augmentation in nanofluid concentration and thermal profile is found to increase the porosity strength of the medium. Nevertheless, an attenuation in the velocity profile is noted under these conditions. The application of this simulation work can be found in optimizing the designs of heat exchanger, in biomedical application and in renewable energy processes.

Advanced Computational Framework to Analyze the Stability of Non-Newtonian Fluid Flow through a Wedge with Non-Linear Thermal Radiation and Chemical Reactions

The main idea of this investigation is to introduce an integrated intelligence approach that investigates the chemically reacting flow of non-Newtonian fluid with a backpropagation neural network (LMS-BPNN). The AI-based LMS-BPNN approach is utilized to obtain the optimal solution of an MHD flow of Eyring–Powell over a porous shrinking wedge with a heat source and nonlinear thermal radiation (Rd). The partial differential equations (PDEs) that define flow problems are transformed into a system of ordinary differential equations (ODEs) through efficient similarity variables. The reference solution is obtained with the bvp4c function by changing parameters as displayed in Scenarios 1–7. The label data are divided into three portions, i.e., 80% for training, 10% for testing, and 10% for validation. The label data are used to obtain the approximate solution using the activation function in LMS-BPNN within the MATLAB built-in command ‘nftool’. The consistency and uniformity of LMS-BPNN are supported by fitness curves based on the MSE, correlation index (R), regression analysis, and function fit. The best validation performance of LMS-BPNN is obtained at 462, 369, 642, 542, 215, 209, and 286 epochs with MSE values of 8.67 × 10−10, 1.64 × 10−9, 1.03 × 10−9, 302 9.35 × 10−10, 8.56 × 10−10, 1.08 × 10−9, and 6.97 × 10−10, respectively. It is noted that f′(η), θ(η), and ϕ(η) satisfy the boundary conditions asymptotically for Scenarios 1–7 with LMS-BPNN. The dual solutions for flow performance outcomes (Cfx, Nux, and Shx) are investigated with LMS-BPNN. It is concluded that when the magnetohydrodynamics increase (M=0.01, 0.05, 0.1), then the solution bifurcates at different critical values, i.e., λc=−1.06329,−1.097,−1.17694. The stability analysis is conducted using an LMS-BPNN approximation, involving the computation of eigenvalues for the flow problem. The deduction drawn is that the upper (first) branch solution remains stable, while the lower branch solution causes a disturbance in the flow and leads to instability. It is observed that the boundary layer thickness for the lower branch (second) solution is greater than the first solution. A comparison of numerical results and predicted solutions with LMS-BPNN is provided and they are found to be in good agreement.