Wyniki wyszukiwania - BazTech

Ograniczanie wyników

1 Computer Methods in Materials Science

1 2006

Znaleziono wyników: 1

Liczba wyników na stronie

Wyniki wyszukiwania

Sortuj według:

Ogranicz wyniki do:

Dimensional reduction for a lattice-like mass-spring polymer model using hp-adaptivity

Jhurani Ch., Demkowicz L.

Computer Methods in Materials Science

2006

Vol. 6, No. 3-4

161-170

We present a method for dimensional reduction of the solution of the static equilibrium of a lattice-like network of point masses and harmonic springs. The dimensional reduction is achieved by energy-driven or goal-oriented hp-adaptivity. A discrete version of the fully-automatic hp-adaptive algorithm is used to constrain the positions of individual masses to conform to a piecewise polynomial. The lattice is divided into elements that consist of adjacent masses and springs. Element refinement is done by reducing the number of masses in an element (h-refinement) or increasing the polynomial degree (p-refinement). The necessity of a refinement and its optimal type depends on a local error estimate. The elements chosen for further refinement are those that give maximum reduction in error per number of invested degrees of freedom. The method generalizes the well-known Quasicontinuum method for molecular simulations (Tadmor, Ortiz, and Phillips 1996). The presented numerical results confirm the exponential reduction in errors measured in the energy norm or in the quantity of interest. The method captures the variation in material constants by choosing large elements in regions of smooth variation of material data resulting in a significant dimensional reduction. The method fits into the "finite element" variational framework for boundary value problems. This enables utilization of an existing hp code for PDE problems with minimal changes.

Artykuł przedstawia metodę redukcji wymiaru rozwiązania problemu statki cząsteczkowej dla sieci punktów materialnych połączonych sprężynami harmonicznymi. Redukcja rozmiaru uzyskana jest poprzez zastosowanie hp adaptacji kierowanej za pomocą normy energetycznej lub funkcjonału celu. Dyskretna wersja algorytmu automatycznej hp adaptacji użyta została w celu aproksymacji położeń poszczególnych cząstek poprzez wielomian kawałkami ciągły. Siatka obliczeniowa podzielona jest na elementy składające się z sąsiadujący cząstek i sprężyn. Elementy adaptowane są poprzez redukcje ilości cząstek w elemencie (h adaptacja) lub poprzez zwiększanie wielomianowego stopnia aproksymacji (p adaptacja), cyzje o konieczności adaptacji oraz ich rodzaju zależą od loki nych oszacowań błędów. Do dalszych adaptacji wybierane są elementy, dla których uzyskamy maksymalną redukcje względem ilości dodanych stopni swobody. Metoda ta uogólnia dobrze znana metodę quasi-ciągłych symulacji molekularnych (Tadmor i in., 1996). Przedstawiona metoda numeryczna osiąga eksponencjal zbieżność błędu numerycznego mierzonego w normie energetycznej lub w normie zdefiniowanej przez funkcjonał Metoda ta automatycznie konstruuje duże elementy w obszarze o małej zmienności stałych materiałowych, i dzięki temu prowadzi do dużej redukcji rozmiaru problemu obliczeniowego. Metodę zaliczyć można do klasy metod elementów skończonych rozwiązujących sformułowanie wariacyjne dla problemu brzegowego. Dzięki temu możliwe było zastosowanie istniejących kodów hp adaptacji dla równań różniczkowych całkowych, z minimalną wymaganą ilością zmian w kodzie.