Spelling suggestions: "subject:"[een] STABILITY ANALYSIS"" "subject:"[enn] STABILITY ANALYSIS""
221 |
Partialkoefficienter för stabilitetsanalys av betongdammarWängberg, Alexander January 2015 (has links)
I Sverige finns det idag ett stort antal dammar och nära hälften av Sveriges elproduktion kommer ifrån vattenkraft. Höga säkerhetskrav ställs på dammarna då konsekvenserna av ett dammbrott kan orsaka stor ekonomisk skada både i form av minskad produktion och som uppbyggnadskostnader. Dammbrott kan även utgöra risk för skador på människor och omkringliggande miljö samt byggnader. Det finns flera dimensioneringsmetoder för att uppfylla gällande säkerhetskrav på konstruktioner. De tillämpningsvägledningar som används vid stabilitetsanalys av betongdammar baseras på gamla deterministiska metoder med säkerhetsfaktorer. Det arbetas med att ta fram nya tillämpningsvägledningar baserat på sannolikhetsbaserade metoder. De sannolikhetsbaserade metoderna tar hänsyn till osäkerheterna i enskilda variabler, vilket förväntas ge effektivare konstruktioner. Stabilitetsanalys med sannolikhetsbaserade metoder i varje enskilt fall är förmodligen det bästa ur säkerhetssynpunkt, men det kan vara väldigt tidskrävande. Ett sett att underlätta stabilitetsberäkningarna, i vanligt förekommande konstruktioner, är användandet av partialkoefficienter. Partialkoefficienter är en semi-probabilistisk metod som kan kalibreras med hjälp av tillförlitlighetsanalysen och appliceras på systemets olika variabler och parametrar. På så vis beaktas osäkerheten i enskilda variabler och parametrar bättre än i deterministiska metoder med säkerhetsfaktorer. Målet med det här examensarbetet var att undersöka om stabilitetsanalys av betongdammar med partialkoefficienter kan vara ett alternativ till de metoder som redan finns. Rapporten kan delas in i tre delar. Den första delen av rapporten beskriver teorin bakom kalibrering av partialkoefficienter med hjälp av tillförlitlighetsteori, FORM. I den andra delen beskrivs metoder och tillämpningsvägledning för stabilitetsanalys av betongdammar med probabilistiska metoder. I den tredje delen används kunskaperna från de två första delarna för att beräkna tillförlitligheten på 15 utav Sveriges dammar. Beräkningarna används sedan för att kalibrera partialkoefficienter. En del kraftiga avgränsningar har gjorts i arbetet, bland annat har beräkningarna utförts för ett statiskt lastfall och bara fokuserat på stabilitetsvillkoret för glidning. Resultatet visar att det inte är rekommenderat att använda sig av partialkoefficienter vid stabilitetsanalys, åtminstone inte med den information och kunskap som idag finns tillgänglig. Osäkerheterna kring vissa modeller och variabler behöver minskas för att partialkoefficienter skall vara ett alternativ. Teorin kring partialkoefficienter kräver även en viss likhet i konstruktioner och konstruktionselement. De dammar som användes i denna studie kan ha varit för olika för att erhålla tillfredsställande resultat med avseende på likartade partialkoefficienter. / In Sweden today there are a large number of dams and nearly half of Sweden's electricity is produced from hydropower. The safety requirements on the dams are high due to the consequences that a dam failure can cause. There are several design methods to achieve the expected safety requirements. The design guidelines used in the stability analysis of concrete dams in Sweden is based on the deterministic methods with safety factors. However, a new proposition for design guidelines based on probabilistic methods is being developed. Compared to the deterministic approach the probabilistic method takes into account the uncertainties in individual variables, which are expected to provide more efficient structures. One problem with stability analysis using probabilistic methods is that it can be very time consuming. Another method which combines the simplicity of the deterministic approach with the effectiveness of the probabilistic method is the use of partial factors. The use of partial factors is a semi-probabilistic method that can be calibrated from the reliability analysis in the probabilistic method and applied to the individual variables and parameters in the system similar to the safety factor. The aim of this thesis was to investigate if stability analysis of concrete dams with partial factors can be an alternative to the methods already available. The report can be divided into three parts. The first part of the report describes the theory behind the calibration process of partial factors using the first order reliability method, FORM. The second part describes existing Swedish methods and application guidelines for stability analysis of concrete dams. The knowledge from the first two parts is then used in order to calculate the reliability of 15 of Sweden's dams. It should be observed that the work contains some limitations, for instance only the sliding stability is studied using one static load case. The results from the calculation of this load case are then used to calibrate the partial factors. The result shows that it is not recommended to use the partial factors for stability analysis of concrete dams, at least not with the information and knowledge available today in the field. The uncertainties surrounding certain models and variables need to be reduced in order for partial factors to be an option. The theory behind partial factors requires a certain degree of uniformness when it comes to the structure or elements used. The difference among the dams analyzed in this study may have impacted the results negatively.
|
222 |
Damping of sub-synchronous control interactions with a STATCOM : Wind farms & series compensated power linesAlvarez Urrutia, Leonardo January 2022 (has links)
The power converter is one of the key components in power system applications such as high voltage direct current (HVDC) systems and the grid connection of intermittent sources such as wind and solar power. However, the increased penetration of converter-based generation introduces challenges, such as sub-synchronous interaction between the converter control system and the grid. These control interactions are characterized by fast-growing, subsynchronous oscillations (SSO). This thesis deals with the analysis of sub-synchronous control interactions (SSCI) between doubly-fed induction generator (DFIG)-based wind farms and series compensated transmission lines. Moreover, the thesis aims to identify a method for mitigating the sub-synchronous oscillations using a static synchronous compensator (STATCOM), with a supplementary damping controller. The study is based on work in PSCAD/EMTDC and uses a system based on the IEEE first benchmark model, acting as a grid, and the scaled power output of a DFIG turbine model, modeling a wind farm. Initial impedance-based analysis in the frequency domain shows that the DFIG wind farm has a negative resistance throughout the sub-synchronous frequency range. A negative resistance may result in negative damping of the system and further introduce the risk of instability. The wind farm resistance and, in turn, system stability is affected by the current control loop of the DFIG-converter. The transmission line compensation factor largely impacts the system stability, while the power output has a minor effect. A time-domain analysis is performed to verify the result of the frequency domain analysis. Further on, a grid-forming STATCOM is added to the system for VAr compensation. Additional stability analysis shows that even though improvingthe stability, the STATCOM alone is not adequate to mitigate the SSCI. The proposed damping strategy is based on modifying the STATCOM voltage reference andcan be divided into three steps: detecting the SSO, estimating the sub-synchronous component, and modifying the extracted signal. The detection algorithm is based on a half-cycle comparator, while the modification is done through a proportional gain. When estimating the sub-synchronous components, two methods are proposed and compared. The first estimation method is based on a conventional power system stabilizer (PSS) method, and the second is afilter-less method.
|
223 |
A computational framework for multidimensional parameter space screening of reaction-diffusion models in biologySolomatina, Anastasia 16 March 2022 (has links)
Reaction-diffusion models have been widely successful in explaining a large variety of patterning phenomena in biology ranging from embryonic development to cancer growth and angiogenesis. Firstly proposed by Alan Turing in 1952 and applied to a simple two-component system, reaction-diffusion models describe spontaneous spatial pattern formation, driven purely by interactions of the system components and their diffusion in space. Today, access to unprecedented amounts of quantitative biological data allows us to build and test biochemically accurate reaction-diffusion models of intracellular processes. However, any increase in model complexity increases the number of unknown parameters and thus the computational cost of model analysis. To efficiently characterize the behavior and robustness of models with many unknown parameters is, therefore, a key challenge in systems biology. Here, we propose a novel computational framework for efficient high-dimensional parameter space characterization of reaction-diffusion models. The method leverages the $L_p$-Adaptation algorithm, an adaptive-proposal statistical method for approximate high-dimensional design centering and robustness estimation. Our approach is based on an oracle function, which describes for each point in parameter space whether the corresponding model fulfills given specifications. We propose specific oracles to estimate four parameter-space characteristics: bistability, instability, capability of spontaneous pattern formation, and capability of pattern maintenance. We benchmark the method and demonstrate that it allows exploring the ability of a model to undergo pattern-forming instabilities and to quantify model robustness for model selection in polynomial time with dimensionality. We present an application of the framework to reconstituted membrane domains bearing the small GTPase Rab5 and propose molecular mechanisms that potentially drive pattern formation.
|
224 |
DYNAMICS OF LARGE ARRAY MICRO/NANO RESONATORSBorra, Chaitanya 15 July 2020 (has links)
No description available.
|
225 |
Reliability – based analysis of embankment damsLiu, Xiaochen January 2019 (has links)
Embankment dams are widely used as water-retaining structures through the world. Two principal structural failure modes of embankment dams are internal erosion and slope instability. The objective of this thesis is to investigate the feasibility of applying probabilistic methods in embankment dam assessment. Compared to deterministic methods, probabilistic methods considers the variability in material properties and the results are more assessable. Meanwhile, the derived sensitivity factors of the variables can be used to find the most influential one.Through a literature review, it is found that the difficulties in defining a limit state function for the full process of internal erosion for embankment dams hold back the application of reliability-based methods. However, reliability-based methods have been recognized as suitable for slope stability assessment.In this thesis, both deterministic and probabilistic calculations have been performed for a homogenous embankment dam. The deterministic calculation is carried out in Slope/W to first find out the most critical slip surface. The failure probability of this slip surface is estimated by both Monte Carlo simulations in Slope/W and FORM analysis in Comrel. Then, the 34 most critical slip surfaces from Slope/W are chosen to calculate the system reliability with simple bounds and applied integration. The results from the deterministic and probabilistic calculations showed that the studied dam has a higher risk of failure than allowed. In the system reliability calculation, simple bounds gives a wide range, especially the upper bounds. It was found that classification of different slip surfaces into different groups with respect to geometry and material was useful to simplify the calculations. Not much improvement in accuracy of the system probability of failure is found when increasing the numbers (from 1 to 5) of the most critical slip surfaces from each group. Even the combination of 1 slip surface from each group only gives an error of 6.5%. Categorizing slip surfaces before performing a system reliability analysis is a good way for simplification of the calculations. Due to this, the criteria used for categorization of the slip surfaces are of interest. In this work, it was found that the slip geometry in combination with the materials included in the slip surface constitute a possible way forward to do this. / Fyllningsdammar används ofta som dämmande konstruktioner runt om i världen. Två potentiella brottmoder för fyllningsdammar är inre erosion och släntinstabilitet.Syftet med detta examensarbete är att undersöka möjligheten att tillämpa probabilistiska metoder vid utvärdering av fyllningsdammars säkerhet för dessa brottmoder. Jämfört med deterministiska metoder betraktar probabilistiska metoder även variationen i materialegenskaper. Erhållna känslighetsfaktorer för ingående variabler kan även användas för att hitta de som har störst inverkan på det studerade problemet.I litteraturstudien framkom att det är svårt att definiera en gränsfunktion som beskriver inre erosion för fyllnadsdammar, vilket begränsar dess användbarhet med probabilistiska metoder. För analys av släntstabilitet är dock probabilistiska metoder användbara. Både deterministiska och probabilistiska beräkningar av släntstabiliteten för en fyllnadsdamm har genomförts inom ramen för detta arbete. De deterministiska beräkningarna har utförts med programmet Slope/W, där den mest kritiska glidytan först har analyserats. Brottsannolikheten för dammen har beräknats med både Monte Carlo-simulering i Slope/W och FORM-analys i Comrel. En begränsning med att enbart beräkna brottsannolikheten för den mest kritiska glidytan är att brottsannolikheten kan underskattas då inverkan från andra potentiella glidytor inte beaktas. De potentiella glidytorna kan betraktas som ett seriesystem, där glidytorna är delvis korrelerade med varandra. För att studera denna fråga valdes de 34 mest kritiska glidytorna från Slope/W ut för att beräkna systemets tillförlitlighet. Resultaten från de deterministiska och probabilistiska beräkningarna visade att den studerade dammen hade en oacceptabel brottsannolikhet. I beräkningarna av gränserna för systemets brottsannolikhet visade resultaten att brottsannolikheten kan variera stort beroende på korrelationen mellan glidytorna. Genom att klassificera glidytorna i olika grupper med hänsyn till geometri och ingående material kan beräkningarna för systemets brottsannolikhet förenklas. Resultaten visade att noggrannheten i beräkningarna för brottsannolikheten av systemet inte förbättrades i någon större utsträckning när antalet valda glidytor (från 1 till 5) från varje grupp av likartade glidytor inkluderades i beräkningen. Kombinationen av 1 glidyta från varje grupp gav endast ett fel på 6,5% jämfört med om alla 34 glidytorna inkluderades. Kategorisering av glidytor av samma karaktär i olika grupper rekommenderas därför att genomföras innan systemets tillförlitlighet analyseras. I detta arbete framkom att glidytan geometri i kombination med ingående material i glidytan utgör ett lämpligt kriterium för att göra denna indelning.
|
226 |
Reliability – based analysis of embankment damsLiu, Xiaochen January 2019 (has links)
Embankment dams are widely used as water-retaining structures through the world. Two principal structural failure modes of embankment dams are internal erosion and slope instability. The objective of this thesis is to investigate the feasibility of applying probabilistic methods in embankment dam assessment. Compared to deterministic methods, probabilistic methods considers the variability in material properties and the results are more assessable. Meanwhile, the derived sensitivity factors of the variables can be used to find the most influential one.Through a literature review, it is found that the difficulties in defining a limit state function for the full process of internal erosion for embankment dams hold back the application of reliability-based methods. However, reliability-based methods have been recognized as suitable for slope stability assessment.In this thesis, both deterministic and probabilistic calculations have been performed for a homogenous embankment dam. The deterministic calculation is carried out in Slope/W to first find out the most critical slip surface. The failure probability of this slip surface is estimated by both Monte Carlo simulations in Slope/W and FORM analysis in Comrel. Then, the 34 most critical slip surfaces from Slope/W are chosen to calculate the system reliability with simple bounds and applied integration. The results from the deterministic and probabilistic calculations showed that the studied dam has a higher risk of failure than allowed. In the system reliability calculation, simple bounds gives a wide range, especially the upper bounds. It was found that classification of different slip surfaces into different groups with respect to geometry and material was useful to simplify the calculations. Not much improvement in accuracy of the system probability of failure is found when increasing the numbers (from 1 to 5) of the most critical slip surfaces from each group. Even the combination of 1 slip surface from each group only gives an error of 6.5%. Categorizing slip surfaces before performing a system reliability analysis is a good way for simplification of the calculations. Due to this, the criteria used for categorization of the slip surfaces are of interest. In this work, it was found that the slip geometry in combination with the materials included in the slip surface constitute a possible way forward to do this.
|
227 |
Grid stability evaluation and techno-economic analysis of shared solar applications within a Swedish energy community.Holsti Heijbel, Hugo January 2023 (has links)
The joint European electricity production in 2022 was 2.9 PWh where 34.3% of the electricity came from sources of natural gas and coal. However, to be able to achieve the goals of the Paris agreement, greenhouse gas emissions must peak before 2025 and thus call for a rapid decarbonatization. For this, using less energy originating from fossil sources and increasing the use of renewable energy is crucial. One way to decentralize and implement renewable energy is using energy communities. This is also a way for actors to participate in energy trading and to enable collective ownership of local production and storage technology. However, since the electricity grid in Sweden is largely outdated, this poses as a problem for continued electrification. To evaluate if decentralisation of the electricity distribution and production can minimise the load on the national grid, The research institute of Sweden (RISE) has initiated pilot studies to assess the possibility and impact of introducing energy communities in Sweden with the help of KTH for simulations. As part of such a pilot study this thesis focuses on a district in Örebro named Tamarinden, where it is planned to implement a local DC-network to allow sharing of locally PV produced electricity. Through using the Python package Panda Power as well as quantitative analyses in Excel, the thesis aim was to propose a configuration of the standard local distribution grid and analyse the impact from PVs, battery storage as well as the DC-network focusing on grid stability and solar penetration. In addition, an economic feasibility analysis was made. The results showed that an energy community greatly impact the solar penetration and that grid stability improvements could not be determined. This could be due to the choice of production technology, solar PV, as these have peak production during low demand season. Furthermore, a proposition for the local grid configuration was constructed. Nevertheless, the economic feasibility showed that even though the impact on PV penetration is great, the cost of the technology is too high in comparison to savings made in terms of lowered electricity demand. In conclusion, the implementation of such an energy community in Sweden is currently unbeneficial from an economic point of view. / Den europeiska elproduktionen var år 2022 2.9 PWh var av 34.3% härstammade från kol och naturgas. För att uppnå målen i Parisavtalet måste växthusgasutsläppen sluta öka senast 2025 därför krävs snabb avkarbonisering. För att åstadkomma detta måste användandet av fossila energikällor minska och användandet av förnyelsebar energi öka. Ett sätt att möjliggöra denna omställning är genom energigemenskaper. Detta är även ett sätt för aktörer att delta i energimarknaden samt att möjliggöra kollektivt ägarskap av lokal produktions och lagrings teknologier. Då elnätet i Sverige till stora delar är utdaterat skapar detta problematik, för en fortsatt elektrifiering. För att utvärdera huruvida decentralisering av eldistributionen och elproduktionen kan minska belastningen på det nationella distributionsnätet, driver Research institute of Sweden (RISE) tillsammans med KTH pilotstudier för att undersöka möjligheterna och påverkan av att introducera energigemenskaper i Sverige. Som en del av en av pilotstudierna fokuserar denna uppsats på en stadsdel i Örebro som heter Tamarinden, där det är planerat att implementera ett lokalt DC-nätverk för att möjliggöra delning av lokalt producerad el. Genom modellering i Pythons paket Panda Power och kvantitativa analyser i Excel, var uppsatsens mål att presentera en möjlig konfiguration av lokalnätet samt att analysera påverkan från solceller, batteri lagring samt DC-nätverket på nätstabiliteten och solpenetrationen. Ytterligare gjordes en ekonomisk genomförbarhetsanalys. Resultaten visade att en energigemenskap kraftigt påverkar solpenetrationen och att förbättringar av nätstabiliteten inte kunde fastställas. Detta kan bero på valet av produktionsteknik, solceller, eftersom dessa har topproduktion under lågsäsong. Dessutom konstruerades ett förslag för den lokala nätkonfigurationen. Den ekonomiska genomförbarheten visade dock att även om effekten på solcellspenetrationen är stor är kostnaden för tekniken för hög i jämförelse med besparingar i form av sänkt efterfrågan på el. Sammanfattningsvis är genomförandet av en sådan energigemenskap i Sverige för närvarande ofördelaktigt ur ekonomisk synvinkel.
|
228 |
Nonlinear Analysis and Control of Standalone, Parallel DC-DC, and Parallel Multi-Phase PWM ConvertersMazumder, Sudip K. 17 August 2001 (has links)
Applications of distributed-power systems are on the rise. They are already used in telecommunication power supplies, aircraft and shipboard power-distribution systems, motor drives, plasma applications, and they are being considered for numerous other applications. The successful operation of these multi-converter systems relies heavily on a stable design. Conventional analyses of power converters are based on averaged models, which ignore the fast-scale instability and analyze the stability on a reduced-order manifold. As such, validity of the averaged models varies with the switching frequency even for the same topological structure.
The prevalent procedure for analyzing the stability of switching converters is based on linearized smooth averaged (small-signal) models. Yet there are systems (in active use) that yield a non-smooth averaged model. Even for systems for which smooth averaged models are realizable, small-signal analyses of the nominal solution/orbit do not provide anything about three important characteristics: region of attraction of the nominal solution, dependence of the converter dynamics on the initial conditions of the states, and the post-instability dynamics. As such, converters designed based on small-signal analyses may be conservative. In addition, linear controllers based on such analysis may not be robust and optimal. Clearly, there is a need to analyze the stability of power converters from a different perspective and design nonlinear controllers for such hybrid systems.
In this Dissertation, using bifurcation analysis and Lyapunov's method, we analyze the stability and dynamics of some of the building blocks of distributed-power systems, namely standalone, integrated, and parallel converters. Using analytical and experimental results, we show some of the differences between the conventional and new approaches for stability analyses of switching converters and demonstrate the shortcomings of some of the existing results. Furthermore, using nonlinear analyses we attempt to answer three fundamental questions: when does an instability occur, what is the mechanism of the instability, and what happens after the instability?
Subsequently, we develop nonlinear controllers to stabilize parallel dc-dc and parallel multi-phase converters. The proposed controllers for parallel dc-dc converters combine the concepts of multiple-sliding-surface and integral-variable-structure control. They are easy to design, robust, and have good transient and steady-state performances. Furthermore, they achieve a constant switching frequency within the boundary layer and hence can be operated in interleaving or synchronicity modes. The controllers developed for parallel multi-phase converters retain many of the above features. In addition, they do not require any communication between the modules; as such, they have high redundancy. One of these control schemes combines space-vector modulation and variable-structure control. It achieves constant switching frequency within the boundary layer and a good compromise between the transient and steady-state performances. / Ph. D.
|
229 |
Subsystem Design in Aircraft Power Distribution Systems using OptimizationChandrasekaran, Sriram 26 June 2000 (has links)
The research reported in this dissertation focuses on the development of optimization tools for the design of subsystems in a modern aircraft power distribution system. The baseline power distribution system is built around a 270V DC bus. One of the distinguishing features of this power distribution system is the presence of regenerative power from the electrically driven flight control actuators and structurally integrated smart actuators back to the DC bus. The key electrical components of the power distribution system are bidirectional switching power converters, which convert, control and condition electrical power between the sources and the loads. The dissertation is divided into three parts.
Part I deals with the formulation of an optimization problem for a sample system consisting of a regulated DC-DC buck converter preceded by an input filter. The individual subsystems are optimized first followed by the integrated optimization of the sample system. It is shown that the integrated optimization provides better results than that obtained by integrating the individually optimized systems.
Part II presents a detailed study of piezoelectric actuators. This study includes modeling, optimization of the drive amplifier and the development of a current control law for piezoelectric actuators coupled to a simple mechanical structure.
Linear and nonlinear methods to study subsystem interaction and stability are studied in Part III. A multivariable impedance ratio criterion applicable to three phase systems is proposed. Bifurcation methods are used to obtain global stability characteristics of interconnected systems. The application of a nonlinear design methodology, widely used in power systems, to incrementally improve the robustness of a system to Hopf bifurcation instability is discussed. / Ph. D.
|
230 |
[pt] ANÁLISE DE ESTABILIDADE DE ESCOAMENTOS VISCOSOS E VISCOELÁSTICOS / [en] LINEAR STABILITY ANALYSIS OF VISCOUS AND VISCOELASTIC FLOWSJULIANA VIANNA VALERIO 04 June 2007 (has links)
[pt] As informações sobre a sensibilidade da solução de um dado
escoamento
mediante a perturbações infinitesimais é importante para o
seu completo
entendimento. A análise de estabilidade de escoamentos
pode ser utilizada
na otimização de processos industriais. Na indústria de
revestimento o controle
da estabilidade é fundamental, uma vez que o escoamento na
região
de aplicação da camada de líquido sobre o substrato, de um
modo geral,
tem que ser laminar, bidimensional e em regime permanente.
O objetive é
determinar, dentro do espaço de parâmetros de operação, a
região onde o escoamento
é estável e conseqüêntemente a camada a ser revestida
uniforme.
Porém, por ser uma análise complexa, só é usada na
indústria em estudos
mais apurados. O sistema linear que descreve a
estabilidade vai ser discretizado
com o método de Galerkin / elementos finitos, dando origem
a um
problema de autovalor generalizado.Tanto para escoamentos
com líquidos
newtonianos como para escoamentos com líquidos
viscoelásticos, uma das
matrizes do problema de autovalor generalizado é singular
e alguns autovalores
se encontram no infinito. No escoamento com líquidos
viscoelásticos
parte do espectro é contínuo, aumentando o grau de
dificuldade da análise
numérica para encontrá-lo. Nesse trabalho, vamos
apresentar um método
baseado em transformações lineares tirando vantagem das
estruturas matriciais
e transformando-as em um problema de autovalor clássico
com dimens
são, pelo menos, três vezes menor que o original. O método
elimina os
autovalores infinitos do problema com um baixo custo
computacional. A
estabilidade de um escoamento de Couette unidimensional de
líquido newtoniano
é analisada como um primeiro exemplo. Depois, o início do
estudo
da estabilidade em um escoamento de Couette bidimensional
e também um
escoamento pistonado com o mesmo líquido. Generaliza-se o
método para
o escoamento de Couette de um líquido viscoelástico, os
resultados para
o escoamento de um líquido cujo comportamento mecânico é
descrito pelo
modelo de Maxwell são apresentados e comparados com a
solução analítica
de Gorodtsov & Leonov, 1967. A relação entre os
autovetores do problema
transformado e do original é apresentada. / [en] Steady state,two-dimensional flows may become unstable
under two and
three-dimensional disturbances, if the flow parameters
exceed some critical
values. In many practical situations, determining the
parameters at
which the flow becomes unstable is essential. The complete
understanding
of viscous and viscoelastic flows requires not only the
steady state solution
of the governing equations, but also its sensitivity to
small perturbations.
Linear stability analysis leads to a generalized
eigenvalue problem, GEVP.
Solving the GEVP is challenging, even for Newtonian
liquids, because the
incompressibility constraint creates singularities that
lead to nonphysical
eigenvalues at infinity. For viscoelastic flows, the
difficulties are even higher
because of the continuous spectrum of eigenmodes
associated with differential
constitutive equations. The complexity and high
computational cost
of solving the GEVP have probably discouraged the use of
linear stability
analysis of incompressible flows as a general engineering
tool for design and
optimization. The Couette flow of UCM liquids has been
used as a classical
problem to address some of the important issues related to
stability analysis
of viscoelastic flows. The spectrum consists of two
discrete eigenvalues
and a continuous segment of eigenvalues with real part
equal to -1/We (We
is the Weissenberg number). Most of the numerical
approximation of the
spectrum of viscoelastic Couette flow presented in the
literature were obtained
using spectral expansions. The eigenvalues close to the
continuous
part of the spectrum show very slow convergence. In this
work, the linear
stability of Couette flow of a Newtonian and UCM liquids
were studied
using finite element method, which makes it easier to
extend the analysis
to complex flows. A new procedure to eliminate the
eigenvalues at infinity
from the GEVP that come from differential equations is
also proposed.
The procedure takes advantage of the structure of the
matrices involved
and avoids the computational effort of common mapping
techniques. With
the proposed procedure, the GEVP is transformed into a
smaller simple
EVP, making the computations more effcient. Reducing the
computational
memory and time. The relation between the eigenvector from
the original
problem and the reduced one is also presented.
|
Page generated in 0.0411 seconds