Skip to content

Commit

Permalink
Fix extinction temperature definition in Python example
Browse files Browse the repository at this point in the history
Assign extinction temperature based on inlet temperatures, for consistency with the use of np.isclose() for the extinction criteria.
  • Loading branch information
d-e-t authored Jun 30, 2021
1 parent 74293ee commit 7d2e0e4
Showing 1 changed file with 2 additions and 1 deletion.
Original file line number Diff line number Diff line change
Expand Up @@ -53,7 +53,7 @@

# Define a limit for the maximum temperature below which the flame is
# considered as extinguished and the computation is aborted
temperature_limit_extinction = 500 # K
temperature_limit_extinction = max(f.oxidizer_inlet.T, f.fuel_inlet.T)

# Initialize and solve
print('Creating the initial solution')
Expand Down Expand Up @@ -114,6 +114,7 @@
strain_factor = alpha[-1] / alpha[n_last_burning]
# Create an initial guess based on the previous solution
# Update grid
# Note that grid scaling changes the diffusion flame width
f.flame.grid *= strain_factor ** exp_d_a
normalized_grid = f.grid / (f.grid[-1] - f.grid[0])
# Update mass fluxes
Expand Down

0 comments on commit 7d2e0e4

Please sign in to comment.