Mathematical modeling of fractional derivatives for magnetohydrodynamic fluid flow between two parallel plates by the radial basis function method

Investigations into the magnetohydrodynamics of viscous fluids have become more important in recent years, owing to their practical significance and numerous applications in astro-physical and geo-physical phenomena. In this paper, the radial base function was utilized to answer fractional equation associated with fluid flow passing through two parallel flat plates with a magnetic field. The magnetohydrodynamics coupled stress fluid flows between two parallel plates, with the bottom plate being stationary and the top plate moving at a persistent velocity. We compared the radial basis function approach to the numerical method (fourth-order Range-Kutta) in order to verify its validity. The findings demonstrated that the discrepancy between these two techniques is quite negligible, indicating that this method is very reliable. The impact of the magnetic field parameter and Reynolds number on the velocity distribution perpendicular to the fluid flow direction is illustrated. Eventually, the velocity parameter is compared for diverse conditions α, Reynolds and position (y), the maximum of which occurs at α = 0.4. Also, the maximum velocity values occur in α=0.4 and Re=1000 and the concavity of the graph is less for α=0.8.     

Methodical investigation of free convection from vertical and horizontal plates

This paper attempts of review the literature on free convection from flat surfaces transferring heat in an unlimited space. Data from the references have been presented as final criterial relations both in tabular form and in graphic system, separately for the vertical plates and the horizontal ones. About 45 results of investigations on free convection have been examined and the following mean criterial relations have been calculated:Nu=0.550·(Ra)^0.250 for vertical plates,Nu=0.310·(Ra)^0.290 for horizontal plates. It has been observed for the vertical plates that divergence of these results relatively to the mean values is constant within the whole variability range of the Rayleigh number and amounts to about ± 15%. In the case of horizontal plates the divergencies are more considerable and vary within ± 50% for the laminar, and within ± 25% for the turbulent range. In the next step the attempt was made to determine the causes of these divergencies and errors. The work was performed because of two main reasons: imperfection of measuring methods and a lack of uniformity in the choice of the characteristic linear dimension.     

Junction of thin plates

The elasticity problem of a thin plate with an edge is considered using asymptotic methods. The small parameter ε describes the relative thickness of the plate. In the case when the elasticity coefficients are everywhere of the same order of magnitude, the asymptotic behaviour of the plate is such that the angle of the edge remains constant under the deformation (the junction is called `rigid'). We also consider a junction mode of a narrow filet (the order of its width is O(ε)) of a `soft' elastic material, the elasticity coefficients being O(ε) with respect to those of the plates. In this case (called `elastic junction'), the asymptotic modelling contains an energy bilinear form associated with the filet which involves the variation of the angle and the sliding along the junction.     

Bifurcation analysis of a tri-neuron neural network model in the frequency domain

In this paper, a class of neural network models with three neurons is considered. By applying the frequency domain approach and analyzing the associated characteristic equation, the existence of the bifurcation parameter point is determined. If the coefficient μ is chosen as a bifurcation parameter, it is found that Hopf bifurcation occurs when the parameter μ passes through a critical value. The direction and the stability of Hopf bifurcation periodic solutions are determined by the Nyquist criterion and the graphical Hopf bifurcation theorem. Some numerical simulations for justifying the theoretical analysis are also provided.     

Couette–Poiseuille flow of Bingham fluids between two porous parallel plates with slip conditions

Analytical solutions of Couette–Poiseuille flow of Bingham fluids between two porous parallel plates are derived. This study extends the work of Tsangaris et al. [S. Tsangaris, C. Nikas, G. Tsangaris, P. Neofytou, Couette flow of a Bingham plastic in a channel with equally porous parallel walls, J. Non-Newtonian Fluid Mech. 144 (2007) 42–48] to a general situation where the slip effect at the porous walls is considered. It is found that the form of the flow inside the channel depends not only on the Bingham number Bn, the Couette number Co (related to the moving wall) and the transverse Reynolds number Re, but also on the slip parameter Cs at the porous walls. In both the Co–Re diagram and the Co–Bn diagram, the region where plug flow appears enlarges as the slip effect increases, especially in the case where Co is negative. In the case where plug flow and double shear flow coexist, the transverse position of the plug flow and the shear rate at the boundaries exhibit two opposite behaviors when Cs increases, depending on the value of the other three dimensionless numbers. In other cases, slippage always weakens the shearing deformation of the flow.     

Homogenization of Stiff Plates and Two-Dimensional High-Viscosity Stokes Equations

The paper deals with the homogenization of stiff heterogeneous plates. Assuming that the coefficients are equi-bounded in L ¹, we prove that the limit of a sequence of plate equations remains a plate equation which involves a strongly local linear operator acting on the second gradients. This compactness result is based on a div-curl lemma for fourth-order equations. On the other hand, using an intermediate stream function we deduce from the plates case a similar result for high-viscosity Stokes equations in dimension two, so that the nature of the Stokes equation is preserved in the homogenization process. Finally, we show that the L ¹-boundedness assumption cannot be relaxed. Indeed, in the case of the Stokes equation the concentration of one very rigid strip on a line induces the appearance of second gradient terms in the limit problem, which violates the compactness result obtained under the L ¹-boundedness condition.     

Transient free convection flow between long vertical parallel plates with ramped wall temperature at one boundary in the presence of thermal radiation and constant mass diffusion

The unsteady laminar free convection flow between two long vertical parallel plates with ramped wall temperature at one boundary has been investigated in the presence of thermal radiation and chemical species concentration. The exact solutions of the momentum, energy and concentration equations have been obtained using the Laplace transform technique. The velocity and temperature profiles, skin-friction and Nusselt number variations are shown graphically and the numerical values of the volume flow rate, the total heat rate and species rate added to the fluid are presented in a table. The influence of different system parameters such as the radiation parameter (R), buoyancy ratio parameter (N), Schmidt number (Sc) and time (t) has been analyzed carefully. A critical analysis of the coupled heat and mass transfer phenomena is provided. The free convective flow due to ramped wall temperature has also been compared with the baseline case of flow due to constant wall temperature.     

Some passages to the limit in elasticity theory and “Sapondzhyan’s paradox”

We consider bending of thin plates with polygonal and curvilinear edges and indicate analogies and differences between the boundary conditions and boundary value problems arising in these two cases if the polygon is inscribed in the curvilinear contour and the number k of vertices of the polygon tends to infinity.We believe that the so-called Sapondzhyan paradox that arises when solving the boundary value problems for supported plates with a curvilinear contour and a k-gonal contour inscribed in it as k → ∞ can be called a paradox only by misunderstanding. Sapondzhyan’s paradox was studied in several papers briefly surveyed in the monograph [1]. Apparently, the interpretation of “paradoxes” and the results proposed in the present paper are published for the first time.Sapondzhyan’s paradox can be generalized to the case of bending of the so-called sliding-fixed plates (i.e., the generalized shear force and the rotation angle are zero on the plate contour) with a curvilinear contour and a k-gonal contour inscribed in it as k → ∞.In the case of three-dimensional elasticity problems, we present boundary conditions and boundary value problems similar to those listed above and consider the situations resulting in “paradoxes” similar to those arising in plate bending. We give the corresponding explanations and interpretations.     

Effective characteristics of corrugated plates

Corrugated plates are widely used in modern constructions and structures, because they, in contrast to plane plates, possess greater rigidity. In many cases, such a plate can be modeled by a homogeneous anisotropic plate with certain effective flexural and tensional rigidities. Depending on the geometry of corrugations and their location, the equivalent homogeneous plate can also have rigidities of mutual influence. These rigidities allow one to take into account the influence of bending moments on the strain in the midplane and, conversely, the influence of longitudinal strains on the plate bending [1]. The behavior of the corrugated plate under the action of a load normal to the midsurface is described by equations of the theory of flexible plates with initial deflection. These equations form a coupled system of nonlinear partial differential equations with variable coefficients [2]. The dependence of the coefficients on the coordinates is determined by the corrugation geometry. In the case of a plate with periodic corrugation, the coefficients significantly vary within one typical element and depend on the values of local variables determined in each of the typical elements. There is a connection between the local and global variables, and therefore, the functions of local coordinates are simultaneously functions of global coordinates, which are sometimes called rapidly oscillating functions [3].One of the methods for solving the equations with rapidly oscillating coefficients is the asymptotic method of small geometric parameter. The standard procedure of this method usually includes preparatory stages. At the first stage, as a rule, a rectangular periodicity cell is distinguished to be a typical element. At the second stage, the scale of global coordinates is changed so that the rectangular structure periodicity cells became square cells of size l × l. The third stage consists in passing to the dimensionless global coordinates relative to the plate characteristic dimension L. As a result, the dependence between the new local variables and the new scaled dimensionless variables is such that the factor 1/α, where α=l/L ? 1 is a small geometric parameter, appears in differentiating any function of the local coordinate with respect to the global coordinate. After this, the solution of the problem in new coordinates is sought as an asymptotic expansion in the small geometric parameter [1], [4–10].We note that, in the small geometric parameter method, the asymptotic series simultaneously have the form of expansions in the gradients of functions depending only on the global coordinates. This averaging procedure can be applied to linear and nonlinear boundary value problems for differential equations with variable periodic coefficients for which the periodicity cell can be affinely transformed into the periodicity cube. In the case of an arbitrary dependence of the coefficients on the coordinates (including periodic dependence), another averaging technique can be used in linear problems. This technique is based on the possibility of the integral representation of the solution of the original problem for the linear equation with variables coefficients in terms of the solution of the same problem for an equation of the same type but with constant coefficients [11–13]. The integral representation implies that the solution of the original problem can be represented in the form of the series in the gradients of the solution of the problem for the equation with constant coefficients [13].The aim of the present paper is to develop methods for calculating effective characteristics of corrugated plates. To this end, we first write out the equilibrium equations for a flexible anisotropic plate, which is inhomogeneous in the thickness direction and in the horizontal projection, with an initial deflection. We write these equations in matrix form, which allows one to significantly reduce the length of the expressions and to simplify further calculations. After this, we average the initial matrix equations with variable coefficients. The averaging procedure implies the statement of problems such that, after solving them, we can calculate the desired effective characteristics. By way of example, we consider the case of a corrugated plate made of a homogeneous isotropic material whose corrugations are hexagonal in the horizontal projection. In this case, we obtain approximate expressions for the components of the effective tensors of flexural rigidity and longitudinal compliance and expressions for the effective plate thickness.     

A Bending-Gradient model for thick plates,Part II: Closed-form solutions for cylindrical bending of laminates

In the first part (Lebée and Sab, 2010a) of this two-part paper we have presented a new plate theory for out-of-plane loaded thick plates where the static unknowns are those of the Kirchhoff–Love theory (3 in-plane stresses and 3 bending moments), to which six components are added representing the gradient of the bending moment. The new theory, called Bending-Gradient plate theory is an extension to arbitrarily layered plates of the Reissner–Mindlin plate theory which appears as a special case when the plate is homogeneous. Moreover, we demonstrated that, in the general case, the Bending-Gradient model cannot be reduced to a Reissner–Mindlin model. In this paper, the Bending-Gradient theory is applied to laminated plates and its predictions are compared to those of Reissner–Mindlin theory and to full 3D (Pagano, 1969) exact solutions. The main conclusion is that the Bending-Gradient gives good predictions of deflection, shear stress distributions and in-plane displacement distributions in any material configuration. Moreover, under some symmetry conditions, the Bending-Gradient model coincides with the second-order approximation of the exact solution as the slenderness ratio L/h goes to infinity.     

Residual mass and flow regimes for the immiscible liquid–liquid displacement in a plane channel

The motion of two immiscible liquids in a plane channel is analyzed for the case in which the flow conditions and the interactions between the liquids and the solid surface maintain the displaced fluid attached to the wall. The Galerkin Finite Element Method is used to compute the velocity field and the configuration of the interface between the two fluids. We compare the residual mass fraction left on the wall with its two counterparts in capillary tubes, namely residual mass fraction and dimensionless layer thickness of the displaced fluid. The main result of this comparison was that although there is a qualitative similarity concerning the layer thickness between the two cases, the residual fraction of mass presented an important difference, showing that when the aspect ratio of the capillary passage is large there is an increase in the displacement efficiency. The thickness of the displaced liquid film attached to the channel walls is a function of the capillary number (Ca) and the viscosity ratio (N_μ). A map of streamlines in the Cartesian space (Ca, N_μ) with the different flow regimes of the problem is presented. We also showed that we can adapt the available analytical results obtained for gas-displacement in capillary tubes to the plane channel case, for low values of Ca.     

Étude théorique et numérique du comportement d'un assemblage de plaquesTheoretical and numerical study of the behaviour of bonded plates

We consider the flexural static behavior of two coplanar Kirchhoff–Love plates bonded by a thin third one, softer and narrower. In the case of perfect adhesion or unilateral conditions between the adherents and the adhesive, we study the asymptotic behavior of the assembly when the stiffness and the narrowness of the intermediate plate tend to zero. To cite this article: F. Zaittouni et al., C. R. Mecanique 330 (2002) 359–364.     

Regularizing Piecewise Smooth Differential Systems: Co-Dimension 2 Discontinuity Surface

In this paper, we are concerned with numerical solution of piecewise smooth initial value problems. Of specific interest is the case when the discontinuities occur on a smooth manifold of co-dimension 2, intersection of two co-dimension 1 singularity surfaces, and which is nodally attractive for nearby dynamics. In this case of a co-dimension 2 attracting sliding surface, we will give some results relative to two prototypical time and space regularizations. We will show that, unlike the case of co-dimension 1 discontinuity surface, in the case of co-dimension 2 discontinuity surface the behavior of the regularized problems is strikingly different. On the one hand, the time regularization approach will not select a unique sliding mode on the discontinuity surface, thus maintaining the general ambiguity of how to select a Filippov vector field in this case. On the other hand, the proposed space regularization approach is not ambiguous, and there will always be a unique solution associated to the regularized vector field, which will remain close to the original co-dimension 2 surface. We will further clarify the limiting behavior (as the regularization parameter goes to 0) of the proposed space regularization to the solution associated to the sliding vector field of Dieci and Lopez (Numer Math 117:779–811, 2011). Numerical examples will be given to illustrate the different cases and to provide some preliminary exploration in the case of co-dimension 3 discontinuity surface.     

A Bending-Gradient model for thick plates. Part I: Theory

This is the first part of a two-part paper dedicated to a new plate theory for out-of-plane loaded thick plates where the static unknowns are those of the Kirchhoff–Love theory (3 in-plane stresses and 3 bending moments), to which six components are added representing the gradient of the bending moment. The new theory, called the Bending-Gradient plate theory is described in the present paper. It is an extension to arbitrarily layered plates of the Reissner–Mindlin plate theory which appears as a special case of the Bending-Gradient plate theory when the plate is homogeneous. However, we demonstrate also that, in the general case, the Bending-Gradient model cannot be reduced to a Reissner–Mindlin model. In part two (Lebée and Sab, 2011), the Bending-Gradient theory is applied to multilayered plates and its predictions are compared to those of the Reissner–Mindlin theory and to full 3D Pagano’s exact solutions. The main conclusion of the second part is that the Bending-Gradient gives good predictions of both deflection and shear stress distributions in any material configuration. Moreover, under some symmetry conditions, the Bending-Gradient model coincides with the second-order approximation of the exact solution as the slenderness ratio L/h goes to infinity.     

Resonance Tongues and Spectral Gaps in Quasi-Periodic Schrödinger Operators with One or More Frequencies. A Numerical Exploration

In this article we investigate numerically the spectrum of some representative examples of discrete one-dimensional Schrödinger operators with quasi-periodic potential in terms of a perturbative constant b and the spectral parameter a. Our examples include the well-known Almost Mathieu model, other trigonometric potentials with a single quasi-periodic frequency and generalisations with two and three frequencies. We computed numerically the rotation number and the Lyapunov exponent to detect open and collapsed gaps, resonance tongues and the measure of the spectrum. We found that the case with one frequency was significantly different from the case of several frequencies because the latter has all gaps collapsed for a sufficiently large value of the perturbative constant and thus the spectrum is a single spectral band with positive Lyapunov exponent. In contrast, in the cases with one frequency considered, gaps are always dense in the spectrum, although some gaps may collapse either for a single value of the perturbative constant or for a range of values. In all cases we found that there is a curve in the (a, b)-plane which separates the regions where the Lyapunov exponent is zero in the spectrum and where it is positive. Along this curve, which is b = 2 in the Almost Mathieu case, the measure of the spectrum is zero.     

An exact solution for a model of pressure-dependent plasticity in an un-steady plane strain process

An incompressible material obeying a pressure-dependent yield condition is confined between two planar plates which are inclined at an angle 2α. The plates intersect in a hinged line and the angle α slowly decreases from an initial value. An initial/boundary value problem for the flow of the material is formulated and solved for the stress and the velocity fields, the solution being in closed form. The material is assumed to obey a special case of the double slip and rotation model, which generalises the classical plastic potential model and is also a variant of the double shearing model. The solution for the velocity field may exhibit sliding or sticking at the plates. Solutions which exhibit sticking may have a rigidly rotating zone in the region adjacent to the plates. It is shown that sliding occurs when the value of α is less than a certain critical value α_c ; that sticking occurs without a rigid zone if α exceeds or equals α_c but is less than a second critical value α₀; and that sticking with a rigid zone adjacent to the plates occurs if α exceeds α₀. The values of α_c and α₀ coincide for a certain range of model parameters. Solutions which exhibit sliding are singular. Qualitative features of the solution found are compared with those of the solution for the classical plastic potential model.     

Thrust enhancement due to flexible trailing-edge of plunging foils

Drag reduction for hydrofoils is studied through thrust generation on foils plunging at low Strouhal numbers in order to simulate the action of the ocean waves. Force, deformation and flow field measurements are presented for a partially flexible plunging foil in water tunnel experiments. The foil is predominantly rigid with a short flexible trailing-edge plate of length: L=0.1c, 0.2c, or 0.3c. Using flexible plates, whose natural structural frequency is much higher than the frequency of the plunge oscillations, increases thrust compared to the rigid case. Flexibility is generally more effective for larger lengths of the flexible plate and smaller plunge amplitudes. The maximum observed is therefore for the largest length and smallest amplitude studied: L=0.3c and a=0.1c and equates to 28% more thrust than the rigid case. Optima are observed in the non-dimensional rigidity (λ) versus flap angle amplitude (δ, which is a measure of the relative deformation) parameter space. These occur at λ≈2 and δ≈7–13° for a wide range of flexible plate length and plunge amplitude. Whilst a satisfactory explanation of why there is an optimal flap amplitude remains unavailable, the case of optimal flap angle amplitude results in increased trailing-edge vortex circulation, giving a stronger reverse Kármán vortex street and thus a stronger time-averaged jet.     

Exact characteristic equations for some of classical boundary conditions of vibrating moderately thick rectangular plates

The dimensionless equations of motion are derived based on the Mindlin plate theory to study the transverse vibration of thick rectangular plates without further usage of any approximate method. The exact closed form characteristic equations are given within the validity of the Mindlin plate theory for plates having two opposite sides simply supported. The six distinct cases considered involve all possible combinations of classical boundary conditions at the other two sides of rectangular plates. Accurate eigenfrequency parameters are presented for a wide range of aspect ratio η and thickness ratio δ for each case. The three dimensional deformed mode shapes together with their associated contour plots obtained from the exact closed form eigenfunctions are also presented. Finally, the effect of boundary conditions, aspect ratios and thickness ratios on the eigenfrequency parameters and vibratory behavior of each distinct cases are studied in detail. It is believed that in the present work, the exact closed form characteristic equations and their associated eigenfunctions, except for the plates with four edges simply supported, for the rest of considered six cases are obtained for the first time.     

Steady and transient sliding under rate-and-state friction

The physics of dry friction is often modelled by assuming that static and kinetic frictional forces can be represented by a pair of coefficients usually referred to as μ_s and μ_k, respectively. In this paper we re-examine this discontinuous dichotomy and relate it quantitatively to the more general, and smooth, framework of rate-and-state friction. This is important because it enables us to link the ideas behind the widely used static and dynamic coefficients to the more complex concepts that lie behind the rate-and-state framework. Further, we introduce a generic framework for rate-and-state friction that unifies different approaches found in the literature.We consider specific dynamical models for the motion of a rigid block sliding on an inclined surface. In the Coulomb model with constant dynamic friction coefficient, sliding at constant velocity is not possible. In the rate-and-state formalism steady sliding states exist, and analysing their existence and stability enables us to show that the static friction coefficient μ_s should be interpreted as the local maximum at very small slip rates of the steady state rate-and-state friction law.Next, we revisit the often-cited experiments of Rabinowicz (J. Appl. Phys., 22:1373–1379, 1951). Rabinowicz further developed the idea of static and kinetic friction by proposing that the friction coefficient maintains its higher and static value μ_s over a persistence length before dropping to the value μ_k. We show that there is a natural identification of the persistence length with the distance that the block slips as measured along the stable manifold of the saddle point equilibrium in the phase space of the rate-and-state dynamics. This enables us explicitly to define μ_s in terms of the rate-and-state variables and hence link Rabinowicz's ideas to rate-and-state friction laws.This stable manifold naturally separates two basins of attraction in the phase space: initial conditions in the first one lead to the block eventually stopping, while in the second basin of attraction the sliding motion continues indefinitely. We show that a second definition of μ_s is possible, compatible with the first one, as the weighted average of the rate-and-state friction coefficient over the time the block is in motion.     

Testing the accuracy of the overlap criterion

Here we investigate the accuracy of the overlap criterion when applied to a simple near-integrable model in both its 2D and 3D versions. To this end, we consider, respectively, two and three quartic oscillators as the unperturbed system, and couple the degrees of freedom by a cubic, non-integrable perturbation. For both systems we compute the unperturbed resonances up to order O(ε²), and model each resonance by means of the pendulum approximation in order to estimate the theoretical critical value of the perturbation parameter for a global transition to chaos. We perform several surface of sections for the bi-dimensional case to derive an empirical value to be compared to our theoretical estimation. Although both values are of the same order of magnitude, there is a significant difference between them. For the 3D case a numerical estimate is attained that we observe matches quite well the critical value resulting from theoretical means. This confirms once again that calculating resonances up to O(ε²) suffices in order the overlap criterion to work out.