Skip to content

Update: comments #1063

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
AlbertoCuadra merged 2 commits into from
Jul 4, 2025
Merged

Update: comments #1063

Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
a624cae
Update: comments
AlbertoCuadra Jul 4, 2025
ee0d179
Merge pull request #1062 from CombustionToolbox/shockTurbulence
AlbertoCuadra Jul 4, 2025
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
44 changes: 43 additions & 1 deletion +combustiontoolbox/+shockdetonation/@JumpConditionsSolver/JumpConditionsSolver.m
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -1,4 +1,46 @@
classdef JumpConditionsSolver < handle
% The :mat:func:`JumpConditionsSolver` class computes the Rankine–Hugoniot
% jump conditions and associated amplification parameters for gas mixtures
% across a shock wave. It supports calorically perfect, thermally perfect, and
% calorically imperfect models. In addition to post-shock thermodynamic properties,
% the class evaluates the dimensionless slopes of the Hugoniot curve—``Gammas2``
% (``\Gamma``), ``Gammas1`` (``\Gamma_\rho``), and ``Gammas3`` (``\Gamma_p``)—which
% quantify the local sensitivity of the post-shock density to changes in downstream
% pressure (``\Gamma``), and the influence of upstream density and pressure
% perturbations on the post-shock state (``\Gamma_\rho``, ``\Gamma_p``),
% as formalized in Ref. [1, 2].
%
% The :mat:func:`JumpConditionsSolver` object can be initialized as follows: ::
%
% solver = JumpConditionsSolver('equilibriumSolver', eqSolver, ...)
%
% Optional name-value arguments can be specified to customize tolerances,
% interpolation settings, plotting options, and solver behavior. Default
% objects for :mat:func:`EquilibriumSolver` and :mat:func:`ShockSolver`
% are created automatically if not provided.
%
% Flags:
% * ``FLAG_FAST``: Use previous state as initial guess for equilibrium solver
% * ``FLAG_INTERPOLATE``: Interpolate results on a reduced Mach number grid
% * ``FLAG_PLOT``: Plot results
% * ``FLAG_PAPER``: Apply the normalization and notation conventions described in Cuadra2025a
% * ``FLAG_RESULTS``: Print computed values to console
% * ``FLAG_TIME``: Print elapsed computation time
% * ``FLAG_REPORT``: Generate predefined plots for diagnostic analysis
%
% See also: :mat:func:`ShockSolver`, :mat:func:`EquilibriumSolver`,
% :mat:func:`solve`, :mat:func:`report`, :mat:func:`plot`, :mat:func:`plotGammas`
%
% References:
%
% [1] Cuadra, A., Di Renzo, M., Hoste, J. J. O., Williams, C. T., Vera, M., & Huete, C. (2025).
% Review of shock-turbulence interaction with a focus on hypersonic flow. Physics of Fluids, 37(4).
% DOI: 10.1063/5.0255816.
%
% [2] Cuadra, A., Williams, C. T., Di Renzo, M. & Huete, C. (2024). Compressibility
% and vibrational-excitation effects in hypersonic shock-turbulence interaction.
% Tech. Rep. Summer Program Proceedings, Center for Turbulence Research,
% Stanford University.

properties
equilibriumSolver % EquilibriumSolver object
Expand All @@ -12,7 +54,7 @@
FLAG_FAST = false % Flag to use fast equilibrium solver (guess from previous state)
FLAG_INTERPOLATE = false % Flag to interpolate results
FLAG_PLOT = false % Flag to plot results
FLAG_PAPER = false % Flag to plot results in paper format
FLAG_PAPER = false % Flag to apply the normalization and notation conventions described in Cuadra2025a
FLAG_RESULTS = true % Flag to print results
FLAG_TIME = true % Flag to print elapsed time
FLAG_REPORT = false % Flag to print predefined plots
Expand Down