Test Problems: JeansInstabilityComoving

Caution

Please do not edit this file (page) manually since the workflow will overwrite your changes.

This file (page) is automatically generated by the workflow Update test problem wiki page using the script tool/wiki/sync_test_problem_pages.py.

The workflow is triggered by push changes to any of example/test_problem/*/*/README.md and tool/wiki/sync_test_problem_pages.py.

configure.py options

• Must enable
  • --model=ELBDM
  • --gravity
  • --comoving
• Must disable
  • --particle
  • --unsplit_gravity
• Available options
  • Miscellaneous Options

Default setup

• 32^3 uniform resolution

  • No grid refinement is adopted in this test

• DT__MAX_DELTA_A is set to an extremely large number so that the time-step will be controlled by the fluid solver

• The default setup in Input__Parameter and Input__TestProb gives an UNSTABLE solution, in which the solution grows exponentially. To test a STABLE solution, one can try the following parameters:

  • Input__Parameter

    Parameter name	Value
    BOX_SIZE	0.5
    A_INIT	1.0e-5
    OUTPUT_DT	2.0e-6

  • Input__TestProb

    Parameter name	Value
    Jeans_RealAmp0	4.0e-8

Caution

This setup requires further verification

Note

• Analytical solution reference: Woo, T. & Chiueh, T. 2009, ApJ, 697, 850

Note

The imaginary part I grows much faster than the real part R. Consequently, the accuracy of linear prediction will be deteriorated when the assumption (2R >> I^2) starts to break down, especially in the higher-resolution tests, and hence the 2nd-order accuracy may no longer hold.

Caution

2nd-order accuracy has NOT been verified for the UNSTABLE solution

Links

Main page of Test Problems