-
Notifications
You must be signed in to change notification settings - Fork 49
Commit
This commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
Merge branch 'master' of https://github.com/lanl/LaGriT
- Loading branch information
Showing
8 changed files
with
307 additions
and
68 deletions.
There are no files selected for viewing
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,45 @@ | ||
| define / MO_H_FIELD / mo_poi_h_field | ||
| define / H_ATTRIBUTE / h_field_att | ||
|
|
||
| # Y = Ax + B | ||
| # | ||
| # Slope values to try | ||
| # 0.1 | ||
| # 0.2 | ||
| # 0.25 | ||
| # 0.29 | ||
| # | ||
| define / SLOPE / 0.29 | ||
| # Intersect | ||
| define / INTERCEPT / 0.0 | ||
| # | ||
| # Create some points and then use the distance from the points to define the user | ||
| # defined resolution field h(x,y) | ||
| # | ||
| cmo / create / mo_pts / / / tet | ||
| createpts/xyz/ 1,1,1 /0.1 0.1 0.0/0.1 0.1 0.0/1,1,1/ | ||
| createpts/xyz/ 1,1,1 /-0.5 -0.7 0.0/-0.5 -0.7 0.0/1,1,1/ | ||
|
|
||
| compute / distance_field / MO_H_FIELD / mo_pts / dfield | ||
|
|
||
| math/multiply/MO_H_FIELD/H_ATTRIBUTE/1,0,0/MO_H_FIELD/dfield/SLOPE/ | ||
| math/add/MO_H_FIELD/H_ATTRIBUTE/1,0,0/MO_H_FIELD/H_ATTRIBUTE/INTERCEPT/ | ||
|
|
||
| math / floor / MO_H_FIELD / H_ATTRIBUTE / 1 0 0 / & | ||
| MO_H_FIELD / H_ATTRIBUTE / H_SCALE | ||
| math / ceiling / MO_H_FIELD / H_ATTRIBUTE / 1 0 0 / & | ||
| MO_H_FIELD / H_ATTRIBUTE / 10H_SCALE | ||
|
|
||
| cmo / printatt / MO_H_FIELD / H_ATTRIBUTE / minmax | ||
| # | ||
| # In order to visualize the h(x,y) field, the user can output | ||
| # mo_poi_h_field and look at the attribute h_field_att. | ||
| # Output will slow things down so output should be commented out | ||
| # in production runs. | ||
| # | ||
| # This is the only place where the user has access to out MO_H_FIELD because it | ||
| # will be deleted when createpts / poisson_disk exits. | ||
| # | ||
| dump / avs / output_h_of_x_y_1.inp / MO_H_FIELD | ||
| cmo / delete / mo_pts | ||
| finish |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,49 @@ | ||
| define / MO_H_FIELD / mo_poi_h_field | ||
| define / H_ATTRIBUTE / h_field_att | ||
|
|
||
| # Y = Ax + B | ||
| # | ||
| # Slope values to try | ||
| # 0.2 works | ||
| # 0.25 works | ||
| # 0.29 works | ||
| # | ||
| define / SLOPE / 0.1 | ||
| # Intersect | ||
| define / INTERCEPT / 0.0 | ||
| # | ||
| # Create some points and then use the distance from the points to define the user | ||
| # defined resolution field h(x,y) | ||
| # | ||
| # use the sin() function to create a set of points in a wave line | ||
| cmo / create / mo_pts / / / tet | ||
| createpts / line / npoints / / / xmin,ymin,zmin / xmax,ymax,zmax / iiz,ijz,ikz / | ||
| createpts / line / 400 / / / -3.95 0.0 0.0 / 3.95 0.0 0.0 / 1 1 1 | ||
| math / sin / mo_pts / yic / 1 0 0 / mo_pts / xic | ||
| math / multiply / mo_pts / yic / 1 0 0 / mo_pts / yic / 0.25 | ||
| math / multiply / mo_pts / xic / 1 0 0 / mo_pts / xic / 0.25 | ||
|
|
||
| compute / distance_field / MO_H_FIELD / mo_pts / dfield | ||
|
|
||
| math/multiply/MO_H_FIELD/H_ATTRIBUTE/1,0,0/MO_H_FIELD/dfield/SLOPE/ | ||
| math/add/MO_H_FIELD/H_ATTRIBUTE/1,0,0/MO_H_FIELD/H_ATTRIBUTE/INTERCEPT/ | ||
|
|
||
| math / floor / MO_H_FIELD / H_ATTRIBUTE / 1 0 0 / & | ||
| MO_H_FIELD / H_ATTRIBUTE / H_SCALE | ||
| math / ceiling / MO_H_FIELD / H_ATTRIBUTE / 1 0 0 / & | ||
| MO_H_FIELD / H_ATTRIBUTE / 10H_SCALE | ||
|
|
||
| cmo / printatt / MO_H_FIELD / H_ATTRIBUTE / minmax | ||
| # | ||
| # In order to visualize the h(x,y) field, the user can output | ||
| # mo_poi_h_field and look at the attribute h_field_att. | ||
| # Output will slow things down so output should be commented out | ||
| # in production runs. | ||
| # | ||
| # | ||
| # This is the only place where the user has access to out MO_H_FIELD because it | ||
| # will be deleted when createpts / poisson_disk exits. | ||
| # | ||
| dump / avs / output_h_of_x_y_2.inp / MO_H_FIELD | ||
| cmo / delete / mo_pts | ||
| finish |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,55 @@ | ||
|
|
||
| Download this Tutorial | ||
|
|
||
| # Tutorial LaGriT Introduction | ||
|
|
||
| LaGriT consists of mesh **commands** and options that are used to create, modify, optimize, and write files. LaGriT can be run on the command line interactively so commands are entered on at a time at the prompt, or with one or more files containing the LaGriT commands. By convention we use ".lgi" for lagrit input file names and ".mlgi" for macro files called from main command file. | ||
|
|
||
| Run lagrit from the command line interactively, or put commands in a file and redirect. | ||
| ``` | ||
| lagrit | ||
| ``` | ||
| or | ||
| ``` | ||
| lagrit < input.lgi | ||
| ``` | ||
| **Important Note:** the **`finish`** command must be used to exit interactive session and must be at end of command file. | ||
|
|
||
|
|
||
|
|
||
| ## Common terms and definitions | ||
|
|
||
| ### General | ||
|
|
||
| Mesh = Grid = elements defining geometric shapes allowing numerical computations for simulations | ||
| Node = Vertice = point = mesh object nodes that can be used to connect into elements. | ||
| Element = cell = nodes are connected into element types such as line, triangle, or tetrahedral (see page) | ||
| structured mesh = elements are arranged in a regular pattern with parallel sides and and orthogonal elements. | ||
| unstructured mesh = elements arranged to fit complex geometries | ||
| Delaunay = tri or tet elements arranged such that mesh is optimized for compuational accuracy | ||
| Voronoi = the vertices of the voronoi volumes correspond to the circumcenters of the Delaunay tri and tets | ||
| control volume (CV) method = ensures conservation of mass and energy in fluid flow and heat transfer equations | ||
|
|
||
| ### LaGriT | ||
|
|
||
| Current Mesh Object = cmo = the mesh object that actions will act on by default. | ||
| imt = default name for node array with positive integer values indicating color or material | ||
| itetclr = default name for element array with positive integer values indicating color or material | ||
| itp = boundary tags | ||
| geometry = created with **`surface`** and **`region`** commands and is used to set mesh materials | ||
|
|
||
|
|
||
| ## Step 1. Create a Mesh Object | ||
|
|
||
| ## Step 2. Check and View the Mesh | ||
|
|
||
| ## Step 3. Assign materials | ||
|
|
||
| ## Step 4. Set boundary nodes | ||
|
|
||
| ## Step 5. Check and View the Mesh | ||
|
|
||
| ## Step 6. Write Mesh and FEHM Setup Files | ||
|
|
||
| ## Step 7. Create a zone list for FEHM | ||
|
|
Oops, something went wrong.