Need to make sure simulation_state components can be properly modified by the solvers, such as the v_inner_solver in the future implementation

[DeerWhale]

Is your feature request related to a problem? Please describe.
In order to solve the v_inner_boundary iteratively, the solver should have the ability to update the properties in simulation_state (including geometry and other components that depends on geometry, such as radiation_field_state.) This is roughly done in PR #2504, but might be more efficient way. Need to double check if the parts in the model/parse_input.py are only setting up initial values (such t_rad and dilution factor initial values) and the properties that the solvers need to be update live outside the parser.

[DeerWhale]

Also check if shrinking the first cell is better or just use the original shell grid is better.