nicolas
/
gaspype
mirror of https://github.com/DLR-Institute-of-Future-Fuels/gaspype.git
1
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

added weighting coefficient for better stability

This commit is contained in:
Nicolas Kruse 2025-06-24 10:55:37 +02:00
parent 9cf901e039
commit df877dcf50
1 changed files with 7 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
src/gaspype/_solver.py
View File

@ -88,6 +88,8 @@ def equilibrium_eq(fs: fluid_system, element_composition: FloatArray, t: float,
# global count # global count
# count = 0 # count = 0
weighting = 100
def residuals(logn: FloatArray) -> tuple[FloatArray, FloatArray]: def residuals(logn: FloatArray) -> tuple[FloatArray, FloatArray]:
# global count # global count
# count += 1 # count += 1
@ -103,11 +105,13 @@ def equilibrium_eq(fs: fluid_system, element_composition: FloatArray, t: float,
j_eq = a - a_sum * n / n_sum j_eq = a - a_sum * n / n_sum
# Residuals from elemental balance: # Residuals from elemental balance:
el_error = np.dot(el_matrix, n) el_sum = np.dot(el_matrix, n)
resid_ab = np.log(el_error) - element_norm_log resid_ab = weighting * (np.log(el_sum) - element_norm_log)
# print(el_sum, element_norm)
# Jacobian # Jacobian
j_ab = el_matrix * n / el_error[:, np.newaxis] j_ab = weighting * el_matrix * n / el_sum[:, np.newaxis]
return (np.hstack([resid_eq, resid_ab]), np.concatenate([j_eq, j_ab], axis=0)) return (np.hstack([resid_eq, resid_ab]), np.concatenate([j_eq, j_ab], axis=0))