Sculpt Smoothing

Sculpt options for specifying how the mesh will be smoothed following mesh generation.

Sculpt includes a tiered approach to smoothing to improve element quality. It starts by applying smoothing to all nodes in the mesh and progressively restricts the smoothing operations to only those nodes that fall below a user-defined scaled Jacobian threshold. Default numbers of iterations and thresholds for each smoothing phase have been tuned for general use, however it may be worthwhile to adjust these parameters. The three smoothing phases include:

Laplacian Smoothing: Applied to all elements. Very inexpensive fast approach to improve quality, but can result in degraded element quality if applied to excess. A fixed default of 2 iterations is applied to all hexes. Increasing the num_laplace parameter can improve some cases, especially convex shapes.

Optimization Smoothing: Applied only to elements who's scaled Jacobian falls below the opt_threshold parameter (default 0.6) and their surrounding elements. This approach uses a more expensive optimization technique to improve regions of elements simultaneously. The max_opt_iters parameter can control the maximum number of iterations applied (default is 5). Iterations will terminate, however, if no further improvement is detected. Because this method optimizes node locations simultaneously, neighboring nodes with competing optimum can sometimes limit mesh quality.

Spot Optimization: Also known as parallel coloring, is applied only to elements who's element quality falls below the pcol_threshold parameter (default 0.2). This technique is the most expensive of the techniques, but focusses only on nodes that are immediately adjacent to poor quality hexes. Each node is smoothed independently of its neighbors, and may require a high number of iterations using the max_pcol_iters to achieve desired results. Increasing the pcol_threshold and max_pcol_iters may yield improved results.

Smoothing 		--smoothing 	-smo
  --smooth,		-S <args>	Smoothing method
  --csmooth,		-CS <args>	Curve smoothing method
  --laplacian_iters,	-LI <args>	Number of Laplacian smoothing iterations
  --max_opt_iters,	-OI <args>	Max. number of parallel Jacobi opt. iters.
  --opt_threshold,	-OT <args>	Stopping criteria for Jacobi opt. smoothing
  --curve_opt_thresh,	-COT <args>	Min metric at which curves won't be honored
  --max_pcol_iters,	-CI <args>	Max. number of parallel coloring smooth iters.
  --pcol_threshold,	-CT <args>	Stopping criteria for parallel color smooth

Sculpt Command Summary

Command: smooth     Smoothing method

Long Name:     --smooth 
Short Name:     -S 
Argument Type:   integer (0, 1, 2, 3) 
Input arguments: off (0)
                 default (1)
                 on (1)
                 fixed_bbox (2)
                 no_surface_projections (3)

Automatic adjustment of node locations following meshing to improve element quality. Controls the combined Laplacian and optimization smoothing procedures applied to volume and surface nodes (see csmooth for curve smoothing options) Uses the laplacian_iters, max_opt_iters, opt_threshold, max_pcol_iters, and pcol_threshold arguments to control the sensitivity and agressiveness of the smoothing operations. In most cases, the default options for these parameters are sufficient, however increasing iterations or threshold values, while potentially causing longer run times, may result in improved mesh quality. Smoothing will adjust the location of nodes on surfaces, projecting them to an approximated surface representation defined by interface reconstruction from volume fractions. In addition to turning smoothing on and off, the surface projection characteristics can be adjusted using the bbox_fixed and no_surface_projections options.

off (0): No volume and surface smoothing is performed.

on/default (1): (Default) Combined Laplacian/Optimization (Hybrid) smoothing both surface and volumes. Automatic boundary buffer layer improvement is performed at interior surfaces intersecting the domain boundary.

fixed_bbox (2): Uses standard hybrid smoothing procedure, however nodes at the the domain boundaries will be projected to one of the six planes of the bounding box. This option turns off the automatic boundary buffer improvement.

no_surface_projections (3): Uses the fixed_bbox method, however interior surfaces are not projected. This can result in smoother interior surface representations for microstructures models. This is effective in smoothing noisy surface data, but can potentially reduce overall volume. This method is default for microstructures file formats.

Command: csmooth     Curve smoothing method

Long Name:     --csmooth 
Short Name:     -CS 
Argument Type:   integer (0, 1, 2, ...6) 
Input arguments: off (0)
                 circle (1)
                 hermite (2)
                 average_tangent (3)
                 neighbor_surface_normal (4)
                 vfrac (5)
                 linear (6)

The csmooth option controls the smoothing method used on curves. In most cases the default should be sufficient, however it may be useful to experiment with different options. The default curve smoothing option is vfrac (5). The following curve smoothing options are available:

off (0): No curve smoothing will be performed.

circle (1): Nodes projected to a fitted circle defined current node and its two neighbors.

hermite (2): Nodes projected based on Hermite interpolation. Note that this method can only be used in serial (-j 1)

average_tangent (3): Nodes projected based on average tangent of neighbors. Note that this method may not be parallel serial consistent.

neighbor_surface_normal (4): Nodes projected based on neighboring surface normals and the resulting intersecting planes.

vfrac (5): (Default) Nodes projected to initial curve interface defined from the original volume fraction data.

linear (6): Nodes projected to the linear segment defined by the node and its two immediate neighbors.

Command: laplacian_iters     Number of Laplacian smoothing iterations

Long Name:     --laplacian_iters 
Short Name:     -LI 
Argument Type:   integer >= 0 

Number of Laplacian smoothing iterations performed when Hybrid smoothing option is used. Default value is 2.

Command: max_opt_iters     Max. number of parallel Jacobi opt. iters.

Long Name:     --max_opt_iters 
Short Name:     -OI 
Argument Type:   integer >= 0 

Indicates the maximum number of iterations of optimization-based smoothing to perform. May complete sooner if no further improvement can be made. Default is 5

Command: opt_threshold     Stopping criteria for Jacobi opt. smoothing

Long Name:     --opt_threshold 
Short Name:     -OT 
Argument Type:   floating point value (-1.0 -> 1.0) 

Indicates the value for scaled Jacobian where Optimization smoothing will be performed. Elements with scaled Jacobian less than opt_threshold and their neighbors will be smoothed. Default value is 0.6

Command: curve_opt_thresh     Min metric at which curves won't be honored

Long Name:     --curve_opt_thresh 
Short Name:     -COT 
Argument Type:   floating point value (-1.0 -> 1.0) 

Indicates the value for scaled Jacobian where if a node that falls on a curve has neighboring quads less than this value, then the smoothing will no longer honor the curve definition. Instead the optimization smoother will attempt to place the node to optimize the neighboring mesh quality, wthout regard for its placement on its owning curve. Normally this value should be set close to zero to avoid too many nodes from floating off of their owning curves, however, if mesh quality is constrained by curve geometry, setting this value higher can help to avoid bad or poor quality elements. Default for this value is 0.1.

Command: max_pcol_iters     Max. number of parallel coloring smooth iters.

Long Name:     --max_pcol_iters 
Short Name:     -CI 
Argument Type:   integer >= 0 

Maximum number of spot smoothing (also known as parallel coloring) iterations to perform. May complete sooner if no further improvement can be made. Default is 100. See also pcol_threshold.

Command: pcol_threshold     Stopping criteria for parallel color smooth

Long Name:     --pcol_threshold 
Short Name:     -CT 
Argument Type:   floating point value (-1.0 -> 1.0) 

Indicates scaled Jacobian threshold for spot smoothing (also known as parallel coloring). A parallel coloring algorithm is used to uniquely identify and isolate nodes to be improved using optimization. Default is 0.2.