OpenSn User Guide (Draft)

• 1. Overview
• 2. Installation
• 3. OpenSn Basics
  • 3.1. The Core Model
  • 3.2. The Typical Workflow
  • 3.3. What Usually Goes Where
  • 3.4. The Three Most Important Practical Rules
  • 3.5. Good Default Habits
  • 3.6. What To Read Next
• 4. Geometry and Mesh
  • 4.1. Overview
  • 4.2. MeshContinuum
  • 4.3. Mesh Generators
  • 4.4. Coordinate Systems
  • 4.5. Partitioning
  • 4.6. Replicated Meshes
  • 4.7. Block IDs
  • 4.8. Boundary Naming and Boundary IDs
  • 4.9. Logical Volumes
  • 4.10. Practical Mesh Workflows
  • 4.11. Inspecting and Exporting Meshes
  • 4.12. Practical Guidance
• 5. Materials and Cross-Section Data
  • 5.1. Overview
  • 5.2. Materials, Block IDs, and xs_map
  • 5.3. Creating a Simple One-Group Cross Section
  • 5.4. Inspecting a MultiGroupXS Object
  • 5.5. Loading OpenSn Cross-Section Files
  • 5.6. Loading OpenMC MGXS Files
  • 5.7. Custom Cross Sections
  • 5.8. Combining Cross Sections
  • 5.9. Scaling Cross Sections
  • 5.10. Choosing Between OpenSn and OpenMC Inputs
  • 5.11. Practical Example
  • 5.12. Cautions and Best Practices
• 6. Angular Quadratures
  • 6.1. Overview
  • 6.2. How to Choose a Quadrature
• 7. Product Quadratures
  • 7.1. Selecting Enough Angular Resolution
  • 7.2. Angle Aggregation Compatibility
  • 7.3. Examples
  • 7.4. Inspecting a Quadrature
  • 7.5. Recommendations
• 8. Groupsets
  • 8.1. Overview
  • 8.2. What a Groupset Does
  • 8.3. Required Core Fields
  • 8.4. Writing a Minimal Groupset
  • 8.5. Common Optional Groupset Controls
  • 8.6. Single Groupset vs Multiple Groupsets
  • 8.7. Groupsets and AGS
  • 8.8. Example Patterns
  • 8.9. Practical Guidance
• 9. Iterative Methods
  • 9.1. Overview
  • 9.2. Where the Iteration Controls Live
  • 9.3. Inner Groupset Methods
  • 9.4. Convergence Controls
  • 9.5. Across-Groupset Iteration
  • 9.6. DSA Options
  • 9.7. Solver-Specific Outer Iteration
  • 9.8. Choosing a Method
  • 9.9. Examples
  • 9.10. Recommendations
• 10. Boundary Conditions and Sources
  • 10.1. Overview
  • 10.2. Boundary Conditions
  • 10.3. Point Sources
  • 10.4. Volumetric Sources
  • 10.5. Attaching Sources at Problem Construction
  • 10.6. Replacing Sources After Construction
  • 10.7. Transient Source Behavior
  • 10.8. Boundary Names and Mesh Labels
  • 10.9. Best Practices
• 11. Discrete-Ordinates Problems
  • 11.1. Overview
  • 11.2. Constructor Summary
  • 11.3. Constructor Inputs
  • 11.4. Shared LBS Interface
  • 11.5. Problem Options
  • 11.6. sweep_type
  • 11.7. Problem Modes
  • 11.8. Field-Function Interface
  • 11.9. Angular-Flux Access
  • 11.10. Balance and Leakage
  • 11.11. Writing and Reading Transport State
  • 11.12. Updating the Problem In Place
  • 11.13. DiscreteOrdinatesCurvilinearProblem
  • 11.14. Typical Construction Patterns
  • 11.15. Practical Guidance
• 12. Solvers
• 13. Overview
• 14. Base Classes
• 15. LBS Problem
• 16. Discrete Ordinates Problem
  • 16.1. Time-Dependent and Steady-State Modes
  • 16.2. Angular-Flux-Specific Features
  • 16.3. DiscreteOrdinatesCurvilinearProblem
• 17. Steady-State Source Solver
  • 17.1. Constructor
  • 17.2. GPU support
  • 17.3. How It Operates
  • 17.4. Adjoint Steady-State Solves
• 18. Transient Solver
  • 18.1. GPU support
  • 18.2. Constructor
  • 18.3. Using Execute
  • 18.4. Using an Explicit Python Timestep Loop
  • 18.5. Balance Reporting
• 19. Power Iteration k-Eigen Solver
  • 19.1. Constructor
  • 19.2. GPU support
  • 19.3. How It Operates
• 20. Nonlinear k-Eigen Solver
  • 20.1. Constructor
  • 20.2. GPU support
• 21. Initialization Order and Common Patterns
• 22. Field Functions and Solver Output
• 23. Example Problems
  • 23.1. Common Input Structure
  • 23.2. Overview
  • 23.3. Common Building Blocks
  • 23.4. Example 1: Source-Driven Steady-State Problem
  • 23.5. Example 2: Multi-Material Fixed-Source Problem
  • 23.6. Example 3: Boundary-Driven Problem
  • 23.7. Example 4: Time-Dependent Source Problem
  • 23.8. Example 5: Explicit Python Time Loop
  • 23.9. Example 6: Power-Iteration k-Eigenvalue Problem
  • 23.10. Example 7: Nonlinear k-Eigenvalue Problem
  • 23.11. Example 8: Curvilinear Problem
  • 23.12. Example 9: Updating a Problem In Place
  • 23.13. Example 10: Field-Function Output
  • 23.14. Example 11: Balance and Leakage Checks
  • 23.15. Example 12: Reusing Flux Data
  • 23.16. Example 13: Writing and Reading Restart Dumps
  • 23.17. Choosing a Template
• 24. Post Processors
  • 24.1. Overview
  • 24.2. Updating Existing Field Functions
  • 24.3. Scalar Flux Field Functions
  • 24.4. Derived Field Functions
  • 24.5. Angular Flux Field Functions
  • 24.6. Exporting Field Functions
  • 24.7. Field-Function Interpolation
  • 24.8. Volume Postprocessor
  • 24.9. Other Useful Post-Processing Paths
  • 24.10. Transient Workflows
  • 24.11. Practical Guidance
• 25. Visualization
  • 25.1. Overview
  • 25.2. Visualizing the Mesh
  • 25.3. Visualizing Scalar Flux
  • 25.4. Visualizing Power
  • 25.5. Exporting Several Fields Together
  • 25.6. Transient Visualization
  • 25.7. Practical Advice
• 26. Running Problems
  • 26.1. Basic command-line usage
  • 26.2. Running under MPI
  • 26.3. Running with GPUs
  • 26.4. Typical execution pattern
  • 26.5. Inputs, logs, and outputs
  • 26.6. Restart and repeated solves
  • 26.7. Practical advice
• 27. Iterative Best Practices
  • 27.1. Overview
  • 27.2. Core Rules
  • 27.3. Upscatter Across Groupsets Must Be Converged
  • 27.4. Single Groupset First, Then Split Deliberately
  • 27.5. Choose the Inner Method Pragmatically
  • 27.6. Use DSA When the Physics Suggests It
  • 27.7. Watch Maximum Iteration Counts
  • 27.8. Be Careful with Nested Solver Expectations
  • 27.9. Transient Best Practices
  • 27.10. Eigenvalue Best Practices
  • 27.11. What to Check First When Convergence Looks Bad
  • 27.12. Recommended Starting Habits
• 28. Troubleshooting
  • 28.1. Start with the simplest question
  • 28.2. If the code fails immediately
  • 28.3. If the solution does not converge
  • 28.4. If the physics looks wrong
  • 28.5. If field functions are confusing
  • 28.6. If parallel behavior looks odd
  • 28.7. If imported cross sections behave unexpectedly
  • 28.8. General debugging workflow
  • 28.9. When to look elsewhere in the manual
• 29. FAQ
  • 29.1. Why is my problem split into a problem object and a solver object?
  • 29.2. Should I start with multiple groupsets?
  • 29.3. Why is my field function not updating?
  • 29.4. How do I get normalized power or XS-weighted outputs?
  • 29.5. When should I use Execute() and when should I use Advance()?
  • 29.6. What quadrature order should I use?
  • 29.7. Is single-angle aggregation always necessary when using unstructured meshes?
  • 29.8. What is the safest default sweep type?
  • 29.9. Do I need to export field functions during every run?