FluxTubes: Alternative approximation - #4798
Open
AHaumer wants to merge 6 commits into
Open
Conversation
AHaumer
requested review from
HansOlsson,
MartinOtter,
ThomasBoedrich,
casella,
christiankral and
maltelenz
August 6, 2026 13:47
Contributor
Author
|
One of the original authors, Thomas Boedrich, wrote to me: |
Contributor
Author
|
I have a set of functions to determine the approximation parameters from raw data. |
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
Add this suggestion to a batch that can be applied as a single commit.This suggestion is invalid because no changes were made to the code.Suggestions cannot be applied while the pull request is closed.Suggestions cannot be applied while viewing a subset of changes.Only one suggestion per line can be applied in a batch.Add this suggestion to a batch that can be applied as a single commit.Applying suggestions on deleted lines is not supported.You must change the existing code in this line in order to create a valid suggestion.Outdated suggestions cannot be applied.This suggestion has been applied or marked resolved.Suggestions cannot be applied from pending reviews.Suggestions cannot be applied on multi-line comments.Suggestions cannot be applied while the pull request is queued to merge.Suggestion cannot be applied right now. Please check back later.
When working with FluxTubes, I found it nearly impossible to implement new materials using the approximation formula of Roschke. Additionally, when the model drives the material very high into saturation ("overfluxed region"), the behaviour is not always as expected (e.g. polarization J doesnt show a horizontal asymptote but decreases with increasing field strength). Unfortunately the original authors are no longer available.
So I decided to implement an alternative approximation: Smoothing Splines with Exponential Extrapolation.
Everything is explained in the documentation, and there will be a contribution to the Modelica 2027 conference.
Additionally, VoestAlpine allowed me to use their measured data of fully processed electric sheet (as used in rotating electric machines). So I was able to add 10 material records. The new parameter records contain the used raw data, so the approximation is traceable. It should be easier for users to implement their own material.
Additonally, 4 new examples have been implemented to demonstrate and test the new approximation thoroughly.
You might call it an enhancement, but it solves the users' problem of implementing new material records.
The solution is 100% backwards-compatible, and has been tested with (in alphabetical order) Dymola, OpenModelica and WolframSystemModeler.