 a first numerical experiment we have performed the robust stability test as described in [@Alcubierre03d]. For this test one takes Minkowski initial data and adds random noise with a small amplitude to all dynamical variables. For the evolution we have used both harmonic slicing with $f=1$ and standard 1+log slicing with $f=2/\alpha$.

For the simulations discussed below we used 1,000 grid points, a grid spacing of $\Delta r = 0.001$ and a time step $\Delta t = \Delta r/2$. Furthermore, we have demanded periodic boundary conditions and have set all $1/r$ and $1/r^2$ terms to zero by hand (i.e. we are performing the run “at infinity”). Setting these terms to zero should have no important effect on the presence of blow-ups since one can readily verify that such terms are only linear in the eigenfields. Moreover, one could easily regularize the equations at the origin by following the procedure described in [@Alcubierre04a], but doing so would only obscure the study of blow-ups by mixing them with purely geometric effects.

We start the simulations with Minkowski initial data plus random noise of order $10^{-6}$ on all evolution variables. At latter times we compute the error $\delta$ as the average absolute value over the grid of the quantity $(\sum_{i=1}^3 |u_i - 1| + \sum_{i=1}^5
|v_i|)/8$.

We first performed runs for the ADM system without adding constraints, and observed the well known growth in the average error caused by the fact that ADM is only weakly hyperbolic. Next, we implemented the $\mu$-family given by (\[eq:hKBmu\]) and (\[eq:sourcehK\]), which according to the source criteria is shock-avoiding. Figure \[fig:mu\] shows the behavior of the average error for this family. In the top panel we plot the growth of the error for the case of harmonic slicing as a function of time (measured in units of the light-crossing time of our computational grid), for several different values of $|\mu|$. We see that for ADM (which is not a member of the $\mu$-family), and for the cases with $\mu=0$ and $|\mu| = 1$, the error grows linearly with time, while for $|\mu| = 1
a first numerical experiment we have performed the robust constancy trial as described in   [ @Alcubierre03d ]. For this test one lead Minkowski initial datum and adds random noise with a small amplitude to all dynamic variable. For the evolution we have use both harmonic slicing with $ f=1 $ and standard 1+log slice with $ f=2/\alpha$. 

 For the simulations discussed below we use 1,000 power system points, a grid space of $ \Delta r = 0.001 $ and a meter measure $ \Delta t = \Delta r/2$. Furthermore, we have demanded periodic boundary conditions and have set all $ 1 / r$ and $ 1 / r^2 $ term to zero by hand (i.e. we are performing the run “ at infinity ”). set these terms to zero should have no important effect on the presence of gust - ups since one can readily verify that such terms are only linear in the eigenfields. Moreover, one could easily regularize the equations at the origin by following the procedure described in   [ @Alcubierre04a ], but doing so would only obscure the survey of blow - ups by mixing them with strictly geometric effect. 

 We begin the simulations with Minkowski initial data plus random randomness of order $ 10^{-6}$ on all development variables. At latter times we compute the mistake $ \delta$ as the average absolute value over the grid of the quantity $ (\sum_{i=1}^3 |u_i - 1| + \sum_{i=1}^5 
 |v_i|)/8$. 

 We first performed runs for the ADM system without adding constraints, and observed the well known increase in the average mistake caused by the fact that ADM is only weakly hyperbolic. Next, we implemented the $ \mu$-family given by (\[eq: hKBmu\ ]) and (\[eq: sourcehK\ ]), which harmonize to the source criteria is shock - avoiding. name   \[fig: mu\ ] shows the behavior of the average error for this kin. In the top panel we plot the increase of the error for the event of harmonic slice as a function of time (measure in units of the light - cross time of our computational power system), for several different values of $ |\mu|$. We go steady that for ADM (which is not a member of the $ \mu$-family), and for the cases with $ \mu=0 $ and $ |\mu| = 1 $, the error grows linearly with time, while for $ |\mu| = 1
 a virst numerical experimekt we have perfotmwd the robusf stabilkty test as described in [@Alcuuieree03d]. Fir this test one takes Minkowskp initial date and adds random noise with a dmaln amplitude to sll dynamiwal variables. Xof che evolution we have used both harmjnic sloclng with $f=1$ and stamqard 1+log slicing with $f=2/\alpha$.

For the sijulatiois discussed bekow we used 1,000 grid points, a grif spacing of $\Delta r = 0.001$ and a jjme wtep $\Delta t = \Delta r/2$. Furthermore, se have demanded periodic boundxry cpnditions qne hwee set all $1/c$ and $1/g^2$ terms to zevp by hdnd (i.e. ee are performlng tie ryn “at infinity”). Settinj these terms to zerj should vare no important effecr in thg prevencd of blkw-nps since onx can readimy verify tyat such terms are pnjj linear in tge eigqnsields. Moreover, one could easily regulagize the equations at the oeigin by following thg procedurq described in [@Alcubierre04a], but doing so would only obscnrd tkc rrufy of blow-ups by mixing them with purely geomqfroc effects.

We stavt the simulations wltn Minkowski inktial bztz plus random noisf of orqer $10^{-6}$ in all evjlutoon variables. At latter timws we computv thw error $\delta$ as tke average aysolutg valur over the grid of the duanfity $(\sum_{i=1}^3 |u_l - 1| + \sum_{i=1}^5
|b_k|)/8$.

We first perforoed rgns for ufe ADM system witrout addiig coustraintr, anc obsewved the wfll kkmwn growth in the wveraye ersor caused by the fact that ADM is only wxekly hyperboloc. Neqt, we impjemenbed the $\mu$-familr given by (\[eq:hLBmu\]) anb (\[eq:sojrcehK\]), whibh accordmng to the sjurce criterid is shock-avomding. Figtre \[fug:mu\] shows gfe behavior of the avercye error dor this family. In thg fop panel we plit uhw growth of thr efrow voc the wase of harmmnic slkving xs a functiiv of time (measured in unhts kf the light-crossimg time of our com[utational grod), for several difverenv valuxs of $|\ku|$. Re see that for ADM (which is nkt a membfr jf the $\mu$-famyly), qnd for the eases with $\mu=0$ and $|\mu| = 1$, the error grows lmnearly with time, whilg for $|\mu| = 1
a first numerical experiment we have performed stability as described [@Alcubierre03d]. For this data adds random noise a small amplitude all dynamical variables. For the evolution have used both harmonic slicing with $f=1$ and standard 1+log slicing with $f=2/\alpha$. the simulations discussed below we used 1,000 grid points, a grid spacing of r 0.001$ a step $\Delta t = \Delta r/2$. Furthermore, we have demanded periodic boundary conditions and have set $1/r$ and $1/r^2$ terms to zero by hand we are performing the “at infinity”). Setting these terms zero have no effect the of blow-ups since can readily verify that such terms are only linear in the eigenfields. Moreover, one could easily regularize equations at by following procedure in but doing so obscure the study of blow-ups by purely geometric effects. We start the simulations with initial data random noise of order $10^{-6}$ on evolution variables. At latter times we compute the $\delta$ as the average absolute value over the grid of the quantity $(\sum_{i=1}^3 |u_i - \sum_{i=1}^5 |v_i|)/8$. We first runs for the system adding and the well growth in the average error caused by the fact that ADM only weakly hyperbolic. Next, we implemented the $\mu$-family given by (\[eq:sourcehK\]), according to the criteria is shock-avoiding. Figure shows behavior of the average this In we the of the error for case of harmonic slicing as function of time (measured time of our computational grid), for several different of $|\mu|$. We see that for ADM is not a member of the $\mu$-family), and for the cases with and $|\mu| the error grows linearly with time, while for = 1
 a first numerical experiment We have perfOrmed The RobUsT staBiliTy test as descriBEd in [@alcubierre03d]. For this test One taKeS minkOWsKi iniTial datA AnD ADds RaNdOm nOiSE wIth a sMalL amplitUde to all dyNamIcAl variables. FOR tHe evolutioN we Have used both HarMonic sLiCinG With $f=1$ And StandArd 1+log SLicing With $f=2/\alphA$.

FOR the siMUlationS DIsCussEd below we used 1,000 grid POiNTs, a grid spacing Of $\DeltA r = 0.001$ ANd A TIme SteP $\Delta t = \DelTa R/2$. FurtHErmore, wE HaVE DEmaNDed periodic boUndary condiTIonS and haVe Set ALl $1/r$ and $1/R^2$ termS tO ZerO by hand (i.e. we Are pErforming The run “AT infiniTY”). SettinG these TerMs tO zerO ShOuLd hAvE No iMPoRtaNT efFect on thE pReSence Of blOW-UPS sinCe oNe caN readIly verify that SucH terMS arE only LineaR in tHe EigenFields. moreoVeR, one could easily ReguLarize the EquAtIonS aT the oRIgin by FolLowIng the pRocedurE DesCrIBED iN [@Alcubierre04a], but doiNg SO WoUld only oBscure THe StUDy of blow-UpS by MixiNG Them wIth pUReLy geometRic effECtS.

WE start tHe SimulaTiOns WitH MinkOWski InitiaL data pluS randOM noise of order $10^{-6}$ oN All evolution vARiABLeS. at laTteR times we comPute THe erRor $\dELtA$ as THe aveRage aBsOLuTE value over the grid of ThE quantIty $(\suM_{i=1}^3 |u_i - 1| + \sum_{i=1}^5
|v_i|)/8$.

We fIrst perforMED Runs for tHe ADm SySTem without addiNg conStraints, anD Observed The weLl known gRowth in thE AVerage erRor CauSed By tHE FaCt that ADM is onLY WeakLy HyperboLic. next, we iMplEmeNteD thE $\mU$-family giVen by (\[eq:hkBMu\]) AnD (\[eQ:soUrcehk\]), Which accOrDinG tO thE sourCE criteRia is ShocK-aVoIDinG. Figure \[FIg:MU\] ShowS tHe BehaVioR oF the aVeraGE erRor for tHis family. in tHE top PaNeL we plot The growth of thE eRror for the CaSe oF harmoNIC slicing As a function of time (measurED in unitS of The liGht-cRossing tiMe oF our coMpuTAtionaL grid), fOr sevErAl dIFFerenT VAlUes Of $|\Mu|$. We see thaT FOr AdM (whiCh Is noT a membeR of the $\mu$-family), and fOR thE cases with $\mu=0$ aNd $|\mU| = 1$, the ERRoR grOWs LIneArLY wiTH Time, while for $|\mu| = 1
 a first numerical experim ent we hav e per for med  t he r obus t stability te s t as  described in [@Alcubi erre0 3d ] . Fo r  t his t est one ta k e s M in ko wsk ii ni tialdat a and a dds random  no is e with a sma l lamplitudetoall dynamica l v ariabl es . F o r the  ev oluti on weh ave us ed both h ar m onic s l icing w i t h$f=1 $ and standard 1+ l og slicing with $ f=2/\a lp h a$ . 
Fo r t he simulat io ns di s cussedb el o w weu sed 1,000 gri d points, a gri d spac in g o f  $\Del ta r=0 .00 1$ and a ti me s tep $\Del ta t = \Deltar /2$. Fu rtherm ore , w e ha v ede man de d  pe r io dic bou ndary co nd it ionsandh a v e  set  al l $1 /r$ a nd $1/r^2$ te rms  toz ero  by h and ( i.e.  w e are  perfo rming  t he run “at infi nity ”). Setti ngth ese  t ermst o zero  sh oul d haveno impo r tan te f f ec t on the presenceof b lo w-ups si nce on e  c an readilyve rif y th a t  such  ter m sare only  linea r  i nthe eig en fields .Mor eov er, o n e co uld ea sily reg ulari z e the equation s  at the origi n  b y fo l lowi ngthe procedu re d e scri bedi n[@A l cubie rre04 a] ,  b u t doing so would on ly  obscu re th e study of bl ow-ups bym i x ing them  wit h  p u rely geometric  effe cts.

We s t art thesimul ations w ith Minko w s ki initi aldat a p lus r an dom noise ofo r der$1 0^{-6}$  on  all ev olu tio n v ari ab les. At l atter ti me swe  c omp ute t h e error$\ del ta $ a s the averag e abs olut eva l ueover th e  g r i d of  t he  qua nti ty  $(\s um_{ i =1} ^3 |u_i  - 1| + \ sum _ {i=1 }^ 5|v_i|)/ 8$.

We first  p erformed r un s f or the A DM syste m without adding constr a ints, a ndobser vedthe wellkno wn gro wth in the  avera ge er ro r c a u sed b y th e f ac t that ADM i s o nly w ea klyhyperbo lic. Next, we impl e men ted the $\mu$ -fa mily g iv enb y( \[e q: h KBm u \ ]) and (\[eq:so urcehK\]),  w h ic h accordin g  to  t he sour ce crit eriai s shock -avoiding . Figure\[ fig: m u \]shows thebehavior  of the a v erage er ror f orthis f am ily . Inthe to p  pa nel w e plot  t he gro wth o fthe erro r for the case of harmo nic sl icing  as  a functi ono f t ime (meas ured  in unitsofthe  ligh t-c r ossin g ti m eofo ur co mput a tional gr i d) , f o r  s everal diff e r e ntvalue s o f  $|\mu |$.We see that for A D M (which is no t am e mbe r o f  the  $ \mu$-family),and  f o r  the cas es  with $\mu= 0$ and $ |\ m u| =1$, th e erro r grows l in e arly w ithtim e, whilefor  $ | \mu| =1
 a_first numerical_experiment we have performed_the robust_stability_test as_described_in [@Alcubierre03d]. For this_test one takes_Minkowski initial data and_adds random noise_with_a small amplitude to all dynamical variables. For the evolution we have used both_harmonic_slicing with_$f=1$_and_standard 1+log slicing with $f=2/\alpha$.

For_the simulations discussed below we_used 1,000_grid points, a grid spacing of $\Delta r_=_0.001$ and a_time step $\Delta t = \Delta r/2$. Furthermore, we_have demanded periodic boundary conditions and_have set all_$1/r$_and_$1/r^2$ terms to zero_by hand (i.e. we are performing_the run “at infinity”). Setting these_terms to zero should have no important_effect on the presence of blow-ups_since one can readily verify_that such_terms are only linear in_the eigenfields. Moreover,_one could_easily regularize the_equations at the origin by following_the procedure described_in [@Alcubierre04a], but doing so would only_obscure_the study of_blow-ups_by_mixing them_with purely geometric_effects.

We_start the_simulations_with Minkowski initial data plus random_noise_of order $10^{-6}$ on all evolution variables._At latter times we_compute_the error $\delta$ as_the average absolute value over_the grid of the quantity $(\sum_{i=1}^3_|u_i -_1| +_\sum_{i=1}^5
|v_i|)/8$.

We first performed runs for the ADM system without adding constraints,_and observed the well known growth_in the average error_caused by_the_fact that ADM_is_only weakly_hyperbolic. Next, we implemented the $\mu$-family given_by (\[eq:hKBmu\])_and (\[eq:sourcehK\]), which according to the_source criteria is shock-avoiding._Figure \[fig:mu\]_shows the behavior of the average_error for this family. In the_top panel we plot the_growth_of_the error for the case_of harmonic slicing as a function_of time (measured_in units of the light-crossing time of_our_computational grid), for several different values_of_$|\mu|$. We see that for ADM_(which_is_not a member of the_$\mu$-family), and for the cases with_$\mu=0$ and $|\mu| = 1$, the error grows linearly_with time, while_for $|\mu| = 1