How do I solve a PDE with space derivatives of order higher than two? For example, the equation
for a function
Introduce names for the second derivatives of u, say
Then the equation can be written
You can now solve the following equivalent system of PDEs for the variables u, P, and Q with COMSOL Multiphysics:
This system can be entered as a General form PDE for the variables u, P, and Q, with the equations:
Px, Py + Qy
For boundary conditions, consider the following examples:
u, uxx, and uyy are given on the boundary. This can be implemented by using Dirichlet conditions for u, P, and Q.
u and its normal derivative du/dn are both given on the boundary. This implies that the derivative of u in the tangential direction can also be computed. Hence, expressions for ux and uy are known on the boundary. These boundary conditions can be implemented by using a Dirichlet condition for u, and Neumann conditions for P and Q:
The supplied example model solves this system of equations with f = 1 and the following boundary conditions:
Boundaries 1 and 2: u = 0, uxx = uyy = 0
Boundary 3: u = x, uxx = 0, uyy = -x
Boundary 4: u = sin(y), ux = sin(y), uy = cos(y)
Other equations of order higher than two can be treated similarly, by introducing names for certain derivatives of u. Use as high an order of Lagrange elements as you can afford, to make the higher order derivatives of u as smooth as possible.
COMSOL makes every reasonable effort to verify the information you view on this page. Resources and documents are provided for your information only, and COMSOL makes no explicit or implied claims to their validity. COMSOL does not assume any legal liability for the accuracy of the data disclosed. Any trademarks referenced in this document are the property of their respective owners. Consult your product manuals for complete trademark details.