Demo problem: Turek & Hron's FSI benchmark problem - Oomph-lib
Demo problem: Turek & Hron's FSI benchmark problem - Oomph-lib
Demo problem: Turek & Hron's FSI benchmark problem - Oomph-lib
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
8 <strong>Demo</strong> <strong>problem</strong>: <strong>Turek</strong> & Hron’s <strong>FSI</strong> <strong>benchmark</strong> <strong>problem</strong><br />
}<br />
// <strong>FSI</strong> parameter<br />
Q=1.429e-6;<br />
// Timestep -- aiming for about 40 steps per period<br />
Dt=0.1;<br />
// Density ratio<br />
Density_ratio=1.0;<br />
// Gravity<br />
Gravity=0.0;<br />
// Max. flux<br />
Max_flux=1.0;<br />
// Ignore fluid<br />
Ignore_fluid_loading=false;<br />
// Compute dependent parameters<br />
// Timescale ratio for solid<br />
Lambda_sq=Re*Q*Density_ratio*St*St;<br />
In the interest of brevity we omit the listings of the assignments for the other cases. Finally, we select the<br />
length of the time interval over which the influx is ramped up from zero to its maximum value to be equal<br />
to 20 timesteps, create a constitutive equation for the solid, and document the parameter values used in the<br />
simulation:<br />
// Ramp period (20 timesteps)<br />
Ramp_period=Dt*20.0;<br />
// "Big G" Linear constitutive equations:<br />
Constitutive_law_pt = new GeneralisedHookean(&Nu,&E);<br />
// Doc<br />
oomph_info