Beta

Latest release notes for the MESHFREE beta executables

Short URL to this page: https://meshfree.pages.fraunhofer.de/docu/beta

WARNING: beta versions -- reduced test suite only

If you would like to be notified about new beta versions of MESHFREE, please send an email to meshfree-beta-subscribe@itwm.fraunhofer.de (or ask our Support team for help).

beta2025.01.0

New:

A missing geometry file now causes an error termination of MESHFREE. For backwards compatibility one can set the environment variable 'MESHFREE_IGNORE_MISSING_GEOMETRY="true"' to restore the deprecated behavior.

Tickets:

#109044: Rebuild SAMG structures less frequently to avoid potential MPI deadlocks.

Bugfixes:

The failing usage of the environment variable MESHFREE_Ucv_DEFAULT is now recovered with MESHFREE_UCV_DEFAULT.

beta2024.11.0

New:

Implemented apparent heat capacity method (ApparentHeatCapacityMethod) for inclusion of latent heat in phase change simulations.
Optional arguments added for MassSpringDashpot to directly provide a spring and damping force.
Force-mapping keyword has been changed for the RIGIDBODY solver. New keyword for RIGIDBODY_pressureToApplyOnBody is RIGIDBODY_ForceToApplyOnBody_PC.
A new force-mapping keyword RIGIDBODY_ForceToApplyOnBody_BE is added. The user provides the force components directly. RIGIDBODY_ForceToApplyOnBody_PC takes priority if both are defined for the same MOVE flag.
Deprecate COMP_SortBEintoBoxes_Version = 1 for 3D, only applicable for 2D (nue = 2 is no longer fully supported).
Deprecate ORGANIZE_ActivateBNDpoints_Version = 2.
Deprecate ORGANIZE_DistanceToBoundary_Version = 2.
Implemented a new algorithm for the version 3 of the NEIGHBOR_FilterMethod based on BEs. This algorithm filters the invalid neighbors based on the boundary elements.

Tickets:

#106304: Allow a linear solver to fail for a prescribed number of TS. Afterwards, all options will be tested
#107536: Fix error in BEmap.

Warning: RIGIDBODY_ForceToApplyOnBody_PC is the new keyword for PC-based force-mapping instead of RIGIDBODY_pressureToApplyOnBody, existing setups must be updated with this change. Known Issues:

LINEQN_solver_ScalarSystems cannot be set to 'AUTO' for both temperature and TURB in the same simulation.

beta2024.09.1

New:

Switched licensing over to FlexNet Publisher (FlexLM). Please check InstallationGuide for instructions and contact Support if you require assistance.

beta2024.09.0

New:

Windows version available (experimental).
SAVE_BE_filter: speed optimized version.
SAMG: More robust handling of singular linear systems
timeAvg() : implement the Fourier transform in order to check frequencies of acoustic sound waves.
approxY() :: setting default value 1.3 for factor_allowed_overshoot, aiming for more stable results. If extrapolation is desired, manually set it to a higher value or 0 (for unlimited).
Co-Simulation with LS-DYNA :
- Force-mapping keywords have been changed for co-simulation with LS-DYNA. New keyword for ForceToApplyOnBody is LSDYNA_ForceToApplyOnBody_BE and for PressureToApplyOnBody is LSDYNA_ForceToApplyOnBody_PC.
- A different shell thickening value can be provided per alias (or PID) for shell elements. The feature is specific to the interface with the co-simulation tool LSDYNA (experimental).
- The keyword LSDYNA_ScaleMinNodeDistance is added to merge nodepoints on MPI-borders. This value allows the user to decide the tolerance to be used for merging nodepoints on the interface with the co-simulation tool LS-DYNA.
- Solid stress and turbulent stress keywords have been changed for co-simulation with LSDYNA. New keyword for SolidStressWeight is LSDYNA_SolidStressWeight and for TurbulentForceWeight is LSDYNA_TurbulentForceWeight.
- LSDYNA_ForceToApplyOnBody_PC uses the PC index, ind_Sn to map pointcloud data on BEs, if the third optional argument > 0.0.
approxY() :: setting default value 1.3 for factor_allowed_overshoot , aiming for more stable results. If extrapolation is desired, manually set it to 0.

Tickets:

#106330: Fixed deadlock in SAMG
#107220: Allow wildcards in BErand()

Bugfixes:

SAVE_BE_filter / SAVE_filter: depending on the order of the begin_save{ environments, the filters were not properly initialized.

Warning: In the co-simulation with LS-DYNA, the keywords for force-mappings have been modified. LSDYNA_ForceToApplyOnBody_BE is the new keyword for BE-based force mapping, whereas LSDYNA_ForceToApplyOnBody_PC is the new keyword for point-cloud based force mapping. Existing setups must be updated with this change.

beta2024.07.0

New:

Enhance the functionality of BC_WettingAngle by the definition of a dynamic contribution to the contact angle of free surfaces in contact with walls.
bubble pressure can now be given as a boundary condition to the free surface EITHER to the hydrostatic OR to the dynamic pressures. Thus, one might switch off calculation of hydrostatic pressure and being able to only live with the dynmic part. See BubbleCorrectionPressure
Random numbers: DEBUG_GeneralParameter(76) is no longer used to switch between PRNG generators.

Tickets:

#96537 and #104477: new option ReadToPCStandby for EnSight ReadInPointCloud .
#105306 Version beta2024.05.0 contained a bug for RIGIDBODY .
#105577 MF will be aborted if the MKS solver makes higher number of rigid-body sub-iterations then the upper-limit specified with the parameters in RIGIDBODY_TimeIntegrationVersion.
#105774 (poor performance at startup due to extensive loops including singleton boundary elements (BND_tria &whatever& x1 y1 z1 x2 y2 z2 x3 y3 z3) ::
- the loops are skipped during read-in in case they are not needed (such as scanning of aliases does not need the loops in the begin/end_boundary_element{} - environment -> save time for interpretation of non-dedicated lines
- the BE_decouple-algorithm is called at the end of reading the singleton boundary elements -> saves a lot of computation time
- DEBUG_GeneralParameter(94)=0 will switch this improvement off (just for performance-reference)

beta2024.05.0

New:

Encryption: The algorithm changed from RSA to AES. Old RSA-encrypted files can no longer be read in.
SAMG: Include Intel MKL
OpenMPI: Updated to version 4.1.6
Printouts of elapsed wall time added in the LSDYNA.Initialization to identify bottlenecks.
Extension of the Coordinate Transformation feature when no MOVE statement is defined for MF-LSDYNA coupling.
DEBUG_GeneralParameter(89) added to improve the performance of the BE.Repair Timecheck. Active only when coupled with LSDYNA.
DEBUG_GeneralParameter(90) added to use a constant length for all BEnormals in ENSIGHT6 format.
RIGIDBODY_TimeIntegrationVersion, RIGIDBODY_TimeIntegrationPPI, RIGIDBODY_TimeIntegrationDamping, RIGIDBODY_UseCollisionModel can now also be set in USER_common_variables.
New CommandLine line option --pedantic
Introduce Ucv_Procedures for generalistic UCV tasks.
The real() options %SAVE_FreeUnit% and %SAVE_FreeUnit100% have been deprecated.

Tickets:

#96537: Improper access to an internal array prevented.
#103974: The implementation of the math function atan2() was corrected and verified.
#104477: A first memory optimization when reading in a STANDBY Pointcloud.

Bugfixes:

Big TIMECHECK output files were generated when MF-LSDYNA coupling was used.
Timecheck for the LSDYNA.Extern produced incorrect results when used with the CouplingRate.
Introduction: fix hyperlinks in linked documents, fix standard solution of Task 19 of the Training Setup.
RIGIDBODY_TimeIntegrationVersion :: bug in the reading of the parameter list: if not all optional parameters were given, the defaults for the optional parameters were overwritten.

beta2024.03.1

Bugfixes:

BEmonApprox(): this function was no longer properly called due to an incorrect warning for unknown acronyms.

beta2024.03.0

New:

ORGANIZE_DistanceToBoundary_Version :: enhancement by an additional option -> actively suppress touch-down of irreducible points (see MONITORPOINTS)
CODIV :: now a prescaling matrix $ \mathbf{B}$ can be used as well, which makes modelling a bit more easy.
additional prefactor matrix $ \mathbf{B}$ in the general formulation of CODIV .
MPI optimized version of saving to ERFHDF5, DEBUG_GeneralParameter(88)
%EQN_default% as well as %EQN_d% :: make the handling of approxY() more easy.
approxY() :: expert mode for defining weights of suppressed points (only in case of two-digit-mode of approximation order)
Non-zero COMP_TimeCheck variable will now print information on elapsed wall time in the MESHFREE startup phase.
increased the number of items, that can be used in BEmonApprox() from 5 to 10.

Tickets:

#86046 :: the two digit-mode of approximation order in approxY() was verified.

Bugfixes:

BDF2 for velocity ( V:BDF2 ) :: did not work safely before due to wrong allocation of BDF2-indices storing the velocity values of previous time steps
Compatibility between LSDYNA and USER_h_PC_refinement
for the heat production due to viscous friction, ONLY use the laminar viscosity %ind_ETA% (resulting in %ind_diss%). Before, %ind_ETA_eff% was used wrongly, taking into account also the turbulent partition of the viscosity. Maybe a minor, problem, but just for correctness.
The default of the second component of COMP_SkipHighVelocities was not always properly set.
add the heatsource term $ q$ to the %BND_ROBIN% boundary conditions, which was missing there in the formulation. Problems did not appear previously, if the control thickness was set to zero (default). Problems did also no appear, if the neat source was zero (which is the case in most applications).

beta2024.02.0

New:

%ACTIVE_init% is now obeyed during the refinement phase of USER_h_PC_refinement.
A new identifier %TOUCH_always_force% is added to ensure unconditional activation of wallpoints.
LINEQN_solver_ScalarSystems :: new option for linear systems of turbulence modelling TURB:aaaa:bbbb:nnnn
BDF2-time-integration :: finished for all items primary items (velocity, temperature, CODI, turbulence) and in the context LAGRANGE, EULERIMPL and EULER. See KepsilonAlgorithm , TurbulenceModel, TimeIntegration, CODI_Integration .

Tickets:

#99356: Activation of boundary points on moving geometry parts for EULERIMPL fixed.
#102484: reading-in of very big restart files :: in the past, having many MPI-processes together with many saved MONITORPOINTS , the reading of a restart file could have failed. Capturing round-offs in integer divisions now fixes the problem.
#102082: now, initialize undefined acronyms ($...$-variables) with a negative index. Thus, avoid consistency problems in case of very many acronym definitions (>400) in the input file

Bugfixes:

Revert default to GEOTREE2_EstablishCON_Version = ( 3, 3 ) to fix problems in some stable tests. This was last changed in beta2023.08.0.
bugfix of a hidden bug, that rarely appeared:: InDom( "Alias" ) could have lead to segemntation faults, if the index of "Alias" was bigger than 2500.
Fixed bug for %MONITORPOINTS_CREATION_ByCluster% that could lead to creation of cluster monitor points even though the creation condition is not satisfied for all points of the cluster.
Fixed bug in DROPLETPHASE_LockedPoints_Collision. There were no locked points created.

beta2024.01.0

New:

BLENDER: Add support for SAVE_CoordinateSystem (regarding positions only).

beta2023.12.1

Bugfixes:

RIGIDBODY and USER_h_PC_refinement could not properly be used together. RIGIDBODY dynamics are now fully skipped during USER_h_PC_refinement iteration.
Fixed errors in reading in input files that could cause incorrect simulation behaviour.

beta2023.12.0

Known Issues:

Experimental case insensitive reading was enabled for Ucv reading. Unfortunately, we observed negative side effects and do not recommend using this version.

New:

Additional option in DropletSource to apply a RHS expression such as propability density functions for the initial velocities
The cv variable 'STRESSTENSOR_NumberSubcycles' is deprecated.

Tickets:

#101715: events triggering restart, stop, abort, resume are no longer effective during the PC refinement phase.
#101594: Limit TurbLengthScale Not Working. If a negative sign in DEBUG_GeneralParameter(54) was used (such as DEBUG_GeneralParameter(54)=-1.5E-7), then MESHFREE steps into a specialized KEPSILON-repair function. Here, the restriction of the turbulent length scale (see COMP_development(1) digits 6...8) was not yet implemented. This now is fixed.

Bugfixes:

COMP_DoOrganizeOnlyAfterHowManyCycles: parameter was not correctly overwritten in USER_h_PC_refinement refinement phase.

beta2023.11.0

New:

New optional parameter for symmetryfaceByCluster{%KeepUserValues%} to be able to keep user defined values.
New parameter compute_FS_InitialFilling to enable free surface computations during initial filling.
ReadInPointCloud: Extended functionality so that the point cloud can also be used as initialisation for a LIQUID phase simulation.

Tickets:

#101264: Representative mass loss at chamfers.
#101268: bug aero kepsilon linear solver does not converge :: the current changes do not strictly avoid the problems of ticket #101268, however some improvements have been made. again: a physical limit/restriction for the turbulent layer thickness %ind_turbulentWallLayer% seems to be given, see section 2.2. and especially equation (2.12) in the document DOCUMATH_NumericalIntegrationOfTurbulence.pdf
#101308: Deprecate COMP_DoOrganizeOnlyAfterHowManyCycles_STANDBY and automatically skip organization steps for the STANDBY point cloud if possible (ensuring compatibility).
#101352: timeAvg() did not work in the 4-Parameter-version (beta v23.10.0).
#101369: Reproducible Random number not properly set. Use DEBUG_GeneralParameter(76) = 0.0 to use old version of random number generation.
#101428: SAVE_format BLENDER now supports ONLY:PARTICLES, GEOMETRY:OFF and GEOMETRY:ONCE options

Bugfixes:

added CV variable SAVE_ERF_retryWrite. If set to 1, MESHFREE will continuously try to write the ERF file even if another program is using the result file instead of crashing.
There was a leak of representative masses %ind_mi_rep% (see RepresentativeMassAlgorithm) in the case when points were removed due to wall distance, but have been isolated at the same time (Y%ind_nbRegularNeighbors%=1). See also ticket #101264.
ticket #101352 was clearly a bug, that is now fixed. In beta2023.10.0 , the 4-parameter-version ( timeAvg( t, u(t), index, Dt ) ) did not work. Only the 5-parameter-version worked ( timeAvg( t, u(t), index, t-dt, t) ) and can be used as workaround.
Interoperability of COMP_DoOrganizeOnlyAfterHowManyCycles and USER_h_PC_refinement: pointcloud organization was previously also skipped in the refinement phase.
Bugfix in SAVE_format BLENDER, when writing out boundary elements

beta2023.10.0

New:

ord_convection: reimplemented option for EULER. Enables users to go back to the linear-order-implementations of versions before 2022, in order to compare solutions.
Trigger option for CODI_eq (third parameter, optional), to better handle expensive evaluations.
Trigger option CODI_executeNow in general for CODI definitions.
smoothing of the representative density DefinitionRepresentativeDensity . Avoid buckling of free surfaces due to VOLUME_correction_local .
see DEBUG_GeneralParameter(82) for controlling the computation of Y%ind_Sn(1)% ... Y%ind_Sn(3)%
Upgraded SAMG to Release 2022.8
new function BErand() for generating point samples on boundary parts.
TurbulenceModels: change the default value of DEBUG_GeneralParameter(57) = 2 (previously 1) -> that pushes k and epsilon in such a way, that the turbulent viscosity will remain nearly constant in the following time cycle.
Feature USER_h_PC_refinement has been refined. (turning off rigid body solver, timecheck statistic writeout)
COMP_development(1)=91001010 (default) providing now also the option to cut the turbulent length scale with respect to the local SMOOTH_LENGTH
Add the four-quadrant inverse tangens atan2() as mathematical function to be used in user defined equations.
timeAvg() - new function in MF_Function to compute time averages of scalar properties
LINEQN_BundleSolvers - possibility to bundle linear solver items for the persomance analysis in the context of LINEQN_solver_ScalarSystems = ... AUTO: ...
SAVE_BE_NODE_ITEM - interoperability with begin_save{

Tickets:

#97289: Dropletsource for boundary elements, see docu of BErand() for application.
#98663: Save_results when velocity solver does not converge: additional option in COMP_SkipHighVelocities.
#100136: exploding value for %ind_tauW% in case of negative/huge values of %ind_TurbulentWallLayer% . %ind_TurbulentWallLayer% is now truncated to the interval (0,1] .
#100647: values not written to timestep file before the 10th timestep.
#100793: confusing/wrong output/behavior for k-epsilon concerning the linear solver automatism, see MESHFREE.AAA.InputFiles.USER_common_variables.__Parameters__.LINEQN_solver_ScalarSystems
#100796: Enhancement of solver automatism for linear systems of (hydrostatic, correction, dynamic) pressure, see LINEQN_BundleSolvers
#101264: rep. mass loss at chamfers
#100819: inconsistent computation of representative masses at %IDENT_inflow% , see BugFixes

Bugfixes:

For *.msh files exported from Gmsh in format Version 4.1 (ASCII) not all nodes were properly read.
inconsistent computation of representative masses at %IDENT_inflow% boundaries fixed. Now also take into account the points, which carry Y%ind_kob% = %BND_edge%
BEmonApprox(): changing the caching strategy for better interoperability with ERFHDF5. We have observed performance issues when using multiple different BEmonApprox()-statements with ERFHDF5.
Fix SAVE_filter for BLENDER SAVE_format.
MPI Shared Memory: Increase maximum number of boundary elements.
There was a leak of representative masses %ind_mi_rep% (see RepresentativeMassAlgorithm ) in the case when points were removed due to wall distance, but have been isolated at the same time (Y%ind_nbRegularNeighbors%=1). See also ticket #101264

beta2023.08.0

New:

New parameter USER_h_PC_refinement for iterating the pointcloud before time integration starts.
Random number generation: we have internally completely switched to a xoshiro256** algorithm for random number generation working towards platform independence.
New parameter %RepMass_DeletedByEvent% to track representative mass deleted by event statements as an addition to %RepMass_DeletedAtMetaplanes%.
new default: GEOTREE2_EstablishCON_Version = ( 3, 5 )
New CommandLine options --model-check, -mc (see ModelCheck) and --limit-timesteps, -lt.

Tickets:

#98788: Revision of AbaqusInterpolation (feature is not maintained) documentation, fix problems during writing of interpolation output.
#99205: Prevent SAMG from returning Not-a-Number if a singular filesystem is encountered in a small (<~100 points) decoupled subsystem.

Bugfixes:

SAVE_BE_ITEM: Time dependent statements were not evaluated correctly for ERFHDF5.

beta2023.07.0

New:

New EnvironmentVariable MESHFREE_SIM_ID for setting the 16-digit MESHFREE SIM ID used in the log files (also see iFPM_process_ID).
The cv variable 'who_am_I' is deprecated.
new experimental function BEmonApprox()
New options for turbulence modelling, see KepsilonAlgorithm .
Introduce DROPLETPHASE_nbTSforParceling to allow skipping of auto parceling

Tickets:

#96615: Remove upper bound for number of arguments in equations.
#99123: ReadInPointcloud: Stationary fields could not be loaded in beta2023.05.0 (if you require a workaround contact support.)
#96909: DROPLETPHASE solver LiquidLayer: Use analytical integration for the viscous forces (decceleration) in normal direction to the boundary layer

Bugfixes:

BEmon_MPIcomm_Version = 2 or 23 does not yield results when shared memory is turned off by COMP_SharedMemoryForBE = false. In this case, we switch now to the old default BEmon_MPIcomm_Version = 1.
CUMU_AVERAGE delivered wrong values, when the boundary elements were initially hit by multiple monitor points.
CUMU_MINIMUM and CUMU_MAXIMUM delivered wrong values due to a wrong initialization.

beta2023.05.0

New:

Output Level: New EnvironmentVariables MESHFREE_OUTPUT_LEVEL and CommandLine option --output-level, -ol. Assigning all existing printouts explicit levels is work in progress.
Deprecated N_addvar and usage of %ind_addvar()%. Use UserDefinedIndices instead.
Changed log folder name from 'FPM_log_ ...' to 'log___SIM_ID_'.
Changed value of internal variable %IDENT_locked% to 1.
New cumulation features for SAVE_BE_MONITOR_ITEM: %CUMU_MAXIMUM%, %CUMU_MINIMUM%, %CUMU_AVERAGE%.
SAVE_BE_MONITOR_ITEM, %CUMU_INTERVAL%: reset the cumulation interval by formula, keyword %CUMU_USER_INTERVAL%.
New cv Parameter BEmon_DefaultValue to manually set the reset value (experimental).
New option for cv Parameter BEmon_MPIcomm_Version = 3 (experimental).
Default value of BEmon_MPIcomm_Version changes from 1 to 2.
Default value of BEmap_MPIcomm_Version changes from 2 to 3.
Adaptive substepping algorithm for LSDYNA and MESHFREE coupling (experimental).

Tickets:

#95001: restart_toberemoved can now also be used with wildcards.
#95981: Fix determination of consistent surfaces (%ind_dA%) and allocation for boundary filling.
#96443: Reactivate ADMIN_TIME_INTEG.SAVE.UDPmonPrep, see NamesOfStopWatches.
#97953: Clarify the documentation of ReadInPointCloud. Only data from EnSight files can be used for STANDBY point clouds, whereas data from ASCII files trigger a point cloud quality check and subsequent termination of the simulation.
#98063: Introduce COMP_DoOrganizeOnlyAfterHowManyCycles_STANDBY (experimental) to enable skipping of STANDBY point cloud organization.
#97468: Fixed Bugs in DROPLETPHASE_tools: Locked Points were not detaching properly.

Bugfixes:

SAVE_format: TIMEACC now allows reading two digit numbers (up to a maximum of 92).

beta2023.03.0

New:

Deprecate DIFFOP_laplace = %DIFFOP_laplace_simplex%.
DROPLETPHASE: Additional parceling functionality DROPLETPHASE_ParcelingDiameterClasses.

Tickets:

#91750: Fixed unclean MPI exit.
#96042: Add option to define IDENT flags for alias definitions as a function of e.g. time, boundary condition flag, midpoint of boundary element (experimental).
#96064: BubbleAlgorithm: Keep irregular bubbles for a specified number of time steps with nb_TStoRescueInvalidBubbles
#96550: SAMG: Fixed treatment of small components where the pressure is only determined up to a constant.
#96652: OBJ Reader will by default split quads into trias, and the normal orientation is flipped.
#96790: In beta2023.02.0, the InDom() function did not work properly. This has been fixed now.

Warning: The orientation of geometry parts read in from OBJ files has been reverted compared to previous versions. It is now consistent with other geometry file formats. For existing setups, please fix the orientation with revOrient{}. To check any geometry orientation, use CONTROL_StopAfterReadingGeometry with SAVE_format code 'NN--' or SimCut.

beta2023.02.0

New:

ReadInPointCloud: add restart and options for skipping and repeating time frames (experimental).
New event handle %EVENT_DeleteResume%.
Y%ind_OrganizeDTB%: New common variable ORGANIZE_ToleranceForIsInGap for detection of %ORGANIZE_IsInGap% and new value %ORGANIZE_IsInNegativeGap%.

Tickets:

#96165: there were two bugs related to multisave ERF causing aborts, both now fixed. Varkeys for the triangulation were duplicated and a corrupt memory in combination with SAVE_BE_filter occured.

Bugfixes:

Fix license check with CommandLine option --check-license, when running with multiple MPI processes.
LSDYNA: Fix license check for floating license model when running with multiple MPI processes.

beta2023.01.1

New:

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

beta2023.01.0

New:

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.

Bugfixes:

Multiple smaller bug fixes.

beta2022.12.0

New:

New index %ind_neighbor_filter_skip%, meant purely for debugging.

Tickets:

#94890: 1) Exclude %IDENT_void% points from clustering of point cloud activated by SCAN_ClustersOfConnectivity. 2) Add missing MPI communication which lead to a simulation abort when simultaneously using %IDENT_void% and BEmap().
#94961: Fixed abort while reading Tables.
#86932: Add experimental support for SAVE_filter when SAVE_format is ASCII.

beta2022.11.3

New:

New experimental feature for automatic parceling and deparceling with COMP_DropletphaseParcels. Additional control by the user with DROPLETPHASE_MinNeighboursForParceling, DROPLETPHASE_MaxNeighboursForDeParceling, DROPLETPHASE_MaxEquivalentParcelDiameter, DROPLETPHASE_MaxRelativeParcelVelocity and DROPLETPHASE_AddParcelingCond.
interoperability of ERFHDF5 and begin_save{.

Tickets:

#94753: Allow for more INTEGRATION statements when saving in ERF format.
#94805: Fixed bug in Encryption and decryption feature.
#94820: Fixed bug in the RIGIDBODY solver (that occured only if one of the two joint-partners was not %MOVE_rigid% and not %MOVE_MassSpringDashpot%).

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.

beta2022.11.2

New:

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

beta2022.11.1

New:

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.

Tickets:

#91711: Fixed bug that sometimes caused a crash during allocating smoothlength arrays.

beta2022.11.0

New:

%MOVE_MassSpringDashpot% :: postprocess spring foces, compression state and compression rate by the function FCOG(), xCOG(), vCOG()

Tickets:

#94338: Fixed bug that produced write out 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.

Bugfixes:

RepresentativeMass_iData :: fix reading bug. Now, also parameter number 18, 19, and 20 (lower soft limit, upper soft limit, upper hard limit) are read-in correctly
RIGIDBODY_TimeIntegrationVersion :: bug correction that lead to very small explicit time steps. In RIGIDBODY_interaction , $ m_{rep}^{ij}$ was wrongly computed with the max-function (now corrected to min-function)

beta2022.10.0

Tickets:

#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% )

Bugfixes:

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

beta2022.09.0

New:

Overall code cleanup, refactoring and improvements.
Fixed typo in ORGANIZE_DeveloperCheck_Version. (Old spelling ORGANIZE_DevelopperCheck_Version still works.)
New index %ind_triggerFreeSurfaceCheck%. In order to trigger free surface check for some point, use this index and set it to 1. Previously, free surface check was also triggered automatically, if Y%ind_event_GeometricalFunctionManipulation% > 0. This is not the case anymore. In order to re-apply the previous behavior, add the following line at the very end of your USER_common_variables-file EVENT(...) = ( [1], %EVENT_FunctionManipulation%, %ind_triggerFreeSurfaceCheck%, [Y%ind_event_GeometricalFunctionManipulation%] ) Sorry for any inconvenience.
RepresentativeMassAlgorithm :: introduce parameters for lowerSoftBound, upperSoftBound, and upperHardLimit, see RepresentativeMass_iData.

Tickets:

#93146 (Bug movement of BE since beta2022.06.0) :: caused due to race-condition if using COMP_SharedMemoryForBE = true together with IDENT%IDENT_slip% (previously %BND_slip%). This race condition is fixed.
#93636: COMP_IsolatedParticles_MinNbOfNeigh and COMP_IsolatedParticles_MinNbOfInteriorNeigh can be used also optionally in USER_common_variables.

Bugfixes:

Fix sign mistake in definition of %MOVE_MassChangeRate% and %MOVE_VolumeChangeRate% and make them work with RepresentativeMass_iData (experimental).
Fix bug in SAVE_filter, when used with equation instead of index for depvar.

beta2022.08.1

New:

DEBUG_GeneralParameter(52): new option to use second order least-squares-transport operator for BC_v($...$) = ( %BND_slip%, ... ) in case of EULERIMPL or EULER.
SAVE_filter: available also for MONITORPOINTS in ENSIGHT6.

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.

beta2022.08.0

New:

ParallelReadingOfSignalFile now works without boost library and hence does not require a separate pCS executable any more.

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).

beta2022.06.0

New:

DELT_dt_variable = 2 :: new option to avoid Runge-Kutta-based time adaptivity in EULERIMPL, because this can lead (temporarily) to small time steps and overshoots in pressure.
BC_v = %BND_slip% :: Implement transport terms in the boundary condition BC_v($BC_name$) = ( %BND_slip%, ... ). See the terms $ \frac{\partial}{\partial t}$ in section 1.6 and 1.7 of DOCUMATH_NumericalIntegrationOfTurbulence.pdf.
BC_k and BC_eps :: implement a control thickness for BC_k and BC_eps (same motivation as the control thickness of BC_v=%BND_free% or BC_v=%BND_slip% ). See sections 1.4 ff in DOCUMATH_NumericalIntegrationOfTurbulence.pdf.
ORGANIZE_CheckPointsAtFS_PerformPreCheck :: additional option to define the precheck-search-radius.
filterBE{ :: new option to further reduce the list of boundary elements in case of boundary manipulation by manipulate{.
Y0() :: new functionalities to compute the y-star of turbulent boundary layer and the solid angle of local tetrahedrization.
new functionalities DEBUG_GeneralParameter(54) until DEBUG_GeneralParameter(58) -> mostly adressing EULERIMPL as well as K-EPSILON formulations

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.

Bugfixes:

COMP_DoOrganizePointsUntil and INTEGRATION :: now integrations carried out correctly if COMP_DoOrganizePointsUntil is used.

beta2022.04.1

Tickets:

#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.

beta2022.04.0

New:

New function indBA() which provides the functionality that was provided by binA(0) before beta2022.02.0. The previously suggested replacement by BEindA(Y%ind_BE1%) is now discouraged, since it can lead to unexpected results for interior points.
New option %VOLUME_DeletedByEvent% for real() function. Volume reduction due to deletion events is no longer included in %VOLUME_DeletedAtMetaplanes%.
New options %RepMass_CreatedByMEMORIZERead%, %RepMass_DeletedByMEMORIZEWrite%, %VOLUME_AddedByMEMORIZERead%, %VOLUME_DeletedByMEMORIZEWrite% for real() function.
Add sub-stopwatches for ADMIN_TIME_INTEG.SAVE, see NamesOfStopWatches.
New function pi() which returns the mathematical constant pi.
KOP(i) = ... PDYN:NONE ... will explicitly switch off the computation of the consistent dynamic pressure (see DynamicPressureAlgorithm ). Instead, the dynamic pressure is derived in the classical Chorin sense from the CorrectionPressureAlgorithm .
Computation of $ \tau_W$ can now be triggered by BCON($BCindex$,%ind_tauW%) = ( ... ). See documentation of boundary conditions: %ind_tauW%
New consistent way to define the turbulent wall layer thickness in %ind_tauW%

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.

beta2022.02.0

New:

%ind_tauWv(1)% ... %ind_tauWv(3)%: vectorial representation of wall stress $ \mathbf{\tau}_W$
In the case of SPM_Regularization_Type = 1, the same regularization technique is performed for SAMG as for other linear solvers.
New entries of DEBUG_GeneralParameter:
- DEBUG_GeneralParameter(37): now consistent formulation of reduction of Y%ind_betaDarcy% for stabilization in case of v--
- DEBUG_GeneralParameter(38): for testing/debugging, one might set the minimum allowed order for Laplacian operators
- DEBUG_GeneralParameter(41): number off smoothing loops for PSI and THETA in DynamicPressureAlgorithm
- DEBUG_GeneralParameter(42): select equilibrium or non-equilibrium ansatz for turbulent logarithmic wall layer formulation (see %ind_tauWv(1)%)

Bugfixes:

Fixed a bug which lead to BE_MAP not taking advantage of COMP_SharedMemoryForBE.

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():
- binA() and BEbinA() no longer support zero or negative arguments. The functionality has been moved to the new function BEindA().
- binA() and BEbinA() now support lists of aliases

beta2022.01.2

New:

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

Bugfixes:

Fixed bug for removeBEonCondition{} which lead to incorrect geometry manipulation in case of shared memory.

beta2022.01.1

New:

Added new function options real(%USER_h_min%) real(%USER_h_max%), see also time dependent USER_h_min and USER_h_max.

Bugfixes:

Small fixes for MPI, fixes compilation on HLRS Hawk.

beta2022.01.0

New:

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

Tickets:

#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.

Bugfixes:

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 occurring during the first time cycles, due to a missing MPI-communication.

beta2021.12.1

Bugfixes:

Small bugfixes and code refactorisation.

beta2021.12.0

New:

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

Tickets:

#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 ) .

Bugfixes:

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")

beta2021.11.0

New:

New INTEGRATION keyword: %SUMMATION_BE%. ATTENTION!!! Renaming of other INTEGRATION keywords: %INTEGRATION_BE_DIRECT%, %INTEGRATION_BE_DIRECT_TIME%, %SUMMATION_BENP%.
MONITORPOINTS: make ind_cluster also available for postprocessing only when the first entry in SCAN_ClustersOfConnectivity < 0) as well as saving + MONITORPOINTS organization is done at the end of the time step (SAVE_atEndOfTimestep = 1).

Tickets:

#85638: Added GeometryRestrictions support for OBJ files.
#86151: LINEQN_solver_ScalarSystems: Added choice of solver for pressure / Temperature / CODI computations (Note: So far this only works in USER_common_variables)

beta2021.10.2

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".

beta2021.10.1

Bugfixes:

Fixed bug in MPI communication that could cause segmentation faults.

beta2021.10.0

New:

Defaults for COMP_SharedMemoryForBE and COMP_SharedMemoryForGT2 changed to true.
New expert feature in cross() function.
For newly created points, the initial position (stored in %ind_x0(1)%, %ind_x0(2)%, %ind_x0(3)%) is now set to the initial position of the parent point. This also influences the results of cross().
%ind_mult% is now taken into account in the LiquidLayer computations.
CODI_eq and CODI_min_max: these two functions are now extendable to all inactive points (Y%ind_act%<=0) of a given material. That helps in defining the smoothing length with USER_h_funct = 'ADTV' in the case of BUBBLE_DoTheManagement > 0. That also solves ticket #86653 .
New environment begin_table{ with corresponding functions table() and table_size().
New function drand().
New experimental parcel injection method for DropletSource
Adaptation to quality parameters for laplace operator.
search tree setup for InDom() and DtDom() computations coupled to classical search tree setup
New functionality %EVENT_PublicTrigger% for point-independent EVENT statements.
Support for aliases and line breaks in EventMessage

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()).

Bugfixes:

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.

beta2021.08.1

New:

New geometry manipulation refineGeometry{.
New MOVE options for free surfaces: %MOVE_MassChangeRate% and %MOVE_VolumeChangeRate% (replacing deprecated %MOVE_foam%).

Bugfixes:

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

beta2021.08.0

New:

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.

beta2021.07.1

Tickets:

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

beta2021.07.0

New:

Performance improvements for the build of the BND-search tree (experimental). Settings to try:
- negative value for GEOTREE2_BND_FinalBoxDimension,
- relatively big value (such as 50) for GEOTREE2_BND_FinalBoxSize,
- relatively big value (such as 1.0 or bigger) for GEOTREE2_BND_FinalBoxDimension.
Default of COMP_SortBEintoBoxes_Version = 21.
Reduced memory usage of BE_MAP in combination with COMP_SharedMemoryForBE.
CODI: Add EULERIMPL approximation for solving the transport terms and second order for %CODI_explicit% (see CODI_Integration).
Arithmetic expressions can now be used in the first column of Curves.
Integration type %INTEGRATION_FLUX_DROPLETPHASE% is now deprecated. Use the more general %INTEGRATION_FLUX% instead.
New integration type %INTEGRATION_FLUX_IN% which generalizes %MASSFLOW_DROPLETPHASE%. %MASSFLOW_DROPLETPHASE% is now deprecated.
New index %ind_DropSrc%.
New index %ind_mult% for DROPLETPHASE points. The values of this index are used in several INTEGRATION statements (e.g. %INTEGRATION_INT%).

Tickets:

#84690: USER_h_min and USER_h_max can now depend on time, and be defined chamberwise.

Bugfixes:

Fixed the interaction of SAVE_atEndOfTimestep = 1, EVENT-triggered saving of results and EVENT-triggered aborting/stopping of a simulation.
Correct filling of %ind_proc% for points created by DropletSource.

beta2021.06.1

Tickets:

#85421: Improve point cloud organisation for EULER case with moving boundaries.

Bugfixes:

Fixed ERFHDF5 crashing with 0 FPM Points.
Fixed ERFHDF5 not working with ExchangeBEonRestart Feature.

beta2021.06.0

New:

Compatibility between begin_save{ and SAVE_BE_ITEM.
New ACTIVE statements during simulation: %ACTIVE_function% and %ACTIVE_interval%.
New options for STRESSTENSOR_Variante.

Tickets:

#84997: Closing of geometry during simulation.

beta2021.05.1

New:

Changing Default of FLIQUID_ConsistentPressure_Version to 1327.

Tickets:

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

beta2021.05.0

New:

New MEMORIZE functionality MEMORIZE_Duplicate.
Speedup of startup phase for large geometries.
approxY(): changing the default for the approximation order iOrder from ord_gradient(iChamber) (default = 3) to 2 to have a more stable approximation behavior in the default.
KindOfProblem : see "Advanced settings", define as time-dependent (interpretable) string or as alias index by a typical RightHandSideExpression .
Interoperability of temporal ACTIVE statement and BubbleAlgorithm.
SAVE_BE_ITEM and parameter SAVE_DefaultAttributesForBE for ERFHDF5.

Tickets:

#84684: Error in %EVENT_FunctionManipulation% : Problem when defining EVENT with numbers, e.g. EVENT(10).
#84713: %MOVE_VirtualRotation% now also works with EULERIMPL and EULER. Take care of the IDENT flag, see %MOVE_TranslationRotation%.
#84766: MOVE_TranslationRotation: %MOVE_InvokeDataCaching% not effective, if not at position 10 on the right hand side.
#84819: Fixed a bug in the search tree for the geometry elements.

Bugfixes:

Bugfix regarding ERFHDF5 and SAVE_CoordinateSystem.
Bugfix with missing geometry in ERFHDF5 in some cases.

beta2021.04.0

New:

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.

Tickets:

#84411: Fixed a bug in BE_MAP which lead to the mapped values being zero everywhere.

Bugfixes:

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.

beta2021.03.0

New:

New parameter SAVE_DefaultAttributesForBE, allows turning off saving of default BE attributes (currently for Ensight6, only).
New BE functions: BEchamber(), BEsymmetryface(), BElayer(), BEident(), BEcluster(), BEshape(), BEident() that can e.g. be used in SAVE_BE_ITEM.
By default there will no longer be an upper bound for the number of submessages an MPI message can be split into (see MPI_AllreduceMaxNbPackages).
Functionality CODIV (COnvection DIffusion of Vector valued (and matrix valued) entities. Applications for fiber orientation, multiple stress tensors (viscoelastic behavior) and similar. TODO: transport in case of EULER, EULERIMPL etc.
ORGANIZE_BE_ClusterNodesPoints_Version = 3 is the new default. Much faster read-in of STL-files. moreover, the clustering algorithm for geometry node points was improved. Now, the coming step in development is to shared memory parallelize the node clustering.
Do organization of MONITORPOINTS also at the end of the time step if SAVE_atEndOfTimestep = 1. Can be switched off by optional second entry.
Speed up the renew SMOOTH_LENGTH algorithm.
Add FAQ section in documentation.

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.

beta2021.02.1

New:

new function strcmp() for comparison of two strings.

Tickets:

#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.

Bugfixes:

Encryption: fixed issues regarding loops, encrypted common_variables.dat, include_Ucv{} for encrypted files.

beta2021.02.0

New:

DROPLETPHASE: New functionality to prevent thin film points from getting stuck in sharp edges and corners. see SlidingThroughCorners
EVENT: The MessageCode defining an EventMessage can now be named with a $$-acronym.
TemporarySwitch_InteriorToBoundary: feature to temporarily declare interior points (Y%ind_kob%=%BND_none%) as boundary points and provide appropriate boundary conditions.
COMP_SortBEintoBoxes_Version: new version number 22, better search-tree-construction performance if the geometry model becomes big (many millions of triangles)
new version of filling interior points (version number 22), to be triggered by ORGANIZE_PSTOneRefillStep3_UseFromWhichTimeStep = -22
new option for filtering neighbors regarding cluster index, see Version=3+4

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{.

Bugfixes:

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.

beta2021.01.0

New:

BCON_CNTCT: Add option to implement own repair strategies to control penetration of coupled phases, see %BND_slip_InContact% and %BND_slip_InContact_Explicit%.
Reactivate option to repair phase boundary via MOVE statement %MOVE_TowardsPhaseBoundary%.
Allow %H_spherical% also around BND_plane.

Tickets:

#81833: allow choice of begin_save environments for %EVENT_SaveResults%
#82683: Fix penetration of two different liquid phases, if the phases are explicitly coupled via %BND_slip_InContact_Explicit%.

beta2020.12.2

Bugfixes:

Further bugfix regarding the Darcy terms in boundary conditions.

beta2020.12.1

New:

Changed default for addition intervals for interior and boundary points from 3 to 1, see int_part_add, int_BND_part_add, COMP_AddInteriorParticles, COMP_AddBoundaryParticles!!!

Tickets:

#81727: Fix multiple definition of DELT_dt and DELT_dt_start in USER_common_variables (also important in case of using a __DEFAULT_configuration_file__, i.e. Ucv_DEFAULT.dat).
#82331: Fix clustering of node points and repairing of boundary elements for stl clean up, especially important for restart_additionalBE.

Bugfixes:

The Darcy terms are now correctly considered in the dynamic pressure boundary conditions according to %ind_DarcyVersion% (see %BND_wall%, %BND_slip%, %BND_inflow%)

beta2020.12.0

New:

New SAVE_type Standby which only saves standby point cloud data. Standby point clouds are no longer included in SAVE_type='PointCloud'.
New cv-parameters MPI_AllreduceMaxNbPackages and MPI_AllreducePackageSize.
New cv-parameter BE_SmoothLength_RefreshAfterHowManyCycles.
New function BEbinA() available in equations which are evaluated for boundary elements
Enhancement of SAVE_format = ASCII . Now putting headers into the first lines of all the result files, additionally write the results of SAVE_BE_ITEM and SAVE_BE_NODE_ITEM definitions into the appropriate files.
Add new integration types with respect to non-irreducible monitor points: %MAXIMUM_MP%, etc.

Tickets:

#82356: allow sequence of creation statements in MONITORPOINTS_CREATION. Index %ind_MP_type% can be used to differentiate between monitor point types.

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.

beta2020.11.0

New:

Feature 'GEOMETRY:ONCE' and 'GEOMETRY:OFF' in SAVE_format (experimental)
use negative numbers for SYMMETRYFACE in order to suppress filling of interior points from selected partitions of the boundary, see SYMMETRYFACE
The ChannelWithFilter example for DROPLETPHASE has been updated
TIMECHECK: New stopwatch for scalar reduce operations to separate these from the general reduce operations.
New parameters DROPLETPHASE_SubcyclingMPICommunication and DROPLETPHASE_SubcyclingReadGravity.
new feature in volume correction: ArtificialPointDisplacement, which only displaces points in order to gain volume correction, it does not affect pressure or velocity. See VolumeCorrection and all items below. Moreover see VOLUME_correction_local, VOLUME_correction, VOLUME_correction_FreeSurface
deprecated parameter UseBoxSystemVersion.
add new friction boundary condition %BND_FRICTION% for the stress tensor

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.

Bugfixes:

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.

beta2020.10.2

Tickets:

#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.

Bugfixes:

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.

beta2020.10.1

Tickets:

#74310: powerflow data can now be read in without problems. Reading EnsightGold format (such as begin_pointcloud{} include{...} format{ensight} ... ) caused problems if data was stored in "scalar per element" or "vector per element" format. That was clearly a bug and is fixed.
#81200: MONITORPOINTS_CREATION($Material$) = ( %MONITORPOINTS_CREATION_PenetrationOfBlindAndEmptyBoundary%, ... ) caused performance issues, especially if many MESHFREE-points are in the simulation. Also that was clearly a bug and is now fixed.

beta2020.10.0

New:

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.

Tickets:

#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 with respect to RESTART!

Bugfixes:

MEMORIZE: Fix writing of memorize data in case of RESTART or resuming (see checkpoint).

beta2020.09.0 = Stable version R2020A

Bugfixes:

EULERIMPL: Improvements and bugfixes in weakly compressible cases.
Minor bugfixes and improvements.