Book a Demo!
CoCalc Logo Icon
StoreFeaturesDocsShareSupportNewsAboutPoliciesSign UpSign In
Download
602 views
unlisted
ubuntu2404
a

7"	i�M�@s�e��dd�ddlZddlTddlZddlmZddlZ	dZ
dZdZej
ded	�ej
d
e
d�ej
d
ed�ej
d
e
d�ej
de
d�ej
de
d�ej
d
dd�ej
ddd�dejd<e�d�Ze�d�Ze�d�Ze�d�Ze	�d�Ze	jddd�Ze	jddd�Ze	�d �Zd!d"�Zd#ed#d$ed%d&ed%ee�ee�ed'ed(d%ed%d%Zed#ee�eeee�d(d)ed(d%ed%d%Zed%ed%ee�ee�ed'ed(d%ed%d%ZeZed%ZeZ d&ed#d*ed%d%ed%d(ee�ed)ed(d%ed%d%Z!eed&eed#d(ee�d)d(d%ed%d%Z"e"ed%Z#e"ed%Z$ed%d+ed,d-edd.ed/d0ed1d2ed#d3ed%d1d4ed%ded/d5ed1d6ed#d7ed%d&ed/ee�d#d8ed9ed,d:edded/d;ed1d<ed#d=ed%d>ed#ee�d&ed?ed5d@ed,dAeddBed/dCed1dDed#dEed%dFee�d&dGed7dHed5dIed,dJeddKed/dLed1dMed#dNed%dOee�d%dPed%dQd%ed%d&dRed%d&ed%eee�d%dRed&d*ed&ed#ee�d&ed%d(ee�d&edQZ%ed%dSed8ded$dTed*dUed&d,ed8ed/ed%d(ed%ee�d&d&ed#d8ed1ded#dVed%d5eee�d%d4eddWed/dXed1dYed#dZed%d,ee�d)ed%ed&d%dRed%d&ed%eee�d%dRed&d*ed&ed#ee�dQZ&ed%d[edd\ed/d]ed1d^ed#d_ed%d`d4ed%d1ed1ed#d#ed%d(ed/ee�d#d8ed,eddaed/dbed1dced#dded%d4ed#ee�d&deedfdgedhdied8djed$dked*dled&dmeee�d&dned5d[ed,doeddped/dqed1dred#dsed%dtee�d%dPed%dQd%ed%d&dRed%d&ed%eee�d%dRed&d*ed&ed#ee�d&ed%d(ee�d&edQZ'e'Z(e'duZ)e'duZ*ed%dveddwed/dxed1dyed#dzed%d{d8ed%d8ed1d$ed%d%ed/ee�d#d&ed|edd}ed/d~ed1ded#d�ed%d�ed#ee�d&d�edfd�edhd�ed8d�ed$d�ed*d�ed&d�eee�d�ed5d�ed,d�edd�ed/d�ed1d�ed#d�ed%d�ee�d%d)ed%ed&d%dRed%d&ed%eee�d%dRed&d*ed&ed#ee�d&ed%d(ee�d&edQZ+ed%d4edd�ed/d�ed1d�ed#dXed%d�d8ed#d&ed#d%ed%d(ed/ee�d#d&ed�edd�ed/d�ed1d~ed#d�ed%d*ed#ee�d&ed|ed,d�edd�ed/d�ed1dTed#d>ed%d,ee�d�ed5d�ed,d�edd�ed/d�ed1d�ed#d�ed%d1ee�d%d)ed%ed&d%dRed%d&ed%eee�d%dRed&d*ed&ed#ee�d&ed%d(ee�d&edQZ,ed%d[edd\ed/d]ed1d^ed#d_ed%d`d4ed%d1ed1ed#d#ed%d(ed/ee�d#d8ed,eddaed/dbed1dced#dded%d4ed#ee�d&deedfdgedhdied8djed$dked*dled&dmeee�d&dned5d[ed,doeddped/dqed1dred#dsed%dtee�d%d�ed%d%ed%d(d&dRed%d&ed%eee�d%dRed&d*ed&ed#ee�d&ed%d(ee�d&edQZ-e-Z.d�d�d�d�d��Z/d�d��d�d��d��Z0d�d��Z1d%e	j2Z3d%e	j2ed(ed(dQZ4d%e	j2d%ed%d(ed%ed%d(d�d(Z5d%e	j2ed(d%ed(d�Z6d#e	j2d%ed%d&e	�7ed%d(ed%ed%d(�d/ed%d1d#ee	�7ed%d(ed%ed%d(�e	�7ed%d(ed%ed%d(�d#d#ed%d(ed&e	�7ed%d(ed%ed%d(�ed#e	�7ed%d(ed%ed%d(�d7ed#ed%d%d%ed%d(d%e	�7ed%d(ed%ed%d(�dQZ8e	j2ed%d(ed(d�Z9e	j2d�ed%d(ed%ed%d(d�ed%ed(d%d#ed�ed%d(dQZ:e	j2ed(d%ed(d�Z;e	j2ed%d(d/ed#ed&ed%ed(d#d/ed%e	�7ed%d(ed%ed%d(�d#ed%d(e	�7ed%d(ed%ed%d(�d(ed�ed%d(dRZ<e1d�e/�\Z=Z>Z?Z@ZAZBZCZDZEZFZGZHZIZJZKZLZMZNZOZPe0d�e/d��ZQed%ed%d(eee3e=e4e>e5e?e6e@e8eAe9eCe:eDe;eEe<eF�Reed%d(d�ee/d�eeQi��S�ZTe0d�e/d��ZQed%ed%d(eee3eGe4eHe5eIe6eJe8eKe9eMe:eNe;eOe<eP�Reed%d(d�ee/d�eeQi��S�ZUe1d�e/�\Z=Z>Z?Z@ZAZBZCZDZEZFZGZHZIZJZKZLZMZNZOZPe0d�e/d��ZQed%ed%d(eee3e=e4e>e5e?e6e@e8eAe9eCe:eDe;eEe<eF�Reed%d(d�ee/d�eeQi��S�ZVe0d�e/d��ZQed%ed%d(eee3eGe4eHe5eIe6eJe8eKe9eMe:eNe;eOe<eP�Reed%d(d�ee/d�eeQi��S�ZWe�Xe�Yd�d�d��e�Yd(d#d��f�ZZe�XeZe�Yd#d�d��f�ZZgZ[gZ\gZ]gZ^e_e`eZ��D]~ZaeZeaZbeT�Reebi�Zce[�debecg�eU�Reebi�Zce\�debecg�eV�Reebi�Zce]�debecg�eW�Reebi�Zce^�debecg��qe�ee[�Z[e�ee\�Z\e�ee]�Z]e�ee^�Z^ejfd%d%ddd�d��\ZgZhdQe/d�d(d�e�ie/d�d(�Zjehd�jke[dd�dfd(e[dd�dfe[dd�d(fd(e[dd�dfd%ejd�d��ehd�jledd�dfedd�d(fd�d�d�dehd��md�dġehd��nd�dơehd�jod/d�dȍehd�jp�qdɡehd�jrd�d$d*dˍe/d�d(d�e�ie/d�d(�Zjehd�jke\dd�dfd(e\dd�dfe\dd�d(fd(e\dd�dfd%ejd�d��ehd�jledd�dfedd�d(fd�d�d�dehd�jp�qdΡehd��md�dϡehd�jrd�d$d*dˍehd�jod/d�dȍe/d�d%Zjehd�jke]dd�dfd(e]dd�dfe]dd�d(fd(e]dd�dfd%ejd�d��ehd�jledd�dfedd�d(fd�d�d�dehd�jp�qdԡehd�jod/d�dȍehd��md�d֡ehd�jrd�d$d*dˍe/d�d(d�Zjehd�jke^dd�dfd(e^dd�dfe^dd�d(fd(e^dd�dfd%ejd�d��ehd�jledd�dfedd�d(fd�d�d�dehd�jp�qd١ehd�jod/d�dȍehd��md�dۡehd�jrd�d�d*dˍe_d%�D]Zaehd(eafjsd�d*dfdߍ�q�dS)��
matplotlib�inline�N)�*�
���font)�size�axes)Z	titlesize)�	labelsize�xtick�ytick�legend)�fontsize�)�labelpad�textT)�usetexi,z
figure.dpizRod_particle_ine_drift.txtzDisk_particle_ine_drift.txtzRod_fluid_ine_drift.txtzDisk_fluid_ine_drift.txtz\kappa�xi)�positivez
\overline{xi}�CcCsdt�|d|d�S)N��?�)�sym�log)�val�r�d/Users/davide/Downloads/Joint_paper_figures/jupyter_notebook_submission/figure_3/Dabade2016_drift.py�invcothLsr�����r�(�ix����iFi��iA.�i5$i������	i i4ii�i�i 
ilil6i�Ri�7iAi�
i�ii�i�i�i�i�i�
iQ�:i�����������i�����s�Z��W��!�4izi!i{i�iA	iV���i�i4����
iD�i�i�(iiO�,���mi�i�ii7�g@�Qi�i2i�i��6�l���E���i����ii�i=i��O�$i�i�i�i
��iic�����D�N�I�#iviiii^���i*ii@ii��a�xg)\���(�?gq���h�?g�������?)�pro_pine�obl_pine�pro_fine�obl_finecCs||dddS�Nr!rrr�Zxi0rrr�<lambda>��r\cCs|ddd|SrZrr[rrrr\�r])�prolate�oblatec+Cs

|dk�r�|d}t�ttdddt|i���}t�ttdddt|i���}t�ttdddt|i���}t�ttdddt|i���}t�ttdddt|i���}t	�ttdddt|i���}t
�ttdddt|i���}	t�ttdddt|i���}
t�ttdddt|i���}t
�ttdddt|i���}|d}t�ttdddt|i���}
t�ttdddt|i���}t�ttdddt|i���}t�ttdddt|i���}t�ttdddt|i���}t�ttdddt|i���}t�ttdddt|i���}t�ttdddt|i���}t�ttdddt|i���}t�ttdddt|i���}�n|d}ttd�ttdddttjtdddi�td}t�|�t|i����}ttd�ttdddttjtdddi�td}t�|�t|i����}ttd�ttdddttjtdddi�td}t�|�t|i����}ttd�ttdddttjtdddi�td}t�|�t|i����}ttd�ttdddttjtdddi�td}t�|�t|i����}t	td�ttdddttjtdddi�td}t�|�t|i����}t
td�ttdddttjtdddi�td}t�|�t|i����}	ttd�ttdddttjtdddi�td}t�|�t|i����}
ttd�ttdddttjtdddi�td}t�|�t|i����}t
td�ttdddttjtdddi�td} t�| �t|i����}|d}ttd�ttdddttjtdddi�td}!t�|!�t|i����}
ttd�ttdddttjtdddi�td}"t�|"�t|i����}ttd�ttdddttjtdddi�td}#t�|#�t|i����}ttd�ttdddttjtdddi�td}$t�|$�t|i����}ttd�ttdddttjtdddi�td}%t�|%�t|i����}ttd�ttdddttjtdddi�td}&t�|&�t|i����}ttd�ttdddttjtdddi�td}'t�|'�t|i����}ttd�ttdddttjtdddi�td}(t�|(�t|i����}ttd�ttdddttjtdddi�td})t�|)�t|i����}ttd�ttdddttjtdddi�td}*t�|*�t|i����}|||||||	|
|||
|||||||||fS)	Nr^rVr!rrrXrWrY)�F1p_def�subs�xibr�evalf�F2p_def�F3p_def�F4p_def�F5p_def�F6p_def�G1p_def�G2p_def�G3p_def�G4p_def�F1f_def�F2f_def�F3f_def�F4f_def�F5f_def�F6f_def�G1f_def�G2f_def�G3f_def�G4f_defr�I�re)+�key�eeZe0ZF1p_effZF2p_effZF3p_effZF4p_effZF5p_effZF6p_effZG1p_effZG2p_effZG3p_effZG4p_effZF1f_effZF2f_effZF3f_effZF4f_effZF5f_effZF6f_effZG1f_effZG2f_effZG3f_effZG4f_effZF1p_tmpZF2p_tmpZF3p_tmpZF4p_tmpZF5p_tmpZF6p_tmpZG1p_tmpZG2p_tmpZG3p_tmpZG4p_tmpZF1f_tmpZF2f_tmpZF3f_tmpZF4f_tmpZF5f_tmpZF6f_tmpZG1f_tmpZG2f_tmpZG3f_tmpZG4f_tmprrr�shape_coeffs�s�
"""""""""""""""""""&@@@@@@@@@@@@@@@@@@@@�r{g�����������rr^rVrXr_rWrYg�h㈵��>g�������?�d)r r)�constrained_layout�sharex�figsizeg�?)rrz$\mathit{\xi_0 = 1.01}$)�label�kz$$\mathit{Dabade \ et \ al., \ 2016}$z--)�cr��lsg{�G�z��g�������?gg�?zlower center)r�locz($\mathit{Prolate, \ particle \ inertia}$zw$\mathit{ \frac{  {St}^{-1} \left( \Delta C_p / (C^2+1)  \right) }{ \left( (\xi_0^2-1)^{3/2} \log(\xi_0-1) \right)}  }$)rr)rrz$\mathit{\xi_0 = 1.1}$z%$\mathit{Prolate, \ fluid \ inertia}$g�������?zy$\mathit{\frac{{Re_p}^{-1} \left( \Delta C_p / (C^2+1)  \right) }{ \left( (\xi_0-1)^{0.5} \log (\xi_0-1) \right)^{-1}} }$zupper center)rrz$\mathit{\xi_0 = 1.0001}$z'$\mathit{Oblate, \ particle \ inertia}$g)\���(��g{�G�z�?zB$\mathit{{St}^{-1} \left( \Delta C_p / (C^2+1)  \right) \xi_0^2 }$)rrz$$\mathit{Oblate, \ fluid \ inertia}$g{�G�zĿg{�G�z�?zL$\mathit{{Re_p}^{-1} \left( \Delta C_f / (C^2+1)  \right) (\xi_0-1)^{0.5} }$g@z$\mathit{C/(C+1)}$)rr)t�get_ipython�run_line_magic�numpy�np�pylab�matplotlib.pyplot�pyplot�pltZsympyrZ
SMALL_SIZEZMEDIUM_SIZEZBIGGER_SIZE�rc�mpl�rcParams�loadtxtZrod_pineZ	disk_pineZrod_fineZ	disk_fine�symbolsr�rrbrrr`rdrerfrgrhrirjrkrlrmrnrorprqrrrsrtrurvrz�rrr{�pi�I1�I2�I3�I4�sqrt�I5I6�J1�J2�J3�J4�F1p�F2p�F3p�F4p�F5p�F6p�G1p�G2p�G3p�G4p�F1f�F2f�F3f�F4f�F5f�F6f�G1f�G2f�G3f�G4f�rrarcZPINE_proZFINE_pro�PINE_obl�FINE_obl�concatenate�linspace�CrangeZCDC_pro_pineZCDC_pro_fineZCDC_obl_pineZCDC_obl_fine�range�len�jZCVAL�data�append�array�subplots�fig�axsrZftr�scatter�plot�set_ylim�set_xlimr�title�set_text�
set_ylabel�
set_xlabelrrrr�<module>s�





	

		T@B
R>RR�J��
��<
�
����$
�
�
�
�
��
��F&����<*�B�F�>&���FB�(
�
���
��"
�
�
����
�
�
�
�
��
��X����

F:�(
�
���
��"
�
�
����V
��>&����H:�>��
��N��.
�
�
��
��>&����FB�Z�N��H
��
��F&����	��	Q
6"L,�$�@�D��8�`j8���	2@��	@��	2@��	@��	"



	&\0$\0\0\0