Latest release notes for the MESHFREE beta executables

Short URL to this page:

WARNING: beta versions -- reduced test suite only

If you would like to be notified about new beta versions of MESHFREE, please ask our Support team to add you to the MESHFREE-Beta mailing list.


  • Allow user provided point weights for MPI bisection (domain decomposition). See MPI_BisectionWeight .


  • DROPLETPHASE: first complete version of locked points for snow-cluster modelling, see DROPLETPHASE_tools. Especially see the convenience function for a simple DEM-implementation in BuiltIn_DEM_method.
  • Updated compilation procedure to allow support for more operating systems.
  • Multiple smaller bug fixes.


New: Tickets:


New: Tickets: Bugfixes:
  • if CONTROL_StopAfterReadingGeometry = 3, this does not work together with COMP_SharedMemoryForGT2 = true. The source of this problem is still unclear. For the time being, MESHFREE puts COMP_SharedMemoryForGT2 = false at the end of initialization in case of CONTROL_StopAfterReadingGeometry = 3.


  • Renamed variables %FPM_<name>% to %<name>%. Old version will still work. See %FPM_variablename%.


  • Previously hidden log folder is now visible as FPM_log___FPM_ID*.
  • Removed warnings about wildcard matched alias names. The alias values used can be found in FPM_log___FPM_ID*/List_of_Aliases.log.
  • #91711: Fixed bug that sometimes caused a crash during allocating smoothlength arrays.


New: Tickets:
  • #94338: Fixed bug that produced writeout strongly increasing with simulation time in case of activated time measurements (see COMP_TimeCheck) in combination with skipping of point cloud organization (see COMP_DoOrganizeOnlyAfterHowManyCycles).
  • #94444: Introduced improved handling of boundary points on nearly singular boundary triangles.
  • #94455: Removed distance check to boundary elements for DropletCollisions.
  • #94462: Fixed bug that lead to deadlocks in SAMG.
  • #94492: Extended InDom()-function to also accept the definition of a reference length for the inside check which is mandatory in case this function is used for the definition of the initial point cloud.


  • #93051: Fixed %EVENT_PerformAfterHowManyTimeCycles%-functionality such that it can be used to execute a single EVENT independent of the point cloud organization controlled by COMP_DoOrganizeOnlyAfterHowManyCycles.
  • #93788: Fixed MPI crash when mapping Monitorpoints to a large geometry.
  • #94102: Fixed crash after "reading GBR". File format GBR no longer supported.
  • #93146: Fixed (geometry faults/wrinkles generated by geometry movement) :: problem accoured due to race-condition in MPI-sharedMemoryModus (only boundary elements defined as %IDENT_slip% )
  • Fix computation of normals for isolated points which caused non-physical behavior close to walls in some cases.
  • removeIsolatedClusters{ now also works in shared memory modus.
  • bugfix on input-file-parsing: now loop constructions using $...$-variables are possible, such as
    begin_loop{"iii", 1,1,&N&} INTEGRATION($INT_&iii&_$) = ( RightHandSideExpression based on counting variable &iii& ) end_loop


New: Tickets: Bugfixes:


New: Tickets:
  • #91246: Fixed bug in the boundary conditions of isolated (or nearly isolated) free surface points.
  • #93583: Fixed bug in SAVE_BE_filter for monitor points that lead to files unreadable by ParaView.


New: Tickets:
  • #88692: Increased the maximum number of characters lines in InputFiles can have to 131072.
  • #91509: Fixed the current number of monitor points %MONITOR_NbParticles% after RESTART as well as the list of indices for monitor points that is used for saving and reading the restart file.
  • #91913: Fixed incorrect triangulation of boundary elements in ERFHDF5.
  • #92111/#92475: Fixed problems with GeometryManipulations in shared memory case (see COMP_SharedMemoryForBE).


New: Tickets:
  • #89713: Removed the upper limit for the number of DropletSource() statements that can be defined (was 100).
  • #90009: SAMG builds a new setup if the global ordering of the particles has changed.
  • #91450: resolve exploding negative values of %ind_tauW% if pressure dependent boundary layer formulation is given (k-epsilon model)
  • #91509: Fix ind_ID for interior points that are projected onto a boundary as well as for slip boundary points that are released from the boundary. Fix merging of ind_ID and indC_... in case of irreducible points created by %MONITORPOINTS_CREATION_IrreducibleFPMpoint% (ind_ID and indC_... values of irreducible points should always be dominant during merging).
  • #91805: Fix extrude{ ... } for BND_quad.


  • #91250: Restartbug in version v22.04.0 :: This was a severe bug, that is not only restricted to restart. It came in with the stopwatch "ADMIN_TIME_INTEG.SAVE.ResetInfo "
  • #91177: Fixed ERF BE ID issues.
  • #91298: Fixed ERF BE ID issues and the resulting crash on HAWK.


New: Tickets:
  • #90008: Pamcrash files now identified by .inc or .pc suffixes of the file names. PART names can now start in the fifth column.
  • #90571: Add %ind_DropSrc% to indices that are saved in the RESTART file.


New: Bugfixes: Tickets:
  • #84401: Fix bugs, improve performance, and add simple examples for ExchangeBEOnRestart.
  • #85969/#87028: Increase the number of time steps that can be saved to 99999 for ENSIGHT6 and ASCII.
  • #88637: Allow comment lines before solid definitions in .stl files, see STL.
  • #88692: Changes to binA() and BEbinA():


  • MESHFREE will no longer stop execution once an undefined alias is found inside a binA() or BEbinA() statement. Instead, a warning is issued and the function will return zero.
  • DEBUG_GeneralParameter(37) :: option to control the stability of the Darcy parameters ( DarcyConstant , ForchheimerConstant ) in case of v-- .
  • SAVE_CoordinateSystem :: additional option to define the camera coordinate system
  • Fixed bug for removeBEonCondition{} which lead to incorrect geometry manipulation in case of shared memory.


New: Bugfixes:
  • Small fixes for MPI, fixes compilation on HLRS Hawk.


  • experimental function to compute the solid angle formed by the surrounding tetras of a point in the point cloud, see %EQN_SolidAngle%
  • experimental reduction of the number of ansatz functions for the neumann operator, see ord_laplace
  • removeIdenticalTriangles{ new feature in order to clean-up geometry after read-in
  • #86994: Fixed sorting of boundary points after read in of point cloud with begin_pointcloud{}.
  • #88353: Fixed bug in preparations for SAMG in case of multi-phase simulations with different IntegrationType (v-- and vp-).
  • #88255: Make sure that the initial filling of DROPLETPHASE points does not lead to overlaps (w.r.t. ind_d30). Currently, this can only be ensured on a single MPI process.
  • #89360: Fixed handling of points on %IDENT_void% boundary elements.
  • Fixed bug for BND_cube, BND_cylinder, BND_disc if index numbers of geometry nodes (STL-file, pc-file, ...) goes beyond 9999999.
  • Fixed bug in case of EULER(IMPL). Spurious point movements could have been occuring during the first time cycles, due to a missing MPI-communication.


  • Small bugfixes and code refactorisation.


  • more detailed documentation on the algorithms and meaning of variables in the context of the k-epsilon turbulence model. See especially KepsilonAlgorithm or %ind_tauW%
  • RIGIDBODY_UseCollisionModel: extend this functionality to RIGIDBODY_CollisionsModel = 1 or 2 or 3 . Also allow exclusion of collisions for dedicated rigid bodies
  • #88526: Fixed a bug which could lead to segmentation faults if large point clouds are saved using ERFHDF5 as SAVE_format.
  • #83844: Made SelectBySwitchOffFunctional more efficient.
  • #86396: Efficient reading of PAMCRASH files. The function works now much more efficient.
  • #87210: water pressing through geometry: some missing ray-intersection checks were added. Moreover, the parameter dist_LayerThickness was slightly enhanced: with the second parameter-entry, one can now limit the distance correction due ro representative masses ( RepresentativeMassAlgorithm ) .
  • valA("alias") now much faster as the numerical values (stored in "alias") are cached, such that the character-string of "alias" does not have to be processed each time valA("alias") is called
  • fix a bug that caused problems for rigid bodies in shared memory mode (racing condition, that did not lead to seg-faults, only SOMETIMES it led to strange rigid body movements)
  • restart for rigid bodies now works for RESTART_useSTREAMfile = true AND RESTART_useSTREAMfile = false (before, it worked only for "true")


New: Tickets:


  • #86992: Fixed bug caused by loops appearing after line continuations.
  • #88070: Fixed bug causing segmentation fault in case of using DtDom() or InDom() with an "undefinedAlias".


  • Fixed bug in MPI communication that could cause segmentation faults.


New: Tickets:
  • #85830: New default settings for improvements to pointcloud organization.
  • #86608: New log file List_of_GeometryNames.log
  • #86653: see New: . CODI_eq and CODI_min_max now with extension to inactive points.
  • #87005: MPI BUG -> it was not primarily an MPI bug, but the problems occur due to Y%ind_OrdApprox(1)%==1. In this case, all gradient approximations are ZERO.
  • #87162: New function erf()
  • #87419: reader identification for geometry files improved
  • #87500: Fixed erroneous values returned by max(...,nrand()).
  • Fixed a bug where equations including a cross() function could return an erroneous result
  • Fixed a bug where convolutions of Functions including nrand() could return erroneous values. See Ticket#87500.


New: Bugfixes:
  • A wrongly issued warning regarding "KOP: Set default solver" has been removed.


  • New (and updated) System Variables for time and date with different separators: @DATEb@, @DATE.@, @TIMEb@, @TIME.@, @TIMESTAMP@. @DATE@ now uses "-" as separator, according to the ISO 8601 extended calendar date format, @TIME@ is unchanged with ":" separator.
  • SAMG: Performance improvements for small point cloud components.


  • #85931: Fixed MPI bug in EULERIMPL solver.
  • #86095: Fixed timestep being incorrect in the first timestep.


New: Tickets: Bugfixes:


  • #85421: Improve point cloud organisation for EULER case with moving boundaries.
  • Fixed ERFHDF5 crashing with 0 FPM Points.
  • Fixed ERFHDF5 not working with ExchangeBEonRestart Feature.


New: Tickets:
  • #84997: Closing of geometry during simulation.


  • Changing Default of FLIQUID_ConsistentPressure_Version to 1327.
  • #85109: reallocating SmoothLengthHelp Array if not big enough instead of program exit. Handle with care! Can lead to performance issues.


New: Tickets: Bugfixes:


  • New version 3 of BEmap_PointFilter_Version which incorporates a comparison of boundary normals.
  • Performance optimizations in distance to boundary computation.
  • Update of training material in Introduction.
  • #84411: Fixed a bug in BE_MAP which lead to the mapped values being zero everywhere.
  • Fixed a bug which would lead to program termination if SimCut or SimCutBoundary were used in USER_common_variables.
  • If too many nodes were contained in an obj file, then the reallocation for nodes was not done properly.


New: Tickets:
  • #82465: Defaults for BEs in Ensight6.
  • #83430: Improve point cloud organization. Avoid filling interior points in the vicinity of intersecting geometries.
  • #83564: The syntax of BE_MAP has been changed: The second argument, previously a chamber index, is now a more general boundary element filter expression (cf. "BE-Filter" in Ticket).
  • #83564: The command BEisCloseToPC() has been introduced and may be used in the boundary element filter of BE_MAP (cf. "BND_wet" in Ticket).
  • #83564: If multiple definitions of BE_MAP exist, the neighborhood information is no longer recomputed for every definition (cf. "Stencil caching"in Ticket).
  • #83564: BE_MAP now uses an improved MPI communication, which was still an experimental feature in previous versions.
  • #83564: The PointFilterExpression in BE_MAP no longer has a reference to the boundary element, i.e. constructs like BEbinA(0)=binA(0) are not possible any more. Instead, use the hardcoded filters supplied via BEmap_PointFilter_Version which should further improve performance.
  • #83830: Fixed a bug in the DROPLETPHASE solver where the angle of boundary normal perturbation for %BND_Collision% would not be read correctly from the UCV.
  • #83975: Fixed a bug in ReadInPointCloud.
  • #84015: Fixed a bug in index management, which caused abort during point filling.
  • #84095: Optimization of memory usage in ERFHDF5 write.
  • #84104: Fixed a bug regarding TemporarySwitch_InteriorToBoundary at RESTART.
  • #84173: SimCut is now available as Ucv parameter.


  • new function strcmp() for comparison of two strings.
  • #83464: Fix the plot command in ComputationalSteering.
  • #83499: Bugfixes and improvements for SAMG solver
  • #83738: Fix segfault when using many aliasses -> memory issue when reallocating the dynamic size array.
  • Encryption: fixed issues regarding loops, encrypted common_variables.dat, include_Ucv{} for encrypted files.


New: Tickets:
  • #79561: Fixed errors in the import (via begin_pointcloud{) of Ensight Gold files which use optional element IDs.
  • #81078: Keywords for IDENT are now prefixed with IDENT_
  • #82757: Fixed a bug which caused the creation of unwanted files when using begin_save{
  • #83025: Fixed effect of IGNORE in alias declaration for corresponding boundary elements created by manipulate{...}duplicate{...} command.
  • #83202: Fixed a bug for the BEmap() function in combination with using shared memory on multiple compute nodes.
  • #83227: Fixed a bug which caused SAVE_CoordinateSystem statements to have no effect
  • #82852: (Euler Aero Explosion II): most probably, the problems mentioned in this ticket can be solved by using the new feature TemporarySwitch_InteriorToBoundary .
  • #82852: Add stabilization strategy for EULERIMPL by reducing to Upwind, if points are in the neighborhood with %ind_OrdApprox(2)% < 1.5.
  • #82996: (SCAN_ClustersOfConnectivity): Fixed a performance issue for the computation of %ind_cluster%. Also, it will be called only if point cloud organization has been executed.
  • #83342: Fixed an error where %EVENT_SaveResults% might not work properly when no begin_save{} environments are used
  • #83559: Fixed a bug which caused the code to exit when there is no definition of SAVE_format in USER_common_variables, even if valid definitions are present in an input file which is included via include_Ucv{.
  • Fixed an error where the restart path would not be set to the outer definition of SAVE_path. This occurred when using begin_save{} environments that have an interior definition of SAVE_path and precede the outer definition.
  • Fixed encryption issues - current implementation is slow and will be fixed soon
Note: As of this version, we no longer deliver executables compiled under centOS 6. This operating system has a serious security issue, for which no fix is available due to end-of-life.


New: Tickets:


  • Further bugfix regarding the Darcy terms in boundary conditions.


New: Tickets: Bugfixes:


New: Tickets: Bugfixes:
  • in all previous versions, Chorin projection ( see CorrectionPressureAlgorithm ) is suppressed, if there is a critical point, for which the velocity update  -\frac{\Delta t_{virt}^{{}}}{\left( 1+\Delta t_{\beta }^{{}}\cdot \beta  \right)}\frac{1}{\rho }\nabla c becomes too big. There was no warning, when this stability constraint took action. With this version, this constraint is taken off the code.


New: Tickets:
  • #78735: Fixed manipulations and visualizations for BND_plane
  • #81681: Fixed abort for CONTROL_StopAfterReadingGeometry with STANDBY pointcloud.
  • #81877: Fixed an error in the calculation regarding CONSTRUCT and many PIDs.
  • Three stabilizing methods for the computation of differential operators are switched off, which are not based on any scientific justification. This has improved the results in many applications. If there are still instabilities in an application, please let us know. There is still the option to turn it back on.


  • #81428: Fixed an error which would lead to a segfault when using USER_h_func = "ADTV" / "ADDS" AND approxY from a chamber AND one MPI process does not contain any points of this chamber.
  • Fixed an error where the mapping triggered by the BE_MAP command would be executed in every time step, even if SAVE_interval was chosen bigger than one. Now, the mapping is only done in time steps which also write out results.




  • New index %ind_ID% providing a particle identifier (kept over time. Experimental! Does not yet work in combination with restart feature)
  • New UCV syntax BE_MAP($$)=RHS and corresponding function BEmap().
  • DROPLETPHASE: Changed behavior of Subcycling (see COMP_DropletphaseSubcycles). Global timestep size is now only influenced by DROPLETPHASE through COEFF_dt or if negative COMP_DropletphaseSubcycles is given.
  • DROPLETPHASE: Added new CV/UCV parameter COEFF_dt_coll which implements a time step restriction for collisions defined via ParticleInteraction and %BND_COLLISION%.
  • DROPLETPHASE: If the coefficient of restitution is different in a particle-particle-collision, now the minimum of both values is taken instead of the previous heuristic choice.
  • #81153: Allow the simultaneous use of MEMORIZE_Write and MEMORIZE_Read statements with differing indices. Due to necessary changes in the restart file structure, backwards compatibility to previous MESHFREE versions is not provided for the MEMORIZE feature wrt RESTART!

beta2020.09.0 = Stable version R2020A

  • EULERIMPL: Improvements and bugfixes in weakly compressible cases.
  • Minor bugfixes and improvements.
List of members:
All Complete release notes for the MESHFREE beta executables