{"id":85,"date":"2017-04-05T17:54:00","date_gmt":"2017-04-05T17:54:00","guid":{"rendered":"https:\/\/therenaissanceraven.com\/?p=85"},"modified":"2017-04-05T17:54:00","modified_gmt":"2017-04-05T17:54:00","slug":"solving-laplaces-equation-with-java","status":"publish","type":"post","link":"https:\/\/therenaissanceraven.com\/?p=85","title":{"rendered":"Solving Laplace’s Equation with Java"},"content":{"rendered":"

To solve Laplace’s equation in Java we are going to use the finite difference method. A very good explanation on how it is derived and the math behind it can be found online [1].<\/p>\n

To start this process, we are going to create a pseudo 2D array in Java. This can be envisioned as discretizing our field into certain x-step and y-step elements (called Blocks in this code). To accomplish this, we will create an object called ValueBox where we cans store its value. Then we simply create a 2D array of the objects set to our desired size:<\/p>\n

\n\n\n\n
\n
 1
 2
 3
 4
 5
 6
 7
 8
 9
10<\/pre>\n<\/td>\n
\n
public<\/span> class<\/span> ValueBox<\/span> {<\/span>

    public<\/span> float<\/span> value =<\/span> 0<\/span>f;<\/span>
    public<\/span> boolean<\/span> boundaryCondition =<\/span> false<\/span>;<\/span>

    public<\/span> ValueBox<\/span>(<\/span>float<\/span> value){<\/span>
        this<\/span>.<\/span>value<\/span> =<\/span> value;<\/span>
    }<\/span>

}<\/span>
<\/pre>\n<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<\/div>\n
<\/p>\n
\n\n\n\n
\n
1<\/pre>\n<\/td>\n
\n
    private<\/span> ValueBox[][]<\/span> valueBoxArray =<\/span> new<\/span> valueBox[<\/span>xBlocks][<\/span>yBlocks];<\/span>
<\/pre>\n<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<\/div>\n
Then we can set our boundary conditions for the system. For Dirichlet boundary conditions, the process is quite simple, we only need set the initial values (boundary conditions) of the boxes before we solve the iteration. We also need to flag the boxes as boundary conditions so that their value remains fixed. The following code snippet will set the bottom of the box, y = 0, to BC1 and the top of the box, y = yBlocks to BC2
<\/p>\n
\n\n\n\n
\n
 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22<\/pre>\n<\/td>\n
\n
for<\/span>(<\/span>int<\/span> x =<\/span> 0<\/span>;<\/span> x <<\/span> xBlocks;<\/span> x++)<\/span> {<\/span>
            for<\/span>(<\/span>int<\/span> y =<\/span> 0<\/span>;<\/span> y <<\/span> yBlocks;<\/span> y++)<\/span> {<\/span>
                \/\/Sets all blocks to an initial value<\/span>
                valueBoxArray[<\/span>x][<\/span>y]<\/span> =<\/span> new<\/span> PathBox(<\/span>0<\/span>f);<\/span>

                \/\/at y=0 boundary condition<\/span>
                if<\/span>(<\/span>y==<\/span>0<\/span>){<\/span>
                    \/\/sets it as a boundary condition<\/span>
                    valueBoxArray[<\/span>x][<\/span>y].<\/span>boundaryCondition<\/span> =<\/span> true<\/span>;<\/span>
                    \/\/sets these blocks to boundary condition 1 (which can be any integer)<\/span>
                    valueBoxArray[<\/span>x][<\/span>y].<\/span>value<\/span> =<\/span> boundaryCondition1;<\/span>
                }<\/span>

                \/\/at y = yBlocks boundary condition<\/span>
                if<\/span>(<\/span>y==<\/span>yBlocks-1<\/b><\/span><\/span>){<\/span>
                    \/\/sets it as a boundary condition<\/span>
                    valueBoxArray[<\/span>x][<\/span>y].<\/span>boundaryCondition<\/span> =<\/span> true<\/span>;<\/span>
                    \/\/sets these blocks to boundary condition 2 (which can be any integer)<\/span>
                    valueBoxArray[<\/span>x][<\/span>y].<\/span>value<\/span> =<\/span> boundaryCondition2;<\/span>
                }<\/span>
            }<\/span>
        }<\/span>
<\/pre>\n<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<\/div>\n
Now that we have are blocks initialized and boundary conditions set, here is the for general for-loop design that you will need to iterate over:
<\/p>\n
\n\n\n\n
\n
 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57<\/pre>\n<\/td>\n
\n
for<\/span> (<\/span>int<\/span> iteration =<\/span> 0<\/span>;<\/span> iteration <<\/span> 500<\/span>;<\/span> iteration++)<\/span> {<\/span>
            for<\/span> (<\/span>int<\/span> x =<\/span> 0<\/span>;<\/span> x <<\/span> xBlocks;<\/span> x++)<\/span> {<\/span>
                for<\/span> (<\/span>int<\/span> y =<\/span> 0<\/span>;<\/span> y <<\/span> yBlocks;<\/span> y++)<\/span> {<\/span>
                    
                    \/\/If the block is not a boundary condition<\/span>
                    if<\/span> (!<\/span>valueBoxArray[<\/span>x][<\/span>y].<\/span>boundaryCondtion<\/span>)<\/span> {<\/span>
                        
                        \/\/If we are at the bottom row, but neither of the corners<\/span>
                        if<\/span> (<\/span>y ==<\/span> 0<\/span> &&<\/span> x !=<\/span> 0<\/span> &&<\/span> x !=<\/span> yBlocks -<\/span> 1<\/span>)<\/span> {<\/span>
                            valueBoxArray[<\/span>x][<\/span>y].<\/span>value<\/span> =<\/span> (<\/span>0.3333f<\/span>)<\/span> *<\/span> ((<\/span>valueBoxArray[<\/span>x +<\/span> 1<\/span>][<\/span>y].<\/span>value<\/span> +<\/span> valueBoxArray[<\/span>x -<\/span> 1<\/span>][<\/span>y].<\/span>value<\/span> +<\/span> valueBoxArray[<\/span>x][<\/span>y +<\/span> 1<\/span>].<\/span>value<\/span>));<\/span>
                        }<\/span>

                        \/\/If we are at the top row, but neither of the corners<\/span>
                        else<\/span> if<\/span> (<\/span>y ==<\/span> yBlocks -<\/span> 1<\/span> &&<\/span> x !=<\/span> 0<\/span> &&<\/span> x !=<\/span> xBlocks -<\/span> 1<\/span>)<\/span> {<\/span>
                            valueBoxArray[<\/span>x][<\/span>y].<\/span>value<\/span> =<\/span> (<\/span>0.3333f<\/span>)<\/span> *<\/span> ((<\/span>valueBoxArray[<\/span>x +<\/span> 1<\/span>][<\/span>y].<\/span>value<\/span> +<\/span> valueBoxArray[<\/span>x -<\/span> 1<\/span>][<\/span>y].<\/span>value<\/span> +<\/span> valueBoxArray[<\/span>x][<\/span>y -<\/span> 1<\/span>].<\/span>value<\/span>));<\/span>
                        }<\/span>

                        \/\/If we are at the right wall but at neither corners<\/span>
                        else<\/span> if<\/span> (<\/span>x ==<\/span> xBlocks -<\/span> 1<\/span> &&<\/span> y !=<\/span> 0<\/span> &&<\/span> y !=<\/span> yBlocks -<\/span> 1<\/span>)<\/span> {<\/span>
                            valueBoxArray[<\/span>x][<\/span>y].<\/span>value<\/span> =<\/span> (<\/span>0.3333f<\/span>)<\/span> *<\/span> ((<\/span>valueBoxArray[<\/span>x -<\/span> 1<\/span>][<\/span>y].<\/span>value<\/span> +<\/span> valueBoxArray[<\/span>x][<\/span>y -<\/span> 1<\/span>].<\/span>value<\/span> +<\/span> valueBoxArray[<\/span>x][<\/span>y +<\/span> 1<\/span>].<\/span>value<\/span>));<\/span>
                        }<\/span>

                        \/\/If we are at the left wall but at neither corners<\/span>
                        else<\/span> if<\/span> (<\/span>x ==<\/span> 0<\/span> &&<\/span> y !=<\/span> 0<\/span> &&<\/span> y !=<\/span> yBlocks -<\/span> 1<\/span>)<\/span> {<\/span>
                            valueBoxArray[<\/span>x][<\/span>y].<\/span>value<\/span> =<\/span> (<\/span>0.3333f<\/span>)<\/span> *<\/span> ((<\/span>valueBoxArray[<\/span>x +<\/span> 1<\/span>][<\/span>y].<\/span>value<\/span> +<\/span> valueBoxArray[<\/span>x][<\/span>y -<\/span> 1<\/span>].<\/span>value<\/span> +<\/span> valueBoxArray[<\/span>x][<\/span>y +<\/span> 1<\/span>].<\/span>value<\/span>));<\/span>
                        }<\/span>

                        \/\/If we are at the origin<\/span>
                        else<\/span> if<\/span> (<\/span>x ==<\/span> 0<\/span> &&<\/span> y ==<\/span> 0<\/span>)<\/span> {<\/span>
                            valueBoxArray[<\/span>x][<\/span>y].<\/span>value<\/span> =<\/span> (<\/span>0.5f<\/span>)<\/span> *<\/span> (<\/span>valueBoxArray[<\/span>x +<\/span> 1<\/span>][<\/span>y].<\/span>value<\/span> +<\/span> valueBoxArray[<\/span>x][<\/span>y +<\/span> 1<\/span>].<\/span>value<\/span>);<\/span>
                        }<\/span>

                        \/\/If we are at the bottom right corner<\/span>
                        else<\/span> if<\/span> (<\/span>x ==<\/span> xBlocks -<\/span> 1<\/span> &&<\/span> y ==<\/span> 0<\/span>)<\/span> {<\/span>
                            valueBoxArray[<\/span>x][<\/span>y].<\/span>value<\/span> =<\/span> (<\/span>0.5f<\/span>)<\/span> *<\/span> (<\/span>valueBoxArray[<\/span>x -<\/span> 1<\/span>][<\/span>y].<\/span>value<\/span> +<\/span> valueBoxArray[<\/span>x][<\/span>y +<\/span> 1<\/span>].<\/span>value<\/span>);<\/span>
                        }<\/span>

                        \/\/If we are at the top left corner<\/span>
                        else<\/span> if<\/span> (<\/span>x ==<\/span> 0<\/span> &&<\/span> y ==<\/span> yBlocks -<\/span> 1<\/span>)<\/span> {<\/span>
                            valueBoxArray[<\/span>x][<\/span>y].<\/span>value<\/span> =<\/span> (<\/span>0.5f<\/span>)<\/span> *<\/span> (<\/span>valueBoxArray[<\/span>x +<\/span> 1<\/span>][<\/span>y].<\/span>value<\/span> +<\/span> valueBoxArray[<\/span>x][<\/span>y -<\/span> 1<\/span>].<\/span>value<\/span>);<\/span>
                        }<\/span>


                        \/\/If we are at the top right corner<\/span>
                        else<\/span> if<\/span> (<\/span>x ==<\/span> xBlocks -<\/span> 1<\/span> &&<\/span> y ==<\/span> yBlocks -<\/span> 1<\/span>)<\/span> {<\/span>
                            valueBoxArray[<\/span>x][<\/span>y].<\/span>value<\/span> =<\/span> (<\/span>0.5f<\/span>)<\/span> *<\/span> (<\/span>valueBoxArray[<\/span>x -<\/span> 1<\/span>][<\/span>y].<\/span>value<\/span> +<\/span> valueBoxArray[<\/span>x][<\/span>y -<\/span> 1<\/span>].<\/span>value<\/span>);<\/span>
                        }<\/span>

                        \/\/If we are anywhere in the middle<\/span>
                        else<\/span> {<\/span>
                            valueBoxArray[<\/span>x][<\/span>y].<\/span>value<\/span> =<\/span> (<\/span>0.25f<\/span>)<\/span> *<\/span> ((<\/span>valueBoxArray[<\/span>x +<\/span> 1<\/span>][<\/span>y].<\/span>value<\/span> +<\/span> valueBoxArray[<\/span>x -<\/span> 1<\/span>][<\/span>y].<\/span>value<\/span> +<\/span> valueBoxArray[<\/span>x][<\/span>y -<\/span> 1<\/span>].<\/span>value<\/span> +<\/span> valueBoxArray[<\/span>x][<\/span>y +<\/span> 1<\/span>].<\/span>value<\/span>));<\/span>
                        }<\/span>
                    }<\/span>
                }<\/span>
            }<\/span>
        }<\/span>
    }<\/span>
<\/pre>\n<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<\/div>\n
Once the iterations are complete, you may simply print out the values or graph them from what is the valueBoxArray. 500 iterations is plenty, for smaller systems where the x-step and y-step elements are less than 100. However, if you are doing very large systems, >500 step elements, you may wish to check convergence by comparing the previous iteration value of a certain box after each loop.<\/p>\n
Here are some examples of solutions where they are color coded based on value. Green = BC1 and Blue = BC2.<\/p>\n\n\n\n\n
<\/a><\/td>\n<\/tr>\n
At Y = 0, Value = BC1 and at Y = yBlocks, Value = BC2<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n
<\/p>\n\n\n\n\n
<\/a><\/td>\n<\/tr>\n
At Y = 0 & Y = yBlocks, Value = BC1 and at X=0 & X=xBlocks, Value = BC2<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n
Citations:
http:\/\/www.public.iastate.edu\/~akmitra\/aero361\/design_web\/Laplace.pdf<\/p>\n","protected":false},"excerpt":{"rendered":"
To solve Laplace’s equation in Java we are going to use the finite difference method. A very good explanation on how it is derived and the math behind it can be found online [1]. To start this process, we are going to create a pseudo 2D array in Java. This can be envisioned as discretizing…<\/p>\n","protected":false},"author":2,"featured_media":0,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[5],"tags":[],"class_list":["post-85","post","type-post","status-publish","format-standard","hentry","category-misc"],"_links":{"self":[{"href":"https:\/\/therenaissanceraven.com\/index.php?rest_route=\/wp\/v2\/posts\/85","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/therenaissanceraven.com\/index.php?rest_route=\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/therenaissanceraven.com\/index.php?rest_route=\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/therenaissanceraven.com\/index.php?rest_route=\/wp\/v2\/users\/2"}],"replies":[{"embeddable":true,"href":"https:\/\/therenaissanceraven.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcomments&post=85"}],"version-history":[{"count":0,"href":"https:\/\/therenaissanceraven.com\/index.php?rest_route=\/wp\/v2\/posts\/85\/revisions"}],"wp:attachment":[{"href":"https:\/\/therenaissanceraven.com\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=85"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/therenaissanceraven.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=85"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/therenaissanceraven.com\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=85"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}