~fluidity-core/fluidity/shallow-water-dev

to get a darcy equation the time, stress and advection terms are not included in a bousinessq formulation (to avoid needing to specify a density). Ideally it would be better to actually have a darcy momentum option to simplify the input.

the absorption term is included in the pressure correction (being a necessity as it is the only term in the matrix equation)</comment>

</simulation_name>

<problem_type>

<string_value lines="1">fluids</string_value>

</problem_type>

<integer_value rank="0">3</integer_value>

</dimension>

<from_file file_name="cube">

<stat>

<include_in_stat/>

</stat>

</from_file>

</mesh>

<from_mesh>

<mesh_shape>

<polynomial_degree>

<integer_value rank="0">0</integer_value>

</polynomial_degree>

</mesh_shape>

<mesh_continuity>

<string_value>discontinuous</string_value>

</mesh_continuity>

<stat>

<exclude_from_stat/>

</stat>

</from_mesh>

</mesh>

<from_mesh>

<mesh_shape>

<polynomial_degree>

<integer_value rank="0">1</integer_value>

</polynomial_degree>

</mesh_shape>

<stat>

<exclude_from_stat/>

</stat>

</from_mesh>

</mesh>

<from_mesh>

<mesh_shape>

<polynomial_degree>

<integer_value rank="0">0</integer_value>

</polynomial_degree>

</mesh_shape>

<mesh_continuity>

<string_value>discontinuous</string_value>

</mesh_continuity>

<stat>

<exclude_from_stat/>

</stat>

</from_mesh>

</mesh>

<from_mesh>

<mesh_continuity>

<string_value>discontinuous</string_value>

</mesh_continuity>

<stat>

<exclude_from_stat/>

</stat>

</from_mesh>

</mesh>

<integer_value rank="0">3</integer_value>

</degree>

</quadrature>

</geometry>

<io>

<dump_format>

<string_value>vtk</string_value>

</dump_format>

<dump_period_in_timesteps>

100

<integer_value rank="0">0</integer_value>

101

</constant>

102

</dump_period_in_timesteps>

103

<output_mesh name="OutputMesh"/>

104

<stat/>

105

</io>

106

107

<current_time>

108

<real_value rank="0">0.0</real_value>

109

</current_time>

110

111

<real_value rank="0">0.1</real_value>

112

</timestep>

113

<finish_time>

114

<real_value rank="0">0.5</real_value>

115

</finish_time>

116

<nonlinear_iterations>

117

<integer_value rank="0">2</integer_value>

118

</nonlinear_iterations>

119

</timestepping>

120

<material_phase name="phase1">

121

<scalar_field name="Pressure" rank="0">

122

123

124

<spatial_discretisation>

125

<continuous_galerkin>

126

<test_continuity_with_cv_dual/>

127

</continuous_galerkin>

128

</spatial_discretisation>

129

130

<poisson_pressure_solution>

131

<string_value lines="1">never</string_value>

132

</poisson_pressure_solution>

133

<use_projection_method/>

134

</scheme>

135

136

<iterative_method name="gmres">

137

138

<integer_value rank="0">30</integer_value>

139

</restart>

140

</iterative_method>

141

142

<relative_error>

143

<real_value rank="0">1.0e-10</real_value>

144

</relative_error>

145

<max_iterations>

146

<integer_value rank="0">1000</integer_value>

147

</max_iterations>

148

<never_ignore_solver_failures/>

149

150

151

</diagnostics>

152

</solver>

153

<boundary_conditions name="left_inflow">

154

<surface_ids>

155

<integer_value shape="1" rank="1">3</integer_value>

156

</surface_ids>

157

158

159

<real_value rank="0">2.0e+05</real_value>

160

</constant>

161

</type>

162

</boundary_conditions>

163

<boundary_conditions name="right_outflow">

164

<surface_ids>

165

<integer_value shape="1" rank="1">4</integer_value>

166

</surface_ids>

167

168

169

<real_value rank="0">0.0</real_value>

170

</constant>

171

</type>

172

</boundary_conditions>

173

174

<stat/>

175

176

<include_in_convergence/>

177

</convergence>

178

179

<exclude_from_detectors/>

180

</detectors>

181

<steady_state>

182

<include_in_steady_state/>

183

</steady_state>

184

<no_interpolation/>

185

</prognostic>

186

</scalar_field>

187

<vector_field name="Velocity" rank="1">

188

189

190

191

<spatial_discretisation>

192

<discontinuous_galerkin>

193

<mass_terms>

194

<exclude_mass_terms/>

195

</mass_terms>

196

<viscosity_scheme>

197

<compact_discontinuous_galerkin/>

198

</viscosity_scheme>

199

<advection_scheme>

200

<none/>

201

<integrate_advection_by_parts>

202

203

</integrate_advection_by_parts>

204

</advection_scheme>

205

</discontinuous_galerkin>

206

<conservative_advection>

207

<real_value rank="0">0.0</real_value>

208

</conservative_advection>

209

</spatial_discretisation>

210

<temporal_discretisation>

211

<theta>

212

<real_value rank="0">1.0</real_value>

213

</theta>

214

215

<real_value rank="0">1.0</real_value>

216

</relaxation>

217

</temporal_discretisation>

218

219

<iterative_method name="gmres">

220

221

<integer_value rank="0">30</integer_value>

222

</restart>

223

</iterative_method>

224

225

<relative_error>

226

<real_value rank="0">1.0e-10</real_value>

227

</relative_error>

228

<max_iterations>

229

<integer_value rank="0">1000</integer_value>

230

</max_iterations>

231

<never_ignore_solver_failures/>

232

233

234

</diagnostics>

235

</solver>

236

<initial_condition name="WholeMesh">

237

238

<real_value shape="3" dim1="dim" rank="1">0.0 0.0 0.0</real_value>

239

</constant>

240

</initial_condition>

241

<boundary_conditions name="no_outflow_top_bottom">

242

<surface_ids>

243

<integer_value shape="2" rank="1">1 2</integer_value>

244

</surface_ids>

245

246

<align_bc_with_cartesian>

247

<z_component>

248

249

<real_value rank="0">0.0</real_value>

250

</constant>

251

</z_component>

252

</align_bc_with_cartesian>

253

</type>

254

</boundary_conditions>

255

<boundary_conditions name="no_outflow_front_back">

256

<surface_ids>

257

<integer_value shape="2" rank="1">5 6</integer_value>

258

</surface_ids>

259

260

<align_bc_with_cartesian>

261

<y_component>

262

263

<real_value rank="0">0.0</real_value>

264

</constant>

265

</y_component>

266

</align_bc_with_cartesian>

267

</type>

268

</boundary_conditions>

269

<vector_field name="Absorption" rank="1">

270

271

272

273

<string_value lines="20" type="python">from numpy import zeros

274

275

# set the dimension

276

dim = 3

277

278

# 1/4

279

quarter = 1.0/4.0

280

281

# set the absorption coeff

282

abs_coeff = 100.0/0.2

283

284

sat = states["phase1"].scalar_fields["PhaseVolumeFraction"]

285

286

for ele in range(field.element_count):

287

288

sat_ele = sat.ele_val(ele)

289

290

# initialise local list of abs coeff

291

sigma_ele = zeros(dim)

292

293

# calc the absorption term for each dim

294

for d in range(dim):

295

296

sigma_ele[d] = abs_coeff * quarter*sum(sat_ele)

297

298

# set the absorption term

299

field.set(field.ele_nodes(ele),sigma_ele)</string_value>

300

</algorithm>

301

302

<stat>

303

<include_in_stat/>

304

</stat>

305

306

<include_in_convergence/>

307

</convergence>

308

309

<include_in_detectors/>

310

</detectors>

311

<steady_state>

312

<include_in_steady_state/>

313

</steady_state>

314

</diagnostic>

315

<include_pressure_correction/>

316

</vector_field>

317

318

<stat>

319

<include_in_stat/>

320

<previous_time_step>

321

<exclude_from_stat/>

322

</previous_time_step>

323

<nonlinear_field>

324

<exclude_from_stat/>

325

</nonlinear_field>

326

</stat>

327

328

<include_in_convergence/>

329

</convergence>

330

331

<include_in_detectors/>

332

</detectors>

333

<steady_state>

334

<include_in_steady_state/>

335

</steady_state>

336

<consistent_interpolation/>

337

</prognostic>

338

</vector_field>

339

<scalar_field name="PhaseVolumeFraction" rank="0">

340

341

342

343

<spatial_discretisation>

344

<control_volumes>

345

<face_value name="FiniteElement">

346

<limit_face_value>

347

348

</limit_face_value>

349

</face_value>

350

<diffusion_scheme name="ElementGradient"/>

351

</control_volumes>

352

<conservative_advection>

353

<real_value rank="0">1.0</real_value>

354

</conservative_advection>

355

</spatial_discretisation>

356

<temporal_discretisation>

357

<theta>

358

<real_value rank="0">1.0</real_value>

359

</theta>

360

<control_volumes>

361

<number_advection_iterations>

362

<integer_value rank="0">3</integer_value>

363

</number_advection_iterations>

364

<limit_theta/>

365

</control_volumes>

366

</temporal_discretisation>

367

368

<iterative_method name="gmres">

369

370

<integer_value rank="0">30</integer_value>

371

</restart>

372

</iterative_method>

373

374

<relative_error>

375

<real_value rank="0">1.0e-10</real_value>

376

</relative_error>

377

<max_iterations>

378

<integer_value rank="0">10000</integer_value>

379

</max_iterations>

380

<never_ignore_solver_failures/>

381

382

383

</diagnostics>

384

</solver>

385

<initial_condition name="WholeMesh">

386

387

<real_value rank="0">0.2</real_value>

388

</constant>

389

</initial_condition>

390

<boundary_conditions name="left_inflow">

391

<surface_ids>

392

<integer_value shape="1" rank="1">3</integer_value>

393

</surface_ids>

394

395

<apply_weakly/>

396

397

<real_value rank="0">0.2</real_value>

398

</constant>

399

</type>

400

</boundary_conditions>

401

<boundary_conditions name="no_outflow_top_bottom">

402

<surface_ids>

403

<integer_value shape="2" rank="1">1 2</integer_value>

404

</surface_ids>

405

406

</boundary_conditions>

407

<boundary_conditions name="no_outflow_front_back">

408

<surface_ids>

409

<integer_value shape="2" rank="1">5 6</integer_value>

410

</surface_ids>

411

412

</boundary_conditions>

413

414

415

<stat/>

416

417

<include_in_convergence/>

418

</convergence>

419

420

<include_in_detectors/>

421

</detectors>

422

<steady_state>

423

<include_in_steady_state/>

424

</steady_state>

425

<consistent_interpolation/>

426

</prognostic>

427

</scalar_field>

428

<scalar_field name="SumVelocityDivergence" rank="0">

429

430

431

432

<test_with_cv_dual/>

433

434

<stat/>

435

436

<include_in_detectors/>

437

</detectors>

438

439

<iterative_method name="cg"/>

440

441

<relative_error>

442

<real_value rank="0">1.0e-10</real_value>

443

</relative_error>

444

<max_iterations>

445

<integer_value rank="0">1000</integer_value>

446

</max_iterations>

447

<never_ignore_solver_failures/>

448

449

450

</diagnostics>

451

</solver>

452

</diagnostic>

453

</scalar_field>

454

<scalar_field name="SumPhaseVolumeFraction" rank="0">

455

456

457

<string_value lines="20" type="python">sat_1 = states["phase1"].scalar_fields["PhaseVolumeFraction"]

458

sat_2 = states["phase2"].scalar_fields["PhaseVolumeFraction"]

459

460

for n in range(field.node_count):

461

462

sat_1_n = sat_1.node_val(n)

463

sat_2_n = sat_2.node_val(n)

464

465

field.set(n,sat_1_n + sat_2_n)</string_value>

466

</algorithm>

467

468

469

<stat/>

470

471

<include_in_convergence/>

472

</convergence>

473

474

<include_in_detectors/>

475

</detectors>

476

<steady_state>

477

<include_in_steady_state/>

478

</steady_state>

479

</diagnostic>

480

</scalar_field>

481

<scalar_field name="ControlVolumeCFLNumber" rank="0">

482

483

484

485

486

<stat/>

487

488

<include_in_convergence/>

489

</convergence>

490

491

<include_in_detectors/>

492

</detectors>

493

<steady_state>

494

<include_in_steady_state/>

495

</steady_state>

496

</diagnostic>

497

</scalar_field>

498

</material_phase>

499

<material_phase name="phase2">

500

<scalar_field name="Pressure" rank="0">

501

502

</scalar_field>

503

<vector_field name="Velocity" rank="1">

504

505

506

507

<spatial_discretisation>

508

<discontinuous_galerkin>

509

<mass_terms>

510

<exclude_mass_terms/>

511

</mass_terms>

512

<viscosity_scheme>

513

<compact_discontinuous_galerkin/>

514

</viscosity_scheme>

515

<advection_scheme>

516

<none/>

517

<integrate_advection_by_parts>

518

519

</integrate_advection_by_parts>

520

</advection_scheme>

521

</discontinuous_galerkin>

522

<conservative_advection>

523

<real_value rank="0">0.0</real_value>

524

</conservative_advection>

525

</spatial_discretisation>

526

<temporal_discretisation>

527

<theta>

528

<real_value rank="0">1.0</real_value>

529

</theta>

530

531

<real_value rank="0">1.0</real_value>

532

</relaxation>

533

</temporal_discretisation>

534

535

<iterative_method name="gmres">

536

537

<integer_value rank="0">30</integer_value>

538

</restart>

539

</iterative_method>

540

541

<relative_error>

542

<real_value rank="0">1.0e-10</real_value>

543

</relative_error>

544

<max_iterations>

545

<integer_value rank="0">1000</integer_value>

546

</max_iterations>

547

<never_ignore_solver_failures/>

548

549

550

</diagnostics>

551

</solver>

552

<initial_condition name="WholeMesh">

553

554

<real_value shape="3" dim1="dim" rank="1">0.0 0.0 0.0</real_value>

555

</constant>

556

</initial_condition>

557

<boundary_conditions name="no_outflow_top_bottom">

558

<surface_ids>

559

<integer_value shape="2" rank="1">1 2</integer_value>

560

</surface_ids>

561

562

<align_bc_with_cartesian>

563

<z_component>

564

565

<real_value rank="0">0.0</real_value>

566

</constant>

567

</z_component>

568

</align_bc_with_cartesian>

569

</type>

570

</boundary_conditions>

571

<boundary_conditions name="no_outflow_front_back">

572

<surface_ids>

573

<integer_value shape="2" rank="1">5 6</integer_value>

574

</surface_ids>

575

576

<align_bc_with_cartesian>

577

<y_component>

578

579

<real_value rank="0">0.0</real_value>

580

</constant>

581

</y_component>

582

</align_bc_with_cartesian>

583

</type>

584

</boundary_conditions>

585

<vector_field name="Absorption" rank="1">

586

587

588

589

<string_value lines="20" type="python">from numpy import zeros

590

591

# set the dimension

592

dim = 3

593

594

# 1/4

595

quarter = 1.0/4.0

596

597

# set the absorption coeff

598

abs_coeff = 100.0/0.8

599

600

sat = states["phase2"].scalar_fields["PhaseVolumeFraction"]

601

602

for ele in range(field.element_count):

603

604

sat_ele = sat.ele_val(ele)

605

606

# initialise local list of abs coeff

607

sigma_ele = zeros(dim)

608

609

# calc the absorption term for each dim

610

for d in range(dim):

611

612

sigma_ele[d] = abs_coeff * quarter*sum(sat_ele)

613

614

# set the absorption term

615

field.set(field.ele_nodes(ele),sigma_ele)</string_value>

616

</algorithm>

617

618

<stat>

619

<include_in_stat/>

620

</stat>

621

622

<include_in_convergence/>

623

</convergence>

624

625

<include_in_detectors/>

626

</detectors>

627

<steady_state>

628

<include_in_steady_state/>

629

</steady_state>

630

</diagnostic>

631

<include_pressure_correction/>

632

</vector_field>

633

634

<stat>

635

<include_in_stat/>

636

<previous_time_step>

637

<exclude_from_stat/>

638

</previous_time_step>

639

<nonlinear_field>

640

<exclude_from_stat/>

641

</nonlinear_field>

642

</stat>

643

644

<include_in_convergence/>

645

</convergence>

646

647

<include_in_detectors/>

648

</detectors>

649

<steady_state>

650

<include_in_steady_state/>

651

</steady_state>

652

<consistent_interpolation/>

653

</prognostic>

654

</vector_field>

655

<scalar_field name="PhaseVolumeFraction" rank="0">

656

657

658

659

660

<stat/>

661

662

<include_in_detectors/>

663

</detectors>

664

</diagnostic>

665

</scalar_field>

666

</material_phase>

667

<porous_media>

668

<scalar_field name="Porosity">

669

670

671

672

673

<real_value rank="0">0.2</real_value>

674

</constant>

675

</value>

676

677

<stat/>

678

679

<exclude_from_detectors/>

680

</detectors>

681

</prescribed>

682

</scalar_field>

683

</porous_media>

684

</fluidity_options>

Older »