From 6ee5bf26ab2e838a5d747cdfbaf5a4c56a1ee419 Mon Sep 17 00:00:00 2001 From: Litvinov Sergey Date: Wed, 22 Feb 2012 10:19:23 +0100 Subject: [PATCH 1/5] Add to gamma-deform-stretching --- .../gamma-deform-stretching.dat | 48 +++++++++++++++++++ 1 file changed, 48 insertions(+) diff --git a/dpd-2d-reverse-polymer/gamma-deform-stretching.dat b/dpd-2d-reverse-polymer/gamma-deform-stretching.dat index 9175264..76fd2b2 100644 --- a/dpd-2d-reverse-polymer/gamma-deform-stretching.dat +++ b/dpd-2d-reverse-polymer/gamma-deform-stretching.dat @@ -94,3 +94,51 @@ 1.000000 50 27.2181987272669 0.24916863105908 1.000000 10 28.0812904848463 0.249102378116986 1.000000 0 28.3331644245877 0.249363882654639 +0.010000 20 0.252445114727354 0.00249793714476833 +0.010000 30 0.216762526242511 0.00246536869450731 +0.010000 40 0.235780955563694 0.00249656400113797 +0.013594 20 0.379480160000615 0.00335585399237701 +0.013594 30 0.383217790907913 0.00334490576185196 +0.013594 40 0.345059981575826 0.00337204865754352 +0.018478 20 0.563036438787973 0.00461109154395204 +0.018478 30 0.508484713332337 0.00462886542472012 +0.018478 40 0.418173436787958 0.00462169867908556 +0.025119 20 0.78130908564268 0.00620651555031252 +0.025119 30 0.738401507663974 0.00615779748348985 +0.025119 40 0.66414524363643 0.00624617438385587 +0.034145 20 0.999793794808231 0.00840288895796307 +0.034145 30 0.915387362214617 0.00846051628837411 +0.034145 40 0.999891760235834 0.00848793687175844 +0.046416 20 1.39491962424244 0.0115302539745124 +0.046416 30 1.2910847030302 0.0114863569160486 +0.046416 40 1.3761933787877 0.0115625110777047 +0.063096 20 1.90584725454534 0.0157250360080374 +0.063096 30 1.77762768484835 0.0156622607095075 +0.063096 40 1.84617730909089 0.0156908505834349 +0.085770 20 2.50281922424203 0.0212951528429003 +0.085770 30 2.43678383636355 0.0213073836921052 +0.085770 40 2.43786583030401 0.0213906659168877 +0.116591 20 3.4043220424238 0.0290853978763829 +0.116591 30 3.29759886666603 0.0289518037485701 +0.116591 40 3.25997795757557 0.0290227531587111 +0.158489 20 4.56097427878803 0.0396112447056648 +0.158489 30 4.54161701818072 0.0394303437839726 +0.158489 40 4.45180503030057 0.0393937587269669 +0.215443 20 6.28499356969619 0.0537733681362065 +0.215443 30 6.05647797575714 0.0534896953753304 +0.215443 40 6.08368792727169 0.0536511100942023 +0.292864 20 8.3519229393931 0.0730334400593969 +0.292864 30 8.007070418181 0.0726043370782388 +0.292864 40 8.14150824848439 0.0727676751587069 +0.398107 20 11.1691059999985 0.0992705353381222 +0.398107 30 11.0259449696953 0.0991783336393523 +0.398107 40 11.0224135757561 0.0984505530554941 +0.541170 20 15.0592363030285 0.134854059233708 +0.541170 30 14.905864060604 0.134936561273912 +0.541170 40 15.1877385454526 0.133980636417747 +0.735642 20 20.4487926666674 0.183326617754556 +0.735642 30 20.2014998181803 0.183633600734454 +0.735642 40 20.1380015757613 0.183525613563737 +1.000000 20 27.7231625454526 0.249366564209995 +1.000000 30 27.3804496969672 0.249557238229797 +1.000000 40 27.3185571515108 0.249836510447018 From b8907d7415953392c8a504c7a47b448aafe9ccc1 Mon Sep 17 00:00:00 2001 From: Litvinov Sergey Date: Wed, 22 Feb 2012 10:20:27 +0100 Subject: [PATCH 2/5] Switching to long simulation --- .../deform-stretching/in.setbonds | 2 +- .../gamma-deform-stretching.dat | 32 +++++++++++++++++++ dpd-2d-reverse-polymer/genid_def.sh | 2 +- .../post-deform-stretching.gp | 8 +++-- 4 files changed, 39 insertions(+), 5 deletions(-) diff --git a/dpd-2d-reverse-polymer/deform-stretching/in.setbonds b/dpd-2d-reverse-polymer/deform-stretching/in.setbonds index 458ce80..e25b4df 100644 --- a/dpd-2d-reverse-polymer/deform-stretching/in.setbonds +++ b/dpd-2d-reverse-polymer/deform-stretching/in.setbonds @@ -3,7 +3,7 @@ special_bonds lj 1 1 1 # FENE parameters variable r0 equal 1.0 variable H equal 100 -variable INF equal 2.0*${xsize} +variable INF equal 10.0*${xsize} variable stfy equal 0.0 variable stfz equal 0.0 bond_coeff 1 ${H} ${r0} ${stfx} ${stfy} ${stfz} -${INF} ${INF} -${INF} ${INF} -${INF} ${INF} diff --git a/dpd-2d-reverse-polymer/gamma-deform-stretching.dat b/dpd-2d-reverse-polymer/gamma-deform-stretching.dat index 76fd2b2..b310844 100644 --- a/dpd-2d-reverse-polymer/gamma-deform-stretching.dat +++ b/dpd-2d-reverse-polymer/gamma-deform-stretching.dat @@ -142,3 +142,35 @@ 1.000000 20 27.7231625454526 0.249366564209995 1.000000 30 27.3804496969672 0.249557238229797 1.000000 40 27.3185571515108 0.249836510447018 +0.010000 1 0.406344920000067 0.00250364334683254 +0.010000 5 0.348317824848613 0.00254396846611246 +0.013594 1 0.540563473333353 0.003440688092407 +0.013594 5 0.423411637575899 0.00339065678012378 +0.018478 1 0.724837998787808 0.004486195883478 +0.018478 5 0.588915413333243 0.00458901161637929 +0.025119 1 0.986850108756034 0.00640982359357683 +0.025119 5 0.887733757720689 0.00617801947326481 +0.034145 1 1.27647630303296 0.00861419058464772 +0.034145 5 1.13772687488871 0.00835324736429371 +0.046416 1 1.68486682424231 0.0115090668576627 +0.046416 5 1.5061245575756 0.0115768442846654 +0.063096 1 2.19490626666754 0.0157299755934559 +0.063096 5 2.12293781818174 0.0158149865597325 +0.085770 1 2.8963003515168 0.0213621387420302 +0.085770 5 2.75717859999994 0.021452995700326 +0.116591 1 3.82192483030212 0.0293065176136166 +0.116591 5 3.69537964242384 0.0290710494562209 +0.158489 1 5.05933147272637 0.0393685568550798 +0.158489 5 4.96591072727077 0.0393513123852328 +0.215443 1 6.69511441211748 0.0536530379741903 +0.215443 5 6.58813273333332 0.0535666740545983 +0.292864 1 8.88676645454288 0.0730077426540564 +0.292864 5 8.8170731212105 0.0730300628742784 +0.398107 1 11.8113883030264 0.0994378762532137 +0.398107 5 11.7246616363655 0.0989616211100995 +0.541170 1 15.6879379394097 0.135207695210306 +0.541170 5 15.6891439394035 0.135161320981083 +0.735642 1 21.007940969693 0.183493879936374 +0.735642 5 20.9717623030262 0.183686588746455 +1.000000 1 28.1846391515133 0.250266467531927 +1.000000 5 28.1654220605909 0.248930021356109 diff --git a/dpd-2d-reverse-polymer/genid_def.sh b/dpd-2d-reverse-polymer/genid_def.sh index fd0f661..8640364 100755 --- a/dpd-2d-reverse-polymer/genid_def.sh +++ b/dpd-2d-reverse-polymer/genid_def.sh @@ -17,7 +17,7 @@ R_list="0.010000 0.541170 \ 0.735642 \ 1.000000" -stfx_list="20 30 40" +stfx_list="1 5" NS="~" diff --git a/dpd-2d-reverse-polymer/post-deform-stretching.gp b/dpd-2d-reverse-polymer/post-deform-stretching.gp index c3fc06b..2b97fcc 100644 --- a/dpd-2d-reverse-polymer/post-deform-stretching.gp +++ b/dpd-2d-reverse-polymer/post-deform-stretching.gp @@ -8,11 +8,13 @@ eta_solvent=9.0*rho plot [][18:50]\ "solvent.dat" u 1:($3/$1) w lp lw 3, \ " Date: Wed, 22 Feb 2012 13:39:09 +0100 Subject: [PATCH 3/5] Run long channel --- dpd-2d-le-fedosov/in.dpd | 1 + dpd-2d-reverse-polymer/genid_def.sh | 22 +++---------------- dpd-2d-reverse-polymer/in.generic | 4 ++-- dpd-2d-reverse-polymer/in.main | 2 +- dpd-2d-reverse-polymer/post-hist.sh | 1 + dpd-2d-reverse-polymer/run.sh | 9 +++----- dpd-2d-reverse-polymer/runone.sh | 9 ++++---- dpd-2d-reverse-polymer/units.mac | 10 ++++++--- .../wall-stretching/in.createforce | 2 +- .../wall-stretching/in.setbonds | 13 ++++++----- .../wall-stretching/in.setdomainsize | 2 +- .../wall-stretching/in.setdump | 19 +++------------- 12 files changed, 34 insertions(+), 60 deletions(-) diff --git a/dpd-2d-le-fedosov/in.dpd b/dpd-2d-le-fedosov/in.dpd index c432ec8..3866ad4 100644 --- a/dpd-2d-le-fedosov/in.dpd +++ b/dpd-2d-le-fedosov/in.dpd @@ -59,6 +59,7 @@ pair_coeff * * ${adpd} ${gamma} ${cutoff} ${m_power} ${cutoff_cons} neighbor 0.5 bin neigh_modify delay 0 every 1 fix idnve all nve + # {erate} value = R # R = engineering strain rate (1/time units) fix 2 all deform 1 xy erate ${R} remap v diff --git a/dpd-2d-reverse-polymer/genid_def.sh b/dpd-2d-reverse-polymer/genid_def.sh index 8640364..e839d72 100755 --- a/dpd-2d-reverse-polymer/genid_def.sh +++ b/dpd-2d-reverse-polymer/genid_def.sh @@ -1,29 +1,13 @@ #! /bin/bash -varlist="R stfx" -R_list="0.010000 - 0.013594 \ - 0.018478 \ - 0.025119 \ - 0.034145 \ - 0.046416 \ - 0.063096 \ - 0.085770 \ - 0.116591 \ - 0.158489 \ - 0.215443 \ - 0.292864 \ - 0.398107 \ - 0.541170 \ - 0.735642 \ - 1.000000" -stfx_list="1 5" +varlist="gx" +gx_list=0.25 NS="~" FS="~" CHECK_VALID_VAL=0 -prefix="Rfx150" +prefix="long" suffix="" diff --git a/dpd-2d-reverse-polymer/in.generic b/dpd-2d-reverse-polymer/in.generic index 68a419f..ee3cabd 100644 --- a/dpd-2d-reverse-polymer/in.generic +++ b/dpd-2d-reverse-polymer/in.generic @@ -4,8 +4,8 @@ #variable simulationtype string "deform-shear" #variable simulationtype string "reverse-poiseuille" #variable simulationtype string "free-polymers" -#variable simulationtype string "wall-stretching" -variable simulationtype string "deform-stretching" +variable simulationtype string "wall-stretching" +#variable simulationtype string "deform-stretching" #variable simulationtype string "deform-relaxation" #variable simulationtype string "wall-stretching" echo both diff --git a/dpd-2d-reverse-polymer/in.main b/dpd-2d-reverse-polymer/in.main index 1c594c7..9a3956c 100644 --- a/dpd-2d-reverse-polymer/in.main +++ b/dpd-2d-reverse-polymer/in.main @@ -29,7 +29,7 @@ variable adpd_sp equal 25.0 variable adpd_pp equal 25.0 variable adpd_sw equal 25.0 variable adpd_pw equal 25.0 -variable ntime equal 360000 +variable ntime equal 36000000 variable savetime equal 10000 #variable gx equal 0.75 # from 0.025 to 0.25 (p.6) diff --git a/dpd-2d-reverse-polymer/post-hist.sh b/dpd-2d-reverse-polymer/post-hist.sh index be5025e..10105b8 100755 --- a/dpd-2d-reverse-polymer/post-hist.sh +++ b/dpd-2d-reverse-polymer/post-hist.sh @@ -13,6 +13,7 @@ for x in $(seq 0.005 0.01 1.0); do echo ${x} ${m} ${vmax} ../scritps/lasvav.sh ${dname}/polymer.ndensyty.hist | awk -v c=${x} '$2==c{print $3, $5}' > ${dname}/polymer.x${x} + ../scritps/lasvav.sh ${dname}/sxy-bond.hist | awk -v c=${x} '$2==c{print $3, $5}' > ${dname}/bond.x${x} done > ${dname}/m.dat echo $* | xargs -n1 ./shift.sh diff --git a/dpd-2d-reverse-polymer/run.sh b/dpd-2d-reverse-polymer/run.sh index 0ea5796..92759eb 100755 --- a/dpd-2d-reverse-polymer/run.sh +++ b/dpd-2d-reverse-polymer/run.sh @@ -3,9 +3,6 @@ set -e set -u -for R in $(./genid.sh list=R); do -for stfx in $(./genid.sh list=stfx); do - echo ${R} - echo ${stfx} -done -done | ~/bin/parallel -N2 -j 8 --verbose ./runone.sh {1} {2} +for gx in $(./genid.sh list=gx); do + echo ${gx} +done | ~/bin/parallel -N1 -j 8 --verbose ./runone.sh {1} diff --git a/dpd-2d-reverse-polymer/runone.sh b/dpd-2d-reverse-polymer/runone.sh index 8a3d187..4443a37 100755 --- a/dpd-2d-reverse-polymer/runone.sh +++ b/dpd-2d-reverse-polymer/runone.sh @@ -3,8 +3,7 @@ set -e set -u -R=$1 -stfx=$2 +gx=$1 Nb=15 if [ -f ".lammps-configs" ]; then @@ -16,10 +15,10 @@ elif [ -f "${HOME}/.lammps-configs" ]; then fi -id=$(./genid.sh R=${R} stfx=${stfx}) +id=$(./genid.sh gx=${gx}) mkdir -p ${id} -vars="-var id ${id} -var ndim 2 -var R ${R} -var stfx ${stfx} -var dpdrandom ${RANDOM}" +vars="-var id ${id} -var ndim 2 -var gx ${gx} -var dpdrandom ${RANDOM}" ${lmp} ${vars} -in in.geninit ../scritps/addpolymer.sh \ @@ -31,4 +30,4 @@ ${lmp} ${vars} -in in.geninit Npoly=full \ addangle=0 -${mpirun} -np 1 ${lmp} ${vars} -in in.main +${mpirun} -np 8 ${lmp} ${vars} -in in.main diff --git a/dpd-2d-reverse-polymer/units.mac b/dpd-2d-reverse-polymer/units.mac index 343c67d..3da3db7 100644 --- a/dpd-2d-reverse-polymer/units.mac +++ b/dpd-2d-reverse-polymer/units.mac @@ -12,11 +12,15 @@ c_fed: c/ (L/T); /* body force */ f_fed: f/ (L/T^2); +/* domain size */ +Lx_fed: Lx / L; + /* set numerical values */ num_env: [dx=2.5e-5, kT=1.38e-14, rho=1e3, m=rho*dx^3]; /* here what you get in Fedosov units */ -float(ev(eta_fed, sol, num_env, eta=1e-4, eval, eval)); -float(ev(c_fed, sol, num_env, c=1e-4, eval, eval)); -float(ev(f_fed, sol, num_env, f=5.0, eval, eval)); +float(ev(eta_fed, sol, num_env, eta=4.5e-5, eval, eval)); +float(ev(c_fed, sol, num_env, c=3e-1, eval, eval)); +float(ev(f_fed, sol, num_env, f=6.12, eval, eval)); +float(ev(Lx_fed, sol, num_env, Lx=1.803e-3, eval, eval)); diff --git a/dpd-2d-reverse-polymer/wall-stretching/in.createforce b/dpd-2d-reverse-polymer/wall-stretching/in.createforce index 375189a..f3b3281 100644 --- a/dpd-2d-reverse-polymer/wall-stretching/in.createforce +++ b/dpd-2d-reverse-polymer/wall-stretching/in.createforce @@ -1,4 +1,4 @@ fix idnve all nve -variable bodyfx atom mass*${gx}*(x>0.75*${xsize}) +variable bodyfx atom mass*${gx} fix body_force flow addforce v_bodyfx 0.0 0.0 diff --git a/dpd-2d-reverse-polymer/wall-stretching/in.setbonds b/dpd-2d-reverse-polymer/wall-stretching/in.setbonds index 3481759..d275ea6 100644 --- a/dpd-2d-reverse-polymer/wall-stretching/in.setbonds +++ b/dpd-2d-reverse-polymer/wall-stretching/in.setbonds @@ -3,9 +3,10 @@ special_bonds lj 1 1 1 # FENE parameters variable r0 equal 1.0 variable H equal 100 -variable INF equal 2.0*${xsize} -variable fx equal 150 -variable fy equal 0.0 -variable fz equal 0.0 -variable xst equal 0.25*${xsize} -bond_coeff 1 ${H} ${r0} ${fx} ${fy} ${fz} 3.0 ${xst} -${INF} ${INF} -${INF} ${INF} +variable INF equal 10.0*${xsize} +variable stfx equal 20 +variable stfy equal 0.0 +variable stfz equal 0.0 +variable xst equal 0.5*${xsize} +variable buff equal 0.45*${xsize} +bond_coeff 1 ${H} ${r0} ${stfx} ${stfy} ${stfz} 3.0 ${xst} -${INF} ${INF} -${INF} ${INF} ${buff} diff --git a/dpd-2d-reverse-polymer/wall-stretching/in.setdomainsize b/dpd-2d-reverse-polymer/wall-stretching/in.setdomainsize index 338feec..19a2f10 100644 --- a/dpd-2d-reverse-polymer/wall-stretching/in.setdomainsize +++ b/dpd-2d-reverse-polymer/wall-stretching/in.setdomainsize @@ -1,5 +1,5 @@ # domain size -variable xsize equal 100 +variable xsize equal 600 # for wall add to 6 extra layers variable ysize equal 25+6 if "${ndim}==3" then & diff --git a/dpd-2d-reverse-polymer/wall-stretching/in.setdump b/dpd-2d-reverse-polymer/wall-stretching/in.setdump index 4c52f8d..0e4d5c2 100644 --- a/dpd-2d-reverse-polymer/wall-stretching/in.setdump +++ b/dpd-2d-reverse-polymer/wall-stretching/in.setdump @@ -18,19 +18,6 @@ fix av_nd_hist polymer ave/spatial 1 10000 10000 & x center 0.01 & y lower 0.01 density/number file ${id}/polymer.ndensyty.hist units reduced ave running -# variable maxi equal 10 -# label loopa -# variable i loop ${maxi} -# # average velocity -# variable x0_rgn equal ${i}*${xsize}/${maxi} -# variable x1_rgn equal (${i}+1)*${xsize}/${maxi} -# region rgn_${i} block ${x0_rgn} ${x1_rgn} INF INF INF INF -# fix av_vx_${i} flow ave/spatial 1 ${Nrepeat} ${Nfreq} & -# y center 0.01 vx file ${id}/vx.av.${i} units box region rgn_${i} -# if $i==${maxi} then "jump SELF break" -# next i -# jump SELF loopa -# label break - - - +fix av_sx_bond_hist flow ave/spatial 1 10000 10000 & + x center 0.01 & + y center 0.01 v_vstress_bond file ${id}/sxy-bond.hist units reduced ave running From bcbf36554decb1331259a6b4d35a33bc1b63a552 Mon Sep 17 00:00:00 2001 From: Litvinov Sergey Date: Wed, 22 Feb 2012 15:01:55 +0100 Subject: [PATCH 4/5] Fixing units.mac --- dpd-2d-reverse-polymer/units.mac | 11 +++++++++-- 1 file changed, 9 insertions(+), 2 deletions(-) diff --git a/dpd-2d-reverse-polymer/units.mac b/dpd-2d-reverse-polymer/units.mac index 3da3db7..3537095 100644 --- a/dpd-2d-reverse-polymer/units.mac +++ b/dpd-2d-reverse-polymer/units.mac @@ -6,7 +6,7 @@ eq[3]: kT = M*L^2/T^2; sol: solve([eq[1], eq[2], eq[3]], [L, M, T])[1]; /* dimensionless form */ -eta_fed: eta/ (L^2/T); +eta_fed: eta/ ( M/(L*T) ); /* sound speed */ c_fed: c/ (L/T); /* body force */ @@ -14,13 +14,20 @@ f_fed: f/ (L/T^2); /* domain size */ Lx_fed: Lx / L; +v_fed: v / (L/T); /* set numerical values */ num_env: [dx=2.5e-5, kT=1.38e-14, rho=1e3, m=rho*dx^3]; /* here what you get in Fedosov units */ -float(ev(eta_fed, sol, num_env, eta=4.5e-5, eval, eval)); +float(ev(eta_fed, sol, num_env, eta=1.50e-2, eval, eval)); float(ev(c_fed, sol, num_env, c=3e-1, eval, eval)); float(ev(f_fed, sol, num_env, f=6.12, eval, eval)); float(ev(Lx_fed, sol, num_env, Lx=1.803e-3, eval, eval)); +float(ev(v_fed, sol, num_env, v=0.041, eval, eval)); + +vmax(f, eta, rho, H):= f/(2*eta/rho) * H/2 * H/2; + +/* Fedosov */ +vmax(0.25, 42, 3, 25); \ No newline at end of file From 374f321f31c05c58e0fa56be621dab9c7b29d2a0 Mon Sep 17 00:00:00 2001 From: Litvinov Sergey Date: Wed, 22 Feb 2012 17:50:12 +0100 Subject: [PATCH 5/5] From kana --- dpd-2d-reverse-polymer/genid_def.sh | 2 +- dpd-2d-reverse-polymer/units.mac | 47 +++++++++++++++++------------ 2 files changed, 28 insertions(+), 21 deletions(-) diff --git a/dpd-2d-reverse-polymer/genid_def.sh b/dpd-2d-reverse-polymer/genid_def.sh index e839d72..01ef835 100755 --- a/dpd-2d-reverse-polymer/genid_def.sh +++ b/dpd-2d-reverse-polymer/genid_def.sh @@ -1,7 +1,7 @@ #! /bin/bash varlist="gx" -gx_list=0.25 +gx_list=0.125 NS="~" diff --git a/dpd-2d-reverse-polymer/units.mac b/dpd-2d-reverse-polymer/units.mac index 3537095..8f1dada 100644 --- a/dpd-2d-reverse-polymer/units.mac +++ b/dpd-2d-reverse-polymer/units.mac @@ -5,29 +5,36 @@ eq[2]: m = M; eq[3]: kT = M*L^2/T^2; sol: solve([eq[1], eq[2], eq[3]], [L, M, T])[1]; -/* dimensionless form */ -eta_fed: eta/ ( M/(L*T) ); -/* sound speed */ -c_fed: c/ (L/T); -/* body force */ -f_fed: f/ (L/T^2); - -/* domain size */ -Lx_fed: Lx / L; -v_fed: v / (L/T); - /* set numerical values */ num_env: [dx=2.5e-5, kT=1.38e-14, rho=1e3, m=rho*dx^3]; -/* here what you get in Fedosov units */ -float(ev(eta_fed, sol, num_env, eta=1.50e-2, eval, eval)); -float(ev(c_fed, sol, num_env, c=3e-1, eval, eval)); -float(ev(f_fed, sol, num_env, f=6.12, eval, eval)); -float(ev(Lx_fed, sol, num_env, Lx=1.803e-3, eval, eval)); -float(ev(v_fed, sol, num_env, v=0.041, eval, eval)); +fed_to_sdpd(units, val):= block([var, aux], + eqlocal: aux=var/units, + float(ev(rhs(solve(eqlocal, var)[1]), sol, num_env, aux=val, eval, eval)) + ); + +sdpd_to_fed(units, val):= block([var, aux], + eqlocal: aux=var*units, + float(ev(rhs(solve(eqlocal, var)[1]), sol, num_env, aux=val, eval, eval)) + ); + +/* sound speed */ +fed_to_sdpd(L/T, 10.0); +/* force */ +fed_to_sdpd(L/T^2, 0.25); +/* velocities */ +fed_to_sdpd(L/T, 1.3); +/* domain size */ +fed_to_sdpd(L, 50); +/* dynamic viscosity */ +fed_to_sdpd(M/(L*T), 42); +/* ks */ +fed_to_sdpd(L^2/T^2*M / L^2, 50); + +/* rmax */ +fed_to_sdpd(L, 1.0); -vmax(f, eta, rho, H):= f/(2*eta/rho) * H/2 * H/2; -/* Fedosov */ -vmax(0.25, 42, 3, 25); \ No newline at end of file +sdpd_to_fed(L, 0.00180281196288426); +sdpd_to_fed(L/T, 0.04);