Theoretical and Computational Fluid Dynamics最新文献

{"title":"Flow dynamics in a model of a dilated thoracic aorta prior to and following prosthetic replacement","authors":"Pier Giuseppe Ledda,&nbsp;Maria Grazia Badas,&nbsp;Gildo Matta,&nbsp;Giorgio Querzoli","doi":"10.1007/s00162-023-00651-4","DOIUrl":"10.1007/s00162-023-00651-4","url":null,"abstract":"<p>We numerically investigate the flow dynamics in a model of a dilated thoracic aorta, and compare the flow features with the case of a prosthetic replacement in its ascending part. The flow is characterized by an inlet jet which impacts the aortic walls and sweeps toward the aortic arch. Secondary flows generated by the transvalvular jet evolve downstream into a helical flow. The small curvature radius at the end of the aortic arch induces flow separation and vortex shedding in the initial part of the descending aorta, during the systole. The implantation of a prosthesis determines several modifications in the global and local flow patterns. An increase of the pulse wave velocity in the aorta leads to larger pressures inside the vessel, due to the geometrical and rigidity modifications. The sweeping jet is more aligned along the axial direction and propagates faster along the aortic arch. Consequently, a stronger separation of the flow downstream of the aortic arch is observed. By also exploiting manifold analysis, we identify regions characterized by near-wall disordered flows which may present intense accumulation and drop of concentration of biochemicals. These regions are localized downstream of the prosthetic replacement, in the aortic arch, and may be more prone to a new emergence of vessel dilation.</p>","PeriodicalId":795,"journal":{"name":"Theoretical and Computational Fluid Dynamics","volume":"37 3","pages":"375 - 396"},"PeriodicalIF":3.4,"publicationDate":"2023-05-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s00162-023-00651-4.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"4398732","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":"Effects of annulation on low Reynolds number flows over an orthocone","authors":"Mitesh Thakor,&nbsp;Kee Horng Seh,&nbsp;Sareta R. Gladson,&nbsp;Martin L. Fernandez,&nbsp;Linda C. Ivany,&nbsp;Melissa Green,&nbsp;Yiyang Sun","doi":"10.1007/s00162-023-00649-y","DOIUrl":"10.1007/s00162-023-00649-y","url":null,"abstract":"<p>This study numerically examines the influences of transverse annulation around a cone surface on the characteristics of a flow over an orthocone. This work is inspired by <i>Spyroceras</i>, a fossilized genus of nautiloid cephalopods from the Paleozoic era, whose method of locomotion is understudied. As a baseline case, a flow over a smooth orthoconic model with a blunt cone end is investigated numerically at Reynolds numbers from 500 to 1500. As Reynolds increases, two different shedding mechanisms—hairpin-vortex wake and spiral-vortex wake—are captured. We notice that an introduction of annulation over the cone surface changes the critical Reynolds number for the transition of the shedding mechanism. The dominant shedding frequency increases with the Reynolds number for the smooth and annulated cone flows. Moreover, the annulation reduces the dominant frequency for the same Reynolds number and increases the time-averaged drag coefficient. Modal decompositions—Proper Orthogonal Decomposition (POD) and Spectral Proper Orthogonal Decomposition (SPOD)—are used to capture the coherent structures and their oscillating frequencies. We have captured modes corresponding to the hairpin-vortex wake and spiral-vortex wake shedding mechanisms. Comparing the leading POD modes for the smooth and the annulated cone flows, we find that the annulation can reduce the twisting effects of the coherent structures in the wake. Additionally, we find that the SPOD analysis can identify modes presenting both hairpin-vortex wake and spiral-vortex wake in one flow condition as leading modes, while the POD leading modes only reveal one shedding mechanism in each flow.</p>","PeriodicalId":795,"journal":{"name":"Theoretical and Computational Fluid Dynamics","volume":"37 3","pages":"357 - 374"},"PeriodicalIF":3.4,"publicationDate":"2023-05-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s00162-023-00649-y.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"4267906","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":"Simulation of droplet impact dynamics on V-shaped walls","authors":"Guoqiang Wu,&nbsp;Sheng Chen","doi":"10.1007/s00162-023-00652-3","DOIUrl":"10.1007/s00162-023-00652-3","url":null,"abstract":"<p>This paper presents the morphological evolution characteristics of a droplet impacting a V-shaped wall by using the lattice Boltzmann method (LBM). Four parameters are investigated comprehensively. The parameters vary over wide ranges: surface wettability (<span>(60^circ le theta ^{eq} le 120^circ )</span>), Weber number (<span>(102.27 le text {We} le 3681.82)</span>), bending angle of the V-shaped wall (90<span>(^circ le theta le 180^circ )</span>), and eccentricity ratio (0 <span>(le b le )</span> 0.5). Two types of collision are observed: deposition and breakage. For breakage, the number of satellite droplets increases against the increment of We. The splashing occurs for a high We. And the lamella ejection is observed on the hydrophilic wall and the neutral wall. The lamella ejection will be slight against the increase of <span>(theta ^{eq})</span>, while it will become obvious against the increment of <span>(theta )</span>. In addition, the nondimensional spreading length, width, and height are measured and analyzed. Regime maps are established based on We, Re, and <span>(theta )</span>.</p>","PeriodicalId":795,"journal":{"name":"Theoretical and Computational Fluid Dynamics","volume":"37 2","pages":"173 - 202"},"PeriodicalIF":3.4,"publicationDate":"2023-05-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"4224402","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":"Vortex dynamics of axisymmetric cones at high angles of attack","authors":"Al Shahriar,&nbsp;Rajan Kumar,&nbsp;Kourosh Shoele","doi":"10.1007/s00162-023-00647-0","DOIUrl":"10.1007/s00162-023-00647-0","url":null,"abstract":"<p>Vortex asymmetry, dynamics, and breakdown in the wake of an axisymmetric cone have been investigated using direct numerical simulation for a wide range of angles of attack. The immersed boundary method is employed with pseudo-body-conformal grids to ensure the accuracy and resolution requirements near the body while being able to account for topology changes near the cone tip. The separated shear layer originated from the surface of the cone swirls into a strong primary vortex. Beneath the primary vortex on the leeward surface of the cone, a well-coherent counter-rotating secondary vorticity is generated. Beyond a particular threshold of swirl, the attached vortex structure breaks and the flow undergoes a chaotic transformation. Depending on the angle of attack, the flow shows different levels of instabilities and the topology of the vortices changes in the wake. In addition to swirl, spiral vortices that revolve around the primary vortex core often merge with the core and play a role in developing the double-helix mode of instability at the onset of the vortex breakdown. At the angle of attack of 60<span>(^circ )</span>, the time-averaged side force becomes asymmetric at the stage where the drag overcomes the lift. At the angle of attack of 75<span>(^circ )</span>, the primary vortex governs the flow asymmetry and the side force. Flow asymmetry is independent of the vortex breakdown. Finally, the contribution of primary vortices to the total forces is quantified using a force partitioning method.</p>","PeriodicalId":795,"journal":{"name":"Theoretical and Computational Fluid Dynamics","volume":"37 3","pages":"337 - 356"},"PeriodicalIF":3.4,"publicationDate":"2023-05-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"4179124","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":"Flow over a hydrofoil subjected to traveling wave-based surface undulation: effect of phase difference between surface waves and wave number","authors":"Sarvesh Shukla,&nbsp;Atul Sharma,&nbsp;Amit Agrawal,&nbsp;Rajneesh Bhardwaj","doi":"10.1007/s00162-023-00646-1","DOIUrl":"10.1007/s00162-023-00646-1","url":null,"abstract":"<p>Flow around a traveling wave-based surface-undulating NACA0012 hydrofoil has been numerically studied. In particular, we determine the effect of the phase speed of the wave, the phase difference between the waves traveling on the top and bottom surfaces, and the wave number on flow dynamics around and behind the hydrofoil and propulsive performance. The flow results in a vortex sheet or a street behind the hydrofoil, where oppositely signed vortices are aligned in either forward or reverse direction. Apart from these, side vortices start forming on either side of the hydrofoil at a higher wave number. The phase difference analysis between the upper and lower surface undulation reveals the configuration better for the hydrofoil’s lateral and longitudinal stability. The hydrofoil can shift from high thrust to high lateral force configuration by changing the phase difference between waves on the top and bottom surfaces. Thrust increases with an increase in the wave number, and a threshold value of phase speed and wave number exists where the drag-to-thrust transition happens. The added mass force-based scaling analysis corroborates with the simulated results.</p>","PeriodicalId":795,"journal":{"name":"Theoretical and Computational Fluid Dynamics","volume":"37 3","pages":"319 - 336"},"PeriodicalIF":3.4,"publicationDate":"2023-04-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s00162-023-00646-1.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"5106989","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":"Passive control of dynamic stall using a flow-driven micro-cavity actuator","authors":"Miguel R. Visbal,&nbsp;Daniel J. Garmann","doi":"10.1007/s00162-023-00645-2","DOIUrl":"10.1007/s00162-023-00645-2","url":null,"abstract":"<p>A novel passive flow control strategy for the mitigation of transient separation and dynamic stall is demonstrated by means of high-fidelity large-eddy simulations. The control technique is based on a properly-sized micro-cavity cut into a wing’s underside near the leading edge, ahead of stagnation. This cavity remains essentially inactive at low incidence. However, as the wing effective angle of attack increases, the stagnation point displaces past the micro-cavity and the accelerating flow grazing the cavity induces a high-frequency resonance phenomenon or so-called Rossiter modes. The self-generated small-scale disturbances are carried around the leading-edge through the boundary layer to the wing’s upper side where the laminar separation bubble (LSB) amplifies these disturbances. This process delays LSB bursting and dynamic stall when the cavity size is selected such that its naturally occurring Rossiter modes are tuned to the receptivity of the LSB. Control effectiveness is explored for a harmonically pitching NACA 0012 wing section with freestream Mach number <span>(M_infty = 0.2)</span>, chord Reynolds numbers <span>(textrm{Re}_textrm{c} = 5 times 10^5)</span>, and maximum angle of attack of <span>(18^circ )</span>. The flow fields are computed employing a validated overset high-order implicit large-eddy simulation (LES) solver based on sixth-order compact schemes for the spatial derivatives augmented with an eighth-order low-pass filter. Despite its simplicity, the micro-cavity resonance is found to be highly effective in preventing the deep dynamic stall experienced by the baseline airfoil. A significant reduction in the cycle-averaged drag and in the force and moment fluctuations is achieved. In addition, the negative (unstable) net-cycle pitch damping found in the baseline cases is eliminated.</p>","PeriodicalId":795,"journal":{"name":"Theoretical and Computational Fluid Dynamics","volume":"37 3","pages":"289 - 303"},"PeriodicalIF":3.4,"publicationDate":"2023-04-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"4998499","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":"Non-isothermal thin-film flow of a viscoplastic material over topography: critical Bingham number for a partial slump","authors":"Miguel Moyers-Gonzalez,&nbsp;James N. Hewett,&nbsp;Dale R. Cusack,&nbsp;Ben M. Kennedy,&nbsp;Mathieu Sellier","doi":"10.1007/s00162-023-00642-5","DOIUrl":"10.1007/s00162-023-00642-5","url":null,"abstract":"<p>This paper considers the non-isothermal flow of a viscoplastic fluid on a horizontal or an inclined surface with a flat, a step-up and a step-down topography. A particular application of interest is the spread of a fixed mass—a block—of material under its own weight. The rheology of the fluid is described by the Bingham model which includes the effect of yield stress, i.e. a threshold stress which must be exceeded before flow can occur. Both the plastic viscosity and the yield stress are modelled with temperature-dependent parameters. The flow is described by a reduced model with a thin-film equation for the height of the block and a depth-averaged energy conservation equation for the heat transfer. Results show that for large values of the yield stress, only the outer fraction of the fluid spreads outward, the inner fraction remaining unyielded, hence the block only partially slumps. Conversely, for small values of the yield stress, the entire block of fluid becomes yielded and therefore slumps. We present an analysis which predicts the critical value of the yield stress for which partial slump occurs and how it depends on temperature.</p>","PeriodicalId":795,"journal":{"name":"Theoretical and Computational Fluid Dynamics","volume":"37 2","pages":"151 - 172"},"PeriodicalIF":3.4,"publicationDate":"2023-03-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s00162-023-00642-5.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"4692654","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":"Direct numerical simulations and flow-pressure acoustic analyses of flapping-uvula-induced flow evolutions within normal and constricted pharynx","authors":"Jinxiang Xi,&nbsp;Junshi Wang,&nbsp;Xiuhua April Si,&nbsp;Haibo Dong","doi":"10.1007/s00162-023-00638-1","DOIUrl":"10.1007/s00162-023-00638-1","url":null,"abstract":"<p>Snoring and obstructive sleep apnea (OSA) are often associated with uvula vibrations and pharynx constrictions. However, successful treatment of snoring or accurate diagnosis of OSA has been proven challenging. This study aimed to identify acoustic indexes that were sensitive to underlying airway structural or kinematic variations. Six physiologically realistic models were developed that consisted of three pharynx constriction levels (M1-3) and two uvula-flapping kinematics (K1-2). Direct numerical simulations (DNS) were performed to resolve spatial and temporal flow dynamics, and an immersed boundary method was used to approximate the uvula vibrations. Time-varying acoustic pressures at six points in the pharynx were analyzed using different algorithms in frequency- or frequency–time domains. Signature flow structures formed near the uvula for different uvula motions and in the pharynx for different pharyngeal constriction levels. The fast Fourier transform showed that the acoustic energy was mainly distributed in four peaks (flapping frequency and three harmonics) with descending magnitudes. Their amplitudes and distribution patterns differed among the six models but were not substantial. The continuous wavelet transforms showed clearly separated acoustic cycles (in both frequency and time) in the uvula-induced flows and revealed a cascading bifurcation pattern in the input–output semblance map. Specifically, the multifractal spectrum was sensitive to uvula flapping kinematics but not pharynx constrictions. By contrast, the input–output cross-correlation and Hilbert phase space showed high sensitivity to pharynx constrictions but low sensitivity to uvula kinematics. The frequency–time analyses of DNS-predicted pressures offered insight into the acoustics signals that were not apparent in original signals and could be used individually or in combination in diagnosis or treatment planning for snoring/OSA patients.</p>","PeriodicalId":795,"journal":{"name":"Theoretical and Computational Fluid Dynamics","volume":"37 2","pages":"131 - 149"},"PeriodicalIF":3.4,"publicationDate":"2023-03-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s00162-023-00638-1.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"4875982","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":"Vortex-induced vibration and flutter of a filament behind a circular cylinder","authors":"Mohd Furquan,&nbsp;Sanjay Mittal","doi":"10.1007/s00162-023-00644-3","DOIUrl":"10.1007/s00162-023-00644-3","url":null,"abstract":"<p>Flow past a flexible filament, a two-dimensional splitter plate with negligible thickness, attached behind a circular cylinder is investigated. The Reynolds number based on the free-stream speed of incoming flow and diameter of the cylinder is <span>(textrm{Re}=100)</span>. At this <span>(textrm{Re})</span>, the flow for a rigid filament is steady. However, a flexible filament undergoes flow-induced vibration for a range of reduced speed, <span>(U^*)</span>, defined as inverse of the first nondimensionalized natural frequency of the filament. Over the wide range of <span>(U^*)</span> considered in this work (<span>(U^*le 240)</span>), it exhibits both flutter and vortex-induced vibration (VIV). Lock-in with various normal modes related to bending of the filament, each in a different regime of reduced speed, is observed during VIV. Interestingly, the fluid–structure system does not lock-in with the first normal mode of bending but with higher modes. The flow is steady for an extended range of reduced speed both before and after the lock-in with second mode. Two patterns of vortex shedding are observed. The <span>(textsf{2P})</span> mode is associated with high-frequency vibration, while the <span>(mathsf {2,S})</span> mode is observed during relatively low-frequency oscillation. A symmetry-breaking pitchfork bifurcation leads to static deflection of the filament during the first steady regime. The filament exhibits flutter response, at large reduced speed, with relatively low amplitude and frequency. No vortex shedding is observed during flutter. The fluid forces that cause flutter arise from asymmetry across the two sides of the filament in the zones of recirculation downstream of the cylinder. Comparison of the space-time patterns of energy transfer at the fluid–filament interface for flutter and vortex-induced vibration reveals that the energy transfer is much smaller during flutter compared to VIV. The point of maximum energy transfer is located close to the root of the filament in case of flutter, while it is near the tip during VIV.\u0000</p>","PeriodicalId":795,"journal":{"name":"Theoretical and Computational Fluid Dynamics","volume":"37 3","pages":"305 - 318"},"PeriodicalIF":3.4,"publicationDate":"2023-03-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"4838864","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":"Backpropagation of neural network dynamical models applied to flow control","authors":"Tarcísio Déda,&nbsp;William R. Wolf,&nbsp;Scott T. M. Dawson","doi":"10.1007/s00162-023-00641-6","DOIUrl":"10.1007/s00162-023-00641-6","url":null,"abstract":"<p>Backpropagation of neural network models (NNMs) is applied to control nonlinear dynamical systems using several different approaches. By leveraging open-loop data, we show the feasibility of building surrogate models with control inputs that are able to learn important features such as types of equilibria, limit cycles and chaos. Two novel approaches are presented and compared to gradient-based model predictive control (MPC): the neural network control (NNC), where an additional neural network is trained as a control law in a recurrent fashion using the nonlinear NNMs, and linear control design, enabled through linearization of the obtained NNMs. The latter is compared with dynamic mode decomposition with control (DMDc), which also relies on a data-driven linearized model. It is shown that the linearized NNMs better approximate the systems’ behavior near an equilibrium point than DMDc, particularly in cases where the data display highly nonlinear characteristics. The proposed control approaches are first tested on low-dimensional nonlinear systems presenting dynamical features such as stable and unstable limit cycles, besides chaos. Then, the NNC is applied to the nonlinear Kuramoto–Sivashinsky equation, exemplifying the control of a chaotic system with higher dimensionality. Finally, the proposed methodologies are tested on the compressible Navier–Stokes equations. In this case, the stabilization of a cylinder vortex shedding is sought using different actuation setups by taking measurements of the lift force with delay coordinates.</p>","PeriodicalId":795,"journal":{"name":"Theoretical and Computational Fluid Dynamics","volume":"37 1","pages":"35 - 59"},"PeriodicalIF":3.4,"publicationDate":"2023-02-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s00162-023-00641-6.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"5331479","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}