Velocity Moments Research Articles

Quasi-linear Analysis of Proton-cyclotron Instability

The proton-cyclotron (PC) instability operates in various space plasma environments. In the literature, the so-called velocity moment-based quasi-linear theory is employed to investigate the physical process of PC instability that takes place after the onset of early linear exponential growth. In this method, the proton velocity distribution function (VDF) is assumed to maintain a bi-Maxwellian form for all time, which substantially simplifies the analysis, but its validity has not been rigorously examined by comparing against the actual solution of the kinetic equation. The present paper relaxes the assumption of the velocity moment-based quasi-linear theory by actually solving for the velocity space diffusion equation under the assumption of separable perpendicular and parallel VDFs, and upon comparison with the simplified velocity moment theory, it demonstrates that the simplified method is largely valid, despite the fact that the method slightly overemphasizes the relaxation of temperature anisotropy when the system is close to the marginally stable state. The overall validation is further confirmed with the results of particle-in-cell and hybrid-code simulations. The present paper thus provides a justification for making use of the velocity moment-based quasi-linear theory as an efficient first-cut theoretical tool for the PC instability.

Chronic Pain Influences Lower Extremity Energetics During Landing Cutting in Patients with Chronic Ankle Instability.

Chronic ankle instability (CAI) patients exhibit altered movement patterns during jump landing/cutting movements. Persistent pain is one of the residual symptoms that may affect movements. Calculating joint energetics affected by chronic pain offers a novel method to understand how chronic pain influences energetics of lower extremity joints in CAI patients. To identify the effects of chronic pain on lower extremity energy dissipation and generation during jump landing and cutting in CAI patients. Cross-sectional Study. Laboratory. Fifteen CAI patients with higher pain (6 males, 9 females; age=22.1±2.1year; height=1.74±0.09m; mass=71.3±10.6kg, pain=66.9±9.4), 15 CAI patients with lower pain (6 males, 9 females; age=22.3±2.1year; height=1.74±0.08m; mass=70.1±10.7kg, pain=89.3±2.6), and 15 healthy controls (6 males, 9 females; age=21.3±1.7year; height=1.73±0.08m; mass=70±10.3kg, pain=100±0). Ground reaction force data were collected during 5 trials of maximal jump landing/cutting tasks. Joint power was defined as the product of angular velocity and joint moment. Energy dissipation and generation by the ankle, knee, and hip joints were calculated by integrating regions of the joint power curve. CAI patients with higher pain displayed less ankle energy dissipation (p=.013 and p=.018) and generation in the ankle (p=.002 and p=.028) than CAI patients with lower pain and healthy controls during the jump landing/cutting phase. CAI patients with higher pain showed more hip energy generation than CAI patients with lower pain (P=.038) and healthy controls (P=.013) during the cutting phase. CAI patients with higher pain changed both energy dissipation and generation in the lower extremities, reducing the burden of the ankle joint during jump landing/cutting and having a hip dominant compensatory strategy during the cutting phase. Our results suggest that chronic pain could be one of the factors that affect motor strategies in the CAI population.

The Role of the Lead Hip in Collegiate Baseball Pitching: Implications for Ball Velocity and Upper-Extremity Joint Moments.

Hip flexibility is an important biomechanical factor for a baseball pitcher. However, there have been limited investigations into the association between upper-extremity joint stresses and ball velocity and hip flexibility, as assessed via motion patterns during the pitch. The purpose of this study was to provide a detailed kinematic description of the lead hip during the pitch and determine the association between lead hip motion and both ball velocity and the elbow varus moment. This study was a secondary analysis of the kinematic and kinetic data previously collected on 99 collegiate-level baseball pitchers using standard optoelectronic motion capture. Significant associations were noted between lead hip internal rotation and both peak ball velocity and the elbow varus moment. The data indicated that for every 10° increase in internal lead hip rotation, ball velocity increased by 0.6m/s (P < .001, r2 = .26), and the elbow varus moment increased by 5N·m (P < .001, r2 = .33). The results of this study suggested that internal hip rotation may be an important means of identifying pitchers that may be at risk for future injury.

The Lower Limit of Dynamical Black Hole Masses Detectable in Virgo Compact Stellar Systems Using the JWST/NIRSpec IFU

Due to observational challenges, the mass function of black holes (BH) at lower masses is poorly constrained in the local universe. Understanding the occupation fraction of BHs in low-mass galaxies is crucial for constraining the origins of supermassive BH seeds. Compact stellar systems (CSSs), including ultracompact dwarf galaxies (UCDs) and compact elliptical galaxies (cEs), are potential intermediate-mass BH hosts. Despite the difficulties posed by their limited spheres of influence, stellar dynamical modeling has been effective in estimating central BH masses in CSSs. Some CSSs may harbor a BH constituting up to 20% of their host stellar mass, while others might not have a central BH. In support of our ongoing efforts to determine the BH masses in select CSSs in the Virgo cluster using JWST/NIRSpec IFU observations and orbit-superposition dynamical models, we create mock kinematic data mimicking the characteristics of observed cEs/UCDs in the Virgo cluster with different BH masses. We then construct a series of dynamical models using the orbit-superposition code FORSTAND and explore the accuracy of recovering the BH mass. We find that the mass of BHs comprising 1% or more of the total host stellar mass can be accurately determined through kinematic maps featuring higher-order velocity moments. We also assess how BH mass measurement is affected by deprojection methods, regularization factors, anisotropy parameters, orbit initial conditions, the absence of higher-order velocity moments, the spatial resolution, and the signal-to-noise ratio.

The effect of stress barriers on unconventional-singularity-driven frictional rupture

Whether or not energy dissipation is localized in the vicinity of the rupture tip, and whether any distal energy dissipation far from the crack tip has a significant influence on rupture dynamics are key questions in the description of frictional ruptures, in particular regarding the application of Linear Elastic Fracture Mechanics (LEFM) to earthquakes. These questions are investigated experimentally using a 40-cm-long experimental frictional interface. Three independent pistons apply a normal load with a fourth piston applying a shear load, enabling the application of a heterogeneous stress state and stress barriers. After loading the frictional interface to a near-critical state, subsequent unloading of one normal-load piston leads to dynamic ruptures which propagate into the heterogeneous stress fields. The ruptures in these experiments are found to be driven by unconventional singularities, characterized by an ever-increasing breakdown work with slip, and as a result do not conform to the assumptions of LEFM. As these experimental stress barriers inhibit slip, they therefore also reduce the breakdown work occurring outside of the cohesive zone. It is shown that this distal weakening, far from the crack tip, must be considered for the accurate prediction of rupture arrest length. These experiments are performed in the context of a proposed stimulation technique for Enhanced Geothermal Systems (EGSs). It has previously been suggested, through theoretical arguments, that stress barriers could be induced through the manipulation of pore pressure such that there is reduced seismic hazard during the shear stimulation of EGSs. This stimulation technique, known as preconditioning, is demonstrated here to reduce the mechanical energy flux to the crack tip, G, while also increasing the fracture energy, Gc. Preconditioning is shown to be capable of arresting seismic rupture and reducing co-seismic slip, slip velocity, and seismic moment at preconditioning stresses which are reasonably achievable in the field. Due to the fully-coupled nature of seismic rupture and fault slip, preconditioning also reduces distal weakening and its contribution to the propagation of induced seismic ruptures. In a similar vein, heterogeneous pore pressure fields associated with some seismic swarms can be used to explain changes in stress drop within the swarm without recourse to material or total-stress heterogeneity.

Robust optimization of geometrical properties of flow diverter stents for treating cerebral aneurysm: A proof-of-concept study

This study presents a novel approach to optimize the design of flow diverter (FD) stents for cerebral aneurysm (CA) treatment. By addressing sources of uncertainty in cardiovascular simulations, including geometrical and physical properties and boundary conditions, we aim to assess the applicability of robust optimization techniques to the FD design, establishing a foundation for acquiring robust FDs that are capable of operating consistently under various real-world scenarios. Blood flow in a simplified 2-dimensional CA and FD model was simulated through computational fluid dynamics (CFD). A design space exploration method, incorporating Latin hypercube sampling and Kriging surrogate models, was employed to obtain the optimal solution. The objective was to maximize the reduction in velocity and vorticity within the CA sac. This study used non-intrusive polynomial chaos expansion (PCE) to quantify and propagate the input uncertainties through the computational model and compute the statistical moments of velocity and vorticity reductions. To assess the effect of uncertain sources on objective functions, a sensitivity analysis method based on Sobol indices was applied. Robust optimization involved simultaneously optimizing the mean and standard deviation of velocity reduction. Additionally, we accounted for patients’ specific conditions and repeated the robust optimization. The results indicate that blood Hematocrit and inlet velocity are the most impactful uncertain sources in FD optimization. Moreover, the obtained Pareto front shows that in robust designs, FD struts are concentrated in the distal region of the CA neck, while optimal designs have more struts in the proximal region. This study proposes an FD that compromises robustness and optimality with a velocity reduction of 72.31 % and a standard deviation of 0.00343.

Hemodynamics in the treatment of pseudoaneurysm caused by extreme constriction of aortic arch with coated stent.

To evaluate the changes in distal vascular morphology and hemodynamics in patients with extremely severe aortic coarctation (CoA) after covered palliative (CP) stent dilation with different surgical strategies. Perioperative computed tomography angiography and digital subtraction angiography were utilized to construct three aortic models with varying stenosis rates and one follow-up model in a patient with extremely severe CoA. The models included: an idealized non-stenosed model (A: 0%), a model post initial stent deployment (B: 28%), a model post balloon expansion (C: 39%), and a model 18 months after post-balloon expansion (D: 39%). Consistent boundary conditions were applied to all models, and hemodynamic simulation was conducted using the pure fluid method. The narrowest and distal diameter of the stent increased by 34.71% and 59.29%, respectively, from model B to C. Additionally, the distal diameter of the stent increased by -13.80% and +43.68% compared to the descending aorta diameter, respectively. Furthermore, the ellipticity of the maximum cross-section of the aneurysm region in model A to D continued to increase. The oscillatory shear index at the stenosis to the region of the aneurysm were found to be higher in Models A and B, and lower in Models C and D. At the moment of maximum flow velocity, the blood flow distribution in models A and B was more uniform in the widest section of the blood vessels at the distal end of the stenosis, whereas models C and D exhibited disturbed blood flow with more than 2 eddy currents. The time-averaged wall shear stress (TAWSS) decreased in the distal and basal aneurysms, while it significantly increased at the step position. The aneurysmal region exhibited an endothelial cell activation potential value lower than 0.4 Pa-1. In patients with extremely severe CoA, it is crucial to ensure that the expanded diameter at both ends of the CP stent does not exceed the native vascular diameter during deployment. Our simulation results demonstrate that overdilation leads to a decrease in the TAWSS above the injured vessel, creating an abnormal hemodynamic environment that may contribute to the development and enlargement of false aneurysms in the early postoperative period. ClinicalTrials.gov, (NCT02917980).

Generalized autocorrelation function in the family of deterministic and stochastic anomalous diffusion processes

We investigate the observables of the one-dimensional model for anomalous transport in semiconductor devices where diffusion arises from scattering at dislocations at fixed random positions, known as the Lévy-Lorentz gas. To gain insight into the microscopic properties of such a stochastically complex system, deterministic dynamics known as the slicer map and fly-and-die dynamics are used. We analytically derive the generalized position autocorrelation function of these dynamics and study the special case, the 3-point position correlation function. For this, we derive single-parameter-dependent scaling and compare it with the numerically estimated 3-point position autocorrelation of the Lévy-Lorentz gas, for which the analytical expression is still an open question. Here we obtained a remarkable agreement between them, irrespective of any functional relationship with time. Moreover, we demonstrate that the position moments and the position autocorrelations of these systems scale in the same fashion, provided the times are large enough and far enough apart. Other observables, such as velocity moments and correlations, are reported to distinguish the systems. Published by the American Physical Society 2024

A Variant of the Local Similarity Theory and Approximations of Vertical Profiles of Turbulent Moments of the Atmospheric Convective Boundary Layer

The approximation of the turbulent moments of the atmospheric convective layer is based on a variant of the local similarity theory using the concepts of the semi-empirical theory of Prandtl turbulence. In the proposed variant of the local similarity theory, the second moment of vertical velocity and the “spectral” Prandtl mixing length are used as basic parameters. This approach allows us to extend Prandtl’s theory to turbulent moments of vertical velocity and buoyancy and additionally offer more than ten new approximations. The comparison of the proposed approximation with other variants of the theory of local similarity is considered. It is shown that the selected basic parameters significantly improve the agreement between the local similarity approximations and experimental data. The approximations are consistent with observations in the turbulent convective layer of the atmosphere, the upper boundary of which nearly corresponds to the lower boundary of the temperature inversion. Analytical approximations of local similarity can find applications in the construction of high-order moment closures in the vortex of resolving numerical turbulence models, as well as in the construction of “mass-flux” parametrization.

Turbulence scaling from deep learning diffusion generative models

Complex spatial and temporal structures are inherent characteristics of turbulent fluid flows and comprehending them poses a major challenge. This comprehension necessitates an understanding of the space of turbulent fluid flow configurations. We employ a diffusion-based generative model to learn the distribution of turbulent vorticity profiles and generate snapshots of turbulent solutions to the incompressible Navier-Stokes equations. We consider the inverse cascade in two spatial dimensions and generate diverse turbulent solutions that differ from those in the training dataset. We analyze the statistical scaling properties of the new turbulent profiles, calculate their structure functions, energy power spectrum, velocity probability distribution function and moments of local energy dissipation. All the learnt scaling exponents are consistent with the expected Kolmogorov scaling This agreement with established turbulence characteristics provides strong evidence of the model's capability to capture essential features of real-world turbulence.

Exact moments for trapped active particles: inertial impact on steady-state properties and re-entrance

In this study, we investigate the behavior of inertial active Brownian particles in a d-dimensional harmonic trap in the presence of translational diffusion. While the solution of the Fokker–Planck equation is generally challenging, it can be utilized to compute the exact time evolution of all time-dependent dynamical moments using a Laplace transform approach. We present the explicit form for several moments of position and velocity in d-dimensions. An interplay of time scales assures that the effective diffusivity and steady-state kinetic temperature depend on both inertia and trap strength, unlike passive systems. The distance from equilibrium, measured by the violation of equilibrium fluctuation-dissipation and the amount of entropy production, decreases with increasing inertia and trap strength. We present detailed ‘phase diagrams’ using kurtosis of velocity and position, showing possibilities of re-entrance to equilibrium.

Parallel expansion of a fuel pellet plasmoid

The problem of the assimilation of a cryogenic fuel pellet injected into a hot plasma is considered. Due to the transparency to ambient particles of the plasmoid, the localised region of high-density plasma created by ionisation of the ablated pellet material, electrons reach a ‘quasiequilibrium’ (QE) state which is characterised by a steady-state on the fastest collisional time scale. The simplified electron kinetic equation of the QE state is solved. Taking a velocity moment of the higher-order electron kinetic equation, which is valid on the expansion time scale, permits a fluid closure, yielding an evolution equation for the macroscopic parameters describing the QE distribution function. In contrast to the Braginskii equations, the closure does not require that electrons have a short mean free path compared with the size of density perturbations, and permits an anisotropic and highly non-Maxwellian distribution function. As the QE distribution function accounts for both trapped and passing electrons, the self-consistent electric potential that causes the expansion can be properly described, in contrast to earlier models of pellet plasmoid expansion with an unbounded potential. The plasmoid expansion is simulated using both a Vlasov model and a cold-fluid model for the ions. During the expansion plasmoid ions and electrons obtain nearly equal amounts of energy; as hot ambient electrons provide this energy in the form of collisional heating of plasmoid electrons, the expansion of a pellet plasmoid is expected to be a potent mechanism for the transfer of energy from electrons to ions on a time scale shorter than that of ion–electron thermalisation.

Microfluidic jet impacts on deep pools transition from capillary-dominated cavity closure to gas pressure-dominated closure at higher Weber numbers.

Studying liquid jet impacts on a liquid pool is crucial for various engineering and environmental applications. During jet impact, the free surface of the pool deforms and a cavity is generated. Simultaneously, the free surface of the cavity extends radially outward and forms a rim. Eventually the cavity collapses by means of gas inertia and surface tension. Our numerical investigation using an axisymmetric model in Basilisk C explores cavity collapse dynamics under different impact velocities and gas densities. We validate our model against theory and experiments across a previously unexplored parameter range. Our results show two distinct regimes in the cavity collapse mechanism. By considering forces pulling along the interface, we derive scaling arguments for the time of closure and maximum radius of the cavity, based on the Weber number. For jets with uniform constant velocity from tip to tail and We ⩽ 150 the cavity closure is capillary dominated and happens below the surface (deep seal). In contrast, for We ⩾ 180 the cavity closure happens above the surface (surface seal) and is dominated by the gas entrainment and the pressure gradient that it causes. Additionally, we monitor gas velocity and pressure throughout the impact process. This analysis reveals three critical moments of maximum gas velocity: before impact, at the instant of cavity collapse, and during droplet ejection following cavity collapse. Our results provide information for understanding pollutant transport during droplet impacts on large bodies of water, and other engineering applications, like additive manufacturing, lithography and needle-free injections.

Presheath-like structures and effusive particle losses for biased probes at and near I–V electron saturation

A theory for presheath-like structures near probes biased at and above the plasma potential is developed for collisionless plasmas with an electron-neutral mean free path on the order of the chamber scale. The theory predicts presheath-like perturbations to the plasma that result from the free streaming of electrons and an effusion loss process from the chamber at the electrode. For these situations, a loss-cone-like velocity distribution function for electrons is predicted where the loss angle of the depletion region corresponds to the angular size of the electrode at a specified distance. The angle of the loss cone becomes 180° at the sheath edge. In comparison with a previous collisional electron presheath model that required electrons satisfy a Bohm criterion at the sheath edge [Scheiner et al., Phys. Plasmas 22, 123520 (2015)], the present work suggests that no such condition is needed for collisionless low pressure plasmas in the ≲10 mTorr range. The theory predicts the generation of a density depletion of roughly 0.5ne and an electron velocity moment of tens of percent of the electron thermal speed by the sheath edge in a presheath with a potential drop of less than Ti/e. The range of this presheath perturbation is determined by the electrode geometry instead of the collisional mean free path. These predictions are tested against previously published particle in cell simulations and are found to be in good agreement.

Simulation of Shear and Tensile Fractures Using Ductile Phase Field Modelling with the Calibration of P Wave Velocity Measurement and Moment Tensor Inversion

The appropriate understanding and formulation of rock discontinuities via FEM is still challenging for rock engineering, as continuous algorithms cannot handle the discontinuities in rock mass. Also, different failure modes of rock samples, containing tensile and shear failure, need to be computed separately. In this study, a novel double-phase field damage model was introduced with two independent phase field damage variables. The construction of the proposed model follows the thermodynamics framework from the overall Helmholtz free energy, with elastic, plastic and surface damage components. The proposed model is calibrated via traditional damage variables, based on ultrasonic wave velocity measurement and acoustic emission monitoring, and both show great consistency between simulation results and laboratory observations. Then the double-phase field damage model is applied to COMSOL software to simulate microcrack propagation in a pre-fractured rock sample. Both lateral and wing cracks are observed in this study, manifested as shear- and tensile-dominated cracks. We also observed different microcracking mechanisms in the proposed numerical models, such as tensile and shear cracking, the influence of plastic strain and the percolation between tensile and shear microcracks. Overall, this study provides valuable insights into the mechanics of microcracking in rocks, and the proposed model shows promising results in simulating crack propagation.

Anisotropy ansatz for the axisymmetric Jeans equations

Abstract The Jeans equations do not form a closed system, and to solve them a parametrization relating the velocity moments is often adopted. For axisymmetric models, a phenomenological choice (the “b-ansatz”) is widely used for the relation between the vertical ($\sigma _z^2$) and radial ($\sigma _R^2$) components of the velocity dispersion tensor, thus breaking their identity present in two-integral systems. However, the way in which the ansatz affects the resulting kinematical fields can be quite complicated, so that the analysis of these fields is usually performed only after numerically computing them. We present here a general procedure to study the properties of the ansatz-dependent fields $\overline{v_{\varphi }^2}$, $\Delta =\overline{v_{\varphi }^2}- \sigma _z^2$ and $\Delta _R= \overline{v_{\varphi }^2}- \sigma _R^2$. Specifically, the effects of the b-ansatz can be determined before solving the Jeans equations once the behaviour over the (R, z)-plane of three easy-to-build ansatz-independent functions is known. The procedure also constrains the ansatz to exclude unphysical results (as a negative $\overline{v_{\varphi }^2}$). The method is illustrated by discussing the cases of three well-known galaxy models: the Miyamoto &amp; Nagai and Satoh disks, and the Binney logarithmic halo, for which the regions and the constraints on the ansatz values can be determined analytically; a two-component (Miyamoto &amp; Nagai plus logarithmic halo) model is also discussed.

Dynamic response of asymmetric suspended structure subjected to random wind excitation

SummarySuspended structures, regarded as the expression of structural beauty, have attracted the focus of many architects and engineers. With the complexity of dynamic responses of the asymmetric suspended structures subjected to the random alongside wind excitation, an innovative method was proposed through combining the finite element method (FEM), the complex mode method (CMM), the pseudo excitation method (PEM), and the quadratic decomposition method (QDM) of frequency response function. The method could obtain unified closed‐form solutions of spectral moment and variance of node displacement and velocity and acceleration of asymmetric structures and its correctness was verified by existing numerical approaches. The proposed method is applicable to the random wind‐induced vibration response analysis of diverse linear complex structures. Notably, its efficacy lies in its circumvention of numerical integration, endowing it with high computational efficiency and accuracy. Extended comparative studies and discussion were also performed on effect of suspended span and comparisons of normal framed structure and symmetry structure. The results showed that the suspended span presented positive relation to horizontal displacement and inverse tendency to horizontal acceleration of the whole structure and vertical and horizontal acceleration of the columns near to suspended part should also be considered in engineering design.

Exact Results for Non-Newtonian Transport Properties in Sheared Granular Suspensions: Inelastic Maxwell Models and BGK-Type Kinetic Model.

The Boltzmann kinetic equation for dilute granular suspensions under simple (or uniform) shear flow (USF) is considered to determine the non-Newtonian transport properties of the system. In contrast to previous attempts based on a coarse-grained description, our suspension model accounts for the real collisions between grains and particles of the surrounding molecular gas. The latter is modeled as a bath (or thermostat) of elastic hard spheres at a given temperature. Two independent but complementary approaches are followed to reach exact expressions for the rheological properties. First, the Boltzmann equation for the so-called inelastic Maxwell models (IMM) is considered. The fact that the collision rate of IMM is independent of the relative velocity of the colliding spheres allows us to exactly compute the collisional moments of the Boltzmann operator without the knowledge of the distribution function. Thanks to this property, the transport properties of the sheared granular suspension can be exactly determined. As a second approach, a Bhatnagar-Gross-Krook (BGK)-type kinetic model adapted to granular suspensions is solved to compute the velocity moments and the velocity distribution function of the system. The theoretical results (which are given in terms of the coefficient of restitution, the reduced shear rate, the reduced background temperature, and the diameter and mass ratios) show, in general, a good agreement with the approximate analytical results derived for inelastic hard spheres (IHS) by means of Grad's moment method and with computer simulations performed in the Brownian limiting case (m/mg→∞, where mg and m are the masses of the particles of the molecular and granular gases, respectively). In addition, as expected, the IMM and BGK results show that the temperature and non-Newtonian viscosity exhibit an S shape in a plane of stress-strain rate (discontinuous shear thickening, DST). The DST effect becomes more pronounced as the mass ratio m/mg increases.

Spectral characteristics of particle preferential concentration in turbulence computed by Eulerian and Lagrangian methods

This study conducts a comprehensive analysis of the Fourier spectra characteristics of particle preferential concentration in turbulent flows obtained by Eulerian and Lagrangian modelling approaches for different Stokes numbers. Particle preferential concentration is characterized by clusters and voids in the spatial particle distributions formed in turbulent flows and can significantly influence processes such as droplet coalescence, evaporation, and gas-particle reactions. The research primarily focuses on comparing the performance of two Eulerian models—one with and one without second-order velocity moments—and one Lagrangian model, which are used to predict particle dispersion. Key findings include the impact of artificial diffusion in Eulerian methods, the superiority of methods incorporating second-order velocity moments contributions at higher Stokes numbers, and the sensitivity of the Lagrangian method to Poisson statistics. Notably, the study reveals minimal variation in mean velocity spectra across different Stokes numbers, as opposed to the marked variations in particle concentration, momentum, and energy spectra. Grid resolution emerges as a crucial factor in enhancing spectral energy predictions in Eulerian methods. The research underscores the nuanced distinctions between Eulerian and Lagrangian methods in modelling preferential concentration, providing a detailed spectral comparison of particle concentration, velocity, momentum, and energy, highlighting the importance of method selection based on specific modelling needs.

Design Method for Impeller of Centrifugal Pump With Guide Vanes Based on Oseen Vortex

Abstract Centrifugal pumps with guide vanes (GVs) are widely used in a pump as turbine, energy storage pump station, and water diversion project. The objective of this study is to establish a design method of centrifugal pumps for matching impeller outlet and diffuser inlet to reduce hydraulic loss. To this end, a theoretical prediction model based on the fluid governing equation and the Oseen vortex model is proposed to predict the velocity moment downstream the impeller of the centrifugal pump. Then, a design method is established to design the impeller of the centrifugal pump with guide vanes. A centrifugal pump with a specific speed of 127 is used to validate the theoretical prediction model, and results of velocity moment show that the deviation between predicted and simulated results is below 0.5% in average. Finally, the design method is applied, and results show that the average efficiency of redesigned pump under the working conditions is 1.04% higher than that of baseline pump, which validates the reliability of proposed design method by theoretical prediction based on Oseen vortex. Analysis on velocity distribution and turbulence eddy dissipation shows that the design method based on Oseen vortex can effectively improve the flow pattern and pump performance.