Version history

ember roughly follows semantic versioning starting at 0.1.0. However, until 1.0.0, minor releases may make breaking changes to the public API without a deprecation period.

0.4.4 (2026-09-30)

  • Add ember.average.mass_band, the part of a structured 2D cut carrying a given band of mass fraction, such as a midspan streamsurface

  • Fix ember.cut.interpolate_to_structured putting nodes in empty space when cutting through a solid

  • Add Reichardt wall function

  • Fix bug with halved interior viscous stresses and heat fluxes

  • Apply the wall-law stress on the wall face itself, not on the face above the wall cell, and stop masking the wall cell’s viscous force

  • Add a prescribed transition property Block.fac_lam, to scale interior turbulent viscosity down and wall function toward laminar

  • Add Solver.sf2_T, a separate weight on the temperature limb of the adaptive smoothing sensor, defaulting to sf2

0.4.3 (2026-09-16)

  • Match periodic patches on a difference in theta, instead of wrapping each face separately, relax tolerance slightly

  • Divergence check now raises on negative pressure

  • Add endwall bands to mixing plane exchange, reducing the mean relaxation approaching the hub and casing by a specified factor, and turning off the harmonic relaxation entirely

0.4.2 (2026-09-12)

  • Add Jameson-style adaptive smoothing

  • Allow cusp patches over partial j ranges

  • Add affordances for faces collapsing to zero area: wall function clipping and NotWallPatch to exempt collapsed nodes from the wall distance calculation

  • Add a soft-start reconfiguration method BaseSolver.soft()

0.4.1 (2026-09-10)

  • Promote cut.signed_distance to public API

  • Made performance improvements to cutting functions

0.4.0 (2026-09-01)

  • Multigrid now applies a uniform coarse correction to the fine nodes, instead of linearly interpolating through the coarse cell centers. This is more stable on clustered grids.

  • The equation-of-state-agnostic fluid base class is now public as Fluid, renamed from _Fluid, and carries the shared from_dict/to_dict and change_* contract.

  • Trimmed the realgas_fit public surface to fit, sample_coolprop, FitResult and FitInfo; legfit2d is now private, and the pure Legendre/coordinate helpers move to the ember._realgas_poly leaf module.

0.3.0 (2026-08-31)

  • Selective frequency damping is available again.

  • Add OutletPatch.P_throttle, the static pressure level a throttled outlet has been adjusted to. It is equal to the initially prescribed pressure plus the mass flow PID correction. Also for consistency renamed the throttle pressure correction from P_throttle to dP_throttle in ConvergenceHistory.

  • Delete Block.xlen_sq_nd and Block.conserved_cell_nd, calculate them inline in the Fortran kernels instead.

  • Add RealFluid real-gas equation of state and fitting module, following Wheeler (2024) polynomial entropy fit.

  • Make Block.cp_nd a nodal array rather than a scalar, pass fields of viscous and thermal transport properties to the source term kernels.

  • Add to_dict and from_dict to fluid classes.

  • Add Fluid.get_kappa thermal conductivity to the fluid API.

  • Add change_visc() to fluid clasess, returning a new instance with viscosity scaled up or down. Perfect gases simply alter their scalar viscosity; real gases carry a scale_visc factor that multiplies the true fitted viscosity.

  • Exchange viscous/thermal halos via six face arrays instead of a full volume array.

  • Run the change limiter (Solver.dampin) before implicit residual smoothing (Solver.sf_resid) rather than after, fused into the kernel already traversing the volume. IRS is linear and the limiter is nonlinear in a global block mean, so the two orderings do not commute: at the default sf_resid=1.0, dampin=25 the residual differs by ~19% of the field scale, growing with sf_resid. The converged answer is unchanged, only the path to it, and the fusion is worth -6% serial and -11% at 100-rank saturation.

  • Block.scratch is now the sole scratch array, and all transient workspace for the Fortan kernels is carved out of it.

  • Various performance improvements to reduce memory footprint.

0.2.0 (2026-08-18)

  • Renamed ember.set_iter to ember.set_iterative.

  • Removed BlockRestart and its helpers.

  • Removed ember.geometry; its functions are now private helpers of their callers.

  • Made ember.struct and the Plot3D I/O module private.

  • Removed Grid.apply_rotation. It assumed a mesh topology — that j=-1 is the casing — of any grid handed to it, which is not ember’s to assume, and its tip_gap and shroud row types were in any case indistinguishable: every wall already takes its block’s Omega, so tip_gap overrode five faces with the value they already had and left the casing turning with the blade. Set Block.Omega for the frame, and add a RotatingPatch only for a wall that is not at that speed.

  • Pruned dead code from ember.util and made its single-caller helpers private.

  • Replaced the reflecting mixing plane with the non-reflecting one.

  • Solver is now a frozen immutable class, derive a variant with dataclasses.replace instead of assigning.

  • Changed Solver defaults: cfl 0.4 to 5.0, dampin None to 25, sf_resid 0.0 to 1.0, and n_levels 0 to 3.

  • Block.flat is now a property, and flattens in Fortran order.

  • Split interp_from into Grid.interp_from_arrays and Grid.interp_from_grid, and transfer primitives rather than conserved variables to prevent datum offset errors.

  • Generalised non-reflecting patches to any surface of revolution.

  • Added OutletPatch.Kp and OutletPatch.Ki, exposing the throttle gains set by set_throttle.

  • Added InletPatch.Po, InletPatch.To, InletPatch.Alpha, InletPatch.Beta and OutletPatch.P, reading back the currently prescribed boundary state in physical units.

  • Added wall_yplus for y+ post-processing.

  • Added Block.memory_usage.

  • Added Block.freeze(), returning an immutable copy.

  • Added block_util.repeat_pitchwise, viewing a periodic passage as a cascade. A periodic solver computes one passage, which is rarely all anyone wants to look at: the passage-to-passage picture is what shows a wake meeting the next leading edge, or a shock crossing a pitch. Returns one copy per passage rather than a concatenated block, and drops patches, a rotated copy being a view of the flow rather than a member of a connected grid.

  • Added util.unwrap_meridional, giving the conformal coordinate m' = int dm/r that a blade-to-blade view is drawn on.

  • Fixed mix_out returning slightly different states for a cut and its k-reversed twin. The Newton loop now iterates to its fixed point instead of stopping at the first iterate inside atol, cutting that difference from ~3e-4 to ~2e-6.

  • Added precompiled wheels for macOS (arm64) and Windows (AMD64), alongside the existing Linux (x86_64) wheels.

  • Fixed mass_average raising on a partial reduction. The averaging already honoured the axes argument, but the zero-flux guard beneath it assumed a scalar, so averaging over the pitch alone and keeping the span raised about an ambiguous truth value instead of returning a profile.

  • Fixed matvec broadcast dispatch order for (ni,1,1) matrices.

  • yaml_util now uses libyaml’s CSafeLoader/CSafeDumper where available for speed.

  • Fixed YAML being read and written in the locale’s encoding rather than UTF-8, so a file carrying a degree sign or a non-ASCII name round-tripped losslessly only by coincidence of environment.

  • Fixed YAML writing a Path raising RepresenterError on Windows.

  • Fixed the sign of the cusp seam correction.

  • Restored the smooth-then-damp order in update_residual.

  • Fixed the mixing-length viscosity returning -inf on flows with no resolved shear when built with gfortran 13.

0.1.0 (2026-08-06)

  • First public release on the Python Package Index.