Perturbation Solution Research Articles

New Numerical Iteration Schemes Based on Perturbation Iteration Algorithms

Perturbation–Iteration algorithms (PIA) have been developed recently to solve differential equations analytically. A continuous solution in terms of closed form functions as an approximation of the original equation can be found using the method. The method has been implemented to algebraic equations, ordinary and partial differential equations successfully. Based on the formalism developed previously, in this work, purely numerical versions of the perturbation–iteration algorithms are proposed for the first time for first-order nonlinear ordinary differential equations. The new algorithms are called as the Numerical Perturbation–Iteration Algorithms (NPIA) to distinguish them from the continuous analytical ones (PIA) and the method in general will be called the Numerical Perturbation–Iteration Method (NPIM). Iteration schemes based on one correction term in the perturbation solution and first-order derivatives in the Taylor series expansions (NPIA(1,1)), and one correction term in the perturbations and second-order derivatives in the series expansions (NPIA(1,2)) are derived. The numerical results are compared with the exact solutions and a very good match is observed. NPIA(1,2) produced slightly better results compared to NPIA(1,1) with total errors being at least [Formula: see text], [Formula: see text] being the step size for both methods. An improvement suggestion for NPIA(1,1) is proposed in the last section to eliminate unrealistic blow-up solutions which are encountered for some specific equations.

Alternating Current Electroosmotic Flow of Maxwell Fluid in a Parallel Plate Microchannel with Sinusoidal Roughness.

The EOF of a viscoelastic Maxwell fluid driven by an alternating pressure gradient and electric field in a parallel plate microchannel with sinusoidal roughness has been investigated within the Debye-Hückel approximation based on boundary perturbation expansion and separation of variables. Perturbation solutions were obtained for the potential distribution, the velocity and the mean velocity, and the relation between the mean velocity and the roughness. There are significant differences in the velocity amplitudes of the Newtonian and Maxwell fluids. It is shown here that the velocity distribution of the viscoelastic fluid is significantly affected by the roughness of the walls, which leads to the appearance of fluctuations in the fluid. Also, the velocity is strongly dependent on the phase difference θ of the roughness of the upper and lower plates. As the oscillation Reynolds number ReΩ increases, the velocity profile and the average velocity um(t) of AC EOF oscillate rapidly but the velocity amplitude decreases. The Deborah number De plays a similar role to ReΩ, which makes the AC EOF velocity profile more likely to oscillate. Meanwhile, phase lag χ (representing the phase difference between the electric field and the mean velocity) decreases when G and θ are increased. However, for larger λ (e.g., λ > 3), it almost has no phase lag χ.

Unification: Emergent universe followed by inflation and dark epochs from multi-field theory

A two scalar field model that incorporates non Riemannian Measures of integration or usually called Two Measures Theory (TMT) is introduced, in order to unify the early and present universe. In the Einstein frame a K-essence is generated and as a consequence for the early universe, we can have a Non Singular Emergent universe followed by Inflation and for the present universe dark epochs with consistent generation of dark energy (DE), dark matter (DM) and stiff matter. The scale invariance is introduced and then is spontaneously broken from the integration of the degrees of freedom associated with the modified measures. The resulting effective potentials and K-essence in the Einstein frame produce three flat regions corresponding to the different epochs mentioned before. For the first flat region we can associate an emergent and an inflationary universe. Here for a parameter-space region this flat plateau possesses a non singular stable emergent universe solution which characterizes an initial epoch of evolution that precedes the inflationary scenario. Also assuming this first plateau, we study the inflation in the framework of the slow-roll approximation. In this scenario under the slow roll approximation we obtain a linear combination that is a constant. The corresponding cosmological perturbations in our model are determined and we also obtain the different constrains on the parameter-space from the Planck data.In the following flat region DE and also the DM, which does not need to be introduced separately, it is instead a result of a K-essence induced by the multi measures, multi field theory. Also stiff matter component is automatically generated from the K-essence theory from two scalar fields. From the perturbative analysis associated to the perturbation solution of background, we find a correlation between the two scalars. Besides, we obtain that our model during the dark epoch has a behavior of tracking freezing model.

Dynamic stiffness matrix perturbation theory for time-varying structural analysis

Matrix perturbation theory (MPT) is widely used in sensitivity analysis and reanalysis of structural eigenvalue problems when small changes are imposed on parameters. However, the conventional MPT based on the finite element method is only applicable to the cases with small parameter variations, and it is difficult to avoid cumulative errors caused by multi-step perturbations. A key issue for the application of MPT is how to evaluate the error of perturbation solution and correct the results during the process of multi-step perturbation solution. By invoking the dynamic stiffness method (DSM), Wittrick-Williams algorithm, and the precise integration method (PIM), in this paper, we propose an improved method with controllable accuracy, self-checking and elimination of accumulated errors, referred to as the dynamic stiffness matrix perturbation theory (DSMPT). In present method, the rationality of the results can be estimated preliminarily according to dynamic stiffness equation, which is not present in the original MPT, can greatly enhance the accuracy of the perturbation results and evaluate the errors. Case studies prove that the DSMPT combined with PIM provides a convenient way to realize accuracy self-inspection and rectification of the perturbation results, and the maximum relative error of eigenvalue and response can be controlled to be smaller than 0.2% and 0.025%, respectively. Besides, DSMPT can also be applied to the analysis of coupled and nonlinear dynamic system with time-varying modal characteristics. Finally, a thermally-induced vibration analysis of tensegrity modules during space deployment is taken as an example to further illustrate the advancements of the DSMPT in solving time-varying dynamic problems. Results show that the DSMPT effectuate an almost sixfold reduction in computational time.

Revisiting coupled CDM-massive neutrino perturbations in diverse cosmological backgrounds

Massive neutrinos are well-known to cause a characteristic suppression in the growth of structures at scales below the neutrino free-streaming length. A detailed understanding of this suppression is essential in the era of precision cosmology we are entering into, enabling us to better constrain the total neutrino mass and possibly probe (beyond)-ΛCDM cosmological model(s). Instead of the usual N-body simulation or Boltzmann solver, in this article we consider a two-fluid framework at the linear scales, where the neutrino fluid perturbations are coupled to the CDM (+ baryon) fluid via gravity at redshifts of interest. Treating the neutrino mass fraction f ν as a perturbative parameter, we find solutions to the system with redshift-dependent neutrino free-streaming length in ΛCDM background via two separate approaches. The perturbative scale-dependent solution is shown to be in excellent agreement with numerical solution of the two-fluid equations valid to all orders in f ν, and also agrees with results from CLASS to a good accuracy. We further generalize the framework to incorporate different evolving dark energy backgrounds and found sub-percent level differences in the suppression, all of which lie within the observational uncertainty of BOSS-like surveys. We also present a brief discussion on the prospects of the current analysis in the context of upcoming missions.

Lubrication approximation of pressure-driven viscoelastic flow in a hyperbolic channel

We investigate theoretically the steady incompressible viscoelastic flow in a hyperbolic contracting channel. The fluid viscoelasticity is modeled using the upper convected Maxwell (UCM), Oldroyd-B, Phan-Thien and Tanner (PTT), Giesekus, and the finite elasticity non-linear elastic dumbbell with the Peterlin approximation (FENE-P) models. We first develop the general governing equations for flow within a non-deformable channel whose cross section varies with the distance from the inlet. We then exploit the classic lubrication approximation, assuming a small aspect ratio of the channel to simplify the original governing equations. The final equations, which we formulate in terms of the stream unction, are then solved analytically using a high-order asymptotic scheme in terms of the Deborah number, De, and the formulas for the average pressure drop are derived up to eight orders in De. The accuracy of the original perturbation solution is enhanced and extended over a wide range of parameters by implementing a convergence acceleration method for truncated series. Furthermore, convergence of the transformed solutions for the average pressure drop is demonstrated. The validity and accuracy of the theoretical results is independently confirmed through comparison with numerical results from simulations performed using high-order finite differences and pseudospectral methods. The results reveal the decrease in the average pressure drop with increasing the Deborah number, the polymer viscosity ratio, and the ratio of the inlet to the outlet height. We also show that the fundamental UCM and Oldroyd-B models can predict the major viscoelastic phenomena for this type of internal and confined lubrication flows, while the effect of the rheological parameters of the PTT, Giesekus, and FENE-P models on the results is minor.

From the string vacuum to FLRW or de Sitter via α' corrections

We first make more precise a recent “Hamiltonian” reformulation of the Hohm-Zwiebach approach to the tree-level, O(d,d)-invariant string cosmology equations at all orders in the α' expansion, and recall how it allows to give a simple characterization of a large class of cosmological scenarios connecting, through a non-singular bounce, two duality-related perturbative solutions at early and late times. We then discuss the effects of adding to the action a non-perturbative, O(d,d)-breaking, dilaton potential V(ϕ). The resulting cosmological solutions, assumed to approach at early times the perturbative string vacuum (with vanishing curvature and string coupling), can stabilize the dilaton at late times and simultaneously approach either a matter-dominated FLRW cosmology or a de-Sitter-like inflationary phase, depending on initial conditions and on the properties of V(ϕ) at moderate-coupling. We also identify a general mechanism for generating isotropic late-time attractors from a large basin of anisotropic initial conditions.

Experimental study on Richtmyer–Meshkov instability at a light–heavy interface over a wide range of Atwood numbers

Richtmyer–Meshkov instability (RMI) at a light–heavy single-mode interface over a wide range of post-shock Atwood numbers $A_1$ is studied systematically through elaborate experiments. The interface generation and $A_1$ variation are achieved by the soap-film technology and gas-layer scheme, respectively. Qualitatively, the nonlinear interface evolution features, including spike, bubble and roll-up structures, are more significant in RMI with higher $A_1$ . Quantitatively, both the impulsive model and an analytical linear model perform well in predicting the linear growth rate under a wide range of $A_1$ conditions. For the weakly nonlinear stage, the significant spike acceleration occurring when $A_1$ is high, which is observed experimentally for the first time, results in the evolution law of RMI with high $A_1$ being different from the counterpart with low or intermediate $A_1$ . None of the considered nonlinear models is found to be applicable for RMI under all $A_1$ conditions, and the predictive capabilities of these models are analysed and summarized. Based on the present experimental results, an empirical nonlinear model is proposed for RMI over a wide range of $A_1$ . Further, modal analysis shows that in RMI with high (low or intermediate) $A_1$ , high-order harmonics evolve rapidly (slowly) and cannot (can) be ignored. Accordingly, for RMI with high (low or intermediate) $A_1$ , the modal model proposed by Zhang & Sohn (Phys. Fluids, vol. 9, 1997, pp. 1106–1124) is less (more) accurate than the one proposed by Vandenboomgaerde et al. (Phys. Fluids, vol. 14, 2002, pp. 1111–1122), since the former ignores perturbation solutions higher than fourth order (the latter retains only terms with the highest power in time).

Maxwell homogenisation methodology for evaluation of effective elastic constants of weakly-nonlinear particulate composites

The far-field methodology, developed by J. C. Maxwell in 1873 for the evaluation of the electrical conductivity of a particle-reinforced material, is utilised to estimate the effective third order elastic constants of composite media containing a random distribution of spherical particles. The results are explicit analytical relationships derived using a perturbation solution, that are in a broad agreement with previous studies employing different homogenisation schemes. In contrast to other homogenisation schemes, Maxwell's methodology relies only on the far-field asymptotics. Therefore, the methodology may be more suitable for the evaluation of effective mechanical and physical properties of composites where the far-field asymptotics are known or can be derived. The main fundamental result of this work is the first demonstration that the Maxwell methodology may be applied to weakly nonlinear problems. Therefore, it can offer a new, simple method to analyse homogenisation problems in many other fields.

Unpolarized QED parton distribution functions in NLO

Perturbative solutions for unpolarized QED parton distribution and fragmentation functions are presented explicitly in the next-to-leading logarithmic approximation. The scheme of iterative solution of QED evolution equations is described in detail. Terms up to are calculated analytically, where is the large logarithm which depends on the factorization energy scale μ F ≫ m e . The results are process independent and relevant for future high-precision experiments.

Macro-Finance Models with Nonlinear Dynamics

We review macro-finance models featuring nonlinear dynamics that have recently been developed in the literature, including models with funding liquidity constraints, market liquidity frictions, and bank run frictions, and discuss the empirical evidence and challenges of this class of models. We also construct an illustrative model featuring financial frictions and nonlinear dynamics for readers who are unfamiliar with the literature. We solve the model using different solution techniques, including both global and perturbation solution methods, and comprehensively compare the accuracy of these solutions. Within this framework, we highlight that local linearization approximations omit important nonlinear dynamics and yield biased impulse responses.

Partial differential equation‐based multiagent systems formation‐containment control with α‐leader input delay

AbstractThis paper investigates the formation‐containment problem for a multiagent system, where the agents are classified into containment followers, formation leaders, anchor leader, and ‐leader. The ‐leader stabilizes formation leaders to desired formation deployment, while formation leaders adopt local state feedback distributed control law without desired formation position and the anchor leader fixes to its desired position. The containment followers adopted consensus distributed control law and converge to the convex hull spanned by the leader agents. Formation leaders' distributed control law parameters are proposed based on the discrete form of the desired partial differential equation (PDE) formation deployment, which is stabilized by the ‐leader agent and transited by communication topology. Based on its neighbor leader agent's delay state, ‐leader control law designed an integral‐type delay‐compensated control law by the backstepping method and the equivalence principle. In order to prove the rightness and existence of delay‐compensated control law, the well‐posedness of kernel functions are given by semigroup perturbation theory and solutions in the sense of distribution. Finally, a simulation example is given to verify the theoretical results.

Quantum perturbative solutions of extended Snyder and Yang models with spontaneous symmetry breaking

We propose ħ-expansions as perturbative solutions of quantum extended Snyder and Yang models, with ħ-independent classical zero-th order terms responsible for the spontaneous breaking of D=4 and D=5 de Sitter symmetries. In such models, with algebraic basis spanned by oˆ(D,1) Lie algebra generators, we relate the vacuum expectation values (VEV) of the spontaneously broken generators with the Abelian set of ten (Snyder, D=4) or fifteen (Yang, D=5) antisymmetric tensorial generalized coordinates, which are also used as zero order input for obtaining the perturbative solutions of quantum extended Snyder and Yang models. In such a way we will attribute to these Abelian generalized coordinates the physical meaning of the order parameters describing spontaneous symmetry breaking (SSB). It appears that the consecutive terms in ħ-power series can be calculated explicitly if we supplement the SSB order parameters by the dual set of tensorial commutative momenta.

Analytical soliton solutions for cubic-quartic perturbations of the Lakshmanan-Porsezian-Daniel equation using the modified extended tanh function method

In this paper, the improved modified extended Tanh function technique is used to find analytical solutions to the cubic-quartic Lakshmanan-Porsezian-Daniel model which describes the pulse propagation in optical fiber. Using the traveling wave theory, the nonlinear partial differential form of the cubic-quartic Lakshmanan-Porsezian-Daniel equation is transformed into an ordinary differential equation. Then applying the proposed methods bright, dark, bright-dark, periodic, singular, and Jacobian elliptic solutions were derived as analytical solutions. Besides, hyperbolic, rational, trigonometric, and exponential functions were also obtained. In order to understand the obtained solutions, graphical simulations were depicted by 3D surfaces and 2D graphs.

A Perturbative Solution for Nonlinear Stratified Upwelling over a Frictional Slope

Abstract A perturbative solution of simplified primitive equations for nonlinear weakly stratified upwelling over a frictional slope is found that resolves the vertical structure of velocity fields and can satisfy Ertel’s potential vorticity conservation in the stratified inviscid interior. The solution uses assumptions consistent with the model proposed by Lentz and Chapman, including a steady-state, constant cross-shore density gradient, no alongshore gradients, laterally inviscid, and consideration of cross-shore advection of alongshore momentum. The solution resolves the vertical structure of velocity fields (including subsurface maxima of compensational flow, not resolved by Lentz and Chapman) and can satisfy Ertel’s potential vorticity conservation in the stratified inviscid interior. The dynamics are similar to Lentz and Chapman; bottom stress balances alongshore wind stress in a homogeneous density ocean and is replaced by nonlinear cross-shore transport of alongshore momentum as the Burger number (S = αN/f, where α, N, and f are the bottom slope, buoyancy frequency, Coriolis frequency, respectively) increases. When the solution uses the empirical relation between cross-shore and vertical density gradients proposed by Lentz and Chapman, vorticity conservation is not satisfied and the nonlinear momentum transport estimated by the solution linearly increases with S, asymptotically matching Lentz and Chapman for S < 1. When the solution conserves interior potential vorticity, the momentum transport is proportional to S2 for S < 1 and is in better agreement with numerical simulations.

Magnetohydrodynamic Shocks Revisited: Magnetically Constraining the Upstream Solar Wind Condition

Magnetohydrodynamic (MHD) shocks are commonly found in various space and astrophysical systems. One often encounters the problem that only the magnetic field data are available when studying interplanetary shocks and planetary bow shocks in the solar system, and not the plasma data due to instrumental and/or telemetry limitations. We aim to reveal the mathematical structure of the MHD Rankine–Hugoniot relation under the condition that magnetic field changes across the shock are known, which can be devised into a diagnosis tool to constrain the shock parameters. The perturbative solution is already known for the Rankine–Hugoniot equation. The solution is combined with the condition of the de Hoffmann–Teller frame for the vanishing convective electric field, and we obtain a two-dimensional matrix equation for the density jump and the Alfvén Mach number. The lesson is that the shock parameters (density jump and Alfvén Mach number) are analytically obtained by inverting the MHD Rankine–Hugoniot equation under the condition that the plasma beta is either a priori set or known. Further numerical studies indicate that the presented shock jump estimation works particularly well for quasi-perpendicular shocks. Moreover, a test case with the Earth bow shock crossing data supports our theoretical development. The presented method opens the door to determining or constraining the shock parameters even if the data are limited to the magnetic field only.

DAMPING OF FINITE AMPLITUDE SOLITARY WAVES IN A FLUME

Solitary wave is a permanent wave when the dissipation is ignored. Many analytical solutions have been developed for finite amplitude solitary waves. In additional to the perturbation solutions, closed form solutions are also available (e.g., McCowan 1891, Clamond and Fructus 2003, being denoted as CF from hereon), which are more accuracy, especially for larger amplitude solitary waves which were discussed in Wang and Liu (2022). In Wang and Liu (2022) ’s laboratory experiments, solitary waves are slowly damped along the wave flume, which can be attributed to the energy dissipation inside the boundary layers on the bottom and sidewalls. Keulegan (1948) first derived an approximate solution to describe the wave damping in the wave flume for small amplitude solitary wave. The basic idea of calculating the damping effect is that the rate of energy dissipation inside the boundary layers must be the same as the rate of wave energy loss. In this work, the same principle of energy balance is adopted for calculating the wave damping for finite amplitude solitary waves in a wave flume. The closed form solutions provided in CF are employed here for analyses.

Morphological/Dynamic Instability of Directional Crystallization in a Finite Domain with Intense Convection

This study is devoted to the morphological/dynamic instability analysis of directional crystallization processes in finite domains with allowance for melt convection. At first, a linear instability theory for steady-state crystallization with a planar solid/liquid interface in the presence of convection was developed. We derived and analyzed a dispersion relation showing the existence of morphological instability over a wide range of wavenumbers. This instability results from perturbations arriving at the solid/liquid interface from the cooled wall through the solid phase. Also, we showed that a planar solid/liquid interface can be unstable when it comes to dynamic perturbations with a zero wavenumber (perturbations in its steady-state velocity). A branch of stable solutions for dynamic perturbations is available too. The crystallizing system can choose one of these branches (unstable or stable) depending of the action of convection. The result of morphological and dynamic instabilities is the appearance of a two-phase (mushy) layer ahead of the planar solid/liquid interface. Therefore, our next step was to analyze the dynamic instability of steady-state crystallization with a mushy layer, which was replaced by a discontinuity interface between the purely solid and liquid phases. This analysis showed the existence of dynamic instability over a wide range of crystallization velocities. This instability appears in the solid material at the cooled wall and propagates to the discontinuity interface, mimicking the properties of a mushy layer. As this takes place, at a certain crystallization velocity, a bifurcation of solutions occurs, leading to the existence of unstable and stable crystallization branches simultaneously. In this case, the system chooses one of them depending of the effect of the convection as before. In general, the crystallizing system may be morphologically/dynamically unstable when it comes to small perturbations arriving at the phase interface due to fluctuations in the heat and mass exchange equipment (e.g., fluctuations in the freezer temperature).

Final sheet thickness variation during the non-isothermal calendering process of rheological Prandtl fluid polymer

Industries widely use the calendering process to create coated fabrics, thermoplastic films, plastic sheeting, and textile fabrics to procure the prerequisite surface smoothness plus texture. This article inspects the non-isothermal calendering phenomena with heat transfer mechanism to report the rheological behavior of Prandtl fluid on the sensitivity of sheet thickness with velocity slip at roll’s surface. Lubrication theory is utilized to simplify the governing expressions before normalization. We obtained the zeroth plus first order perturbative solutions for the velocity profile, pressure gradient, pressure distribution, temperature, and sheet thickness. And other mechanical properties are evaluated with a numerical solution. Through graphs and tabular data, how material parameters influence the various flow and engineering parameters are reported. As it can be seen, the material parameters are the main parameters that control the pressure gradient, leading to different final sheet thickness and separating force. As the material parameters increased, sheet thickness and separating force increased.

Looking at QED with Dyson–Schwinger Equations: Basic Equations, Ward–Takahashi Identities and the Two-Photon-Two-Fermion Irreducible Vertex

A minimal truncated set of the integral Dyson–Schwinger equations, in Minkowski spacetime, that allows to explore QED beyond its perturbative solution is derived for general linear covariant gauges. The minimal set includes the equations for the fermion and photon propagators, the photon-fermion vertex, and the two-photon-two-fermion one-particle-irreducible diagram. If the first three equations are exact, to build a closed set of equations, the two-photon-two-fermion equation is truncated ignoring the contribution of Green functions with large number of external legs. It is shown that the truncated equation for the two-photon-two-fermion vertex reproduces the lowest-order perturbative result in the limit of the small coupling constant. Furthermore, this equation allows to define an iterative procedure to compute higher order corrections in the coupling constant. The Ward–Takahashi identity for the two-photon-two-fermion irreducible vertex is derived and solved in the soft photon limit, where one of the photon momenta vanish, in the low photon momenta limit and for general kinematics. The solution of the Ward–Takahashi identity determines the longitudinal component of the two-photon-two-fermion irreducible vertex, while it is proposed to use the Dyson–Schwinger equation to determine the transverse part of this irreducible diagram. The two-photon-two-fermion DSE is solved in heavy fermion limit, considering a simplified version of the QED vertices. The contribution of this irreducible vertex to a low-energy effective photon-fermion vertex is discussed and the fermionic operators that are generated are computed in terms of the fermion propagator functions.