25.02.2015 Views

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

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!