Class HierarchyAveragedDataManager provides a method of tracking and determining periodic steady states of fields. In particular, this class can be used to determine the average of a field over several instances, for example determining the mean flow of a turbulent flow field. More...
#include <ibtk/HierarchyAveragedDataManager.h>
Private Member Functions | |
| void | getFromRestart () |
| void | commonConstructor (SAMRAI::tbox::Pointer< SAMRAI::tbox::Database > input_db) |
| Registers a scratch variable with the variable database. This data is allocated and deallocated as needed. More... | |
Detailed Description
The average field is defined as \( \bar{u} = \lim_{N\rightarrow\infty}\frac{1}{N} \sum_{i = 0}^{N-1} u(x,t + i*P) \) in which \( P \) is the period. This class computes the average \( \bar{u}_N = \frac{1}{N} \sum_{i = 0}^{N-1} u(x,t + i*P) \) via the relation \( \bar{u}_N = \frac{N - 1}{N} \bar{u}_{N-1} + \frac{1}{N} u(x, t + (N-1)*P) \). The average is considered converged when the L2 norm of \( \frac{1}{N} (\bar{u}_N - u(x, t + (N-1)*P)) \) is less than some specified threshold.
Note when the period is equal to 0.0, the average field reduces to \( \bar{u} = \lim_{T\rightarrow\infty} \frac{1}{T} \int_0^T u(x, t_0 + \tau)d\tau \). In this case, it is important that snapshots be taken at consistent time intervals to ensure the above recurrence relation still holds.
If requested, this class can output visualization files corresponding to the tracked average and convergence criteria (called 'deviation'). Visualization files are written after every call to updateTimeAveragedSnapshot. Side and face centered variables are assumed to represent d-dimensional vector fields, while all other fields output a depth consistent with the provided variable in the constructor.
A common use of this class is to determine an average flow field over the course of a simulation. For example, the below code will take snapshots of a flow field at specified intervals.
Constructor & Destructor Documentation
◆ HierarchyAveragedDataManager() [1/2]
| IBTK::HierarchyAveragedDataManager::HierarchyAveragedDataManager | ( | std::string | object_name, |
| SAMRAI::tbox::Pointer< SAMRAI::hier::Variable< NDIM > > | var, | ||
| SAMRAI::tbox::Pointer< SAMRAI::tbox::Database > | input_db, | ||
| SAMRAI::tbox::Pointer< SAMRAI::hier::GridGeometry< NDIM > > | grid_geom, | ||
| bool | register_for_restart = true |
||
| ) |
The constructor for class HierarchyAveragedDataManager sets some default values and determines data centering. The expected period and number of snapshots must be set in the input database. This class assumes the snapshots will be taken at equidistant points along the periodic interval.
The input database is searched for the following keys:
- 'period_start_time' : Double that represents the beginning of the period
- 'period_end_time' : Double that represents the end of the period. For cases where the period is zero, period_end_time should be set to period_start_time.
- 'threshold' : Double to determine whether a given time point is at a periodic steady state.
- 'num_snapshots' : Integer to determine the number of snapshots taken, equally sampled between t_start and t_end.
- 'enable_logging' : Bool used to determine whether to print convergence information to the log file.
- 'output_data' : Bool used to determine whether to write visualization files for the mean and deviations.
- 'dir_dump_name' : String used to determine which folder to write visualization files to. Defaults to "./".
- 'refine_type' : String to determine which refine operator will be used by default. Defaults to "CONSERVATIVE_LINEAR_REFINE".
◆ HierarchyAveragedDataManager() [2/2]
| IBTK::HierarchyAveragedDataManager::HierarchyAveragedDataManager | ( | std::string | object_name, |
| SAMRAI::tbox::Pointer< SAMRAI::hier::Variable< NDIM > > | var, | ||
| SAMRAI::tbox::Pointer< SAMRAI::tbox::Database > | input_db, | ||
| std::set< double > | snapshot_time_points, | ||
| double | period_start_time, | ||
| double | period_end_time, | ||
| double | threshold, | ||
| SAMRAI::tbox::Pointer< SAMRAI::hier::GridGeometry< NDIM > > | grid_geom, | ||
| bool | register_for_restart = true |
||
| ) |
The constructor for class HierarchyAveragedDataManager sets some default values and determines data centering. In this constructor, the times at which snapshots are taken are set by arguments.
The input database is searched for the following keys:
- 'enable_logging' : Bool used to determine whether to print convergence information to the log file.
- 'output_data' : Bool used to determine whether to write visualization files for the mean and deviations.
- 'dir_dump_name' : String used to determine which folder to write visualization files to. Defaults to "./". Note that dir_dump_name should be unique from other directories where visualization files are written to prevent conflicts in summary files.
- 'refine_type' : String to determine which refine operator will be used by default. Defaults to "CONSERVATIVE_LINEAR_REFINE".
◆ ~HierarchyAveragedDataManager()
|
default |
Member Function Documentation
◆ clearSnapshots()
| void IBTK::HierarchyAveragedDataManager::clearSnapshots | ( | ) |
Clear the snapshots taken with this class.
◆ updateTimeAveragedSnapshot() [1/2]
|
inline |
Update the time averaged snapshot. The update is computed through the relation \( u_avg = u_avg + 1/N * (u - u_avg) \) in which N is the number of updates to the snapshot. If this the first update after clearSnapshots() is called, we simply copy data.
This function returns true if the snapshot is at a steady state, or that 1/N*||u - u_avg|| < threshold. Note that a weight patch index should be supplied to accurately compute the norm.
If specified in the constructor, this function also writes visualization files for both the mean field and the deviation.
◆ updateTimeAveragedSnapshot() [2/2]
| bool IBTK::HierarchyAveragedDataManager::updateTimeAveragedSnapshot | ( | int | u_idx, |
| double | time, | ||
| SAMRAI::tbox::Pointer< SAMRAI::hier::PatchHierarchy< NDIM > > | hierarchy, | ||
| const std::string & | mean_refine_type, | ||
| const int | wgt_idx = IBTK::invalid_index, |
||
| double | tol = 1.0e-8 |
||
| ) |
◆ isAtPeriodicSteadyState() [1/2]
|
inline |
Return true if all the tracked mean fields are at a steady state.
◆ isAtPeriodicSteadyState() [2/2]
|
inline |
Return whether the point at the specified time is at a periodic steady state. If multiple time points are found within the provided tolerance, this function returns the closest time that is less than the requested time.
◆ getSnapshotTimePoints()
|
inline |
Return the time points this class is storing.
◆ getTimePoint()
Returns the exact data point at which points are being stored in the internal map. Note that if a time point is not found within the specified tolerance, an unrecoverable error occurs.
◆ getSnapshotCache()
|
inline |
Get the SnapshotCache object.
◆ putToDatabase()
|
overridevirtual |
Write the averaged data manager to a database
Implements SAMRAI::tbox::Serializable.
◆ getFromRestart()
|
private |
◆ commonConstructor()
|
private |
Member Data Documentation
◆ d_object_name
|
private |
◆ d_var
|
private |
◆ d_scratch_idx
|
private |
◆ d_mean_refine_type
|
private |
◆ d_period_start
|
private |
◆ d_period_end
|
private |
◆ d_period_length
|
private |
◆ d_periodic_thresh
|
private |
◆ d_idx_num_updates_map
◆ d_idx_steady_state_map
◆ d_snapshot_time_pts
|
private |
◆ d_snapshot_cache
|
private |
◆ d_visit_data_writer
|
private |
◆ d_mean_var
|
private |
◆ d_mean_idx
|
private |
◆ d_dev_var
|
private |
◆ d_dev_idx
|
private |
◆ d_visit_ts
|
private |
◆ d_output_data
|
private |
◆ d_enable_logging
|
private |
The documentation for this class was generated from the following file:
- ibtk/include/ibtk/HierarchyAveragedDataManager.h
