idaholab / pika

Phase-field snow microstructure

Pika

The snow micro-structure model, Pika, was developed following the work of Kaempfer and Plapp, 2009 and is a fully-coupled 3D finite element, phase field model capable of tracking the phase transition and capturing the heat and mass transfer at the micro-structure scale in the ice matrix and pore space.

pika's Issues

Add Temporal Scaling to Governing Equations

Due to the large differences in time scales of vapor diffusion and interface migration (~10^-4 and ~10^5 respectively) certain equations should be scaled by xi. This allows us to avoid numerical instabilities while the changes to the end result is negligible (Ref. Kaempfer and Plapp 2009)

Equations Affected:
u --> u / Xi
d'_o --> d'_o / Xi
Beta' --> Beta' / Xi
Lambda --> Xi * Lambda

Requirements for Xi:
10^ -5 < Xi < 1
Time_Diffusion / ( Xi * Interface_Time) << 1

Make test for criteria calculations

Stop using deprecated interface in AppFactory

Migrate to smart pointers in main.C

Time scaling needs its own interface

The time scaling stuff should be a separate interface from CoefficientKernelInterface because some object just need xi and not the other stuff.

Compute Interface Velocity

we need a way to compute the interface velocity

Add outputs to show if sharp interface assumption is violated

Phi = Phi is not informative

Come up with a name for the phase field variable to avoid confusion with which variable to use.

Prepare pika for moose/#6758

rename of HomogenizationHeatConduction to HomogenizedHeatConduction

Testing failing due to TestHarness fix

A second gold file was failing, this wan't being tested before the test harness was fixed in moose.

applyCoupledVar shows up in multiple actions, need to make a PikaAction base class

Clean-up Pika

Use this ticket for general clean-up and documentation at least until things are made public.

Transfer testing failing

phase_evolution/mms.test_double_well_potential does not work in parallel

Add tests for errors in PikaTimeDerivative

• Build RunException tests for the two possible errors in PikaTimeDeriviative
• Should be under tests/check_error directory

CoefficientKernelInterface should use a function pointer

Rather than have an if statement in coefficient on construction a function pointer should be populated that points to the proper call.

Refactor Actions to use smart pointers in MOOSE

heat_equation tests materials not coupling T

The material properties should be coupling in the actual temperature variable, right now they are all using the default value.

Expand Material Properties Testing

Add complete test for lambda, capillary length, interface kinetic coefficient.

How to build the library_path

When I running pika/problems/1D_full/1D_full.i by moose, the error always shows that

*** ERROR ***
A 'ChemicalPotentialIC' is not a registered object.

We loaded objects from the following libraries and still couldn't find your object:
(NONE)

Make sure you have compiled the library and either set the "library_path" variable in your input file or exported "MOOSE_LIBRARY_PATH".

But I can find the ChemicalPotentialIC.lo in pika/src/ics.

It's my first time to use moose, it confused me more than two weeks, I hope some one can help me to solve it.
How can I solve this problem? How to add the library_path to pika?
1D_full.txt

Move VariableGradientComponent to MOOSE

Phase Evolution eq. (33) MMS

Test implementation of the phase evolution equation using method of manufactured solutions.

Create a scaled version of HomogenizedHeatConduction

Make test for PikaInterfaceVelocity

Build convergence plots for MMS tests of each equation.

Temporal scaling test needs review

After fixing PikaTimeDerivative this test is wrong.

Update to ExecFlagEnum

Add in OffDiagonalJacobians to couple kernels

Currently the coupled kernels do not incorporate the off diagonal jacobian. These need to be added to avoid non converging solutions when the problem is fully coupled.

Kernels that need attention:
PikaTimeDerivative
PikaScaledTimeDerivative
PhaseTransition

Input files in the problems directory need to be updated and tested

The problems directory contains a bunch of input files, many of which do not run. These need to be cleaned up and then tests added to make sure they keep working.

ChemicalPotentialBC needs review

Currently: u - u_eq * (1-phi)/2

Should it be: (u - u_eq) * (1-phi)/2

Add moose module

Fixup up test output names that could clash

The new TestHarness patch has identified a few issues in this repo. Fixing before they actually case problems.

MooseEnum Error

When I tried to run the 2D_Bubble problem I get an error similar to ...
attempted to set std::vector MooseEnum "execute on" when it already exists as a MultiMooseEnum.

Update mms tests post PikaTimeDerivative changes

Changed PikaTimeDerivative to actually differentiate the couple variable. This caused a lot of changes with in the tests and they need to be updated to reflect the changes in the time derivative kernel.

Add mms testing for temporally scaled equations

AirProperties depends on ConstantProperties, but does nothing if you don't have it.

AirProperties requires the use of properties defined in ConstantProperties but doesn't check that they exist. This check needs to be implemented in the constructor.

Problem needs to be rescaled spatially

To resolve the interface thickness we run into a petsc error below

Computed Nan differencing parameter h!

Perhaps rescaling spatially will fix the problem.

PhaseFieldProperties needs some clean up

All the calls to getMaterialProperty need to be done in the initialization list

update error macro to use new moose function

PR coming...

phase_evolution/mms.test is slow or not solving

Code Verification of Heat Equation through MMS.

The heat equation from Eqn 34 in Kaempfer and Plapp (2009) needs to be tested using the method of manufacturable solutions.

Code Verification of Mass Transport Eq. using MMS

Use MMS on Eq. 35

Create Python implementation of materials

Develop SymPy calculations of the material properties and compare the results with Pika in a test.

Still using legacy Aux and UO stuff, fix it

Rename CoefficientTimeDerivative to PikaTimeDerviative

Make phi range from 0 to 1 instead of -1 to 1

Forcing Kernels made for MMS have unnecessary parameters

MassTransportSourceMMS and HeatEquationSourceMMS have CoupledVar which was intended for auxVariable phi but a different route was chosen and the code remained. It is neutral code, but it is misleading. MMS code and MMS input files will be the files that need removal of the unnecessary parameter.

heat_equation/mms.test exodiffing

Mass transport MMS test is slow

make it faster

Fix validParams

Fixup some invalid validParams

ACBulk is a template

Fix breakage due to moose/#6198