Integral Momentum Research Articles

Moment of momentum integral analysis of turbulent boundary layers with pressure gradient

Turbulent boundary layers on immersed objects can be significantly altered by the pressure gradients imposed by the flow outside the boundary layer. The interaction of turbulence and pressure gradients can lead to complex phenomena such as relaminarization, history effects and flow separation. The angular momentum integral (AMI) equation (Elnahhas & Johnson, J. Fluid Mech., vol. 940, 2022, A36) is extended and applied to high-fidelity simulation datasets of non-zero pressure gradient turbulent boundary layers. The AMI equation provides an exact mathematical equation for quantifying how turbulence, free-stream pressure gradients and other effects alter the skin friction coefficient relative to a baseline laminar boundary layer solution. The datasets explored include flat-plate boundary layers with nearly constant adverse pressure gradients, a boundary layer over the suction surface of a two-dimensional NACA 4412 airfoil and flow over a two-dimensional Gaussian bump. Application of the AMI equation to these datasets maps out the similarities and differences in how boundary layers interact with favourable and adverse pressure gradients in various scenarios. Further, the fractional contribution of the pressure gradient to skin friction attenuation in adverse-pressure-gradient boundary layers appears in the AMI equation as a new Clauser-like parameter with some advantages for understanding similarities and differences related to upstream history effects. The results highlight the applicability of the integral-based analysis to provide quantitative, interpretable assessments of complex boundary layer physics.

Banana diagrams as functions of geodesic distance

We extend the study of banana diagrams in coordinate representation to the case of curved space-times. If the space is harmonic, the Green functions continue to depend on a single variable – the geodesic distance. But now this dependence can be somewhat non-trivial. We demonstrate that, like in the flat case, the coordinate differential equations for powers of Green functions can still be expressed as determinants of certain operators. Therefore, not-surprisingly, the coordinate equations remain straightforward – while their reformulation in terms of momentum integrals and Picard-Fuchs equations can seem problematic. However we show that the Feynman parameter representation can also be generalized, at least for banana diagrams in simple harmonic spaces, so that the Picard-Fuchs equations retain their Euclidean form with just a minor modification. A separate story is the transfer to the case when the Green function essentially depends on several rather than a single argument. In this case, we provide just one example, that the equations are still there, but conceptual issues in the more general case will be discussed elsewhere.

Artificial Intelligence Empowering the Innovation of Ideological and Political Course Teaching Models in Universities

With the rapid development of information technology, AI has permeated education, showing great potential in innovating teaching models. Ideological and Political Courses in universities face challenges like a disconnect between content and student needs, low interaction, and poor engagement, requiring urgent reforms. This study explores AI-enabled approaches such as personalized teaching, intelligent assessment, and VR/AR technologies. It demonstrates that AI can enhance interaction, improve learning experiences, and support teaching quality. Despite challenges like ethics and data privacy, AI’s continued evolution will lead to broader applications, deeper integration, and new momentum for teaching model improvements in these courses.

Exploring the mechanisms behind dual-peak formation in energy harvesting efficiency of semi-passive flapping airfoils with varied spring stiffness

Flapping airfoil energy harvesting has emerged as a promising paradigm inspired by birds and marine organisms. In this investigation, a semi-passive flapping airfoil device with a predetermined pitching motion was examined using a transient numerical simulation method employing an overlapping grid technique. The impact of the spring stiffness on the energy harvesting characteristics of the flapping airfoil was analyzed through the implementation of the integral momentum theorem and dynamic mode decomposition. The results demonstrated that the efficiency and power curves of the flapping airfoil exhibited a distinctive bimodal M shape in response to variations of the spring stiffness. The initial peak efficiency point, which corresponded to the system's natural frequency in proximity to the pitching frequency, was found to be heavily influenced by the presence of leading-edge vortexes, thereby influencing lift generation. By utilizing the integral momentum theorem, it was determined that the Lamb vector term, fluid acceleration term, and total pressure term predominantly contributed to the lift generation. The emergence of the second peak efficiency point was primarily linked to torque variations induced by trailing-edge vortexes.

Velocity-field based load estimation of an unsteady actuator disc

Traditional methods for determining the load of an object in unsteady cases, such as Control-volume Momentum Integration (CMI) and Noca methods, encounter challenges related to near-body acceleration and pressure estimation. This prompts the need for innovative techniques to overcome these limitations. This study introduces a novel vorticity-based approach to expand classical load determination methods. The numerical results of an unsteady actuator disc (AD) CFD simulation are used to verify and quantify the uncertainty of the three methods. In the AD model, the load is prescribed and used as a source term. The velocity and pressure fields are solved from the Navier-Stokes equations using the open-source CFD package OpenFOAM, which offers a robust validation framework. Systematic examination of the robustness and accuracy of the three methods across various load and motion unsteady scenarios reveals that the proposed vorticity-based method excels in steady-state scenarios, exhibiting the highest level of robustness and accuracy. In contrast, the traditional CMI method proves more robust and accurate in pure motion cases. In scenarios involving static and moving AD with load variations, the new vorticity-based method demonstrates superior robustness and accuracy, mainly when moderate vorticity dynamics. Conversely, in cases with increased vorticity dynamics in the wake, the accuracy of the vorticity-based method decreases, with the CMI method showing superiority.

Precision calculation of the entrance length for laminar flow of Bingham fluid between parallel plates with yield number from 0 to 1000

An estimate of the entrance length for a Bingham fluid flowing between parallel plates is needed in engineering applications. The dimensionless entrance length is a function of yield number. The primary objective of this study was to construct plots of dimensionless entrance length vs. yield number that could be used by practitioners for a quick estimation of the entrance length in a Bingham fluid entering between parallel plates, in a range of yield numbers from 0 to 1000. Three well-known calculation methods were used. The methods were found to yield mutually consistent results. A secondary objective of this study was to examine the sensitivity of the calculated entrance length to the value of a dimensionless parameter somewhat arbitrary chosen during the calculations. This parameter specifies the dimensionless value of the unyielded core velocity that is accepted as the core velocity at the end of the entrance region. In all three calculation methods, the entrance length was found to have a noticeable sensitivity to the variation of this parameter. The least sensitive method was the one based on the momentum integral solution. The findings necessitate further, more fundamental research on the entrance length in non-Newtonian fluids.

Integral methods for friction decomposition and their extensions to rough-wall flows

A primary objective of integral methods, such as the momentum integral method, is to discern the physical processes contributing to skin friction. These methods encompass the momentum, kinetic energy and angular momentum integrals. This paper reformulates existing integrals based on the double-averaged Navier–Stokes equations, and extends their application to flows over rough walls. Our derivation yields distinct decompositions for the bottom-wall viscous friction coefficient, denoted as $C_S$ , and the roughness element drag coefficient $C_R$ . The decompositions comprise three terms: a viscous term, a turbulent term and a roughness (dispersive) term – regardless of the flow configuration, be it channel or boundary layer. Notably, when these integrals are evaluated for laminar flow scenarios, only the viscous term remains significant. In addition, we elucidate the spatial distributions of the terms within these decompositions. To demonstrate the practicality of our formulations, we apply them to analyse data from direct numerical simulations of turbulent half-channel flows. These flows feature aligned and staggered cubical roughness at various packing densities. Our analyses, based on kinetic-energy-oriented decompositions, reveal that when the surface coverage density $\lambda _p$ is small, the dominant terms within the decompositions are the viscous and turbulent terms. With increasing $\lambda _p$ , the viscous dissipation term decreases, while the turbulent production term increases and then decreases. These variations arise from a subdued near-wall cycle and the development of a shear layer at the height of the cubes.

On an Application of the Hill Approach to the General Case of the Three-body Problem

In this work, we use the energy and angular momentum integrals as a resource for applying Hill’s approach to the general three-body problem. As a result, we obtain theorems on the Lagrange stability and Hill stability in the general three-body problem. Also, specific features of the general and restricted three-body problems are discussed.

Predicting differential cross sections of electron scattering from N2 using DCMe method

The e + N2 scattering system is usually a typical testing model in electron scattering from molecule. The accurate resonant and the non-resonant differential cross sections (DCSs), integral cross sections (ICSs) and momentum transfer cross sections (MTCSs) of electron scattering from N2 molecule are studied using limited experimental DCSs and the difference converging method for electron scattering (DCMe) proposed in our previous work. The DCMe results well represent those known experimental data and obtain all unknown ones at given energies. It supplies a theoretical alternative to correctly predict the DCSs those may be experimentally unknown or hard to have quantum mechanically. This study shows that.

Integral Representations over Finite Limits for Quantum Amplitudes

We extend previous research to derive three additional M-1-dimensional integral representations over the interval [0,1]. The prior version covered the interval [0,∞]. This extension applies to products of M Slater orbitals, since they (and wave functions derived from them) appear in quantum transition amplitudes. It enables the magnitudes of coordinate vector differences (square roots of polynomials) |x1−x2|=x12−2x1x2cosθ+x22 to be shifted from disjoint products of functions into a single quadratic form, allowing for the completion of its square. The M-1-dimensional integral representations of M Slater orbitals that both this extension and the prior version introduce provide alternatives to Fourier transforms and are much more compact. The latter introduce a 3M-dimensional momentum integral for M products of Slater orbitals (in M separate denominators), followed in many cases by another set of M-1-dimensional integral representations to combine those denominators into one denominator having a single (momentum) quadratic form. The current and prior methods are also slightly more compact than Gaussian transforms that introduce an M-dimensional integral for products of M Slater orbitals while simultaneously moving them into a single (spatial) quadratic form in a common exponential. One may also use addition theorems for extracting the angular variables or even direct integration at times. Each method has its strengths and weaknesses. We found that these M-1-dimensional integral representations over the interval [0,1] are numerically stable, as was the prior version, having integrals running over the interval [0,∞], and one does not need to test for a sufficiently large upper integration limit as one does for the latter approach. For analytical reductions of integrals arising from any of the three, however, there is the possible drawback for large M of there being fewer tabled integrals over [0,1] than over [0,∞]. In particular, the results of both prior and current representations have integration variables residing within square roots asarguments of Macdonald functions. In a number of cases, these can be converted to Meijer G-functions whose arguments have the form (ax2+bx+c)/x, for which a single tabled integral exists for the integrals from running over the interval [0,∞] of the prior paper, and from which other forms can be found using the techniques given therein. This is not so for integral representations over the interval [0,1]. Finally, we introduce a fourth integral representation that is not easily generalizable to large M but may well provide a bridge for finding the requisite integrals for such Meijer G-functions over [0,1].

Effect of finite volume on thermodynamics of quark-hadron matter

The effects of a finite system volume on thermodynamic quantities, such as the pressure, energy density, specific heat, speed of sound, conserved charge susceptibilities, and correlations, in hot and dense strongly interacting matter are studied within the parity-doublet chiral mean field model. Such an investigation is motivated by relativistic heavy-ion collisions, which create a blob of hot QCD matter of a finite volume, consisting of strongly interacting hadrons and potentially deconfined quarks and gluons. The effect of the finite volume of the system is incorporated by introducing lower momentum cutoffs in the momentum integrals appearing in the model, the numerical value of the momentum cutoff being related to the de Broglie wavelength of the given particle species. It is found that some of these quantities show a significant volume dependence, in particular, those sensitive to pion degrees of freedom, and the crossover transition is generally observed to become smoother in finite volume. These findings are relevant for the effective equation of state used in fluid dynamical simulations of heavy-ion collisions and efforts to extract the freeze-out properties of heavy-ion collisions with susceptibilities involving electric charge and strangeness. Published by the American Physical Society 2024

Positron scattering from structurally related biomolecules.

We report the integral elastic, momentum transfer, and inelastic (positronium formation and ionisation) cross sections for positron scattering from structurally related molecules. The molecules chosen for the current investigation are formamide, formylphosphine, formic acid, N-methylformamide, acetone, acetic acid, and formaldehyde. The cross sections were calculated using the optical potential approach and the complex scattering potential-ionisation contribution method for impact energies between 1 and 5 keV. A sizable repository of data is now available for positron scattering from various atoms and molecules; however, data on the impact of positrons on current targets is still scarce and fragmented. While most cross sections are the first of their kind, we analyze our total cross sections (TCSs) with the previous literature available, which has become attractive to researchers trying to model the tracks of charged particles in matter. TCSs have recently seen a resurgence in popularity thanks to their utility in specifying the mean-free path between the collisions of such simulations. We find good qualitative convergence between experimental and theoretical results below and above the positronium formation threshold. However, around the threshold region, a significant discrepancy is encountered, which can be accounted for due to the experiment's lack of forward angle scattering effect discrimination. This level of agreement evolves to become quantitative at intermediate and higher energies.

A two-dimensional friction-decoupling method based on numerical integration method and momentum theorem

This article introduces a simplified Two-Dimensional Friction-Coupling Model (TDFC) designed for typical surface-to-surface contact friction systems in engineering structures. With the assumption of a constant friction force magnitude within each calculation time step, the alteration of the friction direction of TDFC within each calculation time step can be visualized as the rotation of a circular radius under horizontal bidirectional dynamic loads. This dynamic force analysis approach is referred to as the “Circular making”. It is noteworthy that the friction direction change in each calculation time step demonstrates unidirectional rotatability. Leveraging this rotatability and assuming uniform rotation, the momentum theorem accurately determines the direction of the plane friction acting on the TDFC model at the end of each calculation time step. This allows for the decoupled calculation of plane friction and precise response calculation of the TDFC model. The study verifies the reliability of theoretical assumptions, such as constant friction force magnitude and uniform rotation, with a time step of 0.01 s. Additionally, by introducing an equivalent spring effect to modify the TDFC model, it is observed that the modified TDFC effectively simulates the behavior of the Friction Pendulum Bearing (FPB) under horizontal bidirectional seismic loads. Compared to the ABAQUS finite element model of FPB, the modified TDFC model provides a more efficient approach for rapidly and accurately calculating the response of various FPB types under horizontal bidirectional seismic loads by adjusting parameter values of the modified TDFC model.

Close Binary Stars Modeled by Two Prolate Ellipsoids in Synchronous Rotation

The presence of tidal deformations in close binary stars has already been confirmed by astronomical observations. The present paper aims to simply address an astronomy problem, studying the relative movement of close binaries disturbed by their mutual deformation through some basic concepts and tools of celestial mechanics. For this purpose, the tidal effect is modeled by considering that each star is an elongated revolution ellipsoid in such a way that axes of revolution are coincident, and their largest axes point toward each other along the motion. The potential for mutual attraction is then obtained, resulting in a perturbed Keplerian system with perturbation proportional to the inverse of the cubic distance between the stars, thus being a one-degree-of-freedom problem and, therefore, integrable. The effective potential, the integrals of energy and angular momentum, and the Laplace vector are used to obtain qualitative information about the dynamics before integrating it. The motion describes a rosette-like orbit with periodic osculating elements, or a circle when the energy is a local minimum. Finally, an analytical solution is presented in terms of elliptic functions by using a regularizing and linearizing function.

Factorization of covariant Feynman graphs for the effective action

We prove a neat factorization property of Feynman graphs in covariant perturbation theory. The contribution of the graph to the effective action is written as an integral over Schwinger parameters whose integrand is a product of a massless scalar momentum integral that only depends on the basic graph topology, and a background-field dependent piece that contains all the information of spin, gauge representations, masses etc. We give a closed expression for the momentum integral in terms of four graph polynomials whose properties we derive in some detail. Our results can also be useful for standard (non-covariant) perturbation theory.

The cosmological tree theorem

A number of diagrammatic “cutting rules” have recently been developed for the wavefunction of the Universe which determines cosmological correlation functions. These leverage perturbative unitarity to relate particular “discontinuities” in Feynman-Witten diagrams (with cosmological boundary conditions) to simpler diagrams, in much the same way that the Cutkosky rules relate different scattering amplitudes. In this work, we make use of a further causality condition to derive new cutting rules for Feynman-Witten diagrams on any time-dependent spacetime background. These lead to the cosmological analogue of Feynman’s tree theorem for amplitudes, which can be used to systematically expand any loop diagram in terms of (momentum integrals of) tree-level diagrams. As an application of these new rules, we show that certain singularities in the wavefunction cannot appear in equal-time correlators due to a cancellation between “real” and “virtual” contributions that closely parallels the KLN theorem. Finally, when combined with the Bunch-Davies condition that certain unphysical singularities are absent, these cutting rules completely determine any tree-level exchange diagram in terms of simpler contact diagrams. Altogether, these results remove the need to ever perform nested time integrals when computing cosmological correlators.

Mode transition mechanisms in semi-passive flapping-wing propulsion or energy extraction

To facilitate the smooth transition of operational modes in self-sustaining micro-thrusters, this study focuses on a device employing semi-passive flapping wings. Transient numerical methods are used to analyze the variation of flap energy acquisition and propulsion characteristics with pitching frequency, and the relationship between the dominant vortex structure and the force characteristics is analyzed by using the integral momentum theorem and the dynamic mode decomposition (DMD) method. The results indicate that, under the investigated operational conditions, with an increase in the pitching frequency, distinct wake evolution characteristics were observed. In the energy harvesting operational regime, the wake patterns manifest as 2P + mS, 2S + mS, and mS types. During the transitional phase from energy harvesting to propulsion, the wake patterns shift from 2S + mS to 2S transitional types, eventually leading to the manifestation of a reverse Karman vortex street (2S RBVK). In the propulsion operational regime, the wake patterns consist of a reverse Karman vortex street and asymmetric reverse Karman vortex street phenomena. Simultaneously, it was observed that the transition of flapping-wing performance from energy harvesting to propulsion conditions delays behind the transformation of vortex street structures. This delay is attributed to the necessity for the flapping-wing device to overcome its own resistance before generating a net effective propulsive force. The contributions of unsteady wake to thrust primarily encompass vortex thrust, thrust due to localized fluid acceleration, induced momentum force, and induced pressure force. As the pitching frequency increases, the influence of wake vortices on propulsion also intensifies. The contribution of wake vortices to flapping-wing propulsion is determined by the spatial distribution of Lamb vectors and localized fluid accelerations. The conclusions drawn from the dynamic modal analysis and reconstructed flow field analysis of wake vortices align with the findings of the investigation of wake vortices based on the integral momentum theorem.

Theoretical investigations of positron collisions with phosphorus-containing compounds

A theoretical investigation of positron scattering from phosphorus-containing compounds (viz., PH3, P2H4, PCl3, PF3, PBr3, POF3, POCl3, and H2PO4) is reported in this article. The quantum mechanical potential scattering approach is utilized to calculate integral elastic, excitation, momentum transfer, direct ionization, positronium formation, total ionization, inelastic, differential, and total cross sections on a fine energy grid from 1 to 5000 eV. The ionization contribution in the inelastic scattering is estimated using the complex scattering potential-ionization contribution technique. Prior research on positron collisions with these targets is scarce; as a result, the purpose of this study is to make up, at least in part, for this deficiency in cross-section data. In addition to being pertinent to positron transport analyses, such as Monte Carlo methods, the current results should be useful to benchmark the accuracy and validity of positron molecule collision computations and, more significantly, to compare these calculations with related electron scattering outcomes. Furthermore, the calculated cross sections of PH3 are compared with NH3 and other phosphorus-containing compounds. The analysis makes it abundantly evident that the atoms on the periphery of a molecule have a substantially larger impact on the scattering process than the central atom. To analyze the scattering dynamics of positrons and their anti-particle electrons, a comparative study of cross sections of H2PO4 and H2SO4 is also presented. For most of these targets, positron calculations are carried out for the first time.

An integral representation for quantum amplitudes

The central impediment to reducing multidimensional integrals of transition amplitudes to analytic form, or at least to a fewer number of integral dimensions, is the presence of magnitudes of coordinate vector differences (square roots of polynomials) in disjoint products of functions. Fourier transforms circumvent this by introducing a three-dimensional momentum integral for each of those products, followed in many cases by another set of integral representations to move all of the resulting denominators into a single quadratic form in one denominator whose square my be completed. Gaussian transforms introduce a one-dimensional integral for each such function while squaring the square roots of coordinate vector differences and moving them into a common exponential. Addition theorems may also be used for extracting the angular variables, and sometimes direct integration is even possible. Each method has its strengths and weaknesses. An integral representation is derived herein that stands as an alternative to these four approaches. A number of consequent integrals of Macdonald functions, hypergeometric functions, and Meijer G-functions with complicated arguments are given.

A unified thermodynamic/Lamb-vector-based analysis of the aerodynamic force

The Lamb vector, cross-product of flow vorticity and velocity, is at the basis of different far-field methods developed in the last decades for the aerodynamic force analysis and decomposition, as an alternative to the nowadays well-assessed thermodynamic methods. We here propose a mixed approach, where exact Lamb-vector-based force formulas are used in combination with a thermodynamic-based calculation of the Lamb vector through Crocco's equation. In computational fluid dynamics, this way of calculating the Lamb vector, therefore, inherits from the numerical form of the flow momentum equation and discretely satisfies the local (and integral) momentum balance on which far-field methods rely. The resulting hybrid method, which does not require an explicit vorticity calculation, provides results in far better agreement with regard to near-field force data when compared to standard vorticity-based approaches, especially in the presence of shock waves, where inaccuracies of domain integrals involving the Lamb vector were systematically reported by different authors. In addition, it overcomes the limitations of previous thermodynamic methods, which only compute the drag force.