scholarly journals Filtering Mixing Planes for Low Reduced Frequency Analysis of Turbomachines

Author(s):  
G. Pullan ◽  
J. J. Adamczyk

A class of problems in turbomachinery is characterised by unsteady interactions at low reduced frequencies. These interactions are often the result of perturbations with length-scale on the order of the machine circumference and examples include axial compressors operating with inlet distortion, fans with downstream pylons, and turbine rotors downstream of mid-frame struts. Typically, this unsteadiness is accompanied by higher frequency fluctuations caused by perturbations with a length-scale on the order of a blade pitch. Conventional numerical analysis of this class of problem requires computations with a time step governed by the high frequency content but a greatly reduced run time could be achieved if the time step were dictated solely by the low reduced frequency, long length-scale, interaction of interest. In this paper, a filtering mixing plane technique is proposed that removes unwanted short length-scale perturbations at the interfaces between blade rows. This approach gives the user control over the amount of mixing that occurs at these interfaces with the limits being fully mixed-out to pitchwise uniformity (conventional mixing plane) or no mixing (conventional sliding plane). By choosing to retain only enough harmonics to resolve the low reduced frequency interaction of interest, an order of magnitude reduction in run time can be achieved.

2017 ◽  
Vol 139 (9) ◽  
Author(s):  
G. Pullan ◽  
J. J. Adamczyk

A class of problems in turbomachinery is characterized by unsteady interactions at low reduced frequencies. These interactions are often the result of perturbations with length-scale on the order of the machine circumference and examples include axial compressors operating with inlet distortion, fans with downstream pylons, and turbine rotors downstream of midframe struts. Typically, this unsteadiness is accompanied by higher frequency fluctuations caused by perturbations with a length-scale on the order of a blade pitch. Conventional numerical analysis of this class of problem requires computations with a time step governed by the high-frequency content but a greatly reduced run time could be achieved if the time step was dictated solely by the low reduced frequency, long length-scale, interaction of interest. In this paper, a filtering mixing plane technique is proposed that removes unwanted short length-scale perturbations at the interfaces between blade rows. This approach gives the user control over the amount of mixing that occurs at these interfaces with the limits being fully mixed-out to pitchwise uniformity (conventional mixing plane) or no mixing (conventional sliding plane). By choosing to retain only enough harmonics to resolve the low reduced frequency interaction of interest, an order of magnitude reduction in run time can be achieved.


Author(s):  
A. Grimaldi ◽  
V. Michelassi

This paper discusses the impact of inlet flow distortions on centrifugal compressors based upon a large experimental data base in which the performance of several impellers in a range of corrected flows and corrected speeds have been measured after been coupled with different inlet plenums technologies. The analysis extends to centrifugal compressor inlets including a side stream, typical of liquefied natural gas applications. The detailed measurements allow a thorough characterization of the flow field and associated performance. The results suggest that distortions can alter the head by as much as 3% and efficiency of around 1%. A theoretical analysis allowed to identify the design features that are responsible for this deviation. In particular, an extension of the so-called “reduced-frequency,” a coefficient routinely used in axial compressors and turbine aerodynamics to weigh the unsteadiness generated by upstream to downstream blade rows, allowed to determine a plenum-to-impeller reduced frequency that correlates very well with the measured performance. The theory behind the new coefficient is discussed together with the measurement details and validates the correlation that can be used in the design phase to determine the best compromise between the inlet plenum complexity and impact on the first stage.


Author(s):  
A. Grimaldi ◽  
V. Michelassi

This paper discusses the impact of inlet flow distortions on centrifugal compressors based upon a large experimental data base in which the performance of several impellers in a range of corrected flows and corrected speeds have been measured after been coupled with different inlet plenums technologies. The analysis extends to centrifugal compressor inlets including a side stream, typical of Liquefied Natural Gas (LNG) applications. The detailed measurements allow a thorough characterization of the flow field and associated performance. The results suggest that distortions can alter the head by as much as 3% and efficiency of around 1%. A theoretical analysis allowed to identify the design features that are responsible for this deviation. In particular, an extension of the so-called “reduced-frequency”, a coefficient routinely used in axial compressors and turbine aerodynamics to weigh the unsteadiness generated by upstream to downstream blade rows, allowed to determine a plenum-to-impeller reduced frequency that correlates very well with the measured performance. The theory behind the new coefficient is discussed together with the measurement details, and validate the correlation that can be used in the design phase to determine the best compromise between the inlet plenum complexity and impact on the first stage.


2021 ◽  
Vol 11 (4) ◽  
pp. 1399
Author(s):  
Jure Oder ◽  
Cédric Flageul ◽  
Iztok Tiselj

In this paper, we present uncertainties of statistical quantities of direct numerical simulations (DNS) with small numerical errors. The uncertainties are analysed for channel flow and a flow separation case in a confined backward facing step (BFS) geometry. The infinite channel flow case has two homogeneous directions and this is usually exploited to speed-up the convergence of the results. As we show, such a procedure reduces statistical uncertainties of the results by up to an order of magnitude. This effect is strongest in the near wall regions. In the case of flow over a confined BFS, there are no such directions and thus very long integration times are required. The individual statistical quantities converge with the square root of time integration so, in order to improve the uncertainty by a factor of two, the simulation has to be prolonged by a factor of four. We provide an estimator that can be used to evaluate a priori the DNS relative statistical uncertainties from results obtained with a Reynolds Averaged Navier Stokes simulation. In the DNS, the estimator can be used to predict the averaging time and with it the simulation time required to achieve a certain relative statistical uncertainty of results. For accurate evaluation of averages and their uncertainties, it is not required to use every time step of the DNS. We observe that statistical uncertainty of the results is uninfluenced by reducing the number of samples to the point where the period between two consecutive samples measured in Courant–Friedrichss–Levy (CFL) condition units is below one. Nevertheless, crossing this limit, the estimates of uncertainties start to exhibit significant growth.


2018 ◽  
Vol 64 (247) ◽  
pp. 745-758 ◽  
Author(s):  
E. DE ANDRÉS ◽  
J. OTERO ◽  
F. NAVARRO ◽  
A. PROMIŃSKA ◽  
J. LAPAZARAN ◽  
...  

ABSTRACTWe have developed a two-dimensional coupled glacier–fjord model, which runs automatically using Elmer/Ice and MITgcm software packages, to investigate the magnitude of submarine melting along a vertical glacier front and its potential influence on glacier calving and front position changes. We apply this model to simulate the Hansbreen glacier–Hansbukta proglacial–fjord system, Southwestern Svalbard, during the summer of 2010. The limited size of this system allows us to resolve some of the small-scale processes occurring at the ice–ocean interface in the fjord model, using a 0.5 s time step and a 1 m grid resolution near the glacier front. We use a rich set of field data spanning the period April–August 2010 to constrain, calibrate and validate the model. We adjust circulation patterns in the fjord by tuning subglacial discharge inputs that best match observed temperature while maintaining a compromise with observed salinity, suggesting a convectively driven circulation in Hansbukta. The results of our model simulations suggest that both submarine melting and crevasse hydrofracturing exert important controls on seasonal frontal ablation, with submarine melting alone not being sufficient for reproducing the observed patterns of seasonal retreat. Both submarine melt and calving rates accumulated along the entire simulation period are of the same order of magnitude, ~100 m. The model results also indicate that changes in submarine melting lag meltwater production by 4–5 weeks, which suggests that it may take up to a month for meltwater to traverse the englacial and subglacial drainage network.


2011 ◽  
Vol 689 ◽  
pp. 32-74 ◽  
Author(s):  
C.-K. Kang ◽  
H. Aono ◽  
C. E. S. Cesnik ◽  
W. Shyy

AbstractEffects of chordwise, spanwise, and isotropic flexibility on the force generation and propulsive efficiency of flapping wings are elucidated. For a moving body immersed in viscous fluid, different types of forces, as a function of the Reynolds number, reduced frequency (k), and Strouhal number (St), acting on the moving body are identified based on a scaling argument. In particular, at the Reynolds number regime of $O(1{0}^{3} \ensuremath{-} 1{0}^{4} )$ and the reduced frequency of $O(1)$, the added mass force, related to the acceleration of the wing, is important. Based on the order of magnitude and energy balance arguments, a relationship between the propulsive force and the maximum relative wing-tip deformation parameter ($\gamma $) is established. The parameter depends on the density ratio, St, k, natural and flapping frequency ratio, and flapping amplitude. The lift generation, and the propulsive efficiency can be deduced by the same scaling procedures. It seems that the maximum propulsive force is obtained when flapping near the resonance, whereas the optimal propulsive efficiency is reached when flapping at about half of the natural frequency; both are supported by the reported studies. The established scaling relationships can offer direct guidance for micro air vehicle design and performance analysis.


2016 ◽  
Vol 3 (1) ◽  
pp. 40
Author(s):  
Uriel Goldberg

A turbulence closure based on transport equations for the square-root of the kinetic energy of turbulence, q=k1/2 and the length-scale, , is proposed and tested. The model is topography parameter free (no wall distance needed), uses local wall proximity indicators instead, and is meant to be applicable to both wall-bounded and free shear flows. Solving directly for the turbulence length-scale, invoking Dirichlet boundary conditions for both q and  and the fact that q varies linearly across the viscous sublayer contribute to reduced sensitivity of this model to near-wall grid concentration (as long as the sublayer is resolved) and to less numerical stiffness, hence faster convergence. A variable Cm parameter is featured in this model to account for non-simple shear where mean strain and vorticity rates are different. Several cases, covering a wide variety of flows, are presented to demonstrate the model’s performance. Fluids engineers whose work involves complex 3D topologies, particularly with non-stationary grids which require re-computing wall distance arrays at each time-step (a heavy demand on time and budget) may appreciate the fact that no distance arrays are needed for the q-  model.


2020 ◽  
Vol 495 (4) ◽  
pp. 4306-4313 ◽  
Author(s):  
Michael Y Grudić ◽  
Philip F Hopkins

ABSTRACT We describe a new adaptive time-step criterion for integrating gravitational motion, which uses the tidal tensor to estimate the local dynamical time-scale and scales the time-step proportionally. This provides a better candidate for a truly general-purpose gravitational time-step criterion than the usual prescription derived from the gravitational acceleration, which does not respect the equivalence principle, breaks down when $\boldsymbol {a}=0$, and does not obey the same dimensional scaling as the true time-scale of orbital motion. We implement the tidal time-step criterion in the simulation code gizmo, and examine controlled tests of collisionless galaxy and star cluster models, as well as galaxy merger simulations. The tidal criterion estimates the dynamical time faithfully, and generally provides a more efficient time-stepping scheme compared to an acceleration criterion. Specifically, the tidal criterion achieves order-of-magnitude smaller energy errors for the same number of force evaluations in potentials with inner profiles shallower than ρ ∝ r−1 (i.e. where $\boldsymbol {a}\rightarrow 0$), such as star clusters and cored galaxies. For a given problem these advantages must be weighed against the additional overhead of computing the tidal tensor on-the-fly, but in many cases this overhead is small.


2015 ◽  
Vol 45 (10) ◽  
pp. 2497-2521 ◽  
Author(s):  
Benjamin D. Mater ◽  
Subhas K. Venayagamoorthy ◽  
Louis St. Laurent ◽  
James N. Moum

AbstractOceanic density overturns are commonly used to parameterize the dissipation rate of turbulent kinetic energy. This method assumes a linear scaling between the Thorpe length scale LT and the Ozmidov length scale LO. Historic evidence supporting LT ~ LO has been shown for relatively weak shear-driven turbulence of the thermocline; however, little support for the method exists in regions of turbulence driven by the convective collapse of topographically influenced overturns that are large by open-ocean standards. This study presents a direct comparison of LT and LO, using vertical profiles of temperature and microstructure shear collected in the Luzon Strait—a site characterized by topographically influenced overturns up to O(100) m in scale. The comparison is also done for open-ocean sites in the Brazil basin and North Atlantic where overturns are generally smaller and due to different processes. A key result is that LT/LO increases with overturn size in a fashion similar to that observed in numerical studies of Kelvin–Helmholtz (K–H) instabilities for all sites but is most clear in data from the Luzon Strait. Resultant bias in parameterized dissipation is mitigated by ensemble averaging; however, a positive bias appears when instantaneous observations are depth and time integrated. For a series of profiles taken during a spring tidal period in the Luzon Strait, the integrated value is nearly an order of magnitude larger than that based on the microstructure observations. Physical arguments supporting LT ~ LO are revisited, and conceptual regimes explaining the relationship between LT/LO and a nondimensional overturn size are proposed. In a companion paper, Scotti obtains similar conclusions from energetics arguments and simulations.


2019 ◽  
Vol 622 ◽  
pp. A79 ◽  
Author(s):  
Mika Juvela

Context. Thermal dust emission carries information on physical conditions and dust properties in many astronomical sources. Because observations represent a sum of emission along the line of sight, their interpretation often requires radiative transfer (RT) modelling. Aims. We describe a new RT program, SOC, for computations of dust emission, and examine its performance in simulations of interstellar clouds with external and internal heating. Methods. SOC implements the Monte Carlo RT method as a parallel program for shared-memory computers. It can be used to study dust extinction, scattering, and emission. We tested SOC with realistic cloud models and examined the convergence and noise of the dust-temperature estimates and of the resulting surface-brightness maps. Results. SOC has been demonstrated to produce accurate estimates for dust scattering and for thermal dust emission. It performs well with both CPUs and GPUs, the latter providing a speed-up of processing time by up to an order of magnitude. In the test cases, accelerated lambda iterations (ALIs) improved the convergence rates but was also sensitive to Monte Carlo noise. Run-time refinement of the hierarchical-grid models did not help in reducing the run times required for a given accuracy of solution. The use of a reference field, without ALI, works more robustly, and also allows the run time to be optimised if the number of photon packages is increased only as the iterations progress. Conclusions. The use of GPUs in RT computations should be investigated further.


Sign in / Sign up

Export Citation Format

Share Document