Spelling suggestions: "subject:"summation"" "subject:"summations""
81 |
Parameter Tuning for the NFFT Based Fast Ewald SummationNestler, Franziska 14 September 2016 (has links)
The computation of the Coulomb potentials and forces in charged particle systems under 3d-periodic boundary conditions is possible in an efficient way by utilizing the Ewald summation formulas and applying the fast Fourier transform (FFT). In this paper we consider the particle-particle NFFT (P2NFFT) approach, which is based on the fast Fourier transform for nonequispaced data (NFFT) and compare the error behaviors regarding different window functions, which are used in order to approximate the given continuous charge distribution by a mesh based charge density. Typically B-splines are applied in the scope of particle mesh methods, as for instance within the well-known particle-particle particle-mesh (P3M) algorithm. The publicly available P2NFFT algorithm allows the application of an oversampled FFT as well as the usage of different window functions. We consider for the first time also an approximation by Bessel functions and show how the resulting root mean square errors in the forces can be predicted precisely and efficiently. The results show that, if the parameters are tuned appropriately, the Bessel window function is in many cases even the better choice in terms of computational costs. Moreover, the results indicate that it is often advantageous in terms of efficiency to spend some oversampling within the NFFT while using a window function with a smaller support.
|
82 |
High order summation-by-parts methods in time and spaceLundquist, Tomas January 2016 (has links)
This thesis develops the methodology for solving initial boundary value problems with the use of summation-by-parts discretizations. The combination of high orders of accuracy and a systematic approach to construct provably stable boundary and interface procedures makes this methodology especially suitable for scientific computations with high demands on efficiency and robustness. Most classes of high order methods can be applied in a way that satisfies a summation-by-parts rule. These include, but are not limited to, finite difference, spectral and nodal discontinuous Galerkin methods. In the first part of this thesis, the summation-by-parts methodology is extended to the time domain, enabling fully discrete formulations with superior stability properties. The resulting time discretization technique is closely related to fully implicit Runge-Kutta methods, and may alternatively be formulated as either a global method or as a family of multi-stage methods. Both first and second order derivatives in time are considered. In the latter case also including mixed initial and boundary conditions (i.e. conditions involving derivatives in both space and time). The second part of the thesis deals with summation-by-parts discretizations on multi-block and hybrid meshes. A new formulation of general multi-block couplings in several dimensions is presented and analyzed. It collects all multi-block, multi-element and hybrid summation-by-parts schemes into a single compact framework. The new framework includes a generalized description of non-conforming interfaces based on so called summation-by-parts preserving interpolation operators, for which a new theoretical accuracy result is presented.
|
83 |
Algorithms for Molecular Dynamics SimulationsHedman, Fredrik January 2006 (has links)
<p>Methods for performing large-scale parallel Molecular Dynamics(MD) simulations are investigated. A perspective on the field of parallel MD simulations is given. Hardware and software aspects are characterized and the interplay between the two is briefly discussed. </p><p>A method for performing <i>ab initio </i>MD is described; the method essentially recomputes the interaction potential at each time-step. It has been tested on a system of liquid water by comparing results with other simulation methods and experimental results. Different strategies for parallelization are explored.</p><p>Furthermore, data-parallel methods for short-range and long-range interactions on massively parallel platforms are described and compared. </p><p>Next, a method for treating electrostatic interactions in MD simulations is developed. It combines the traditional Ewald summation technique with the nonuniform Fast Fourier transform---ENUF for short. The method scales as <i>N log N</i>, where <i>N </i>is the number of charges in the system. ENUF has a behavior very similar to Ewald summation and can be easily and efficiently implemented in existing simulation programs.</p><p>Finally, an outlook is given and some directions for further developments are suggested.</p>
|
84 |
Summation By Parts Finite Difference Methods with Simultaneous Approximation Terms for the Heat Equation with Discontinuous CoefficientsKåhlman, Niklas January 2019 (has links)
In this thesis we will investigate how the SBP-SAT finite difference method behave with and without an interface. As model problem, we consider the heat equation with piecewise constant coefficients. The thesis is split in two main parts. In the first part we look at the heat equation in one-dimension, and in the second part we expand the problem to a two-dimensional domain. We show how the SAT-parameters are chosen such that the scheme is dual consistent and stable. Then, we perform numerical experiments, now looking at the static case. In the one-dimensional case we see that the second order SBP-SAT method with an interface converge with an order of two, while the second order SBP-SAT method without an interface converge with an order of one.
|
85 |
Simulações numéricas de Monte Carlo aplicadas no estudo das transições de fase do modelo de Ising dipolar bidimensional / Numerical Monte Carlo simulations applied to study of phase transitions in two-dimensional dipolar Ising modelRizzi, Leandro Gutierrez 24 April 2009 (has links)
O modelo de Ising dipolar bidimensional inclui, além da interação ferromagnética entre os primeiros vizinhos, interações de longo alcance entre os momentos de dipolo magnético dos spins. A presença da interação dipolar muda completamente o sistema, apresentando um rico diagrama de fase, cujas características têm originado inúmeros estudos na literatura. Além disso, a possibilidade de explicar fenômenos observados em filmes magnéticos ultrafinos, os quais possuem diversas aplicações em àreas tecnológicas, também motiva o estudo deste modelo. O estado fundamental ferromagnético do modelo de Ising puro é alterado para uma série de fases do tipo faixas, as quais consistem em domínios ferromagnéticos de largura $h$ com magnetizações opostas. A largura das faixas depende da razao $\\delta$ das intensidades dos acoplamentos ferromagnético e dipolar. Através de simulações de Monte Carlo e técnicas de repesagem em histogramas múltiplos identificamos as temperaturas críticas de tamanho finito para as transições de fase quando $\\delta=2$, o que corresponde a $h=2$. Calculamos o calor específico e a susceptibilidade do parâmetro de ordem, no intervalo de temperaturas onde as transições são observadas, para diferentes tamanhos de rede. As técnicas de repesagem permitem-nos explorar e identificar máximos distintos nessas funções da temperatura e, desse modo, estimar as temperaturas críticas de tamanho finito com grande precisão. Apresentamos evidências numéricas da existência de uma fase nemática de Ising para tamanhos grandes de rede. Em nossas simulações, observamos esta fase para tamanhos de rede a partir de $L=48$. Para verificar o quanto a interação dipolar de longo alcance afeta as estimativas físicas, nós calculamos o tempo de autocorrelação integrado nas séries temporais da energia. Inferimos daí quão severo é o critical slowing down (decaimento lento crítico) para esse sistema próximo às transições de fase termodinâmicas. Os resultados obtidos utilizando um algoritmo de atualização local foram comparados com os resultados obtidos utilizando o algoritmo multicanônico. / Two-dimensional spin model with nearest-neighbor ferromagnetic interaction and long-range dipolar interactions exhibit a rich phase diagram, whose characteristics have been exploited by several studies in the recent literature. Furthermore, the possibility of explain observed phenomena in ultrathin magnetic films, which have many technological applications, also motivates the study of this model. The presence of dipolar interaction term changes the ferromagnetic ground state expected for the pure Ising model to a series of striped phases, which consist of ferromagnetic domains of width $h$ with opposite magnetization. The width of the stripes depends on the ratio $\\delta$ of the ferromagnetic and dipolar couplings. Monte Carlo simulations and reweighting multiple histograms techniques allow us to identify the finite-size critical temperatures of the phase transitions when $\\delta=2$, which corresponds to $h=2$. We calculate, for different lattice sizes, the specific heat and susceptibility of the order parameter around the transition temperatures by means of reweighting techniques. This allows us to identify in these observables, as functions of temperature, the distinct maxima and thereby to estimate the finite-size critical temperatures with high precision. We present numerical evidence of the existence of a Ising nematic phase for large lattice sizes. Our results show that simulations need to be performed for lattice sizes at least as large as $L=48$ to clearly observe the Ising nematic phase. To access how the long-range dipolar interaction may affect physical estimates we also evaluate the integrated autocorrelation time in energy time series. This allows us to infer how severe is the critical slowing down for this system with long-range interaction and nearby thermodynamic phase transitions. The results obtained using a local update algorithm are compared with results obtained using the multicanonical algorithm.
|
86 |
Development of an Instrumented and Powered Exoskeleton for the Rehabilitation of the HandAbolfathi, Peter Puya January 2008 (has links)
Doctor of Philosophy (PhD) / With improvements in actuation technology and sensory systems, it is becoming increasingly feasible to create powered exoskeletal garments that can assist with the movement of human limbs. This class of robotics referred to as human-machine interfaces will one day be used for the rehabilitation of paralysed, damaged or weak upper and lower extremities. The focus of this project was the development of an exoskeletal interface for the rehabilitation of the hands. A novel sensor was designed for use in such a device. The sensor uses simple optical mechanisms centred on a spring to measure force and position simultaneously. In addition, the sensor introduces an elastic element between the actuator and its corresponding hand joint. This will allow series elastic actuation (SEA) to improve control and safely of the system. The Hand Rehabilitation Device requires multiple actuators. To stay within volume and weight constraints, it is therefore imperative to reduce the size, mass and efficiency of each actuator without losing power. A method was devised that allows small efficient actuating subunits to work together and produce a combined collective output. This work summation method was successfully implemented with Shape Memory Alloy (SMA) based actuators. The actuation, sensory, control system and human-machine interface concepts proposed were evaluated together using a single-joint electromechanical harness. This experimental setup was used with volunteer subjects to assess the potentials of a full-hand device to be used for therapy, assessment and function of the hand. The Rehabilitation Glove aims to bring significant new benefits for improving hand function, an important aspect of human independence. Furthermore, the developments in this project may one day be used for other parts of the body helping bring human-machine interface technology into the fields of rehabilitation and therapy.
|
87 |
Algorithms for Molecular Dynamics SimulationsHedman, Fredrik January 2006 (has links)
Methods for performing large-scale parallel Molecular Dynamics(MD) simulations are investigated. A perspective on the field of parallel MD simulations is given. Hardware and software aspects are characterized and the interplay between the two is briefly discussed. A method for performing ab initio MD is described; the method essentially recomputes the interaction potential at each time-step. It has been tested on a system of liquid water by comparing results with other simulation methods and experimental results. Different strategies for parallelization are explored. Furthermore, data-parallel methods for short-range and long-range interactions on massively parallel platforms are described and compared. Next, a method for treating electrostatic interactions in MD simulations is developed. It combines the traditional Ewald summation technique with the nonuniform Fast Fourier transform---ENUF for short. The method scales as N log N, where N is the number of charges in the system. ENUF has a behavior very similar to Ewald summation and can be easily and efficiently implemented in existing simulation programs. Finally, an outlook is given and some directions for further developments are suggested.
|
88 |
Physical modelling of acoustic shallow-water communication channelsSvensson, Elin January 2007 (has links)
Akustiska kanaler för undervattenskommunikation är flervägskanaler där ljudet reflekteras från havets yta och botten och bryts vid ljudhastighetsförändringar. I grunt vatten är impulssvaret långt jämfört med symbolernas tidslängd i ett utskickat meddelande. Detta orsakar intersymbolinterferens, vilket gör det svårt att återskapa meddelandet. Denna avhandling behandlar fysikalisk modellering av kommunikationskanalen. Sådan modellering kan öka insikten om de svårigheter som finns vid design av kommunikationssystem och kan vara till hjälp vid utveckling av lämpliga modulationstekniker och avkodningsalgoritmer. Ljudutbredningen simuleras med en strålgångsmetod med reflektionskoefficienter beräknade för plana vågor. I artikel 4 och 5 utvidgas modellen till en algoritm för gaussisk strålsummation. De viktigaste vetenskapliga bidragen är följande. Artikel 1: Aktuell kunskap om impulssvaret hos kommunikationskanalen gör det betydligt lättare att tolka det mottagna meddelandet. I denna artikel studeras tidsvariabiliteten hos impulssvaret i termer av bitfelssannolikhet när ett gammalt impulssvar används för kanalutjämning. Tidsvariabiliteten visar sig variera avsevärt med mottagarpositionen, inte bara när det gäller avståndet till sändaren, utan även när det gäller placeringen i djupled. Artikel 2: En hybridmetod presenteras, där strålgång i ett avståndsberoende medium kombineras med lokal fullfältsmodellering av interaktionen med havsbottnen. Metoden används för simuleringar av akustisk kommunikation i grunt vatten. Artikel 3: För att kunna göra tillförlitliga simuleringar av ljudutbredning behöver man god kunskap om mediets ljudhastighetsprofil --- information som inte alltid är tillgänglig. I denna artikel används den hybrida strålgångsmetoden från artikel 2 för att skatta ljudhastighetsprofilen från kommunikationsdata. Miljöparametrarna som beräknats genom inversionen minskar avvikelsen mellan simulerade och observerade skattningar av impulssvaret jämfört med avvikelsen då en uppmätt, två dagar gammal ljudhastighetsprofil används vid simuleringen. Miljömodellen används också för en genomgång av alternativa käll- och mottagarpositioner. Artikel 4: Gaussisk strålsummation är en vågutbredningsmodell som liknar strålgång men kan ge korrekta resultat i strålgångens singulära områden, som skuggzoner och kaustikor. I denna artikel diskuteras hur några olika val av den komplexa strålparametern $\epsilon$ fungerar i grunda vågledare. Bäst resultat erhålls om man väljer $\epsilon$ så att strålen blir smal i punkten närmast mottagaren och får en plan vågfront där. Artikel 5: En adaptiv metod för gaussisk strålsummation i grunda vågledare presenteras. Algoritmen ger en noggrannhet som är bättre än eller minst lika bra som strålgångens, även i grunda vågledare med en starkt djupberoende ljudhastighet. / Acoustic underwater communication channels are multipath channels where sound is reflected from the surface and the bottom of the sea and refracted by sound speed variations. In shallow water, the impulse response is typically long compared to the time length of the symbols in a transmitted message. This causes inter-symbol interference, which makes the message difficult to decode. This thesis deals with physical modelling of the communication channel. Such modelling can provide insight into the difficulties of communication system design and may serve as an aid in the development of appropriate modulation techniques and decoding algorithms. The sound propagation is simulated by a ray tracing method with plane-wave reflection coefficients, in papers 4 and 5 expanded to a Gaussian beam summation algorithm. The main scientific contributions are the following. Paper 1: Up-to-date knowledge of the impulse response of the communication channel considerably simplifies the extraction of information from a detected signal. In this paper the time variability of the impulse response is studied in terms of the bit-error rate, when an old impulse response is used for channel equalisation. The time variability is found to vary significantly with the receiver position, not only in range, but also in depth. Paper 2: A hybrid raytrace method is presented, combining ray tracing in a range-dependent water column with local full-field modelling of the seabed interaction. The method is applied to simulations of acoustic communication in shallow water. Paper 3: To be able to make reliable simulations of sound propagation, one needs to know, quite accurately, the sound speed profile of the medium --- information which is not always available. In this paper the hybrid raytrace method from paper 2 is used to estimate the sound speed profile from communication data. The environmental model obtained by the inversion reduces the mismatch between the modelled and the observed impulse response estimates, compared to the mismatch with a two-days-old sound speed profile. The model is also used to investigate alternative source-receiver configurations. Paper 4: Gaussian beam summation is a wave propagation model similar to ray tracing, which can yield correct results in singular regions like shadow zones and caustic points. In this paper some different choices of the complex beam parameter $\epsilon$ are discussed for shallow waveguides. Best results are observed when $\epsilon$ is chosen so that each beam is narrow at the point where it is closest to the receiver and has a plane wavefront there. Paper 5: An adaptive method for Gaussian beam summation in shallow waveguides is presented. The algorithm yields better or at least as good accuracy as ray tracing, even in shallow waveguides with a strongly depth-dependent sound speed profile. / QC 20100819
|
89 |
Development of an Instrumented and Powered Exoskeleton for the Rehabilitation of the HandAbolfathi, Peter Puya January 2008 (has links)
Doctor of Philosophy (PhD) / With improvements in actuation technology and sensory systems, it is becoming increasingly feasible to create powered exoskeletal garments that can assist with the movement of human limbs. This class of robotics referred to as human-machine interfaces will one day be used for the rehabilitation of paralysed, damaged or weak upper and lower extremities. The focus of this project was the development of an exoskeletal interface for the rehabilitation of the hands. A novel sensor was designed for use in such a device. The sensor uses simple optical mechanisms centred on a spring to measure force and position simultaneously. In addition, the sensor introduces an elastic element between the actuator and its corresponding hand joint. This will allow series elastic actuation (SEA) to improve control and safely of the system. The Hand Rehabilitation Device requires multiple actuators. To stay within volume and weight constraints, it is therefore imperative to reduce the size, mass and efficiency of each actuator without losing power. A method was devised that allows small efficient actuating subunits to work together and produce a combined collective output. This work summation method was successfully implemented with Shape Memory Alloy (SMA) based actuators. The actuation, sensory, control system and human-machine interface concepts proposed were evaluated together using a single-joint electromechanical harness. This experimental setup was used with volunteer subjects to assess the potentials of a full-hand device to be used for therapy, assessment and function of the hand. The Rehabilitation Glove aims to bring significant new benefits for improving hand function, an important aspect of human independence. Furthermore, the developments in this project may one day be used for other parts of the body helping bring human-machine interface technology into the fields of rehabilitation and therapy.
|
90 |
Simulações numéricas de Monte Carlo aplicadas no estudo das transições de fase do modelo de Ising dipolar bidimensional / Numerical Monte Carlo simulations applied to study of phase transitions in two-dimensional dipolar Ising modelLeandro Gutierrez Rizzi 24 April 2009 (has links)
O modelo de Ising dipolar bidimensional inclui, além da interação ferromagnética entre os primeiros vizinhos, interações de longo alcance entre os momentos de dipolo magnético dos spins. A presença da interação dipolar muda completamente o sistema, apresentando um rico diagrama de fase, cujas características têm originado inúmeros estudos na literatura. Além disso, a possibilidade de explicar fenômenos observados em filmes magnéticos ultrafinos, os quais possuem diversas aplicações em àreas tecnológicas, também motiva o estudo deste modelo. O estado fundamental ferromagnético do modelo de Ising puro é alterado para uma série de fases do tipo faixas, as quais consistem em domínios ferromagnéticos de largura $h$ com magnetizações opostas. A largura das faixas depende da razao $\\delta$ das intensidades dos acoplamentos ferromagnético e dipolar. Através de simulações de Monte Carlo e técnicas de repesagem em histogramas múltiplos identificamos as temperaturas críticas de tamanho finito para as transições de fase quando $\\delta=2$, o que corresponde a $h=2$. Calculamos o calor específico e a susceptibilidade do parâmetro de ordem, no intervalo de temperaturas onde as transições são observadas, para diferentes tamanhos de rede. As técnicas de repesagem permitem-nos explorar e identificar máximos distintos nessas funções da temperatura e, desse modo, estimar as temperaturas críticas de tamanho finito com grande precisão. Apresentamos evidências numéricas da existência de uma fase nemática de Ising para tamanhos grandes de rede. Em nossas simulações, observamos esta fase para tamanhos de rede a partir de $L=48$. Para verificar o quanto a interação dipolar de longo alcance afeta as estimativas físicas, nós calculamos o tempo de autocorrelação integrado nas séries temporais da energia. Inferimos daí quão severo é o critical slowing down (decaimento lento crítico) para esse sistema próximo às transições de fase termodinâmicas. Os resultados obtidos utilizando um algoritmo de atualização local foram comparados com os resultados obtidos utilizando o algoritmo multicanônico. / Two-dimensional spin model with nearest-neighbor ferromagnetic interaction and long-range dipolar interactions exhibit a rich phase diagram, whose characteristics have been exploited by several studies in the recent literature. Furthermore, the possibility of explain observed phenomena in ultrathin magnetic films, which have many technological applications, also motivates the study of this model. The presence of dipolar interaction term changes the ferromagnetic ground state expected for the pure Ising model to a series of striped phases, which consist of ferromagnetic domains of width $h$ with opposite magnetization. The width of the stripes depends on the ratio $\\delta$ of the ferromagnetic and dipolar couplings. Monte Carlo simulations and reweighting multiple histograms techniques allow us to identify the finite-size critical temperatures of the phase transitions when $\\delta=2$, which corresponds to $h=2$. We calculate, for different lattice sizes, the specific heat and susceptibility of the order parameter around the transition temperatures by means of reweighting techniques. This allows us to identify in these observables, as functions of temperature, the distinct maxima and thereby to estimate the finite-size critical temperatures with high precision. We present numerical evidence of the existence of a Ising nematic phase for large lattice sizes. Our results show that simulations need to be performed for lattice sizes at least as large as $L=48$ to clearly observe the Ising nematic phase. To access how the long-range dipolar interaction may affect physical estimates we also evaluate the integrated autocorrelation time in energy time series. This allows us to infer how severe is the critical slowing down for this system with long-range interaction and nearby thermodynamic phase transitions. The results obtained using a local update algorithm are compared with results obtained using the multicanonical algorithm.
|
Page generated in 0.1069 seconds