Spelling suggestions: "subject:"finitedifference"" "subject:"clinicaldifferences""
201 |
Defektkorrekturverfahren für singulär gestörte Randwertaufgaben / Defect Correction Methods for Singularly Perturbed Boundary Value ProblemsFröhner, Anja 27 December 2002 (has links) (PDF)
Wir untersuchen ein Defektkorrekturverfahren, das ein einfaches Upwind-Differenzenverfahren erster Ordnung mit einem zentralen Differenzenverfahren kombiniert, für ein- und zweidimensionale singulär gestörte Konvektions-Diffusions-Probleme auf einer Klasse von Shishkin-Typ-Gittern. Im eindimensionalen Fall wird nachgewiesen, dass das Verfahren von (fast) zweiter Ordnung, gleichmäßig bezüglich des Diffusionsparameters $\epsilon$ konvergiert. Zur Konvergenzanalyse für das zweidimensionale Modellproblem werden verschiedene Techniken diskutiert. In einem Spezialfall kann auf einem stückweise uniformen Shishkin-Gitter die $\epsilon$-gleichmäßige Konvergenz des Verfahrens von fast zweiter Ordnung gezeigt werden. Ferner sind die bisher bekannten Stabilitätsaussagen und ihre Verwendung zur Konvergenzanalysis der betrachteten Differenzenverfahren sowie Methoden zur Analyse von Defektkorrekturverfahren zusammengestellt. Einige Bemerkungen zu Defektkorrekturverfahren und Finite-Elemente-Methoden schließen die Arbeit ab. Numerische Experimente untermauern die theoretischen Resultate. / We consider a defect correction method that combines a first-order upwinded difference scheme with a second-order central difference scheme for model singularly perturbed convection-diffusion problems in one and two dimensions on a class of Shishkin-Type meshes. In one dimension, the method is shown to be convergent uniformly in the diffusion parameter $\epsilon$ of second order in the discrete maximum norm. To analyze the two-dimensional case, we discuss several proof techniques for defect correction methods. For a special problem with constant coefficients on a piecewise uniform Shishkin-mesh we can show the second order convergence of the considered scheme, uniformly with respect to the diffusion parameter. Moreover the known stability properties and their impact on the convergence analysis of the considered differnce schemes are compiled. Some remarks on defect correction and finite elements conclude the theses. Numerical experiments support our theoretical results.
|
202 |
Signal and power integrity co-simulation using the multi-layer finite difference methodBharath, Krishna 26 March 2009 (has links)
Mixed signal system-on-package (SoP) technology is a key enabler for increasing functional integration, especially in mobile and wireless
systems. Due to the presence of multiple dissimilar modules, each having unique power supply requirements, the design of the power distribution network (PDN) becomes critical.
Typically, this PDN is designed as alternating layers of power and ground planes with signal interconnects routed in between or on top of the planes.
The goal for the simulation of multi-layer power/ground planes, is the following:
Given a stack-up and other geometrical information, it is required to find the
network parameters (S/Y/Z) between port locations.
Commercial packages have extremely complicated stack-ups, and the trend to increasing
integration at the package level only points to increasing complexity. It is computationally
intractable to solve these problems using these existing methods.
The approach proposed in this thesis for obtaining the response of the PDN is the multi-layer finite difference method (M-FDM).
A surface mesh / finite difference based approach is developed, which leads to a system matrix that is
sparse and banded, and can be solved efficiently.
The contributions of this research are the following:
1. The development of a PDN modeler for multi-layer packages and boards called the the multi-layer finite difference method.
2. The enhancement of M-FDM using multi-port connection networks to include the effect of fringe fields and gap coupling.
3. An adaptive triangular mesh based scheme called the multi-layer finite element method (MFEM) to address the limitations of M-FDM
4. The use of modal decomposition for the co-simulation of signal nets with the PDN.
5. The use of a robust GA-based optimizer for the selection and placement of decoupling capacitors in multi-layer geometries.
6. Implementation of these methods in a tool called MSDT 1.
|
203 |
ANÁLISE DINÂMICA DE TORRES ESTAIADAS DE LINHAS DE TRANSMISSÃO SUBMETIDAS À RUPTURA DE CABO / DINAMIC ANALYSIS OF TRANSMISSION LINES GUYED TOWERS SUBMITTED TO BROKEN CABLECarlos, Thiago Brazeiro 10 February 2015 (has links)
Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / The need for transmission lines expansion (TL) is evident in the face of growing electrical energy demand in Brazil as well as the great distances between generation and consumption centers. It is known that numerous failures can interrupt the energy transmission and the towers collapse is among the main reasons, causing huge losses to both electrical energy companies and final users. Thus, it is crucial to understand the mechanical behavior of these structures, considering the different actions that they can act so that the energy transmission can take place in a safe, efficient and economical way. In this context, this paper deals with the static and dynamic response evaluation of guyed lattice metallic towers of TLs submitted to broken conductor and broken lightning-rod cable. For this purpose, numerical models of isolated tower and completed stretch of TL were developed, including all components, towers, insulator chains, lightning-rod cables, conductor cables and prestressed stay cables. The models are subjected to static analysis through "equivalent static loads" application coming from the broken cable and dynamic analysis, in the time domain, using the method of direct integration of equations of motion explicitly, with central finite differences. The static and dynamic responses are given in terms of nodal displacements at the top and normal stresses in the tower bars. Finally, the results of the dynamic analysis in the time domain are compared with those determined by static analysis, employed in the usual design practice, that is, the peak and final values (after the structure stops vibrating) are compared with the static response values in order to verify the validity of the model used in the design practice. / A necessidade da expansão de linhas de transmissão (LTs) fica evidente diante da crescente demanda por energia elétrica no Brasil, além das grandes distâncias entre a geração e os centros consumidores. Sabe-se que inúmeras falhas podem interromper a transmissão da energia, e o colapso de torres está entre as principais, causando enormes prejuízos, tanto as companhias de energia elétrica quanto aos consumidores finais. Dessa forma, torna-se fundamental o conhecimento do comportamento mecânico destas estruturas, frente às diferentes ações que podem atuar, para que o transporte de energia ocorra de maneira segura, eficiente e econômica. Neste contexto, o presente trabalho trata da avaliação da resposta estática e dinâmica de torres metálicas treliçadas estaiadas de LTs submetidas à ruptura de cabo condutor e de cabo para-raios. Para isso, foram desenvolvidos modelos numéricos de torre isolada e de um trecho completo de LT, incluindo todos componentes, torres, cadeias de isoladores, cabos para-raios, cabos condutores e cabos de estais pré-tensionados. Os modelos são submetidos a análises estáticas, através da aplicação de cargas estáticas equivalentes oriundas da ruptura de cabo, e análises dinâmicas, no domínio do tempo, empregando o método de integração direta das equações do movimento, de forma explícita, com diferenças finitas centrais. As respostas estática e dinâmica são dadas em termos de deslocamentos nodais no topo e esforços normais nas barras da torre. Por fim, os resultados da análise dinâmica no domínio do tempo são comparados com os determinados por meio da análise estática, empregada na prática usual de projeto, ou seja, são comparados os valores de pico e os valores finais (após a estrutura parar de vibrar) com os valores da resposta estática, a fim de verificar a validade do modelo empregado na prática de projeto.
|
204 |
Modelo neural por padrões proximais de aprendizagem para automação personalizada de conteúdos didáticosMelo, Francisco Ramos de 25 March 2012 (has links)
This study presents a model for the organization of educational content customized for environments of individual studies. For many students the availability of content in general form can not be efficient. It proposed a multilevel structure of concepts to provide the development of different combinations to show the same content. The work shows that it is possible to customize the content in order to encourage other students with the use of proximal learning standards. These patterns are obtained from the analysis of the action of students with positive results in the individual organization of the content. The formal representation establishes the definition of the student profile, multi-level content, the distribution plan of correction of concepts and teaching career. The structure of the trajectory of student teaching is formally established by the method of finite differences. The system uses artificial intelligence techniques to organize and personalize content reactively. Customization is provided by an artificial neural network that enables the classification of the student profile and assign that profile to a standard proximal learning. To mediate and adjust the contents of a reactive system was inserted into a set of rules from experts in teaching. The experiment showed the applicability and appropriateness of the proposed model. The results indicated the suitability of the approach by automating the organization\'s custom content so adaptive and reactive. The intelligent system to establish the structure of the custom content to be presented was considered efficient, giving the student a better use of the content, with higher and lower final average study time and content presented. / Este trabalho apresenta uma modelagem para a organização personalizada de conteúdos didáticos para ambientes de estudos individuais. Para muitos estudantes a disponibilização do conteúdo em formato generalizado pode não ser eficiente. É proposta uma estrutura multinível de conceitos para proporcionar o desenvolvimento de diferentes combinações para a apresentação do mesmo conteúdo. O trabalho mostra que é possível personalizar o conteúdo de forma a favorecer outros estudantes com o uso de padrões proximais de aprendizagem. Estes padrões são obtidos da análise da ação de estudantes, com resultados positivos na organização individual do conteúdo. A representação formal estabelece a definição do perfil do estudante, o conteúdo multinível, o plano de distribuição dos conceitos e a correção da trajetória didática. A estruturação da trajetória didática do estudante é formalmente estabelecida pelo método das diferenças finitas. O sistema utiliza técnicas de inteligência artificial para organizar e personalizar reativamente o conteúdo. A personalização é proporcionada por uma rede neural artificial que possibilita a classificação do perfil do estudante e associa esse perfil a um padrão proximal de aprendizagem. Para mediar e ajustar o conteúdo de forma reativa foi inserido no sistema um conjunto de regras de especialistas em docência. O experimento realizado mostrou a aplicabilidade e a adequação da modelagem proposta. Os resultados indicaram a adequação da abordagem, automatizando a organização personalizada do conteúdo de forma adaptativa e reativa. O sistema inteligente ao estabelecer a estruturação personalizada do conteúdo a ser apresentado foi considerado eficiente, proporcionando ao estudante um melhor aproveitamento do conteúdo, com maior média final e menor tempo de estudo e conteúdo apresentado. / Doutor em Ciências
|
205 |
Contribution to thermal radiation to dust flame propagation : application to aluminium dust explosions / Étude de la contribution des échanges thermiques radiatifs aux processus de propagation des flammes de poussières : application aux explosions de poussières d’aluminiumBen Moussa, Rim 20 December 2017 (has links)
Ces travaux de thèse sont consacrés à l’examen du rôle du rayonnement thermique dans le processus de propagation des flammes issues de la combustion des particules d’aluminium dans l’air. Le sujet étant complexe et d’un intérêt industriel, il nécessite de prendre en compte le couplage de nombreux phénomènes physico-chimiques afin de prédire finement les conséquences des explosions de poussières. Une analyse bibliographique approfondie est proposée, concernant les mécanismes d’inflammation et de combustion des particules d’aluminium et aussi concernant les connaissances relatives à la propagation des flammes de poussières. La question spécifique de la nature des échanges thermiques et de l’influence du rayonnement thermique est étudiée. La revue bibliographique souligne les approximations et les hypothèses simplificatrices utilisées dans la littérature permettant donc de définir les pistes d’améliorations. Compte tenu des limitations importantes concernant la physique de ces flammes, un outil de simulation de physique numérique nommé « RADIAN », proche de la simulation numérique directe, a été développé proposant un couplage fin entre les différents modes d’échanges thermiques et la combustion pour modéliser la propagation de la flamme dans un nuage de poussières. La méthode des éléments discrets (MED) est utilisée pour modéliser numériquement les échanges radiatifs entre les particules et les échanges conductifs entre gaz et particules. La méthode des différences finies est utilisée pour modéliser numériquement la conduction thermique dans la phase gazeuse et la combustion. Un modèle radiatif est proposé se basant sur la théorie de Mie sur les interactions rayonnement-particules. Les résultats des simulations sont comparés avec des solutions analytiques et des données expérimentales de la littérature. Mais en plus, une étude expérimentale est aussi conduite afin de mesurer la distribution du flux radiatif devant la flamme et la vitesse de combustion laminaire pour des flammes Méthane-Sic, Méthane-Alumine et Al-air. Un bon accord entre les simulations et les expériences est démontré. La loi de Beer-Lambert relative au transfert radiatif devant le front de flamme s’avère inapplicable et une nouvelle solution analytique est proposée. La présence de particules absorbantes du rayonnement promeut la propagation de la flamme. En particulier, il a été montré expérimentalement et confirmé numériquement que les mélanges riches d’AL-air sont susceptibles d’accélérer rapidement. / In this thesis, the role of thermal radiation in aluminum-air flames propagation is studied. The subject being complex and of industrial interest, it requires the coupling of many physiochemical phenomena to accurately predict the consequences of dust explosions. A thorough literature review is proposed about the ignition and the combustion of aluminum particles and about the available theoretical models of dust flames propagation. The specific question of the nature of thermal exchanges and the influence of thermal radiation is studied. The bibliographic review underlines the simplifying assumptions and hypotheses used in the literature making possible the definition of improvement areas. Because of the limited amount of knowledge available to address these questions, a numerical tool “RADIAN” is developed enabling an accurate coupling between the different modes of heat exchange and combustion. The Discrete Element Method (DEM) is used to numerically model the radiative exchanges between particles and the gas-particle thermal conduction. The Finite Difference method is used to numerically model the thermal conduction through the gas phase and combustion. A radiative model based on Mie theory for radiation-particles interactions is incorporated. The results of the simulations are compared with available analytical solutions and experimental data. An original experimental study is also conducted to measure the distribution of irradiance ahead of the flame front and the laminar burning velocity for methane-air-Sic, methane-air-alumina and Al-air flames. A good agreement between numerical simulations and experiments is demonstrated. The Beer-Lambert’s law for radiative transfer in front of the flame front is found to be inapplicable and a new analytical solution is proposed. The presence of absorbing particles may promote the flame propagation. In particular, it is shown experimentally and confirmed theoretically/numerically that Al-air rich mixtures are likely to rapidly accelerate.
|
206 |
Computação paralela aplicada a problemas eletromagneticos utilizando o metodo FDTD / Parallel computing applied to electromagnetic problems using the FDTD methodSantos, Carlos Henrique da Silva 08 May 2005 (has links)
Orientador: Hugo Enrique Hernandez Figueroa / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Eletrica e de Computação / Made available in DSpace on 2018-08-05T08:10:41Z (GMT). No. of bitstreams: 1
Santos_CarlosHenriquedaSilva_M.pdf: 1752834 bytes, checksum: 8ed5b0902bb130762ff802db03187fbb (MD5)
Previous issue date: 2005 / Resumo: Esse trabalho tem por objetivo desenvolver soluções computacionais de alto desempenho a um baixo custo, seguindo as propostas incentivadoras do Governo Federal para adoção de software livre. Essas soluções possibilitam simular, de maneira eficiente, os domínios computacionais de médio e grande porte utilizados no eletromagnetismo computacional. Os bons resultados obtidos nesse trabalho mostram a importância e eficiência da computação massivamente paralela utilizando cluster Beowulf para o processamento do método FDTD aplicado em estruturas complexas, porém a um baixo custo financeiro. O desempenho desse sistema ficou comprovado na realização de experimentos para analisar a SAR na cabeça humana e estudar os efeitos de estruturas metamateriais / Abstract: This work has as objective to develop high performance computational solutions to a low cost, following the stimulated proposals of the Federal Government for adoption of free software. They make possible to simulate, in efficient way, the computational domains of middle and high size useful on the computational electromagnetism. The good results gotten in these work showed the importance and efficiency of the massive parallel computation using the Beowulf cluster for the process the FDTD method applied on complex structures, however to a low financial cost. The performance of this system was proved in the realization of experiments to analyze the SAR on the human head and to study the effects of metamarial structures / Mestrado / Telecomunicações e Telemática / Mestre em Engenharia Elétrica
|
207 |
Modelagem e solução numérica de equações reação-difusão em processos biológicosRodrigues, Daiana Aparecida 29 August 2013 (has links)
Submitted by Renata Lopes (renatasil82@gmail.com) on 2016-04-11T19:27:27Z
No. of bitstreams: 1
daianaaparecidarodrigues.pdf: 8225936 bytes, checksum: 96ec323f343f92c319f4e261145f9c6a (MD5) / Approved for entry into archive by Adriana Oliveira (adriana.oliveira@ufjf.edu.br) on 2016-04-24T03:34:16Z (GMT) No. of bitstreams: 1
daianaaparecidarodrigues.pdf: 8225936 bytes, checksum: 96ec323f343f92c319f4e261145f9c6a (MD5) / Made available in DSpace on 2016-04-24T03:34:16Z (GMT). No. of bitstreams: 1
daianaaparecidarodrigues.pdf: 8225936 bytes, checksum: 96ec323f343f92c319f4e261145f9c6a (MD5)
Previous issue date: 2013-08-29 / CAPES - Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Fenômenos biológicos são todo e qualquer evento que possa ser observado nos seres vivos.
O estudo desses fenômenos permite propor explicações para o seu mecanismo, a m
de entender as causas e efeitos. Pode-se citar como exemplos de fenômenos biológicos o
comportamento das células como respiração, reprodução, metabolismo e morte celular.
Equações de reação-difusão são frequentemente utilizadas para modelar fenômenos bioló-
gicos. Sistemas de reação-difusão podem produzir padrões espaciais estáveis a partir de
uma distribuição inicial uniforme esse fenômeno é conhecido como instabilidade de Turing.
Este trabalho apresenta a análise da instabilidade de Turing bem como resultados
numéricos para a solução de três modelos biológicos, modelo de Schnakenberg, modelo de
glicólise e modelo da coagulação sanguínea. O modelo de Schnakenberg é utilizado para
descrever uma reação química autocatalítica e o modelo de glicólise é relativo ao processo
de degradação metabólica da molécula de glicose para proporcionar energia para o metabolismo
celular, esses dois modelos são frequentemente relatados na literatura. O terceiro
modelo é mais recente e descreve o fenômeno da coagulação sanguínea. Nas soluções
numéricas se utiliza o método das linhas onde a discretização espacial é feita através de
um esquema de diferenças nitas. O sistema de equações diferencias ordinárias resultante
é resolvido por um esquema de integração adaptativo, com a utilização de pacote para
computação cientí ca da linguagem Python, Scipy. / Biological phenomena are all and any event that can be observed in living beings. The
study of these phenomena enables us to propose explanations for its mechanisms in order
to understand causes and e ects. One can cite as examples of biological phenomena
the behavior of cells as respiration, reproduction, metabolism and cell death. Reactiondi
usion equations are often used to model biological phenomena. Reaction-di usion
systems can produce stable spatial patterns from a uniform initial distribution, this phenomenon
is known as Turing instability. This dissertation presents an analysis of the
Turing instability as well as numerical results for the solution of three biological models,
model Schnakenberg, model of glycolysis and model of blood coagulation. The Schnakenberg
model is used to describe an autocatalytic chemical reaction and glycolysis model
refers to the process of metabolic breakdown of the glucose molecule to provide energy for
cellular metabolism, these two models are frequently reported in the literature. The third
model is newer and describes the phenomenon of blood coagulation. The method of lines
is used in the numerical solutions, where the spatial discretization is done through a nite
di erence scheme. The resulting system of ordinary di erential equations is then solved
by an adaptive integration scheme with the use of the package for scienti c computing of
Python language, Scipy.
|
208 |
Different Approaches to Model Cover-Cracking of RC Structures due to CorrosionRoshan, Arman January 2018 (has links)
This thesis presents three different approaches to model corrosion-induced crack propagation in reinforced concrete structures. The first approach is solved numerically using finite differences to model the softening behaviour of concrete in tension. The second approach idealizes the concrete cover as either a brittle elastic or an elastoplastic material so that it may be solved using a closed-form solution. Both approaches are based on a thick-walled cylinder (TWC) analogy and consider rust compressibility and rust diffusion into cracks. The third approach uses finite element modelling to validate the application of the TWC and perform a parametric study. The results obtained using each approach are compared against each other as well as against experimental results. The TWC was found to be an appropriate analogy for the geometries and reinforcement configurations considered. Analytical models were found to provide upper and lower limits to the results based on the numerical model. The experimental data found in the literature showed reasonable agreement with predictions from the numerical and elastoplastic models.
|
209 |
Numerical analyses of segmental tunnel lining under static and dynamic loads / Analyses numériques de revêtement articulé de tunnel sous charges statique et dynamiqueDo, Ngoc Anh 07 July 2014 (has links)
Cette thèse vise à étudier le comportement de revêtement articulé du tunnel en développant une nouvelle approche numérique à la Méthode de Réaction Hyperstatique (HRM) et la production des modèles numériques en deux dimensions et trois dimensions à l'aide de la méthode des différences finies (FDM). L'étude a été traitée d'abord sous charges statiques, puis effectuée sous charges dynamiques. Tout d'abord, une étude bibliographique a été effectuée. Une nouvelle approche numérique appliquée à la méthode HRM a ensuite été développée. En même temps, un modèle numérique en deux dimensions est programmé sur les conditions de charge statique dans le but d'évaluer l'influence des joints, en termes de la distribution et des caractéristiques des joints, sur le comportement du revêtement articulé de tunnel. Après cela, des modèles complets en trois dimensions d'un seul tunnel, de deux tunnels horizontaux et de deux tunnels empilés, dans lesquels le système des joints est simulé, ont été développés. Ces modèles en trois dimensions permettent d'étudier le comportement non seulement du revêtement du tunnel, mais encore le déplacement du sol entourant le tunnel lors de l’excavation. Un modèle numérique en trois dimensions simplifié a ensuite été réalisé afin de valider la nouvelle approche numérique appliquée à la méthode HRM.Dans la dernière partie de ce mémoire, la performance du revêtement articulé du tunnel sous chargements dynamiques est prise en compte par l’analyse quasi-statique et dynamique complète en utilisant le modèle numérique en deux dimensions (FDM). Un modèle HRM a également été développé prenant en compte des charges quasi-statiques. Les différences de comportement de tunnel sous chargements statiques et sismiques sont mises en évidence et expliquées. / This PhD thesis has the aim to study the behaviour of segmental tunnel lining by developing a new numerical approach to the Hyperstatic Reaction Method (HRM) and producing two-dimensional (2D) and three-dimensional (3D) numerical models using the finite difference method (FDM). The study first deals with under static loads, and then performs under dynamic loads. Firstly, a literature review has been conducted. A new numerical approach applied to the HRM has then been developed. At the same time, a 2D numerical model is programmed regarding static loading conditions in order to evaluate the influence of the segmental joints, in terms of both joint distribution and joint stiffness characteristics, on the tunnel lining behaviour. After that, full 3D models of a single tunnel, twin horizontal tunnels and twin tunnels stacked over each other, excavated in close proximity in which the joint pattern is simulated, have been developed. These 3D models allow one to investigate the behaviour of not only the tunnel lining but also the displacement of the ground surrounding the tunnel during the tunnel excavation. A simplified 3D numerical model has then been produced in order to validate the new numerical approach applied to the HRM. In the last part of the manuscript, the performance of the segmental tunnel lining exposed to dynamic loading is taken into consideration through quasi-static and full dynamic analyses using 2D numerical models (FDM). A new HRM model has also been developed considering quasi-static loads. The differences of the tunnel behaviour under static and seismic loadings are highlighted.
|
210 |
Simulations numériques d’écoulements incompressibles interagissant avec un corps déformable : application à la nage des poissons / Numerical simulation of incompressible flows interacting with forced deformable bodies : Application to fish swimmingGhaffari Dehkharghani, Seyed Amin 15 December 2014 (has links)
Une méthode numérique précise et efficace est proposée pour la simulation de corps déformables interagissant avec un écoulement incompressible. Les équations de Navier-Stokes, considérées dans leur formulation vorticité fonction de courant, sont discrétisées temporellement et spatialement à l'aide respectivement d'un schéma d'ordre 4 de Runge-Kutta et par des différences finies compactes. Grâce à l'utilisation d'un maillage uniforme, nous proposons un nouveau solveur direct au quatrième ordre pour l'équation de Poisson, permettant de garantir l'incompressibilité au zéro machine sur une grille optimale. L'introduction d'un corps déformable dans l'écoulement de fluide est réalisée au moyen d'une méthode de pénalisation de volume. La déformation du corps est imposée par l'utilisation d'un maillage lagrangien structuré mobile qui interagit avec le fluide environnant en raison des forces hydrodynamiques et du moment (calculés sur le maillage eulérien de référence). Une loi de contrôle efficace de la courbure d'un poisson anguilliforme nageant vers une cible prescrite est proposée. La méthode numérique développée prouve son efficacité et précision tant dans le cas de la nage du poisson mais aussi plus d'un grand nombre de problèmes d'interactions fluide-structure. / We present an efficient algorithm for simulation of deformable bodies interacting with two-dimensional incompressible flows. The temporal and spatial discretizations of the Navier--Stokes equations in vorticity stream-function formulation are based on classical fourth-order Runge--Kutta and compact finite differences, respectively. Using a uniform Cartesian grid we benefit from the advantage of a new fourth-order direct solver for the Poisson equation to ensure the incompressibility constraint down to machine zero over an optimal grid. For introducing a deformable body in fluid flow, the volume penalization method is used. A Lagrangian structured grid with prescribed motion covers the deformable body which is interacting with the surrounding fluid due to the hydrodynamic forces and the torque calculated on the Eulerian reference grid. An efficient law for controlling the curvature of an anguilliform fish, swimming toward a prescribed goal, is proposed which is based on the geometrically exact theory of nonlinear beams and quaternions. Validation of the developed method shows the efficiency and expected accuracy of the algorithm for fish-like swimming and also for a variety of fluid/solid interaction problems.
|
Page generated in 0.0467 seconds