Stardis Solver Release Notes git repository
- The net flux imposed can be combined with other boundary/connection conditions, i.e. a net flux can be set in addition to convective exchange and radiative transfer.
- Added support for plain text log messages. Until now, log messages were intended to be read by a VT100-like terminal and could therefore contain escape sequences that required post-processing to store them in plain text log files.
- Added support for user-defined signature on the green function. It allows to check that, when reloaded, a green function is the one expected by the user according to its own constraints that the green function cannot check itself such as, for example, that the same deltas are used in conductive random walks.
- Changes the value of the constant
SDIS_VOLUMIC_POWER_NONE. Its previous value of zero caused problems during the evaluation of the propagator: a media with a power density of zero was not registered in the list of media with a volumic power. A volumic power that was not zero was therefore not taken into account during the re-evaluation of the propagator. The constant is now set to
DBL_MAX, which means that the medium has no power density, while a value of 0 is now treated as any valid power density term.
Fixed compilation errors and compilation warnings displayed on some versions of GCC.
Non linear radiative transfer
Uses a new iterative numerical method to estimate radiative transfer. With a recursion level of 1, this is equivalent to a linearization of the radiative transfer but with a reference temperature that can vary in time and space. By using a higher-order recursion, one can converge towards a rigorous estimate that takes into account the non-linearity of the radiative transfer; the higher the recursion order, the better the convergence, but with the counterpart of an increase in calculation time.
Distributed memory parallelism
Uses message passing interface to distribute computation across multiple computers. Stardis-Solver now, uses a mixed parallelism: on one computer (i.e. a node), it uses a shared memory parallelism and relies on the message passing interface to parallelize calculations between several nodes.
Type and state of the random number generator
Adds the member input variable
rng_type to the solve functions. It defines
the type of random number generator to use when no generator is defined. Note
sdis_solve_camera function does not have a random number generator
as an input variable and has therefore been updated to support it.
Reading the source code
Refactoring and deep rewriting of the source code to simplify its reading.
Fix green paths ending in a fluid (transcient computation): The path's end was not correctly registred and the path was later treated as failed.
- Sets the required version of Star-SampPling to 0.12. This version fixes compilation errors with gcc 11 but introduces API breaks.
- Fix warnings detected by gcc 11.
Updates the way numerical issues are handled during a conductive random walk. Previously, a zealous test would report a numerical error and stop the calculations when that error could be handled.
Add the support of thermal contact resistance between two solids: the new
thermal_contact_resistance functor on the data structure
sdis_interface_shader defines the thermal resistance contact in K.m^2.W^-1 at
a given time and at a specific position onto the interface.
- Add support of unsteady green evaluation. The resulting green function can then be used to quickly evaluate the system at the same time but with different limit and initial conditions, volumetric powers and imposed fluxes.
- Add checks on green re-evaluation to ensure that the system remains unchanged regarding its scale factor and its reference temperature.
- Remove the ambient radiative temperature, the reference temperature and the geometry scale factor from the list of arguments submitted to the solve functions. They become scene arguments defined on scene creation.
- Update the
sdis_scene_[2d_]createfunction profile: its data are now grouped into a variable of type
- In green function estimation, the time sent to the user callbacks is no more the elapsed time from the beginning of the realisation: as in a regular computation, it is now the observation time.
- Fix the flux computation for boundaries with an imposed flux: it was
previously ignored. The new
sdis_estimator_get_imposed_fluxfunction returns this estimated flux component.
- Return an error if the flux is computed at a boundary whose temperature is known: this configuration is not currently supported.
- Fix build with the CL compiler.
- Add support of green function [de]serialization. The
sdis_green_function_writefunction serializes the green function into a stream while the
sdis_green_function_create_from_streamfunction deserialize it. Note that the scene used to deserialize the green function must be the same of the one used to estimate it: the media and the interfaces have to be created in the same order, the scene geometry must be the same, etc.
- Add the
sdis_scene_find_closest_pointfunction: search the point onto the scene geometry that is the closest of the submitted position.
- Add the
sdis_compute_powerfunction that evaluates the power of a medium.
- Update the solver: the time of the sampled path is now rewind on solid reinjection.
- Update the API of the solve functions: the parameters of the simulation are now grouped into a unique data structure rather than separately submitted as function arguments. Thank to this structure and its default value, updating input parameters should now affect marginally the calling code.
- Improve the logger. Add a prefix to the printed text to indicate the type of the message (info, error or warning). Add a progress message during simulation.
- Bump the version of the Star-Enclosures <2D|3D> dependencies to 0.5
- Fix an issue when the
sdis_solve_boundary_fluxfunction was invoked on a boundary with radiative transfer: several sampled paths were rejected due to data inconsistencies.
- Fix a memory leak when the scene creation failed.
- Enable parallelism on Star-Enclosure[2D] to improve the performances of the enclosure extraction on the setup of the Stardis-Solver scene.
- Fix a solver issue that led to reject valid sampled paths.
- Bump the version of the Star-Enclosure[2D] libraries to 0.4.2. These versions fix a numerical issue that might led to an infinite loop at the scene creation.
- Drastically improve the robustness of the solver~: far less realisations are now rejected.
- Add the estimation of the time spent per realisation estimate. Add the
sdis_estimator_get_realisation_timefunction that returns this estimate.
- Add the
sdis_estimator_bufferAPI~: it manages a two dimensional array of regular estimators and provides global estimations over the whole estimators saved into the buffer.
- Update the signature of the
sdis_solve_camerafunction~: it now returns a
sdis_estimator_buffer. It now also supports time integration as well as heat paths registration.
Add Green function support
Provide new solve functions that compute and save the Green function, i.e. the propagator used in regular solvers. The resulting Green function can be then evaluated to obtain an estimate of the temperature.
The time spent to compute the Green function is comparable to the computation time of regular solvers; actually, they rely on the same code. However, its evaluation is instantaneous while it still handles the limit conditions, the boundary fluxes and the power term of the media at the moment of the evaluation. This means that one can estimate the Green function of a system only one time and then evaluate it with different limit conditions, boundary fluxes or power terms with negligible computation costs.
Currently, Stardis-Solver assumes that during the Green function estimation, the properties of the system do not depend on time. In addition, it assumes that the boundary fluxes and the volumetric powers are constants in time and space. Anyway, on Green function evaluation, the limit conditions of the system can still vary in time and space; systems in steady state can be simulated with Green functions.
Add heat path registration
int register_paths mask to almost all solve functions to enable the
registration against the returned estimator of the failure and/or successful
random paths used by the solvers. For each path, the registered data are:
- the vertices of the path;
- the type of the path (failed or succeed);
- the type of the path vertices (conductive, convective or radiative);
- the Monte-Carlo weight of each path vertex;
- the current time of each path vertex.
Note that the amount of registered data can be huge if too more paths are registered. Consequently, this functionality should be used with few realisations to obtain a subset of representative paths, or to only register the few paths that failed in order to diagnose what went wrong.
- Add the
sdis_solve_mediumfunction: it estimates the average temperature of a medium.
- Fix the setup of the interfaces: the interface associated to a geometric primitive could not be the right one.
- Bump version of the Star-Enclosures[2D] dependencies: the new versions fix issues in the construction of fluid enclosures.
- Bump version of the Star-<2D|3D> dependencies: the new versions rely on Embree3 rather than on Embree2 for their ray-tracing back-end.
- Add the
sdis_solve_boundaryfunction: it computes the average temperature on a subset of geometric primitives.
- Add flux solvers: the new
sdis_solve_boundary_fluxfunctions estimate the convective and radiative fluxes at a given surface position or for a sub-set of geometric primitives, respectively.
- Add support of time integration: almost all solvers can estimate the average
temperature on a given time range. Only the
sdis_solve_camerafunction does not support time integration, yet.
- Add support of an explicit initial time
t0for the fluid.
- Fix a bug in the estimation of unknown fluid temperatures: the associativity between the internal Stardis-Solver data and the user defined data was wrong.
Add support of fluid enclosure with unknown uniform temperature.
- The convection coefficient of the surfaces surrounding a fluid whose
temperature is unknown can vary in time and space. Anyway, the caller has to
ensure that for each triangle of the fluid enclosure, the convection
coefficient returned by its
struct sdis_interface_shader- at a given position and time - is less than or equal to the
convection_coef_upper_boundparameter of the shader.
Full rewrite of how the volumetric power is taken into account.
- Change the scheme of the random walk "solid re-injection": use a 2D re-injection scheme in order to handle 2D effects. On one hand, this scheme drastically improves the accuracy of the temperature estimation in solid with a volumetric power term. On the other hand it is more sensible to numerical imprecisions. The previous 1D scheme is thus used in situations where the 2D scheme exhibits too numerical issues, i.e. on sharp angles.
- Add the missing volumetric power term on solid re-injection.
- Add a corrective term to fix the bias on the volumetric power introduced when
the random walk progresses at a distance of
deltaof a boundary.
- Add several volumetric power tests.
- Remove the
delta_boundaryparameter of the
struct sdis_solid_shaderdata structure.
- Some interface properties become double sided: the temperature, emissivity
and specular fraction is defined for each side of the interface. Actually,
only the convection coefficient is shared by the 2 sides of the interface.
The per side interface properties are grouped into the new
struct sdis_interface_side_shaderdata structure.
- Add the support of fixed fluxes: the flux is a per side interface property. Currently, the flux is handled only for the interface sides facing a solid medium.
- Add the
sdis_scene_boundary_project_posfunction that computes the parametric coordinates of a world space position projected onto a given primitive with respect to its normal. If the projection lies outside the primitive, its parametric coordinates are wrapped against its boundaries in order to ensure that they are valid coordinates into the primitive. Actually, this function was mainly added to help in the definition of the probe position onto a boundary as expected by the
- Update the default comportment of the interface shader when a function is not set.
- Rename the
- Rename the
enum sdis_side_flagenumerate in
enum sdis_sideand update its values.
- Add the support of volumic power to solid media: add the
volumic_powerfunctor to the
sdis_solid_shaderdata structure that, once defined, should return the volumic power of the solid at a specific position and time. On solve invocation, the conductive random walks take into account this spatio-temporal volumic power in the computation of the solid temperature.
- Add the
sdis_solve_probe_boundaryfunction: it computes the temperature at a given position and time onto a geometric primitive. The probe position is defined by the index of the primitive and a parametric coordinates onto it.
- Add the
sdis_scene_get_boundary_positionfunction: it computes a world space position from the index of a geometric primitive and a parametric coordinate onto it.
- Fix how the
sdis_solve_probewas parallelised. The submitted
threads_hintparameter was not correctly handled.
- Add the support of radiative temperature.
- Add the
sdis_cameraAPI: it defines a pinhole camera into the scene.
- Add the
sdis_accum_bufferAPI: it is a pool of MC accumulators, i.e. a sum of MC weights and square weights.
- Add the
sdis_solve_camerafunction: it relies on a
sdis_accum_bufferto compute the radiative temperature that reaches each pixel of an image whose definition is defined by the caller. Note that actually this function uses the same underlying MC algorithm behind the
First version and implementation of the Stardis-Solver API.
- Support fluid/solid and solid/solid interfaces.
- Only conduction is currently fully supported: convection and radiative temperature are not computed yet. Fluid media can be added to the system but currently, Stardis-Solver assumes that their temperature are known.