Theoretical and Computational Fluid Dynamics最新文献

{"title":"An adjoint-based methodology for calculating manufacturing tolerances for natural laminar flow airfoils susceptible to smooth surface waviness","authors":"Mohammad Moniripiri,&nbsp;Pedro P. C. Brito,&nbsp;André V. G. Cavalieri,&nbsp;Ney R. Sêcco,&nbsp;Ardeshir Hanifi","doi":"10.1007/s00162-023-00681-y","DOIUrl":"10.1007/s00162-023-00681-y","url":null,"abstract":"<p>An adjoint-based method is presented for determining manufacturing tolerances for aerodynamic surfaces with natural laminar flow subjected to wavy excrescences. The growth of convective unstable disturbances is computed by solving Euler, boundary layer, and parabolized stability equations. The gradient of the kinetic energy of disturbances in the boundary layer (<i>E</i>) with respect to surface grid points is calculated by solving adjoints of the governing equations. The accuracy of approximations of <span>(Delta E)</span>, using gradients obtained from adjoint, is investigated for several waviness heights. It is also shown how second-order derivatives increase the accuracy of approximations of <span>(Delta E)</span> when surface deformations are large. Then, for specific flight conditions, using the steepest ascent and the sequential least squares programming methodologies, the waviness profile with minimum <span>(L2-)</span>norm that causes a specific increase in the maximum value of <i>N</i>- factor, <span>(Delta N)</span>, is found. Finally, numerical tests are performed using the NLF(2)-0415 airfoil to specify tolerance levels for <span>(Delta {N})</span> up to 2.0 for different flight conditions. Most simulations are carried out for a Mach number and angle of attack equal to 0.5 and <span>(1.25^{circ })</span>, respectively, and with Reynolds numbers between <span>(9times 10^6)</span> and <span>(15times 10^6)</span> and for waviness profiles with different ranges of wavelengths. Finally, some additional studies are presented for different angles of attack and Mach numbers to show their effects on the computed tolerances.\u0000</p>","PeriodicalId":795,"journal":{"name":"Theoretical and Computational Fluid Dynamics","volume":"38 1","pages":"15 - 37"},"PeriodicalIF":2.2,"publicationDate":"2023-12-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s00162-023-00681-y.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139056735","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Stability of supersonic boundary layer over an unswept wing with a parabolic airfoil","authors":"P. V. Chuvakhov,&nbsp;I. M. Ilyukhin,&nbsp;A. V. Fedorov","doi":"10.1007/s00162-023-00680-z","DOIUrl":"10.1007/s00162-023-00680-z","url":null,"abstract":"<p>Under the low-noise Mach 3 flight conditions for a supersonic passenger aircraft having unswept wings with a thin parabolic airfoil, laminar-turbulent transition is due to amplification of the first mode. Stability of a local self-similar boundary layer over such a wing is investigated both using the <span>(e^{N})</span> method in the framework of linear stability theory and direct numerical simulation (DNS). It is found that the instability amplitude should reach a maximum over the entire spectral range above the profiles of 2.5% and thicker. The locus of maximum appears at the trailing edge and moves to the leading edge as the profile becomes thicker, while the maximum amplitude decreases. The theoretical findings are supported by DNS of the linear wave packets propagating in the boundary layer. Significance of these results to the design of laminar supersonic wings is discussed.</p>","PeriodicalId":795,"journal":{"name":"Theoretical and Computational Fluid Dynamics","volume":"38 1","pages":"1 - 13"},"PeriodicalIF":2.2,"publicationDate":"2023-12-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139056889","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Inviscid modeling of unsteady morphing airfoils using a discrete-vortex method","authors":"Alfonso Martínez-Carmena,&nbsp;Kiran Ramesh","doi":"10.1007/s00162-023-00678-7","DOIUrl":"10.1007/s00162-023-00678-7","url":null,"abstract":"<p>A low-order physics-based model to simulate the unsteady flow response to airfoils undergoing large-amplitude variations of the camber is presented in this paper. Potential-flow theory adapted for unsteady airfoils and numerical methods using discrete-vortex elements are combined to obtain rapid predictions of flow behavior and force evolution. To elude the inherent restriction of thin-airfoil theory to small flow disturbances, a time-varying chord line is proposed in this work over which to satisfy the appropriate boundary condition, enabling large deformations of the camber line to be modeled. Computational fluid dynamics simulations are performed to assess the accuracy of the low-order model for a wide range of dynamic trailing-edge flap deflections. By allowing the chord line to rotate with trailing-edge deflections, aerodynamic loads predictions are greatly enhanced as compared to the classical approach where the chord line is fixed. This is especially evident for large-amplitude deformations.</p>","PeriodicalId":795,"journal":{"name":"Theoretical and Computational Fluid Dynamics","volume":"38 6","pages":"845 - 862"},"PeriodicalIF":2.2,"publicationDate":"2023-12-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s00162-023-00678-7.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138556801","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Investigation of Stokes flow in a grooved channel using the spectral method","authors":"Mainendra Kumar Dewangan","doi":"10.1007/s00162-023-00679-6","DOIUrl":"10.1007/s00162-023-00679-6","url":null,"abstract":"<p>Pressure-driven Newtonian fluid flow between grooved and flat surfaces is analysed with no-slip boundary conditions at walls. The effect of corrugation on the fluid flow is investigated using the mesh-free spectral method. The primary aim of the present work is to develop an asymptotic/semi-analytical theory for confined transverse flows to bridge the gap between the limits of thin and thick channels. The secondary aim is to calculate permeability with reference to the effect of wall corrugation (roughness) without the restriction of pattern amplitude. We performed mathematical modelling and evaluated the analytical solution for hydraulic permeability with respect to the flat channel. The Pad<span>(acute{e})</span> approximate is employed to improve the solution accuracy of an asymptotic model. The results elucidate that permeability always follows a decreasing trend with increasing pattern amplitude using the spectral approach at the long-wave and short-wave limits. The prediction of the spectral model is more accurate than the asymptotic-based model by Stroock et al. (Anal Chem 74(20):5306, 2002) and Pad<span>(acute{e})</span> approximate, regardless of the grooved depth and wavelength of the channel. The finite-element-based numerical simulation is also used to understand the usefulness of theoretical models. A very low computational time is required using the mesh-free spectral model as compared to the numerical study. The agreement between the present model and the fully resolved numerical results is gratifying. Regarding numerical values, we calculated the relative error for different theoretical models such as an asymptotic model, Pad<span>(acute{e})</span> approximate, and a mesh-free spectral model. The spectral model always predicts the maximum relative error as less than <span>(3 %)</span>, regardless of the large pattern amplitude and wavelength. In addition, the results of the molecular dynamic (MD) simulations by Guo et al. (Phys Rev Fluids 1(7):074102, 2016) and the theoretical model by Wang (Phys Fluids 15(5):1121, 2003) are found to be quantitatively compatible with the predictions of effective slip length from the spectral model in the thick channel limit.\u0000</p>","PeriodicalId":795,"journal":{"name":"Theoretical and Computational Fluid Dynamics","volume":"38 1","pages":"39 - 59"},"PeriodicalIF":2.2,"publicationDate":"2023-11-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135321602","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"GPU computing of yield stress fluid flows in narrow gaps","authors":"Ivonne Leonor Medina Lino,&nbsp;Mariana Carrasco-Teja,&nbsp;Ian Frigaard","doi":"10.1007/s00162-023-00674-x","DOIUrl":"10.1007/s00162-023-00674-x","url":null,"abstract":"<p>We present a Graphic Processing Units (GPU) implementation of non-Newtonian Hele-Shaw flow that models the displacement of Herschel-Bulkley fluids along narrow eccentric annuli. This flow is characteristic of many long-thin flows that require extensive calculation due to an inherent nonlinearity in the constitutive law. A common method of dealing with such flows is via an augmented Lagrangian algorithm, which is often painfully slow. Here we show that such algorithms, although involving slow iterations, can often be accelerated via parallel implementation on GPUs. Indeed, such algorithms explicitly solve the nonlinear aspects only locally on each mesh cell (or node), which makes them ideal candidates for GPUs. Combined with other advances, the optimized GPU implementation takes <span>(approx 2.5%)</span> of the time of the original algorithm.</p>","PeriodicalId":795,"journal":{"name":"Theoretical and Computational Fluid Dynamics","volume":"37 5","pages":"661 - 680"},"PeriodicalIF":3.4,"publicationDate":"2023-10-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"71909902","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Application of the lattice Boltzmann method to the study of ultrasound propagation and acoustic streaming in three-dimensional cavities: advantages and limitations","authors":"Jaouad Benhamou,&nbsp;Bjarne Vincent,&nbsp;Sophie Miralles,&nbsp;Mohammed Jami,&nbsp;Daniel Henry,&nbsp;Ahmed Mezrhab,&nbsp;Valéry Botton","doi":"10.1007/s00162-023-00676-9","DOIUrl":"10.1007/s00162-023-00676-9","url":null,"abstract":"<p>The paper presents a three-dimensional numerical study of the acoustic streaming induced by the dissipation of ultrasounds during their propagation in the air. The waves are generated by a circular acoustic source positioned at the center of the left wall of a parallelepipedic cavity. The simulations are performed with the lattice Boltzmann method associated with the D3Q19 multiple relaxation time model. A validation of this model is first performed by comparing the numerical and analytical acoustic intensities along the central axis of the acoustic source. The main objective of this study is to use two different methods to calculate the acoustic streaming flow. The first method is the direct calculation of the mean velocity fields as the mean values of the instantaneous velocities. The second method is an indirect technique, which first calculates the acoustic streaming force and then injects this force into the numerical code to produce the streaming. A comparison between the results obtained by the two methods was carried out and a good agreement was found between them. These different investigations, rather new in three-dimensional configurations, have allowed us to discuss the advantages and limitations of the lattice Boltzmann approach to simulate real situations of wave propagation and acoustic streaming.</p>","PeriodicalId":795,"journal":{"name":"Theoretical and Computational Fluid Dynamics","volume":"37 6","pages":"725 - 753"},"PeriodicalIF":3.4,"publicationDate":"2023-10-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"71909987","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Inverted stochastic lattice Boltzmann-Lagrangian model for identifying indoor particulate pollutant sources","authors":"Jinghong Qin","doi":"10.1007/s00162-023-00675-w","DOIUrl":"10.1007/s00162-023-00675-w","url":null,"abstract":"<p>This paper studies the inverted stochastic lattice Boltzmann-Lagrangian approach for identifying indoor particulate pollutant sources. The dynamics of the fluid (indoor air) as well as the transport of the particles in the Eulerian description are solved using the lattice Boltzmann method. The particles regard as rigid bodies, and the data interactions between lattice fluid and particle movement are implemented by calculating for interaction force and void fraction. Particle-wall collision process is based on the softball model which describes the dynamic characteristics of particles in microscopic state. The results are shown that the particle forward and inverted drifting paths and its mechanisms are investigated clearly than previous methods. Indoor particulate pollutant sources can exactly identify with this approach. This research can offer theoretical relevance to the modeling of multi-phase particle fluid.</p>","PeriodicalId":795,"journal":{"name":"Theoretical and Computational Fluid Dynamics","volume":"37 6","pages":"755 - 765"},"PeriodicalIF":3.4,"publicationDate":"2023-10-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"71909310","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Network-theoretic modeling of fluid–structure interactions","authors":"Aditya G. Nair,&nbsp;Samuel B. Douglass,&nbsp;Nitish Arya","doi":"10.1007/s00162-023-00673-y","DOIUrl":"10.1007/s00162-023-00673-y","url":null,"abstract":"<p>The coupling interactions between deformable structures and unsteady fluid flows occur across a wide range of spatial and temporal scales in many engineering applications. These fluid–structure interactions (FSI) pose significant challenges in accurately predicting flow physics. In the present work, two multi-layer network approaches are proposed that characterize the interactions between the fluid and structural layers for an incompressible laminar flow over a two-dimensional compliant flat plate at a 35<span>(^{circ })</span> angle of attack. In the first approach, the network nodes are formed by wake vortices and bound vortexlets, and the edges of the network are defined by the induced velocity between these elements. In the second approach, coherent structures (fluid modes), contributing to the kinetic energy of the flow, and structural modes, contributing to the kinetic energy of the compliant structure, constitute the network nodes. The energy transfers between the modes are extracted using a perturbation approach. Furthermore, the network structure of the FSI system is simplified using the community detection algorithm in the vortical approach and by selecting dominant modes in the modal approach. Network measures are used to reveal the temporal behavior of the individual nodes within the simplified FSI system. Predictive models are then built using both data-driven and physics-based methods. Overall, this work sets the foundation for network-theoretic reduced-order modeling of fluid–structure interactions, generalizable to other multi-physics systems.</p>","PeriodicalId":795,"journal":{"name":"Theoretical and Computational Fluid Dynamics","volume":"37 6","pages":"707 - 723"},"PeriodicalIF":3.4,"publicationDate":"2023-10-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"71909715","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Aerodynamic and aeroacoustic performance of a pitching foil with trailing edge serrations at a high Reynolds number","authors":"Xueyu Ji,&nbsp;Li Wang,&nbsp;Sridhar Ravi,&nbsp;John Young,&nbsp;Joseph C. S. Lai,&nbsp;Fang-Bao Tian","doi":"10.1007/s00162-023-00677-8","DOIUrl":"10.1007/s00162-023-00677-8","url":null,"abstract":"<p>The aerodynamic and aeroacoustic performance of a low-aspect-ratio (<span>(hbox {AR}=0.2)</span>) pitching foil during dynamic stall are investigated numerically with focus on the effects of trailing edge serrations. A hybrid method coupling an immersed boundary method for incompressible flows with the Ffowcs Williams–Hawkings acoustic analogy is employed. Large eddy simulation and turbulent boundary layer equation wall model are also employed to capture the turbulent effects. A modified NACA0012 foil with a rectangular trailing edge flap attached to the trailing edge (baseline case) undergoing pitching motion is considered. Trailing edge serrations are applied to the trailing edge flap and their effects on the aerodynamic and aeroacoustic performance of the oscillating airfoil are considered by varying the wave amplitude (<span>(2h^*= 0.05, 0.1)</span>, and 0.2) at a Reynolds number of 100,000 and a Mach number of 0.05. It is found that the reduction of the sound pressure level at the dimensionless frequency band <span>(St_{b}in [1.25,4])</span> can be over 4 dB with the presence of the trailing edge serrations (<span>(2h^*=0.1)</span>), while the aerodynamic performance and its fluctuations are not significantly altered except a reduction around 10% in the negative moment coefficient and it fluctuations. This is due to the reduction of the average spanwise coherence function and the average surface pressure with respect to that of the baseline case, suggesting the reduction of the spanwise coherence and the noise source may result in the noise reduction. Analysis of the topology of the near wake coherent structure for <span>(2h^*=0.1)</span> reveals that the alignment of the streamwise-oriented vortex with the serration edge may reduce the surface pressure fluctuation.</p>","PeriodicalId":795,"journal":{"name":"Theoretical and Computational Fluid Dynamics","volume":"38 6","pages":"825 - 844"},"PeriodicalIF":2.2,"publicationDate":"2023-10-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s00162-023-00677-8.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135093791","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Galerkin spectral estimation of vortex-dominated wake flows","authors":"Katherine J. Asztalos,&nbsp;Abdulrahman Almashjary,&nbsp;Scott T. M. Dawson","doi":"10.1007/s00162-023-00670-1","DOIUrl":"10.1007/s00162-023-00670-1","url":null,"abstract":"<p>We propose a technique for performing spectral (in time) analysis of spatially-resolved flowfield data, without needing any temporal resolution or information. This is achieved by combining projection-based reduced-order modeling with spectral proper orthogonal decomposition. In this method, space-only proper orthogonal decomposition is first performed on velocity data to identify a subspace onto which the known equations of motion are projected, following standard Galerkin projection techniques. The resulting reduced-order model is then utilized to generate time-resolved trajectories of data. Spectral proper orthogonal decomposition (SPOD) is then applied to this model-generated data to obtain a prediction of the spectral content of the system, while predicted SPOD modes can be obtained by lifting back to the original velocity field domain. This method is first demonstrated on a forced, randomly generated linear system, before being applied to study and reconstruct the spectral content of two-dimensional flow over two collinear flat plates perpendicular to an oncoming flow. At the range of Reynolds numbers considered, this configuration features an unsteady wake characterized by the formation and interaction of vortical structures in the wake. Depending on the Reynolds number, the wake can be periodic or feature broadband behavior, making it an insightful test case to assess the performance of the proposed method. In particular, we show that this method can accurately recover the spectral content of periodic, quasi-periodic, and broadband flows without utilizing any temporal information in the original data. To emphasize that temporal resolution is not required, we show that the predictive accuracy of the proposed method is robust to using temporally-subsampled data.</p>","PeriodicalId":795,"journal":{"name":"Theoretical and Computational Fluid Dynamics","volume":"38 6","pages":"801 - 823"},"PeriodicalIF":2.2,"publicationDate":"2023-08-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"44195277","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}