1D Discretization Sample Clauses

1D Discretization. Let us consider the 1D discretization. If we take the model problem: ∂ u ∂2u ∂t − µ∂x2 = f (3.8) We replace the variable u representing the solution to the problem, with uh a piecewise polynomial. Let us introduce a mesh. h is a representative dimension of the size of the elements. We can consider uh as a sum of Lagrange piecewise polynomials. uh = Σ Uj(t)φj(x) (3.9) Here φj(x) are the Lagrange piecewise polynomials, meaning φj = 1 at node xj and 0 on all other nodes, and is piecewise linear around the node xj. Now we will look at the weak formulation: ∂t φi + = ∫ b ∂uh ∫ b ∂ uh ∂ φi ∫ b and thus Σ ∫ b dUj(t) ∫ b ∂ φj ∂ φi ∫ b ( dt φjφi + Uj ∂x = Now we can define the matrix M with the entries K, the matrix with the entries ∫ b ∂ φj ∂ φi (3.13) and vector b with entries ∫ Then we are left with the ordinary differential equation system: M + KU = b. (3.15) Under the theory of FDM we now perform a time discretization with a θ method. As stated above the FDM approximates the derivative and solves considering the following general equation: lim ∆x→0 yn+1 − yn = f (y , xn+1 ) (3.16) The theta method considers the following: 