dolfin team mailing list archive

[cutejeff <question109445@xxxxxxxxxxxxxxxxxxxxx>]

Question #109445 on DOLFIN changed:
https://answers.launchpad.net/dolfin/+question/109445

cutejeff posted a new comment:
from dolfin import *

# Create mesh and function space
mesh = Interval(5000,0,1)
V = FunctionSpace(mesh, "CG", 8)

f = Expression("9.0*pi*pi*sin(3.0*pi*x[0])*sin(3.0*pi*x[0])*cos(3.0*pi*x[0])*cos(3.0*pi*x[0])+9.0*pi*pi*sin(3.0*pi*x[0])")
u = Function(V)

# Define variational problem
v  = TestFunction(V)
du = TrialFunction(V)
L  = u*u*u.dx(0)*u.dx(0)*v*dx+u.dx(0)*v.dx(0)*dx - v*f*dx
a  = derivative(L, u, du)

# Sub domain for Dirichlet boundary condition 1
u0 = Constant(0.0)

class LeftDirichletBoundary(SubDomain):
    def inside(self, x, on_boundary):
	tol=1E-14
        return on_boundary and abs(x[0]) < tol

bc_1 = DirichletBC(V, u0, LeftDirichletBoundary())

# Sub domain for Dirichlet boundary condition 2
u1 = Constant(0.0)

class RightDirichletBoundary(SubDomain):
    def inside(self, x, on_boundary):
	tol=1E-14
        return on_boundary and abs(x[0]-1) < tol

bc_2 = DirichletBC(V, u1, RightDirichletBoundary())
bc=[bc_1,bc_2]

# Solve PDE and plot solution
problem = VariationalProblem(a, L, bc, nonlinear=True)
problem.solve(u)

# Save solution to file
file = File("power2.pvd")
file << u

# Plot solution
p=plot(u)
p.write_png("1.png")

-- 
You received this question notification because you are a member of
DOLFIN Team, which is an answer contact for DOLFIN.