Repository navigation
Use 'thermal speed' with parallel temperature in bc constraint forcing - #482
Merged
Merged
Conversation
Where constraints are applied as part of the ion sheath-entrance boundary condition, a 'filter' is applied so that the corrections don't affect the shape function near v_parallel=0 or affect too strongly the large-v_parallel tail. Previously the small-v_parallel filter used a velocity scale based on sqrt(2*T_parallel/m_i) while the large-v_parallel filter used vth=sqrt(2*T/m_i), which is sqrt(2*T_parallel/3/m_i) in 1V simulations. It would be more consistent to use either T_parallel or T for both. As the modifications are to the shape of the correction in v_parallel, T_parallel makes sense, and makes the scales consistent between 1V and 2V simualtions, so this commit changes the large-v_parallel filter to use sqrt(2*T_parallel/m_i). Expected results for a few test results change, but only by very small amounts.
johnomotani
added
enhancement
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
Add this suggestion to a batch that can be applied as a single commit.This suggestion is invalid because no changes were made to the code.Suggestions cannot be applied while the pull request is closed.Suggestions cannot be applied while viewing a subset of changes.Only one suggestion per line can be applied in a batch.Add this suggestion to a batch that can be applied as a single commit.Applying suggestions on deleted lines is not supported.You must change the existing code in this line in order to create a valid suggestion.Outdated suggestions cannot be applied.This suggestion has been applied or marked resolved.Suggestions cannot be applied from pending reviews.Suggestions cannot be applied on multi-line comments.Suggestions cannot be applied while the pull request is queued to merge.Suggestion cannot be applied right now. Please check back later.
Where constraints are applied as part of the ion sheath-entrance boundary condition, a 'filter' is applied so that the corrections don't affect the shape function near v_parallel=0 or affect too strongly the large-v_parallel tail. Previously the small-v_parallel filter used a velocity scale based on sqrt(2*T_parallel/m_i) while the large-v_parallel filter used vth=sqrt(2*T/m_i), which is sqrt(2*T_parallel/3/m_i) in 1V simulations. It would be more consistent to use either T_parallel or T for both. As the modifications are to the shape of the correction in v_parallel, T_parallel makes sense, and makes the scales consistent between 1V and 2V simualtions, so this commit changes the large-v_parallel filter to use sqrt(2*T_parallel/m_i).
Expected results for a few test results change, but only by very small amounts.
Resolves #475.