Description:Brief Description of the Drawings

Generally, the present disclosure relates to systems for analyzing fluid flow. Particularly, the present disclosure relates to a system for fluid flow stability analysis.
Background
The background description includes information that may be useful in understanding the present invention. It is not an admission that any of the information provided herein is prior art or relevant to the presently claimed invention, or that any publication specifically or implicitly referenced is prior art.
In the domain of computational fluid dynamics (CFD), the analysis of fluid flow stability plays a crucial role in various engineering applications. These applications range from aerospace to civil engineering, where the understanding of fluid behavior under different conditions can lead to improvements in design and operational efficiency. Analysis of fluid flow stability involves the study of how infinitesimal disturbances affect the flow field, which can be critical in predicting and mitigating undesirable phenomena such as turbulence or separation.
One common approach to analyzing fluid flow stability is through the derivation of stability equations that model the behavior of fluid under small amplitude disturbances. This process requires the collection and processing of data pertaining to the base flow, which is typically a steady, two-dimensional laminar flow. The derivation of such equations enables researchers and engineers to predict the onset of instabilities and transition to turbulence in fluid flows.
Further, the discretization of the derived stability equations is an essential step in their analysis. The Chebyshev spectral collocation method is a widely used numerical technique for this purpose, offering advantages in terms of accuracy and computational efficiency. This method involves transforming the continuous problem into a discrete one, which can be solved using computational resources. The accurate discretization of stability equations is vital for the subsequent steps in the analysis process.
Moreover, the application of boundary conditions is critical in the study of fluid flow stability. Boundary conditions at the domain inlet, outlet, wall, and free-stream play a significant role in accurately modeling the physical problem. These conditions ensure that the analysis reflects the behavior of the fluid under realistic constraints, thereby enhancing the relevance and applicability of the results.
Furthermore, the transformation of discretized stability equations into a generalized eigenvalue problem facilitates modal analysis. This transformation is crucial for identifying the modes of instability within the fluid flow. Modal analysis enables the examination of the eigenvalues and eigenvectors of the system, providing insight into the stability characteristics of the flow.
Additionally, the interpretation and presentation of results from modal analysis are fundamental in conveying the stability status of the fluid flow. An output module that effectively presents these results can aid in the understanding of complex flow phenomena and inform decisions regarding the design and operation of systems involving fluid flow.
In light of the above discussion, there exists an urgent need for solutions that overcome the problems associated with conventional systems and techniques for analyzing fluid flow stability.
Summary
The following presents a simplified summary of various aspects of this disclosure in order to provide a basic understanding of such aspects. This summary is not an extensive overview of all contemplated aspects, and is intended to neither identify key or critical elements nor delineate the scope of such aspects. Its purpose is to present some concepts of this disclosure in a simplified form as a prelude to the more detailed description that is presented later.
The following paragraphs provide additional support for the claims of the subject application.
The disclosure outlines an advanced system for fluid flow stability analysis, primarily focusing on the study of boundary layer (BL) stability over axisymmetric bodies, a topic of paramount importance in aerodynamics and hydrodynamics. These fields benefit significantly from understanding BL stability due to its direct correlation with drag reduction on streamlined vehicles, which exhibit lower drag forces compared to conventional shapes. This reduction is pivotal for creating faster, quieter vehicles with a broad spectrum of applications. Drag reduction strategies are categorized into passive methods, such as shape modification to shift the BL transition towards the stern, potentially decreasing drag by up to 30%, and active flow control (AFC) techniques like suction and injection that manipulate the flow to stabilize the BL. While passive methods face challenges due to the complexity of manufacturing curved shapes, AFC methods, particularly wall suction, can reduce drag by up to 50%, thereby lowering fuel consumption and emissions by delaying the transition from laminar to turbulent flow and extending the laminar flow region.
Given the high costs and time demands of experimental investigations, numerical simulations have become essential for predicting flow transitions and identifying critical parameters. This system employs a data collection unit for deriving stability equations that account for unsteady disturbances on steady laminar flow. A processor, utilizing a Chebyshev spectral collocation method, discretizes these equations to accurately capture variations in disturbance near the wall and minimize numerical errors. The system accurately represents boundary conditions, including no-slip conditions at the wall, homogeneous Dirichilet conditions at the inlet, and linear extrapolate conditions at the outlet. The system also models Neumann conditions for velocity disturbances in the free-stream, as well as LPPE (linear pressure Poisson equation) conditions at the wall for pressure disturbances. Additionally, it takes into account the effects of suction and injection on the BL disturbances.
Furthermore, an eigenvalue problem solver, in tandem with the processor, facilitates modal analysis for assessing BL stability by evaluating eigenvectors and eigenvalues. An output module interprets these findings, while an energy balance analysis offers deeper insight into the energy transfer mechanisms within the BL. This analysis elucidates the reasons behind the growth or decay of disturbances, thereby enhancing our understanding of BL behavior and the efficiency of control methods. Such a comprehensive approach not only advances aerodynamic design but also contributes significantly to energy efficiency improvements.

Field of the Invention

The features and advantages of the present disclosure would be more clearly understood from the following description taken in conjunction with the accompanying drawings in which:
FIG. 1 illustrates a block diagram of a system (100) for fluid flow stability analysis, in accordance with the embodiments of the present disclosure.
FIG. 2 illustrates a block diagram which indicates a stepwise procedure to perform modal analysis, in accordance with the embodiments of the present disclosure.
FIG. 3 illustrates a block diagram which indicates the stepwise procedure to perform energy balance analysis, in accordance with the embodiments of the present disclosure.
FIG. 4, spanning subfigures 4(a) to 4(b), presents a detailed schematic representation of the boundary layer dynamics around an axisymmetric cylinder subjected to the effects of mass transpiration, in accordance with the embodiments of the present disclosure.
FIG. 5 illustrates two graphs comparing the streamwise velocity component (Ub) and its radial derivative (∂Ub/∂r) in an axisymmetric boundary layer influenced by various profiles of suction and injection, all at a normal wall angle (θ = 90°), in accordance with the embodiments of the present disclosure.
FIG. 6(a) and (b) elucidate the influence of oblique angles on the Ub velocity profiles in an axisymmetric boundary layer, exhibiting how linear suction profiles and quadratic injection profiles affect the flow at a streamwise position of x = 770, in accordance with the embodiments of the present disclosure.
FIG. 7 presents a comparative eigenspectrum for an axisymmetric boundary layer (BL) subject to varying transpiration profiles: none, linear, quadratic, and uniform, in accordance with the embodiments of the present disclosure.
FIG. 8(a) and 8(b) elucidate the dynamics of stability in an axisymmetric boundary layer by depicting the eigenspectrum across a spectrum of oblique angles for linear suction and injection profiles, respectively, in accordance with the embodiments of the present disclosure.
FIG. 9 presents the eigenspectra of two distinctive boundary layer profiles—a flat plate and an axisymmetric layer—under the influence of linear injection at three different Reynolds numbers: 195, 285, and 414, in accordance with the embodiments of the present disclosure.
FIG. 10 presents the component-wise energy balance analysis of the axisymmetric boundary layer with a wall-normal uniform injection profile, in accordance with the embodiments of the present disclosure.
FIG. 11 (FIG. 11(a) and FIG. 11(b)) illustrates a comparison of energy balance integral for different (a) profiles of suction with θ = 300 and (b) oblique angles (θ) with quadratic injection profile, corresponding to Re = 285, I = 0.5% of U∞ and N = 1, in accordance with the embodiments of the present disclosure.
Detailed Description
In the following detailed description of the invention, reference is made to the accompanying drawings that form a part hereof, and in which is shown, by way of illustration, specific embodiments in which the invention may be practiced. In the drawings, like numerals describe substantially similar components throughout the several views. These embodiments are described in sufficient detail to claim those skilled in the art to practice the invention. Other embodiments may be utilized and structural, logical, and electrical changes may be made without departing from the scope of the present invention. The following detailed description is, therefore, not to be taken in a limiting sense, and the scope of the present invention is defined only by the appended claims and equivalents thereof.
The use of the terms “a” and “an” and “the” and “at least one” and similar referents in the context of describing the invention (especially in the context of the following claims) are to be construed to cover both the singular and the plural, unless otherwise indicated herein or clearly contradicted by context. The use of the term “at least one” followed by a list of one or more items (for example, “at least one of A and B”) is to be construed to mean one item selected from the listed items (A or B) or any combination of two or more of the listed items (A and B), unless otherwise indicated herein or clearly contradicted by context. The terms “comprising,” “having,” “including,” and “containing” are to be construed as open-ended terms (i.e., meaning “including, but not limited to,”) unless otherwise noted. Recitation of ranges of values herein are merely intended to serve as a shorthand method of referring individually to each separate value falling within the range, unless otherwise indicated herein, and each separate value is incorporated into the specification as if it were individually recited herein. All methods described herein can be performed in any suitable order unless otherwise indicated herein or otherwise clearly contradicted by context. The use of any and all examples, or exemplary language (e.g., “such as”) provided herein, is intended merely to better illuminate the invention and does not pose a limitation on the scope of the invention unless otherwise claimed. No language in the specification should be construed as indicating any non-claimed element as essential to the practice of the invention.
Pursuant to the "Detailed Description" section herein, whenever an element is explicitly associated with a specific numeral for the first time, such association shall be deemed consistent and applicable throughout the entirety of the "Detailed Description" section, unless otherwise expressly stated or contradicted by the context.
The term “system for fluid flow stability analysis” as used throughout the present disclosure relates to an assembly of components configured to analyze the stability of fluid flow, especially focusing on the impact of infinitesimal amplitude disturbances on a steady, two-dimensional laminar base flow. The system aims to provide insights into the transition of fluid flow from laminar to turbulent states, which is critical in various engineering and scientific applications.
The term “data collection unit” as used throughout the present disclosure relates to a component of the system tasked with deriving stability equations. These equations are based on the application of infinitesimal amplitude disturbances to a steady, two-dimensional laminar base flow. The purpose of the data collection unit is to capture and process initial flow conditions and disturbances to formulate mathematical models that represent the stability characteristics of the fluid flow.
The term “processor” as used throughout the present disclosure relates to a computational component within the system. This processor is adapted to discretize the derived stability equations using a Chebyshev spectral collocation method. By discretizing the stability equations, the processor converts continuous mathematical models into a form suitable for numerical analysis, thereby facilitating the detailed examination of fluid flow stability.
The term “boundary condition applicator” as used throughout the present disclosure relates to a component designed to impose predetermined boundary conditions at various boundaries of the computational domain. These boundaries include domain inlet, outlet, free-stream, wall, and porous where velocity and pressure disturbances are considered. The boundary condition applicator ensures that the numerical simulations accurately reflect the physical constraints and conditions present in real-world fluid flow scenarios.
The term “eigenvalue problem solver” as used throughout the present disclosure relates to a component in communication with the processor. It is configured to transform the discretized stability equations into a generalized eigenvalue problem suitable for modal analysis. The eigenvalue problem solver plays a crucial role in identifying and analyzing the modes of instability within the fluid flow, thereby providing valuable information on the flow’s stability status.
The term “output module” as used throughout the present disclosure relates to a component structured to interpret and present the results from the eigenvalue problem solver. The output module effectively communicates the stability status of the fluid flow, based on the analysis conducted by the eigenvalue problem solver. It presents the results in a format that can be readily understood and utilized by researchers, engineers, and scientists for further analysis or decision-making.
FIG. 1 illustrates a block diagram of a system (100) for fluid flow stability analysis, in accordance with the embodiments of the present disclosure. The system (100) is composed of various interconnected components, each serving a distinct function in the analysis process. A data collection unit (102) is provided for deriving stability equations based on infinitesimal amplitude disturbances upon a steady, two-dimensional laminar base flow. Said data collection unit (102) captures initial conditions and disturbances to establish mathematical models for stability assessment. In communication with the data collection unit (102), a processor (104) is adapted for discretizing the derived stability equations using a Chebyshev spectral collocation method. Said processor (104) facilitates the transformation of continuous stability models into a discrete form, enabling detailed numerical analysis. A boundary condition applicator (106) is included for the imposition of predetermined boundary conditions at domain inlet, outlet, free-stream, wall, and porous boundaries. Such a boundary condition applicator (106) ensures that simulation conditions reflect realistic operational parameters for the velocity and pressure disturbances within the computational domain. An eigenvalue problem solver (108) is in communication with the processor (104), tasked with transforming the discretized stability equations into a generalized eigenvalue problem suitable for modal analysis. Said eigenvalue problem solver (108) plays a crucial role in identifying the modes of instability within the fluid flow. Additionally, an output module (110) is structured to interpret and display results from the eigenvalue problem solver (108). Such an output module (110) presents the findings in a manner indicative of the stability status of the analyzed fluid flow, thereby assisting in the decision-making process for system improvements or further research.
In an embodiment, the system (100) further comprises a data processor to apply numerical solutions to the derived base Navier-Stokes (N-S) equations relevant to the steady two-dimensional laminar base flow. The inclusion of a data processor enhances the system's (100) capability by enabling the application of numerical methods to solve the Navier-Stokes equations, which are fundamental to understanding fluid dynamics and flow stability. The N-S equations describe the motion of fluid substances and are pivotal in the analysis of fluid flow, particularly in predicting the behavior under various conditions. By solving these equations, the system (100) can accurately model the fluid flow, taking into account factors such as viscosity and pressure, thereby providing a comprehensive analysis of the flow's stability. This addition not only extends the system's (100) analytical depth but also augments its utility in complex fluid dynamics studies.
In another embodiment, the data processor within the system (100) is further configured to implement Chebyshev polynomials for the discretization of governing stability equations and to determine the number of collocation points required for said discretization. This configuration enables the system (100) to utilize Chebyshev polynomials, known for their efficiency in dealing with differential equations, thereby enhancing the precision of the discretization process. The determination of collocation points, which are used in the Chebyshev spectral collocation method, is crucial for achieving accurate results in the analysis of fluid flow stability. The method's ability to handle boundary conditions effectively makes it an ideal choice for this application. The adoption of Chebyshev polynomials and the strategic determination of collocation points significantly improve the system's (100) capability to model and analyze fluid dynamics with higher accuracy and less computational effort.
In a further embodiment, the eigenvalue problem solver (108) within the system (100) employs a Krylov subspace method with specific shift values and a predefined number of iterations to solve the generalized eigenvalue problem. The Krylov subspace method, renowned for its efficiency in solving large-scale eigenvalue problems, is particularly suited for the analysis of fluid flow stability. By employing this method with specific shift values and a predefined number of iterations, the eigenvalue problem solver (108) enhances the system's (100) ability to accurately and efficiently determine the stability of the fluid flow. This approach not only speeds up the computation process but also ensures that the analysis captures the essential characteristics of the flow's behavior. The use of the Krylov subspace method signifies a significant advancement in the system's (100) capability to tackle complex fluid dynamics problems, offering a robust tool for researchers and engineers.
In another embodiment, the output module (110) within the system (100) is additionally configured to conduct a modal analysis to extract the least stable eigenmode, as well as its temporal and spatial growth characteristics. This configuration allows the system (100) to not only identify instabilities within the fluid flow but also to understand their evolution over time and space. Modal analysis is a powerful tool in the study of fluid dynamics, providing insights into the mechanisms driving flow instability. By focusing on the least stable eigenmode, the system (100) pinpoints the conditions under which the flow is most susceptible to transition from laminar to turbulent, enabling targeted interventions. The ability to extract detailed information about the temporal and spatial growth characteristics of instabilities greatly enhances the system's (100) utility in designing more efficient and stable fluid flow systems.
In a further embodiment, the boundary condition applicator (106) within the system (100) is adapted to impose boundary conditions consisting of homogeneous Dirichlet conditions at the inlet and porous section, linear extrapolation conditions at the outlet, no-slip conditions on the wall, and Neumann conditions along the free-stream boundaries. This adaptation ensures that the system (100) can accurately simulate real-world fluid flow scenarios by applying relevant boundary conditions. Homogeneous Dirichlet conditions, which specify the value of the function on the boundary, and Neumann conditions, which specify the value of the derivative of the function on the boundary, are critical for accurately modeling the flow near boundaries. The ability to apply these conditions enhances the system's (100) fidelity in replicating actual fluid flow behaviors, thereby improving the reliability of the stability analysis.
In another embodiment, the data processor within the system (100) executes an algorithm for the integration of kinetic energy equations over the computational domain to determine total mechanical energy (TME) and energy production due to Reynolds stress (EPRS). This algorithm allows for a comprehensive analysis of the energy dynamics within the fluid flow, providing insights into the distribution and transformation of energy. By calculating TME and EPRS, the system (100) can assess the efficiency of the flow and identify potential improvements. The integration of kinetic energy equations offers a quantitative measure of the flow's stability, making it an invaluable tool for optimizing fluid flow systems for energy efficiency and performance.
In a further embodiment, the algorithm within the system (100) is configured to evaluate contributions to the kinetic energy from viscous dissipation (VD), pressure work (PW), and streamline curvature (SC) effects. This configuration enables a detailed analysis of the factors influencing the kinetic energy within the fluid flow. Viscous dissipation, pressure work, and streamline curvature are critical elements that determine the flow's energy dynamics. By assessing their contributions, the system (100) provides a nuanced understanding of how energy is generated, transformed, and lost in fluid flow. This insight is crucial for designing flows that are stable, efficient, and optimized for specific applications.
In another embodiment, the eigenvalue problem solver (108) within the system (100) is configured to process the stability equations for average disturbances over a single time and azimuthal mode to discern energy balance components and their influence on fluid flow stability. This configuration allows for a focused analysis of how disturbances affect the flow's energy balance, offering insights into the mechanisms behind flow instability. By concentrating on average disturbances over a single mode, the system (100) can identify specific factors contributing to instability, enabling targeted interventions to enhance flow stability. This approach not only simplifies the analysis but also provides valuable information for the development of strategies to maintain or improve fluid flow stability.
In a further embodiment, the data processor within the system (100) is calibrated to adjust boundary conditions dynamically based on the flow characteristics at the inlet, outlet, wall, and free-stream to maintain the fidelity of the stability analysis. This dynamic adjustment allows the system (100) to adapt to changes in flow conditions in real-time, ensuring that the stability analysis remains accurate under varying operational scenarios. By calibrating the data processor to modify boundary conditions in response to observed flow characteristics, the system (100) enhances its ability to simulate real-world fluid dynamics accurately. This adaptability is crucial for conducting reliable fluid flow stability analyses, as it accounts for the inherent variability in fluid flow systems. Such dynamic adjustment of boundary conditions represents a significant advancement in the system’s (100) capabilities, enabling it to provide more precise and relevant results. This, in turn, aids in the development of more efficient and effective fluid flow systems by ensuring that stability analyses reflect actual conditions as closely as possible. The inclusion of this feature underscores the system's (100) comprehensive approach to fluid flow stability analysis, highlighting its utility in a wide range of applications where accurate modeling of fluid behavior is critical.
FIG. 2 illustrates a block diagram which indicates a stepwise procedure to perform modal analysis, in accordance with the embodiments of the present disclosure. FIG. 2 depicts stepwise procedure to compute the leading bi-global temporal eigenmode and its spatial structure for modal analysis. This procedure is adopted for modal analysis of any BL problem. The first step is to derive governing stability equations from the base and instantaneous equations. Then discretized using spectral method and imposed appropriate stability boundary conditions. The stability equations also contain the base flow solution, which is obtained by solving 2D N-S steady equation using a finite volume code. The discretized stability equations are written in the form of EVP and solved by Arnoldi iterative algorithm and Krylov subspace method. The solution of the EVP contains eigenvalues and vectors, and extracting the only leading eigenmode and corresponding eigenvectors indicate the flow stability and spatial structure of disturbances, respectively.
FIG. 3 illustrates a block diagram which indicates the stepwise procedure to perform energy balance analysis, in accordance with the embodiments of the present disclosure. The first step is to derive the kinetic energy equation for the small disturbance from the momentum equation through algebraic manipulation. This equation is then integrated over the entire computational domain, and averaging over a single time and azimuthal mode yields the different energy balance components. Based on whether the total modal energy (TME) is positive or negative, the flow is determined to be modally unstable or stable, respectively.
FIG. 4, spanning subfigures 4(a) to 4(b), presents a detailed schematic representation of the boundary layer dynamics around an axisymmetric cylinder subjected to the effects of mass transpiration, in accordance with the embodiments of the present disclosure. The diagram delineates the boundary layer as it responds to both uniform and non-uniform transpiration profiles at various oblique angles, captured within the computational domain labeled 'abcd.' Within this domain, 'L' signifies the cylinder's total length, while 'Lx' and 'Lr' define the dimensions of the computational space. The investigation encapsulated in the diagram is a novel exploration into the global stability analysis of the axisymmetric boundary layer, an area previously unexplored in contrast to the local stability analyses predominantly featured in existing literature. This analysis places particular emphasis on the boundary layer's reaction to the non-uniform, oblique suction, and injection of mass across its surface—a factor proven to be pivotal in influencing stability outcomes and effective in boundary layer control. The core objectives that emerge from a synthesis of literature reviews are twofold: firstly, to conduct a comprehensive modal stability analysis of the axisymmetric boundary layer, taking into account the implications of oblique and non-uniform transpiration, and secondly, to perform an energy balance analysis aimed at elucidating the fundamental physical processes that facilitate energy exchange, contributing to boundary layer destabilization. Uniform and non-uniform transpiration profiles—including linear and quadratic variations—as well as oblique transpiration at specified angles ranging from 30 to 150 degrees, are integral to this study. By achieving these objectives, the study aims to advance boundary layer management strategies and reduce the drag forces that impinge upon cylindrical bodies within marine and aerospace environments, optimizing their performance and operational efficiency..
In FIG. 5, two graphs are presented, comparing the streamwise velocity component (Ub) and its radial derivative (∂Ub/∂r) in an axisymmetric boundary layer influenced by various profiles of suction and injection, all at a normal wall angle (θ = 90°), in accordance with the embodiments of the present disclosure. The experiments, set at a Reynolds number of 195 with an intensity of 0.5% of the free-stream velocity (U∞) at a streamwise position of x = 360, reveal how different treatments—ranging from no suction or injection to uniform, linear, and quadratic approaches—affect the velocity within the boundary layer. The first graph (a) displays Ub and shows each transpiration strategy's impact on the velocity profile, crucial for understanding how the wall boundary conditions alter the flow. The second graph (b) focuses on the radial gradient of Ub (∂Ub/∂r), indicating how the velocity changes across the boundary layer's thickness and is indicative of the momentum exchange rates. Notably, a grey rectangular box highlights the 'point of inflection' in these profiles, a critical feature that may signal the transition between laminar and turbulent flow. Such data is instrumental in devising control strategies for boundary layers, potentially leading to significant advancements in fluid dynamics applications, particularly in optimizing energy efficiency by managing flow separation and drag forces.
FIG. 6(a) and (b) elucidate the influence of oblique angles on the Ub velocity profiles in an axisymmetric boundary layer, exhibiting how linear suction profiles and quadratic injection profiles affect the flow at a streamwise position of x = 770, in accordance with the embodiments of the present disclosure. In these figures, the transpiration velocity, resultant from suction and injection processes, is decomposed into two distinct components—tangential (Vw cosθ) and normal (Vw sinθ) to the wall. With the increment of the angle θ from 30 degrees toward 90 degrees, there is an escalation in the wall-normal component, reaching its zenith at 90 degrees, concurrently, the tangential component diminishes, ceasing entirely at θ = 90 degrees. This phenomenon causes a noticeable alteration in the velocity profiles: during suction, the profile gravitates closer to the cylinder surface, while during injection, it shifts away. Contrarily, when the angle θ surpasses 90 degrees, advancing towards 150 degrees, an inverse trend is observed; the wall-normal component dwindles, and the tangential component augments in the reverse direction, mitigating its influence on the base flow due to its relative insignificance—a maximum of 2.5% of the free-stream velocity U∞ in this study. This leads to a discernible similarity in the base flow profiles at angles θ = 30 degrees and 150 degrees, as well as θ = 60 degrees and 120 degrees, with disparities so minute they generally escape unaided visual detection, rendering the base flow profiles virtually indistinguishable due to equivalent wall-normal components of transpiration velocity.
FIG. 7 presents a comparative eigenspectrum for an axisymmetric boundary layer (BL) subject to varying transpiration profiles: none, linear, quadratic, and uniform, in accordance with the embodiments of the present disclosure. In scenarios devoid of transpiration, the boundary layer displays modal stability since the temporal growth rate (ωi) of the primary Tollmien-Schlichting (T-S) wave is negative, indicating no amplification of disturbances under the stipulated parameters. This stability landscape shifts under the influence of suction, as observed in FIG. 7(a), where suction encourages a transition of the T-S branch deeper into the damped region of the eigenspectrum, signifying enhanced stability. This effect arises due to suction attenuating the energy exchange between the base flow and perturbations, consequently starving disturbances of the sustenance required for growth. Conversely, injection manifests an inverse dynamic; the influx of fluid into the BL reduces wall shear stress and expands the layer, magnifying the non-parallel effects and broadening the perturbation profile. As FIG. 7(b) illustrates, this escalation in disturbance energy, coupled with a decrease in viscous energy dissipation, propels the T-S branch toward the upper half of the eigenspectrum, indicative of decreased modal stability. With uniform injection, the BL's modal stability is compromised even at lower Reynolds numbers and transpiration velocities, as some T-S modes traverse into the unstable region. Conversely, linear and quadratic injection profiles maintain stability at subdued velocities and Reynolds numbers but surrender to instability as these parameters increase. Thus, while suction is synonymous with the suppression of instability modes like T-S waves through the extraction of high-energy fluid layers, injection exacerbates these instabilities by bolstering the fluid's energy near the wall. The interplay of these factors outlines a complex framework wherein the BL's stability is highly contingent on the mode of transpiration, the nature of the profiles employed, and the magnitude of the transpiration velocity relative to the free-stream speed. FIG. 8(a) and 8(b) elucidate the dynamics of stability in an axisymmetric boundary layer by depicting the eigenspectrum across a spectrum of oblique angles for linear suction and injection profiles, respectively, in accordance with the embodiments of the present disclosure. As the angle of suction ascends from 30°, progressing to the perpendicular angle of 90°, the Tollmien-Schlichting (T-S) branch is observed to gravitate towards a more damped state within the eigenspectrum, suggesting an increase in the boundary layer's stability. This trajectory implies that as the suction's obliquity heightens, the effectiveness in damping disturbances enhances, correlating to a sturdier flow. In stark contrast, when examining injection profiles, an elevation in the oblique angle marks a migration of the T-S branch towards the unstable region of the spectrum, signaling a decrease in flow stability; this movement denotes that the injected fluid's momentum at greater angles introduces energy into the system, facilitating the growth of instabilities.
As the angular scope extends beyond the perpendicular, from 90° to the complete reversal at 180°, for both suction and injection, an inversion in trends is recorded, with the eigenspectrum revealing a symmetry in stability outcomes. This pattern reflects an invariance in base flow behavior across the angular intervals, with only slight discrepancies emerging due to the relatively low transpiration intensities—capped at 2.5% of the free-stream velocity U∞ in the study at hand. This consistency is also mirrored in other suction and injection profiles, underpinning a general stability characteristic that oblique angles, whether acute or obtuse, influence the boundary layer's stability in a predictable and replicable manner.
FIG. 9 presents the eigenspectra of two distinctive boundary layer profiles—a flat plate and an axisymmetric layer—under the influence of linear injection at three different Reynolds numbers: 195, 285, and 414, in accordance with the embodiments of the present disclosure. The conditions maintained for the comparison include a transpiration angle of 30° and an intensity level at 0.5% of the free-stream velocity U∞, with a specified mode of N = 0 for the axisymmetric boundary layer. In each graph, the Tollmien-Schlichting (T-S) modes are identified, providing a clear visualization of how the stability characteristics of each boundary layer configuration vary with the Reynolds number. As Reynolds number increases, the T-S modes, which represent the primary instability mechanisms within the boundary layers, exhibit changes in their temporal growth rates (ωi). The depicted trends offer insights into the stability of the boundary layers: the axisymmetric layer is shown to have a different stability response when compared to the flat plate, as evidenced by the positioning and spread of the T-S modes in the eigenspectra. These distinctions hold implications for understanding the fluid dynamics and potential transition to turbulence in different flow geometries, particularly in the context of engineering applications where flow control and drag reduction are paramount.
FIG. 10 presents the component-wise energy balance analysis of the axisymmetric boundary layer with a wall-normal uniform injection profile, in accordance with the embodiments of the present disclosure. This analysis is performed corresponding to the leading T-S mode. The results indicate that one of the energy production terms,
𝑃
2
=
∬

(

𝑢
𝑝

𝑣
𝑝

𝜕
𝑈
𝑏
/𝜕𝑟)𝑑𝑟𝑑𝑥, where subscript ‘p’ refers to disturbance quantity), is dominant over the other three terms (P1, P3, and P4), and the viscous dissipation (VD) term is also dominant. Therefore, these two parameters significantly affect the stability results, while the rest of the parameters have negligible values. The total mechanical energy (TME ≈ P2 – VD ≈2ωi) being positive indicates that the flow is modally unstable.
FIG. 11 (FIG. 11(a) and FIG. 11(b)) illustrates a comparison of energy balance integral for different (a) profiles of suction with θ = 300 and (b) oblique angles (θ) with quadratic injection profile, corresponding to Re = 285, I = 0.5% of U∞ and N = 1, in accordance with the embodiments of the present disclosure.FIG. 11(a) examines how different suction profiles influence the energy dynamics within an axisymmetric boundary layer by analyzing the interplay between the energy production term P2, the viscous dissipation (VD), and the total mechanical energy (TME). The findings underscore that TME registers negative values under suction, a hallmark of modal stability characterized by diminished energy production and augmented viscous dissipation. Among these profiles, the uniform suction stands out, yielding the lowest TME and thereby indicating a superior stability compared to its non-uniform counterparts. This suggests that uniform suction is more effect

I/We Claims

A system (100) for fluid flow stability analysis comprising:
a data collection unit (102) configured to derive stability equations based on imposed infinitesimal amplitude disturbances upon a steady, two-dimensional laminar base flow;
a processor (104) adapted to discretize the derived stability equations utilizing a Chebyshev spectral collocation method;
a boundary condition applicator (106) operable to impose predetermined boundary conditions at domain inlet, outlet, free-stream, wall, and porous boundaries of the computational domain for both velocity and pressure disturbances;
an eigenvalue problem solver (108) in communication with said processor (104), configured to transform the discretized stability equations into a generalized eigenvalue problem for modal analysis; and
an output module (110) structured to interpret and present results from the eigenvalue problem solver (108) indicating the stability status of the fluid flow.
The system (100) of claim 1, further comprising a data processor to apply numerical solutions to the derived base Navier-Stokes (N-S) equations relevant to the steady two-dimensional laminar base flow.
The system (100) of claim 2, wherein the data processor is further configured to implement Chebyshev polynomials for the discretization of governing stability equations, and to determine the number of collocation points required for said discretization.
The system (100) of claim 3, wherein said eigenvalue problem solver (108) employs a Krylov subspace method with specific shift values and a predefined number of iterations to solve the generalized eigenvalue problem.
The system (100) of claim 4, wherein the output module (110) is additionally configured to conduct a modal analysis to extract the least stable eigenmode, as well as its temporal and spatial growth characteristics.
The system (100) of claim 5, wherein the boundary condition applicator (106) is adapted to apply boundary conditions that comprises homogeneous Dirichlet conditions at the inlet, linear extrapolate conditions at the outlet, no-slip conditions at the wall, and Neumann conditions at the free-stream boundaries for velocity disturbances, while LPPE conditions at the wall boundary for pressure disturbances.
The system (100) of claim 6, wherein the data processor further executes an algorithm for the integration of kinetic energy equations over the computational domain to determine total mechanical energy (TME) and energy production due to Reynolds stress (EPRS).
The system (100) of claim 7, wherein the algorithm is further configured to evaluate contributions to the kinetic energy from viscous dissipation (VD), pressure work (PW), and streamline curvature (SC) effects.
The system (100) of claim 8, wherein said eigenvalue problem solver (108) is configured to process the stability equations for average disturbances over a single time and azimuthal mode to discern energy balance components and their influence on fluid flow stability.
The system (100) of claim 2, wherein the data processor is further calibrated to adjust boundary conditions dynamically based on the flow characteristics at the inlet, outlet, wall, and free-stream to maintain the fidelity of the stability analysis.

SYSTEM FOR ANALYZING FLUID FLOW STABILITY

The present disclosure provides a system for fluid flow stability analysis. The system comprises a data collection unit configured to derive stability equations based on imposed infinitesimal amplitude disturbances upon a steady, two-dimensional laminar base flow; a processor adapted to discretize the derived stability equations utilizing a Chebyshev spectral collocation method; a boundary condition applicator operable to impose predetermined boundary conditions at domain inlet, outlet, free-stream, wall, and porous boundaries within the computational domain for both velocity and pressure disturbances; an eigenvalue problem solver in communication with the processor, configured to transform the discretized stability equations into a generalized eigenvalue problem for modal analysis; and an output module structured to interpret and present results from the eigenvalue problem solver indicating the stability status of the fluid flow.

, Claims:I/We Claims

A system (100) for fluid flow stability analysis comprising:
a data collection unit (102) configured to derive stability equations based on imposed infinitesimal amplitude disturbances upon a steady, two-dimensional laminar base flow;
a processor (104) adapted to discretize the derived stability equations utilizing a Chebyshev spectral collocation method;
a boundary condition applicator (106) operable to impose predetermined boundary conditions at domain inlet, outlet, free-stream, wall, and porous boundaries of the computational domain for both velocity and pressure disturbances;
an eigenvalue problem solver (108) in communication with said processor (104), configured to transform the discretized stability equations into a generalized eigenvalue problem for modal analysis; and
an output module (110) structured to interpret and present results from the eigenvalue problem solver (108) indicating the stability status of the fluid flow.
The system (100) of claim 1, further comprising a data processor to apply numerical solutions to the derived base Navier-Stokes (N-S) equations relevant to the steady two-dimensional laminar base flow.
The system (100) of claim 2, wherein the data processor is further configured to implement Chebyshev polynomials for the discretization of governing stability equations, and to determine the number of collocation points required for said discretization.
The system (100) of claim 3, wherein said eigenvalue problem solver (108) employs a Krylov subspace method with specific shift values and a predefined number of iterations to solve the generalized eigenvalue problem.
The system (100) of claim 4, wherein the output module (110) is additionally configured to conduct a modal analysis to extract the least stable eigenmode, as well as its temporal and spatial growth characteristics.
The system (100) of claim 5, wherein the boundary condition applicator (106) is adapted to apply boundary conditions that comprises homogeneous Dirichlet conditions at the inlet, linear extrapolate conditions at the outlet, no-slip conditions at the wall, and Neumann conditions at the free-stream boundaries for velocity disturbances, while LPPE conditions at the wall boundary for pressure disturbances.
The system (100) of claim 6, wherein the data processor further executes an algorithm for the integration of kinetic energy equations over the computational domain to determine total mechanical energy (TME) and energy production due to Reynolds stress (EPRS).
The system (100) of claim 7, wherein the algorithm is further configured to evaluate contributions to the kinetic energy from viscous dissipation (VD), pressure work (PW), and streamline curvature (SC) effects.
The system (100) of claim 8, wherein said eigenvalue problem solver (108) is configured to process the stability equations for average disturbances over a single time and azimuthal mode to discern energy balance components and their influence on fluid flow stability.
The system (100) of claim 2, wherein the data processor is further calibrated to adjust boundary conditions dynamically based on the flow characteristics at the inlet, outlet, wall, and free-stream to maintain the fidelity of the stability analysis.

SYSTEM FOR ANALYZING FLUID FLOW STABILITY