# Bursting Bubble Simulation - Default Parameters (Newtonian)
# Format: key=value (one per line, # for comments)
# This file provides a typical configuration for bursting bubble simulations

# ============================================================
# Case Identification
# ============================================================
# Case number for folder naming (4-digit: 1000-9999)
# Output will be created in simulationCases/<CaseNo>/
CaseNo=1000

# ============================================================
# Physical Parameters (Dimensionless Numbers)
# ============================================================
# Ohnesorge number: Oh = mu/sqrt(rho*sigma*R) (viscous / inertial-capillary)
# Controls viscous vs surface tension effects
# Typical values: 1e-3 (water-like) to 1e-1 (viscous)
Oh=1e-2

# Bond number: Bo = rho*g*R^2/sigma (gravity / surface tension)
# Controls gravitational effects on bubble shape
# Note: Currently only Bo=0.0010 geometry available (Bo0.0010.dat)
Bond=1e-3

# Gas/liquid Ohnesorge ratio: the gas-phase Ohnesorge number is
# Oha = OhRatio*Oh, setting the viscosity contrast across the interface
OhRatio=2e-2

# ============================================================
# Mesh Adaptation (space)
# ============================================================
# Maximum refinement level (2^MAXlevel cells at finest resolution)
# Level 10 -> 1024 cells, Level 11 -> 2048 cells, Level 12 -> 4096 cells
# Higher levels give better accuracy but slower computation
MAXlevel=12

# Far-field coarsening floor (the interface is always kept at MAXlevel;
# this only limits how coarse the quiescent bulk is allowed to become)
MINlevel=4

# Initial uniform grid level used before adaptation kicks in
init_grid_level=5

# Wavelet error tolerance on the VOF tracer f (interface location)
fErr=1e-3

# Wavelet error tolerance on the velocity field
VelErr=1e-3

# Wavelet error tolerance on interface curvature
KErr=1e-6

# ============================================================
# Time Resolution (adaptive)
# ============================================================
# dtmax is a CEILING on the timestep, not a fixed step: surface tension
# (tension.h) reduces the real step to the capillary-wave limit each
# iteration, so the effective timestep is set adaptively at run time.

# Advective CFL number
CFL=0.1

# Timestep ceiling; the capillary-wave limit sets the real adaptive step
dtmax=1e-2

# Poisson/viscous solver convergence tolerance
TOLERANCE=1e-4

# Snapshot/restart dump interval (capillary time units)
tsnap=1e-2

# ============================================================
# Time Control
# ============================================================
# Maximum simulation time (dimensionless, based on capillary time scale)
# t_cap = sqrt(rho * R^3 / sigma)
# Typical jet formation occurs at t ~ 0.5-2.0
tmax=1.5

# ============================================================
# Domain Configuration
# ============================================================
# zWall: distance from bubble south pole to bottom wall
# Controls domain size: Ldomain = min(zWall + 6.0, 16.0)
# Typical values:
#   zWall=0.05: Bubble very close to wall (current default)
#   zWall=0.025: Bubble very close to wall
#   zWall=2.0: Domain x in [-4, 4] (matches original setup)
#   zWall=4.0: Domain x in [-6, 4], wall at x=-6
zWall=0.05