Contains functions reponsible for the pebble disk initialisation, evolution due to source terms and pebble accretion. More...

#include "fargo.h"

Include dependency graph for Pebbles.c:

Go to the source code of this file.

Macros
#define	TRAPEZMAX 35

#define	TRAPEZEPS 1.0e-7

Functions
void	InitPebbleArrays ()
	Initialise polar arrays associated with the pebble disk. More...

void	EquilPebbleDisk (PolarGrid Rho, PolarGrid Vrad, PolarGrid *Vtheta)
	Sets up a pebble disk in a coagulation-drift equilibrium. More...

void	RestartPebbleDisk (PolarGrid *Rho, int index)
	Reads arrays necessary to restart the pebble disk. More...

void	InitPebblesViaFlux (PolarGrid Rho, PolarGrid Vrad)
	Imposing the radial mass flux of pebbles in a coagulation-drift equilibrium, initialises their surface density, flow velocity and local Stokes numbers. More...

void	BckpFieldsForBC ()
	Backs up the initial state of the pebble disk to impose damping boundary conditions later. More...

void	EtaPressureSupport (PolarGrid *Vtheta)
	Calculates tha gas rotation parameter eta. More...

void	PebbleStokesNumbers (PolarGrid *Rho)
	Calculates the local Stokes number using the dominant pebble size, local gas density and parametric pebble material density. More...

real	Trapzd (int n, real a, real b, real H2)
	Trapezoidal rule to integrate the column mass of pebbles located within the accretion radius. More...

real	IntegrateColumnMass (real a, real b, real Hpeb)
	Primitive algorithm to find the mass in a vertical column of pebbles overlapping the accretion radius. More...

void	AccretePebblesOntoPlanets (PlanetarySystem sys, PolarGrid Rho, PolarGrid Energy, PolarGrid Vtheta, real dt)
	Finds the amount of pebbles to be transfered from the pebble disk onto the planets. More...

void	CorrectPebblesVtheta (real domega)
	Corrects the azimuthal flow velocity of pebbles to keep up with the frame rotation. More...

void	SourceTermsPebbles (PolarGrid Vrad, PolarGrid Vtheta, real dt)
	Calculates the source terms acting on pebbles. More...

void	SubStep1Pebbles (PolarGrid Vrad, PolarGrid Vtheta, real dt)
	Applies a semi-implicit method to evolve the pebbles dynamically. More...

void	ParticleDiffusion (PolarGrid *Rho)
	Applies the particle diffusion term acting on pebbles. More...

void	EvolvePebbleDisk (real dt)
	Calls the transport routines and applies the boundary conditions for pebbles. More...

void	SynchronizePebbleDisc ()
	Synchronises pebble fluid hydrodynamic quantities among the overlapping grid zones. More...

void	CriticalCharTime (PolarGrid Vrad, PolarGrid Vtheta)
	Restricts the time step using the CFL condition for the pebble fluid. More...

void	WritePebbles (int index)
	Outputs the pebble fluid arrays. More...

boolean	DetectCrashPebbles ()
	Safety check for negative pebble densities. More...

void	ParametricAccretion (PlanetarySystem *sys, real dt)
	Writes filtering factors if needed. More...

Variables
boolean	Restart

boolean	FastTransport

static PolarGrid *	PebbleDensTemp

static PolarGrid *	PebbleAccelrad

static PolarGrid *	PebbleAcceltheta

static PolarGrid *	PebbleSize

static PolarGrid *	EtaFaceCentered

static PolarGrid *	EtaCellCentered

static PolarGrid *	AccretedMass

static real	pebbulkdens

Detailed Description

Contains functions reponsible for the pebble disk initialisation, evolution due to source terms and pebble accretion.

Also controls several outputs.

Author: Ondřej Chrenko chren.nosp@m.ko@s.nosp@m.irrah.nosp@m..tro.nosp@m.ja.mf.nosp@m.f.cu.nosp@m.ni.cz

The pebble disk equilibrium model is inspired by Lambrechts & Johansen (2014). The evolution follows a standard two-fluid approximation with linear drag coupling and particle diffusion.

LICENSE

Definition in file Pebbles.c.

Macro Definition Documentation

#define TRAPEZEPS 1.0e-7

Definition at line 22 of file Pebbles.c.

Referenced by IntegrateColumnMass().

#define TRAPEZMAX 35

Definition at line 21 of file Pebbles.c.

Referenced by IntegrateColumnMass().

Function Documentation

void AccretePebblesOntoPlanets	(	PlanetarySystem *	sys,
		PolarGrid *	Rho,
		PolarGrid *	Energy,
		PolarGrid *	Vtheta,
		real	dt
	)

Finds the amount of pebbles to be transfered from the pebble disk onto the planets.

Definition at line 312 of file Pebbles.c.

References AccretHeating, CellAbscissa, CellOrdinate, DT, polargrid::Field, GetGlobalIFrac(), HEATINGDELAY, heatsrc, heatsrc_index, heatsrc_max, IntegrateColumnMass(), InvRmed, Max_or_active, MPI_Allreduce(), MPI_COMM_WORLD, MPI_DOUBLE, MPI_INT, mpi_make1Dprofile(), MPI_SUM, polargrid::Nrad, polargrid::Nsec, OmegaFrame, OmegaInv, PEBBLEALPHA, PebbleDens, PebbleVrad, PebbleVtheta, PhysicalTime, PI, PLANETARYDENSITY, RHO2CGS, Rinf, Rmed, Rsup, SoundSpeed, SQRT2PI_INV, SQRT_ADIABIND_INV, StokesNumber, Surf, vt1D, VXplanet, VYplanet, Xplanet, Yplanet, and Zero_or_active.

Referenced by AlgoGas().

Here is the call graph for this function:

Here is the caller graph for this function:

void BckpFieldsForBC ( )

Backs up the initial state of the pebble disk to impose damping boundary conditions later.

Definition at line 168 of file Pebbles.c.

References polargrid::Field, polargrid::Nrad, polargrid::Nsec, OmegaFrame, PebbleDens, PebbleVrad, PebbleVtheta, PebDensInit, PebVradInit, PebVthetaInit, and Rmed.

Referenced by RestartPebbleDisk().

Here is the caller graph for this function:

void CorrectPebblesVtheta ( real domega )

Corrects the azimuthal flow velocity of pebbles to keep up with the frame rotation.

Definition at line 516 of file Pebbles.c.

References CorrectVtheta(), and PebbleVtheta.

Referenced by AlgoGas().

Here is the call graph for this function:

Here is the caller graph for this function:

void CriticalCharTime	(	PolarGrid *	Vrad,
		PolarGrid *	Vtheta
	)

Restricts the time step using the CFL condition for the pebble fluid.

Definition at line 702 of file Pebbles.c.

References FastTransport, polargrid::Field, invdtpeb_sq, max2(), Max_or_active, polargrid::Nsec, One_or_active, PebbleVrad, PebbleVtheta, PI, Rinf, Rmed, and Rsup.

Referenced by AlgoGas().

Here is the call graph for this function:

Here is the caller graph for this function:

boolean DetectCrashPebbles ( )

Safety check for negative pebble densities.

Definition at line 753 of file Pebbles.c.

References DetectCrash(), and PebbleDens.

Referenced by AlgoGas().

Here is the call graph for this function:

Here is the caller graph for this function:

void EquilPebbleDisk	(	PolarGrid *	Rho,
		PolarGrid *	Vrad,
		PolarGrid *	Vtheta
	)

Sets up a pebble disk in a coagulation-drift equilibrium.

Definition at line 58 of file Pebbles.c.

References EtaPressureSupport(), and InitPebblesViaFlux().

Referenced by Initialization().

Here is the call graph for this function:

Here is the caller graph for this function:

void EtaPressureSupport ( PolarGrid * Vtheta )

Calculates tha gas rotation parameter eta.

Definition at line 195 of file Pebbles.c.

References polargrid::Field, polargrid::Nrad, OmegaFrame, and Rmed.

Referenced by EquilPebbleDisk().

Here is the caller graph for this function:

void EvolvePebbleDisk ( real dt )

Calls the transport routines and applies the boundary conditions for pebbles.

Definition at line 675 of file Pebbles.c.

References DampPebbles(), PebbleDens, PebbleVrad, PebbleVtheta, and TransportPebbles().

Referenced by AlgoGas().

Here is the call graph for this function:

Here is the caller graph for this function:

void InitPebbleArrays ( )

Initialise polar arrays associated with the pebble disk.

Definition at line 37 of file Pebbles.c.

References CreatePolarGrid(), DragForceRad, DragForceTheta, GasAccelrad, GasAcceltheta, NRAD, NSEC, PEBBLEBULKDENS, PebbleDens, PebbleVrad, PebbleVtheta, pebbulkdens, RHO2CGS, and StokesNumber.

Referenced by InitEuler().

Here is the call graph for this function:

Here is the caller graph for this function:

void InitPebblesViaFlux	(	PolarGrid *	Rho,
		PolarGrid *	Vrad
	)

Imposing the radial mass flux of pebbles in a coagulation-drift equilibrium, initialises their surface density, flow velocity and local Stokes numbers.

Definition at line 80 of file Pebbles.c.

References ADIABIND, CPU_Master, CPU_Number, polargrid::Field, FLUX2CU, Merge, OmegaFrame, OUTPUTDIR, PEBBLECOAGULATION, PebbleDens, PEBBLEFLUX, PebbleVrad, PebbleVtheta, pebbulkdens, PebDensInit, PebVradInit, PebVthetaInit, PI, Rmed, StokesNumber, and WriteDiskPolar().