From cec2146cc2e4507635ef4f31efc0ba81d50d0f9c Mon Sep 17 00:00:00 2001 From: Haoyuan Li Date: Sat, 11 Jul 2026 12:09:56 -0700 Subject: [PATCH] Prescribed condition with particles --- doc/modules/changes/20260711_lhy11009 | 3 + .../particle/property/initial_composition.h | 7 + .../particle/property/initial_composition.cc | 38 ++ ...omposition_initial_condition_particles.prm | 65 +++ .../screen-output | 31 ++ .../solution/solution-00000.0000.gnuplot | 455 ++++++++++++++++++ .../solution/solution-00001.0000.gnuplot | 455 ++++++++++++++++++ 7 files changed, 1054 insertions(+) create mode 100644 doc/modules/changes/20260711_lhy11009 create mode 100644 tests/prescribed_solution_composition_initial_condition_particles.prm create mode 100644 tests/prescribed_solution_composition_initial_condition_particles/screen-output create mode 100644 tests/prescribed_solution_composition_initial_condition_particles/solution/solution-00000.0000.gnuplot create mode 100644 tests/prescribed_solution_composition_initial_condition_particles/solution/solution-00001.0000.gnuplot diff --git a/doc/modules/changes/20260711_lhy11009 b/doc/modules/changes/20260711_lhy11009 new file mode 100644 index 00000000000..f7e05640500 --- /dev/null +++ b/doc/modules/changes/20260711_lhy11009 @@ -0,0 +1,3 @@ +Changed: add support for updating the InitialComposition particle property using the prescribed solution plugin. +
+(Haoyuan Li, 2026/07/11) diff --git a/include/aspect/particle/property/initial_composition.h b/include/aspect/particle/property/initial_composition.h index 9cc4cdb8228..d7b0fdd8510 100644 --- a/include/aspect/particle/property/initial_composition.h +++ b/include/aspect/particle/property/initial_composition.h @@ -57,6 +57,13 @@ namespace aspect initialize_one_particle_property (const Point &position, std::vector &particle_properties) const override; + /** + * @copydoc aspect::Particle::Property::Interface::update_particle_properties() + */ + void + update_particle_properties (const ParticleUpdateInputs &inputs, + typename ParticleHandler::particle_iterator_range &particles) const override; + /** * Returns an enum, which determines how this particle property is * initialized for particles that are created later than the initial diff --git a/source/particle/property/initial_composition.cc b/source/particle/property/initial_composition.cc index 8c452052041..c62090d7e6a 100644 --- a/source/particle/property/initial_composition.cc +++ b/source/particle/property/initial_composition.cc @@ -20,6 +20,7 @@ #include #include +#include namespace aspect { @@ -80,6 +81,43 @@ namespace aspect return property_information; } + + template + void + InitialComposition::update_particle_properties(const ParticleUpdateInputs &/*inputs*/, + typename ParticleHandler::particle_iterator_range &particles) const + { + const aspect::PrescribedSolution::Manager &prescribed_solution_manager = this->get_prescribed_solution(); + const auto &plugin_objects = prescribed_solution_manager.get_active_plugins(); + + + std::vector> evaluation_points(1); + std::vector component_indices(1); + std::vector component_is_constrained(1); + std::vector constrained_component_value(1); + + for (auto &particle : particles) + { + evaluation_points[0] = particle.get_location(); + + for (unsigned int j = 0; j < this->n_compositional_fields(); ++j) + { + component_indices[0] = this->introspection().component_indices.compositional_fields[j]; + component_is_constrained[0] = false; + constrained_component_value[0] = 0.0; + + for (const auto &plugin : plugin_objects) + { + typename DoFHandler::active_cell_iterator cell; + plugin->constrain_solution(cell, evaluation_points, component_indices, component_is_constrained, constrained_component_value); + } + + if (component_is_constrained[0]) + particle.get_properties()[this->data_position+j] = constrained_component_value[0]; + } + } + } + } } } diff --git a/tests/prescribed_solution_composition_initial_condition_particles.prm b/tests/prescribed_solution_composition_initial_condition_particles.prm new file mode 100644 index 00000000000..7154c44b386 --- /dev/null +++ b/tests/prescribed_solution_composition_initial_condition_particles.prm @@ -0,0 +1,65 @@ +######################################################### +# This is modified from the composition_passive.prm +# parameter file by prescribing the compositional solution +# from initial compositional fields and using the particle +# method + +include $ASPECT_SOURCE_DIR/cookbooks/composition_passive/composition_passive.prm + +# Modify the original cookbook to end earlier +set Start time = 0 +set End time = 0.1 +set Use years instead of seconds = false + +# Modify the original cookbook to have coarser meshes +subsection Mesh refinement + set Initial adaptive refinement = 0 + set Initial global refinement = 3 + set Time steps between mesh refinement = 0 +end + +# Modify the compositional fields to use particles +subsection Compositional fields + set Number of fields = 2 + set Compositional field methods = particles, particles + set Mapped particle properties = C_1: initial C_1, C_2: initial C_2 +end + +# Include the particles +subsection Particles + set Minimum particles per cell = 25 + set Maximum particles per cell = 100 + set Load balancing strategy = remove and add particles + set List of particle properties = initial composition, position + set Interpolation scheme = bilinear least squares + set Update ghost particles = true + set Particle generator name = reference cell + subsection Generator + subsection Reference cell + set Number of particles per cell per direction = 7 + end + end +end + +# The next section prescribes the composition inside the domain +# using a spatially varying function. An indicator function selects +# the region where the composition is constrained. +subsection Prescribed solution + set List of model names = initial composition + subsection Initial composition + subsection Indicator function + set Variable names = x, y + set Function expression = (x<1) ? 1:0 + end + end +end + +subsection Postprocess + set List of postprocessors = visualization + + subsection Visualization + set Interpolate output = false + set Time between graphical output = 0.1 + set Output format = gnuplot + end +end diff --git a/tests/prescribed_solution_composition_initial_condition_particles/screen-output b/tests/prescribed_solution_composition_initial_condition_particles/screen-output new file mode 100644 index 00000000000..f8789af7171 --- /dev/null +++ b/tests/prescribed_solution_composition_initial_condition_particles/screen-output @@ -0,0 +1,31 @@ + +Number of active cells: 64 (on 4 levels) +Number of degrees of freedom: 1,526 (578+81+289+289+289) + +*** Timestep 0: t=0 seconds, dt=0 seconds + Solving temperature system... 0 iterations. + Advecting particles... done. + Solving Stokes system (GMG)... 11+0 iterations. + + Postprocessing: + Writing graphical output: output-prescribed_solution_composition_initial_condition_particles/solution/solution-00000 + +*** Timestep 1: t=0.0625 seconds, dt=0.0625 seconds + Solving temperature system... 10 iterations. + Advecting particles... done. + Solving Stokes system (GMG)... 7+0 iterations. + + Postprocessing: + +*** Timestep 2: t=0.1 seconds, dt=0.0375 seconds + Solving temperature system... 10 iterations. + Advecting particles... done. + Solving Stokes system (GMG)... 12+0 iterations. + + Postprocessing: + Writing graphical output: output-prescribed_solution_composition_initial_condition_particles/solution/solution-00001 + +Termination requested by criterion: end time + + + diff --git a/tests/prescribed_solution_composition_initial_condition_particles/solution/solution-00000.0000.gnuplot b/tests/prescribed_solution_composition_initial_condition_particles/solution/solution-00000.0000.gnuplot new file mode 100644 index 00000000000..27f61771c61 --- /dev/null +++ b/tests/prescribed_solution_composition_initial_condition_particles/solution/solution-00000.0000.gnuplot @@ -0,0 +1,455 @@ +# This file was generated by the deal.II library. + + +# +# For a description of the GNUPLOT format see the GNUPLOT manual. +# +#

+0 0 0 0 3397.39 1 1 0 +0.25 0 0.178855 0 3397.17 1 1 0 + +0 0.125 0 -0.099234 2972.9 0.875 1 0 +0.25 0.125 0.179098 -0.071298 2972.66 0.875 1 0 + + +0.25 0 0.178855 0 3397.17 1 1 0 +0.5 0 0.260099 0 3396.64 1 1 0 + +0.25 0.125 0.179098 -0.071298 2972.66 0.875 1 0 +0.5 0.125 0.259728 -0.00659826 2972.09 0.875 1 0 + + +0 0.125 0 -0.099234 2972.9 0.875 1 0 +0.25 0.125 0.179098 -0.071298 2972.66 0.875 1 0 + +0 0.25 0 -0.199109 2548.54 0.75 0 0 +0.25 0.25 0.178238 -0.141739 2548.24 0.75 0 0 + + +0.25 0.125 0.179098 -0.071298 2972.66 0.875 1 0 +0.5 0.125 0.259728 -0.00659826 2972.09 0.875 1 0 + +0.25 0.25 0.178238 -0.141739 2548.24 0.75 0 0 +0.5 0.25 0.255646 -0.0129817 2547.55 0.75 0 0 + + +0.5 0 0.260099 0 3396.64 1 1 0 +0.75 0 0.203464 0 3396.08 1 1 0 + +0.5 0.125 0.259728 -0.00659826 2972.09 0.875 1 0 +0.75 0.125 0.203477 0.0608429 2971.5 0.875 1 0 + + +0.75 0 0.203464 0 3396.08 1 1 0 +1 0 0.0388256 0 3395.79 1 1 0 + +0.75 0.125 0.203477 0.0608429 2971.5 0.875 1 0 +1 0.125 0.038285 0.095923 2971.19 0.875 1 0 + + +0.5 0.125 0.259728 -0.00659826 2972.09 0.875 1 0 +0.75 0.125 0.203477 0.0608429 2971.5 0.875 1 0 + +0.5 0.25 0.255646 -0.0129817 2547.55 0.75 0 0 +0.75 0.25 0.201694 0.1209 2546.83 0.75 0 0 + + +0.75 0.125 0.203477 0.0608429 2971.5 0.875 1 0 +1 0.125 0.038285 0.095923 2971.19 0.875 1 0 + +0.75 0.25 0.201694 0.1209 2546.83 0.75 0 0 +1 0.25 0.0366056 0.193194 2546.45 0.75 0 0 + + +0 0.25 0 -0.199109 2548.54 0.75 0 0 +0.25 0.25 0.178238 -0.141739 2548.24 0.75 0 0 + +0 0.375 0 -0.299157 2124.33 0.625 0 0 +0.25 0.375 0.171121 -0.208281 2123.91 0.625 0 0 + + +0.25 0.25 0.178238 -0.141739 2548.24 0.75 0 0 +0.5 0.25 0.255646 -0.0129817 2547.55 0.75 0 0 + +0.25 0.375 0.171121 -0.208281 2123.91 0.625 0 0 +0.5 0.375 0.237993 -0.0186818 2123.01 0.625 0 0 + + +0 0.375 0 -0.299157 2124.33 0.625 0 0 +0.25 0.375 0.171121 -0.208281 2123.91 0.625 0 0 + +0 0.5 0 -0.396595 1700.36 0.5 0 0 +0.25 0.5 0.147628 -0.262746 1699.69 0.5 0 0 + + +0.25 0.375 0.171121 -0.208281 2123.91 0.625 0 0 +0.5 0.375 0.237993 -0.0186818 2123.01 0.625 0 0 + +0.25 0.5 0.147628 -0.262746 1699.69 0.5 0 0 +0.5 0.5 0.187137 -0.0225651 1698.47 0.5 0 0 + + +0.5 0.25 0.255646 -0.0129817 2547.55 0.75 0 0 +0.75 0.25 0.201694 0.1209 2546.83 0.75 0 0 + +0.5 0.375 0.237993 -0.0186818 2123.01 0.625 0 0 +0.75 0.375 0.192122 0.177108 2122.08 0.625 0 0 + + +0.75 0.25 0.201694 0.1209 2546.83 0.75 0 0 +1 0.25 0.0366056 0.193194 2546.45 0.75 0 0 + +0.75 0.375 0.192122 0.177108 2122.08 0.625 0 0 +1 0.375 0.033547 0.291924 2121.54 0.625 0 0 + + +0.5 0.375 0.237993 -0.0186818 2123.01 0.625 0 0 +0.75 0.375 0.192122 0.177108 2122.08 0.625 0 0 + +0.5 0.5 0.187137 -0.0225651 1698.47 0.5 0 0 +0.75 0.5 0.162803 0.221016 1697.24 0.5 0 0 + + +0.75 0.375 0.192122 0.177108 2122.08 0.625 0 0 +1 0.375 0.033547 0.291924 2121.54 0.625 0 0 + +0.75 0.5 0.162803 0.221016 1697.24 0.5 0 0 +1 0.5 0.0283452 0.389607 1696.39 0.5 0 0 + + +1 0 0.0388256 0 3395.79 1 1 0 +1.25 0 -0.139022 0 3395.94 1 1 0 + +1 0.125 0.038285 0.095923 2971.19 0.875 1 0 +1.25 0.125 -0.140949 0.0729854 2971.34 0.875 1 0 + + +1.25 0 -0.139022 0 3395.94 1 1 0 +1.5 0 -0.222486 0 3396.41 1 1 0 + +1.25 0.125 -0.140949 0.0729854 2971.34 0.875 1 0 +1.5 0.125 -0.224567 0.00659826 2971.86 0.875 1 0 + + +1 0.125 0.038285 0.095923 2971.19 0.875 1 0 +1.25 0.125 -0.140949 0.0729854 2971.34 0.875 1 0 + +1 0.25 0.0366056 0.193194 2546.45 0.75 0 0 +1.25 0.25 -0.145244 0.146464 2546.65 0.75 0 0 + + +1.25 0.125 -0.140949 0.0729854 2971.34 0.875 1 0 +1.5 0.125 -0.224567 0.00659826 2971.86 0.875 1 0 + +1.25 0.25 -0.145244 0.146464 2546.65 0.75 0 0 +1.5 0.25 -0.228033 0.0129817 2547.3 0.75 0 0 + + +1.5 0 -0.222486 0 3396.41 1 1 0 +1.75 0 -0.163631 0 3396.91 1 1 0 + +1.5 0.125 -0.224567 0.00659826 2971.86 0.875 1 0 +1.75 0.125 -0.165328 -0.0625304 2972.41 0.875 1 0 + + +1.75 0 -0.163631 0 3396.91 1 1 0 +2 0 0 0 3397.13 1 1 0 + +1.75 0.125 -0.165328 -0.0625304 2972.41 0.875 1 0 +2 0.125 0 -0.092612 2972.64 0.875 1 0 + + +1.5 0.125 -0.224567 0.00659826 2971.86 0.875 1 0 +1.75 0.125 -0.165328 -0.0625304 2972.41 0.875 1 0 + +1.5 0.25 -0.228033 0.0129817 2547.3 0.75 0 0 +1.75 0.25 -0.1687 -0.125625 2547.97 0.75 0 0 + + +1.75 0.125 -0.165328 -0.0625304 2972.41 0.875 1 0 +2 0.125 0 -0.092612 2972.64 0.875 1 0 + +1.75 0.25 -0.1687 -0.125625 2547.97 0.75 0 0 +2 0.25 0 -0.187279 2548.27 0.75 0 0 + + +1 0.25 0.0366056 0.193194 2546.45 0.75 0 0 +1.25 0.25 -0.145244 0.146464 2546.65 0.75 0 0 + +1 0.375 0.033547 0.291924 2121.54 0.625 0 0 +1.25 0.375 -0.147055 0.218878 2121.84 0.625 0 0 + + +1.25 0.25 -0.145244 0.146464 2546.65 0.75 0 0 +1.5 0.25 -0.228033 0.0129817 2547.3 0.75 0 0 + +1.25 0.375 -0.147055 0.218878 2121.84 0.625 0 0 +1.5 0.375 -0.223447 0.0186818 2122.72 0.625 0 0 + + +1 0.375 0.033547 0.291924 2121.54 0.625 0 0 +1.25 0.375 -0.147055 0.218878 2121.84 0.625 0 0 + +1 0.5 0.0283452 0.389607 1696.39 0.5 0 0 +1.25 0.5 -0.136723 0.283646 1696.87 0.5 0 0 + + +1.25 0.375 -0.147055 0.218878 2121.84 0.625 0 0 +1.5 0.375 -0.223447 0.0186818 2122.72 0.625 0 0 + +1.25 0.5 -0.136723 0.283646 1696.87 0.5 0 0 +1.5 0.5 -0.191257 0.0225651 1698.13 0.5 0 0 + + +1.5 0.25 -0.228033 0.0129817 2547.3 0.75 0 0 +1.75 0.25 -0.1687 -0.125625 2547.97 0.75 0 0 + +1.5 0.375 -0.223447 0.0186818 2122.72 0.625 0 0 +1.75 0.375 -0.168056 -0.187706 2123.64 0.625 0 0 + + +1.75 0.25 -0.1687 -0.125625 2547.97 0.75 0 0 +2 0.25 0 -0.187279 2548.27 0.75 0 0 + +1.75 0.375 -0.168056 -0.187706 2123.64 0.625 0 0 +2 0.375 0 -0.284692 2124.06 0.625 0 0 + + +1.5 0.375 -0.223447 0.0186818 2122.72 0.625 0 0 +1.75 0.375 -0.168056 -0.187706 2123.64 0.625 0 0 + +1.5 0.5 -0.191257 0.0225651 1698.13 0.5 0 0 +1.75 0.5 -0.151898 -0.241915 1699.41 0.5 0 0 + + +1.75 0.375 -0.168056 -0.187706 2123.64 0.625 0 0 +2 0.375 0 -0.284692 2124.06 0.625 0 0 + +1.75 0.5 -0.151898 -0.241915 1699.41 0.5 0 0 +2 0.5 0 -0.38262 1700.09 0.5 0 0 + + +0 0.5 0 -0.396595 1700.36 0.5 0 0 +0.25 0.5 0.147628 -0.262746 1699.69 0.5 0 0 + +0 0.625 0 -0.486039 1276.77 0.375 0 0 +0.25 0.625 0.088844 -0.286361 1275.58 0.375 0 0 + + +0.25 0.5 0.147628 -0.262746 1699.69 0.5 0 0 +0.5 0.5 0.187137 -0.0225651 1698.47 0.5 0 0 + +0.25 0.625 0.088844 -0.286361 1275.58 0.375 0 0 +0.5 0.625 0.0697766 -0.0223763 1273.92 0.375 0 0 + + +0 0.625 0 -0.486039 1276.77 0.375 0 0 +0.25 0.625 0.088844 -0.286361 1275.58 0.375 0 0 + +0 0.75 0 -0.552448 854.151 0.25 0 0.142857 +0.25 0.75 -0.0453658 -0.24535 851.377 0.25 0 0.142857 + + +0.25 0.625 0.088844 -0.286361 1275.58 0.375 0 0 +0.5 0.625 0.0697766 -0.0223763 1273.92 0.375 0 0 + +0.25 0.75 -0.0453658 -0.24535 851.377 0.25 0 0.142857 +0.5 0.75 -0.15993 -0.0149564 849.32 0.25 0 0.142857 + + +0.5 0.5 0.187137 -0.0225651 1698.47 0.5 0 0 +0.75 0.5 0.162803 0.221016 1697.24 0.5 0 0 + +0.5 0.625 0.0697766 -0.0223763 1273.92 0.375 0 0 +0.75 0.625 0.0916352 0.233515 1272.34 0.375 0 0 + + +0.75 0.5 0.162803 0.221016 1697.24 0.5 0 0 +1 0.5 0.0283452 0.389607 1696.39 0.5 0 0 + +0.75 0.625 0.0916352 0.233515 1272.34 0.375 0 0 +1 0.625 0.0183469 0.480676 1270.83 0.375 0 0 + + +0.5 0.625 0.0697766 -0.0223763 1273.92 0.375 0 0 +0.75 0.625 0.0916352 0.233515 1272.34 0.375 0 0 + +0.5 0.75 -0.15993 -0.0149564 849.32 0.25 0 0.142857 +0.75 0.75 -0.0643719 0.182271 847.665 0.25 0 0.142857 + + +0.75 0.625 0.0916352 0.233515 1272.34 0.375 0 0 +1 0.625 0.0183469 0.480676 1270.83 0.375 0 0 + +0.75 0.75 -0.0643719 0.182271 847.665 0.25 0 0.142857 +1 0.75 -0.00742234 0.549509 844.294 0.25 0 0.142857 + + +0 0.75 0 -0.552448 854.151 0.25 0 0.142857 +0.25 0.75 -0.0453658 -0.24535 851.377 0.25 0 0.142857 + +0 0.875 0 -0.542292 433.474 0.125 0 1 +0.25 0.875 -0.354591 -0.0954565 426.662 0.125 0 1 + + +0.25 0.75 -0.0453658 -0.24535 851.377 0.25 0 0.142857 +0.5 0.75 -0.15993 -0.0149564 849.32 0.25 0 0.142857 + +0.25 0.875 -0.354591 -0.0954565 426.662 0.125 0 1 +0.5 0.875 -0.534008 -0.00327196 424.489 0.125 0 1 + + +0 0.875 0 -0.542292 433.474 0.125 0 1 +0.25 0.875 -0.354591 -0.0954565 426.662 0.125 0 1 + +0 1 0 0 21.0699 0 0 1 +0.25 1 -1 0 -2.52827 0 0 1 + + +0.25 0.875 -0.354591 -0.0954565 426.662 0.125 0 1 +0.5 0.875 -0.534008 -0.00327196 424.489 0.125 0 1 + +0.25 1 -1 0 -2.52827 0 0 1 +0.5 1 -1 0 -0.645742 0 0 1 + + +0.5 0.75 -0.15993 -0.0149564 849.32 0.25 0 0.142857 +0.75 0.75 -0.0643719 0.182271 847.665 0.25 0 0.142857 + +0.5 0.875 -0.534008 -0.00327196 424.489 0.125 0 1 +0.75 0.875 -0.393632 0.0404765 423.936 0.125 0 1 + + +0.75 0.75 -0.0643719 0.182271 847.665 0.25 0 0.142857 +1 0.75 -0.00742234 0.549509 844.294 0.25 0 0.142857 + +0.75 0.875 -0.393632 0.0404765 423.936 0.125 0 1 +1 0.875 -0.114817 0.542012 415.782 0.125 0 1 + + +0.5 0.875 -0.534008 -0.00327196 424.489 0.125 0 1 +0.75 0.875 -0.393632 0.0404765 423.936 0.125 0 1 + +0.5 1 -1 0 -0.645742 0 0 1 +0.75 1 -1 0 5.20377 0 0 1 + + +0.75 0.875 -0.393632 0.0404765 423.936 0.125 0 1 +1 0.875 -0.114817 0.542012 415.782 0.125 0 1 + +0.75 1 -1 0 5.20377 0 0 1 +1 1 -1 0 -21.0889 0 0 1 + + +1 0.5 0.0283452 0.389607 1696.39 0.5 0 0 +1.25 0.5 -0.136723 0.283646 1696.87 0.5 0 0 + +1 0.625 0.0183469 0.480676 1270.83 0.375 0 0 +1.25 0.625 -0.0958487 0.322758 1271.71 0.375 0 0 + + +1.25 0.5 -0.136723 0.283646 1696.87 0.5 0 0 +1.5 0.5 -0.191257 0.0225651 1698.13 0.5 0 0 + +1.25 0.625 -0.0958487 0.322758 1271.71 0.375 0 0 +1.5 0.625 -0.0961879 0.0223763 1273.55 0.375 0 0 + + +1 0.625 0.0183469 0.480676 1270.83 0.375 0 0 +1.25 0.625 -0.0958487 0.322758 1271.71 0.375 0 0 + +1 0.75 -0.00742234 0.549509 844.294 0.25 0 0.142857 +1.25 0.75 0.0163829 0.298966 846.489 0.25 0 0.142857 + + +1.25 0.625 -0.0958487 0.322758 1271.71 0.375 0 0 +1.5 0.625 -0.0961879 0.0223763 1273.55 0.375 0 0 + +1.25 0.75 0.0163829 0.298966 846.489 0.25 0 0.142857 +1.5 0.75 0.114901 0.0149564 849.004 0.25 0 0.142857 + + +1.5 0.5 -0.191257 0.0225651 1698.13 0.5 0 0 +1.75 0.5 -0.151898 -0.241915 1699.41 0.5 0 0 + +1.5 0.625 -0.0961879 0.0223763 1273.55 0.375 0 0 +1.75 0.625 -0.0986399 -0.269912 1275.31 0.375 0 0 + + +1.75 0.5 -0.151898 -0.241915 1699.41 0.5 0 0 +2 0.5 0 -0.38262 1700.09 0.5 0 0 + +1.75 0.625 -0.0986399 -0.269912 1275.31 0.375 0 0 +2 0.625 0 -0.475313 1276.51 0.375 0 0 + + +1.5 0.625 -0.0961879 0.0223763 1273.55 0.375 0 0 +1.75 0.625 -0.0986399 -0.269912 1275.31 0.375 0 0 + +1.5 0.75 0.114901 0.0149564 849.004 0.25 0 0.142857 +1.75 0.75 0.035389 -0.235888 851.117 0.25 0 0.142857 + + +1.75 0.625 -0.0986399 -0.269912 1275.31 0.375 0 0 +2 0.625 0 -0.475313 1276.51 0.375 0 0 + +1.75 0.75 0.035389 -0.235888 851.117 0.25 0 0.142857 +2 0.75 0 -0.546571 853.909 0.25 0 0.142857 + + +1 0.75 -0.00742234 0.549509 844.294 0.25 0 0.142857 +1.25 0.75 0.0163829 0.298966 846.489 0.25 0 0.142857 + +1 0.875 -0.114817 0.542012 415.782 0.125 0 1 +1.25 0.875 0.311707 0.145132 421.401 0.125 0 1 + + +1.25 0.75 0.0163829 0.298966 846.489 0.25 0 0.142857 +1.5 0.75 0.114901 0.0149564 849.004 0.25 0 0.142857 + +1.25 0.875 0.311707 0.145132 421.401 0.125 0 1 +1.5 0.875 0.485501 0.00327196 424.679 0.125 0 1 + + +1 0.875 -0.114817 0.542012 415.782 0.125 0 1 +1.25 0.875 0.311707 0.145132 421.401 0.125 0 1 + +1 1 -1 0 -21.0889 0 0 1 +1.25 1 1 0 0.472162 0 0 1 + + +1.25 0.875 0.311707 0.145132 421.401 0.125 0 1 +1.5 0.875 0.485501 0.00327196 424.679 0.125 0 1 + +1.25 1 1 0 0.472162 0 0 1 +1.5 1 1 0 0.645742 0 0 1 + + +1.5 0.75 0.114901 0.0149564 849.004 0.25 0 0.142857 +1.75 0.75 0.035389 -0.235888 851.117 0.25 0 0.142857 + +1.5 0.875 0.485501 0.00327196 424.679 0.125 0 1 +1.75 0.875 0.350748 -0.0901517 426.336 0.125 0 1 + + +1.75 0.75 0.035389 -0.235888 851.117 0.25 0 0.142857 +2 0.75 0 -0.546571 853.909 0.25 0 0.142857 + +1.75 0.875 0.350748 -0.0901517 426.336 0.125 0 1 +2 0.875 0 -0.541731 433.296 0.125 0 1 + + +1.5 0.875 0.485501 0.00327196 424.679 0.125 0 1 +1.75 0.875 0.350748 -0.0901517 426.336 0.125 0 1 + +1.5 1 1 0 0.645742 0 0 1 +1.75 1 1 0 -3.14767 0 0 1 + + +1.75 0.875 0.350748 -0.0901517 426.336 0.125 0 1 +2 0.875 0 -0.541731 433.296 0.125 0 1 + +1.75 1 1 0 -3.14767 0 0 1 +2 1 0 0 21.1079 0 0 1 + + diff --git a/tests/prescribed_solution_composition_initial_condition_particles/solution/solution-00001.0000.gnuplot b/tests/prescribed_solution_composition_initial_condition_particles/solution/solution-00001.0000.gnuplot new file mode 100644 index 00000000000..3b34a31e853 --- /dev/null +++ b/tests/prescribed_solution_composition_initial_condition_particles/solution/solution-00001.0000.gnuplot @@ -0,0 +1,455 @@ +# This file was generated by the deal.II library. + + +# +# For a description of the GNUPLOT format see the GNUPLOT manual. +# +#

+0 0 0 0 3397.85 1 1 0 +0.25 0 0.199038 0 3397.62 1 1 0 + +0 0.125 0 -0.107827 2973.36 0.856873 1 0 +0.25 0.125 0.19717 -0.0832456 2973.12 0.859392 1 0 + + +0.25 0 0.199038 0 3397.62 1 1 0 +0.5 0 0.313234 0 3397.06 1 1 0 + +0.25 0.125 0.19717 -0.0832456 2973.12 0.859392 1 0 +0.5 0.125 0.308523 -0.0273944 2972.52 0.865453 0.944588 0 + + +0 0.125 0 -0.107827 2973.36 0.856873 1 0 +0.25 0.125 0.19717 -0.0832456 2973.12 0.859392 1 0 + +0 0.25 0 -0.214352 2549 0.727018 0 0 +0.25 0.25 0.190329 -0.163392 2548.7 0.733351 0 0 + + +0.25 0.125 0.19717 -0.0832456 2973.12 0.859392 1 0 +0.5 0.125 0.308523 -0.0273944 2972.52 0.865453 0.944588 0 + +0.25 0.25 0.190329 -0.163392 2548.7 0.733351 0 0 +0.5 0.25 0.291532 -0.0525803 2548.01 0.74759 0 0 + + +0.5 0 0.313234 0 3397.06 1 1 0 +0.75 0 0.306111 0 3396.39 1 1 0 + +0.5 0.125 0.308523 -0.0273944 2972.52 0.865453 0.944588 0 +0.75 0.125 0.30133 0.0340602 2971.83 0.87263 1 0 + + +0.75 0 0.306111 0 3396.39 1 1 0 +1 0 0.184922 0 3395.87 1 1 0 + +0.75 0.125 0.30133 0.0340602 2971.83 0.87263 1 0 +1 0.125 0.181689 0.0823554 2971.28 0.876823 0.977343 0 + + +0.5 0.125 0.308523 -0.0273944 2972.52 0.865453 0.944588 0 +0.75 0.125 0.30133 0.0340602 2971.83 0.87263 1 0 + +0.5 0.25 0.291532 -0.0525803 2548.01 0.74759 0 0 +0.75 0.25 0.284474 0.0661604 2547.22 0.760874 0 0 + + +0.75 0.125 0.30133 0.0340602 2971.83 0.87263 1 0 +1 0.125 0.181689 0.0823554 2971.28 0.876823 0.977343 0 + +0.75 0.25 0.284474 0.0661604 2547.22 0.760874 0 0 +1 0.25 0.171571 0.163922 2546.58 0.766825 0 0 + + +0 0.25 0 -0.214352 2549 0.727018 0 0 +0.25 0.25 0.190329 -0.163392 2548.7 0.733351 0 0 + +0 0.375 0 -0.317586 2124.79 0.591315 0 0 +0.25 0.375 0.174516 -0.235405 2124.38 0.600962 0 0 + + +0.25 0.25 0.190329 -0.163392 2548.7 0.733351 0 0 +0.5 0.25 0.291532 -0.0525803 2548.01 0.74759 0 0 + +0.25 0.375 0.174516 -0.235405 2124.38 0.600962 0 0 +0.5 0.375 0.253119 -0.0723694 2123.51 0.621173 0 0 + + +0 0.375 0 -0.317586 2124.79 0.591315 0 0 +0.25 0.375 0.174516 -0.235405 2124.38 0.600962 0 0 + +0 0.5 0 -0.414035 1700.81 0.456465 0 0 +0.25 0.5 0.141765 -0.28986 1700.16 0.470657 0 0 + + +0.25 0.375 0.174516 -0.235405 2124.38 0.600962 0 0 +0.5 0.375 0.253119 -0.0723694 2123.51 0.621173 0 0 + +0.25 0.5 0.141765 -0.28986 1700.16 0.470657 0 0 +0.5 0.5 0.176326 -0.0816306 1699.01 0.496078 0 0 + + +0.5 0.25 0.291532 -0.0525803 2548.01 0.74759 0 0 +0.75 0.25 0.284474 0.0661604 2547.22 0.760874 0 0 + +0.5 0.375 0.253119 -0.0723694 2123.51 0.621173 0 0 +0.75 0.375 0.247399 0.0927002 2122.58 0.640247 0 0 + + +0.75 0.25 0.284474 0.0661604 2547.22 0.760874 0 0 +1 0.25 0.171571 0.163922 2546.58 0.766825 0 0 + +0.75 0.375 0.247399 0.0927002 2122.58 0.640247 0 0 +1 0.375 0.152997 0.241665 2121.75 0.650531 0 0 + + +0.5 0.375 0.253119 -0.0723694 2123.51 0.621173 0 0 +0.75 0.375 0.247399 0.0927002 2122.58 0.640247 0 0 + +0.5 0.5 0.176326 -0.0816306 1699.01 0.496078 0 0 +0.75 0.5 0.17457 0.106891 1697.92 0.519557 0 0 + + +0.75 0.375 0.247399 0.0927002 2122.58 0.640247 0 0 +1 0.375 0.152997 0.241665 2121.75 0.650531 0 0 + +0.75 0.5 0.17457 0.106891 1697.92 0.519557 0 0 +1 0.5 0.121835 0.307562 1696.77 0.535119 0 0 + + +1 0 0.184922 0 3395.87 1 1 0 +1.25 0 0.00644438 0 3395.72 1 1 0 + +1 0.125 0.181689 0.0823554 2971.28 0.876823 0.977343 0 +1.25 0.125 0.00190887 0.0872069 2971.12 0.874561 1 0 + + +1.25 0 0.00644438 0 3395.72 1 1 0 +1.5 0 -0.122484 0 3395.97 1 1 0 + +1.25 0.125 0.00190887 0.0872069 2971.12 0.874561 1 0 +1.5 0.125 -0.129013 0.0348194 2971.4 0.867622 1 0 + + +1 0.125 0.181689 0.0823554 2971.28 0.876823 0.977343 0 +1.25 0.125 0.00190887 0.0872069 2971.12 0.874561 1 0 + +1 0.25 0.171571 0.163922 2546.58 0.766825 0 0 +1.25 0.25 -0.0106018 0.176875 2546.39 0.762626 0 0 + + +1.25 0.125 0.00190887 0.0872069 2971.12 0.874561 1 0 +1.5 0.125 -0.129013 0.0348194 2971.4 0.867622 1 0 + +1.25 0.25 -0.0106018 0.176875 2546.39 0.762626 0 0 +1.5 0.25 -0.146619 0.0702378 2546.78 0.749739 0 0 + + +1.5 0 -0.122484 0 3395.97 1 1 0 +1.75 0 -0.118205 0 3396.36 1 1 0 + +1.5 0.125 -0.129013 0.0348194 2971.4 0.867622 1 0 +1.75 0.125 -0.123298 -0.0379676 2971.84 0.860955 0.9375 0 + + +1.75 0 -0.118205 0 3396.36 1 1 0 +2 0 0 0 3396.55 1 1 0 + +1.75 0.125 -0.123298 -0.0379676 2971.84 0.860955 0.9375 0 +2 0.125 0 -0.0718522 2972.05 0.85815 0.927235 0 + + +1.5 0.125 -0.129013 0.0348194 2971.4 0.867622 1 0 +1.75 0.125 -0.123298 -0.0379676 2971.84 0.860955 0.9375 0 + +1.5 0.25 -0.146619 0.0702378 2546.78 0.749739 0 0 +1.75 0.25 -0.136995 -0.0793218 2547.38 0.734826 0 0 + + +1.75 0.125 -0.123298 -0.0379676 2971.84 0.860955 0.9375 0 +2 0.125 0 -0.0718522 2972.05 0.85815 0.927235 0 + +1.75 0.25 -0.136995 -0.0793218 2547.38 0.734826 0 0 +2 0.25 0 -0.149447 2547.67 0.728174 0 0 + + +1 0.25 0.171571 0.163922 2546.58 0.766825 0 0 +1.25 0.25 -0.0106018 0.176875 2546.39 0.762626 0 0 + +1 0.375 0.152997 0.241665 2121.75 0.650531 0 0 +1.25 0.375 -0.0276367 0.270472 2121.5 0.643837 0 0 + + +1.25 0.25 -0.0106018 0.176875 2546.39 0.762626 0 0 +1.5 0.25 -0.146619 0.0702378 2546.78 0.749739 0 0 + +1.25 0.375 -0.0276367 0.270472 2121.5 0.643837 0 0 +1.5 0.375 -0.168316 0.106059 2122.09 0.624547 0 0 + + +1 0.375 0.152997 0.241665 2121.75 0.650531 0 0 +1.25 0.375 -0.0276367 0.270472 2121.5 0.643837 0 0 + +1 0.5 0.121835 0.307562 1696.77 0.535119 0 0 +1.25 0.5 -0.0430006 0.36691 1696.36 0.52535 0 0 + + +1.25 0.375 -0.0276367 0.270472 2121.5 0.643837 0 0 +1.5 0.375 -0.168316 0.106059 2122.09 0.624547 0 0 + +1.25 0.5 -0.0430006 0.36691 1696.36 0.52535 0 0 +1.5 0.5 -0.178527 0.139427 1697.33 0.500217 0 0 + + +1.5 0.25 -0.146619 0.0702378 2546.78 0.749739 0 0 +1.75 0.25 -0.136995 -0.0793218 2547.38 0.734826 0 0 + +1.5 0.375 -0.168316 0.106059 2122.09 0.624547 0 0 +1.75 0.375 -0.153584 -0.125877 2123.01 0.602917 0 0 + + +1.75 0.25 -0.136995 -0.0793218 2547.38 0.734826 0 0 +2 0.25 0 -0.149447 2547.67 0.728174 0 0 + +1.75 0.375 -0.153584 -0.125877 2123.01 0.602917 0 0 +2 0.375 0 -0.236905 2123.44 0.592756 0 0 + + +1.5 0.375 -0.168316 0.106059 2122.09 0.624547 0 0 +1.75 0.375 -0.153584 -0.125877 2123.01 0.602917 0 0 + +1.5 0.5 -0.178527 0.139427 1697.33 0.500217 0 0 +1.75 0.5 -0.159942 -0.175104 1698.74 0.472517 0 0 + + +1.75 0.375 -0.153584 -0.125877 2123.01 0.602917 0 0 +2 0.375 0 -0.236905 2123.44 0.592756 0 0 + +1.75 0.5 -0.159942 -0.175104 1698.74 0.472517 0 0 +2 0.5 0 -0.334358 1699.46 0.45776 0 0 + + +0 0.5 0 -0.414035 1700.81 0.456465 0 0 +0.25 0.5 0.141765 -0.28986 1700.16 0.470657 0 0 + +0 0.625 0 -0.498847 1277.21 0.323135 0 0 +0.25 0.625 0.0760705 -0.308103 1276.06 0.344419 0 0 + + +0.25 0.5 0.141765 -0.28986 1700.16 0.470657 0 0 +0.5 0.5 0.176326 -0.0816306 1699.01 0.496078 0 0 + +0.25 0.625 0.0760705 -0.308103 1276.06 0.344419 0 0 +0.5 0.625 0.0349932 -0.0740824 1274.45 0.371728 0 0 + + +0 0.625 0 -0.498847 1277.21 0.323135 0 0 +0.25 0.625 0.0760705 -0.308103 1276.06 0.344419 0 0 + +0 0.75 0 -0.559048 854.552 0.193413 0 0.180478 +0.25 0.75 -0.0602112 -0.258061 851.841 0.223728 0 0.19805 + + +0.25 0.625 0.0760705 -0.308103 1276.06 0.344419 0 0 +0.5 0.625 0.0349932 -0.0740824 1274.45 0.371728 0 0 + +0.25 0.75 -0.0602112 -0.258061 851.841 0.223728 0 0.19805 +0.5 0.75 -0.203285 -0.0476832 849.786 0.243145 0 0.164067 + + +0.5 0.5 0.176326 -0.0816306 1699.01 0.496078 0 0 +0.75 0.5 0.17457 0.106891 1697.92 0.519557 0 0 + +0.5 0.625 0.0349932 -0.0740824 1274.45 0.371728 0 0 +0.75 0.625 0.040814 0.0993703 1273.31 0.396331 0 0 + + +0.75 0.5 0.17457 0.106891 1697.92 0.519557 0 0 +1 0.5 0.121835 0.307562 1696.77 0.535119 0 0 + +0.75 0.625 0.040814 0.0993703 1273.31 0.396331 0 0 +1 0.625 0.0666743 0.342908 1271.59 0.419277 0 0 + + +0.5 0.625 0.0349932 -0.0740824 1274.45 0.371728 0 0 +0.75 0.625 0.040814 0.0993703 1273.31 0.396331 0 0 + +0.5 0.75 -0.203285 -0.0476832 849.786 0.243145 0 0.164067 +0.75 0.75 -0.18788 0.064798 848.85 0.26676 0 0.226453 + + +0.75 0.625 0.040814 0.0993703 1273.31 0.396331 0 0 +1 0.625 0.0666743 0.342908 1271.59 0.419277 0 0 + +0.75 0.75 -0.18788 0.064798 848.85 0.26676 0 0.226453 +1 0.75 -0.0473949 0.311424 846.367 0.301418 0 2.77556e-17 + + +0 0.75 0 -0.559048 854.552 0.193413 0 0.180478 +0.25 0.75 -0.0602112 -0.258061 851.841 0.223728 0 0.19805 + +0 0.875 0 -0.545909 433.826 0.0838035 0 1 +0.25 0.875 -0.365977 -0.0957433 427.102 0.132366 0 1 + + +0.25 0.75 -0.0602112 -0.258061 851.841 0.223728 0 0.19805 +0.5 0.75 -0.203285 -0.0476832 849.786 0.243145 0 0.164067 + +0.25 0.875 -0.365977 -0.0957433 427.102 0.132366 0 1 +0.5 0.875 -0.557665 -0.0232964 424.851 0.138601 0 1 + + +0 0.875 0 -0.545909 433.826 0.0838035 0 1 +0.25 0.875 -0.365977 -0.0957433 427.102 0.132366 0 1 + +0 1 0 0 21.7435 0 0 1 +0.25 1 -1 0 -2.76748 0 0 1 + + +0.25 0.875 -0.365977 -0.0957433 427.102 0.132366 0 1 +0.5 0.875 -0.557665 -0.0232964 424.851 0.138601 0 1 + +0.25 1 -1 0 -2.76748 0 0 1 +0.5 1 -1 0 1.4437 0 0 1 + + +0.5 0.75 -0.203285 -0.0476832 849.786 0.243145 0 0.164067 +0.75 0.75 -0.18788 0.064798 848.85 0.26676 0 0.226453 + +0.5 0.875 -0.557665 -0.0232964 424.851 0.138601 0 1 +0.75 0.875 -0.536303 0.0286538 424.839 0.148707 0 1 + + +0.75 0.75 -0.18788 0.064798 848.85 0.26676 0 0.226453 +1 0.75 -0.0473949 0.311424 846.367 0.301418 0 2.77556e-17 + +0.75 0.875 -0.536303 0.0286538 424.839 0.148707 0 1 +1 0.875 -0.331303 0.145945 421.246 0.183053 0 1 + + +0.5 0.875 -0.557665 -0.0232964 424.851 0.138601 0 1 +0.75 0.875 -0.536303 0.0286538 424.839 0.148707 0 1 + +0.5 1 -1 0 1.4437 0 0 1 +0.75 1 -1 0 -0.538797 0 0 1 + + +0.75 0.875 -0.536303 0.0286538 424.839 0.148707 0 1 +1 0.875 -0.331303 0.145945 421.246 0.183053 0 1 + +0.75 1 -1 0 -0.538797 0 0 1 +1 1 -1 0 0.986867 0 0 1 + + +1 0.5 0.121835 0.307562 1696.77 0.535119 0 0 +1.25 0.5 -0.0430006 0.36691 1696.36 0.52535 0 0 + +1 0.625 0.0666743 0.342908 1271.59 0.419277 0 0 +1.25 0.625 -0.0471352 0.461399 1270.82 0.404839 0 0 + + +1.25 0.5 -0.0430006 0.36691 1696.36 0.52535 0 0 +1.5 0.5 -0.178527 0.139427 1697.33 0.500217 0 0 + +1.25 0.625 -0.0471352 0.461399 1270.82 0.404839 0 0 +1.5 0.625 -0.145424 0.158507 1272.45 0.37569 0 0 + + +1 0.625 0.0666743 0.342908 1271.59 0.419277 0 0 +1.25 0.625 -0.0471352 0.461399 1270.82 0.404839 0 0 + +1 0.75 -0.0473949 0.311424 846.367 0.301418 0 2.77556e-17 +1.25 0.75 -0.0232171 0.5375 844.3 0.279567 0 5.55112e-17 + + +1.25 0.625 -0.0471352 0.461399 1270.82 0.404839 0 0 +1.5 0.625 -0.145424 0.158507 1272.45 0.37569 0 0 + +1.25 0.75 -0.0232171 0.5375 844.3 0.279567 0 5.55112e-17 +1.5 0.75 -0.00704215 0.133468 847.699 0.247935 0 0.262413 + + +1.5 0.5 -0.178527 0.139427 1697.33 0.500217 0 0 +1.75 0.5 -0.159942 -0.175104 1698.74 0.472517 0 0 + +1.5 0.625 -0.145424 0.158507 1272.45 0.37569 0 0 +1.75 0.625 -0.128509 -0.212487 1274.66 0.3457 0 0 + + +1.75 0.5 -0.159942 -0.175104 1698.74 0.472517 0 0 +2 0.5 0 -0.334358 1699.46 0.45776 0 0 + +1.75 0.625 -0.128509 -0.212487 1274.66 0.3457 0 0 +2 0.625 0 -0.436174 1275.87 0.324005 0 0 + + +1.5 0.625 -0.145424 0.158507 1272.45 0.37569 0 0 +1.75 0.625 -0.128509 -0.212487 1274.66 0.3457 0 0 + +1.5 0.75 -0.00704215 0.133468 847.699 0.247935 0 0.262413 +1.75 0.75 -0.00315912 -0.200201 850.535 0.224669 0 0.792906 + + +1.75 0.625 -0.128509 -0.212487 1274.66 0.3457 0 0 +2 0.625 0 -0.436174 1275.87 0.324005 0 0 + +1.75 0.75 -0.00315912 -0.200201 850.535 0.224669 0 0.792906 +2 0.75 0 -0.523452 853.278 0.193818 0 0.866006 + + +1 0.75 -0.0473949 0.311424 846.367 0.301418 0 2.77556e-17 +1.25 0.75 -0.0232171 0.5375 844.3 0.279567 0 5.55112e-17 + +1 0.875 -0.331303 0.145945 421.246 0.183053 0 1 +1.25 0.875 0.0954015 0.541379 415.826 0.159831 0 1 + + +1.25 0.75 -0.0232171 0.5375 844.3 0.279567 0 5.55112e-17 +1.5 0.75 -0.00704215 0.133468 847.699 0.247935 0 0.262413 + +1.25 0.875 0.0954015 0.541379 415.826 0.159831 0 1 +1.5 0.875 0.34388 0.0151214 423.659 0.139313 0 1 + + +1 0.875 -0.331303 0.145945 421.246 0.183053 0 1 +1.25 0.875 0.0954015 0.541379 415.826 0.159831 0 1 + +1 1 -1 0 0.986867 0 0 1 +1.25 1 1 0 -21.7192 0 0 1 + + +1.25 0.875 0.0954015 0.541379 415.826 0.159831 0 1 +1.5 0.875 0.34388 0.0151214 423.659 0.139313 0 1 + +1.25 1 1 0 -21.7192 0 0 1 +1.5 1 1 0 6.32223 0 0 1 + + +1.5 0.75 -0.00704215 0.133468 847.699 0.247935 0 0.262413 +1.75 0.75 -0.00315912 -0.200201 850.535 0.224669 0 0.792906 + +1.5 0.875 0.34388 0.0151214 423.659 0.139313 0 1 +1.75 0.875 0.3294 -0.0684911 425.877 0.132141 0 1 + + +1.75 0.75 -0.00315912 -0.200201 850.535 0.224669 0 0.792906 +2 0.75 0 -0.523452 853.278 0.193818 0 0.866006 + +1.75 0.875 0.3294 -0.0684911 425.877 0.132141 0 1 +2 0.875 0 -0.541208 432.707 0.0835731 0 1 + + +1.5 0.875 0.34388 0.0151214 423.659 0.139313 0 1 +1.75 0.875 0.3294 -0.0684911 425.877 0.132141 0 1 + +1.5 1 1 0 6.32223 0 0 1 +1.75 1 1 0 -5.47853 0 0 1 + + +1.75 0.875 0.3294 -0.0684911 425.877 0.132141 0 1 +2 0.875 0 -0.541208 432.707 0.0835731 0 1 + +1.75 1 1 0 -5.47853 0 0 1 +2 1 0 0 21.7589 0 0 1 + +