betaPsi = -N * 2 * Phi / (1 + PhiSq)
A = alpha**2 * c * (1.+ c * beta) * betaPsi
D = alpha**2 * (1.+ c * beta)**2
dxi = phase.getFaceGrad()._take((1, 0), axis = 1) * (-1, 1)
anisotropySource = (A * dxi).getDivergence()
from fipy.terms.transientTerm import TransientTerm
from fipy.terms.explicitDiffusionTerm import ExplicitDiffusionTerm
from fipy.terms.implicitSourceTerm import ImplicitSourceTerm
phaseEq = TransientTerm(tau) == ExplicitDiffusionTerm(D) + \
ImplicitSourceTerm(mVar * ((mVar < 0) - phase)) + \
((mVar > 0.) * mVar * phase + anisotropySource)
from fipy.terms.implicitDiffusionTerm import ImplicitDiffusionTerm
temperatureEq = TransientTerm() == \
ImplicitDiffusionTerm(tempDiffusionCoeff) + \
(phase - phase.getOld()) / timeStepDuration
bench.stop('terms')
phase.updateOld()
temperature.updateOld()
phaseEq.solve(phase, dt=timeStepDuration)
temperatureEq.solve(temperature, dt=timeStepDuration)
steps = 10
bench.start()
for i in range(steps):
phase.updateOld()
temperature.updateOld()
from fipy.meshes.grid1D import Grid1D
mesh = Grid1D(nx = nx, dx = dx)
from fipy.variables.cellVariable import CellVariable
u = CellVariable(name="solution variable",
mesh=mesh,
value=0,hasOld=1)
from fipy.terms.implicitDiffusionTerm import ImplicitDiffusionTerm
from fipy.terms.transientTerm import TransientTerm
a=1
timeStepDuration = 0.01
steps=1000
eq= ImplicitDiffusionTerm(coeff=(a,a))+(u.getOld().getOld() - 2*u.getOld() + u) / timeStepDuration ** 2
from fipy.boundaryConditions.fixedValue import FixedValue
from fipy.boundaryConditions.fixedFlux import FixedFlux
from fipy.boundaryConditions.nthOrderBoundaryCondition import NthOrderBoundaryCondition
BCs = ( NthOrderBoundaryCondition(faces=mesh.getFacesRight(), value=0,order=2),
NthOrderBoundaryCondition(faces=mesh.getFacesRight(), value=1,order=3),
FixedValue(faces=mesh.getFacesLeft(), value=0),
FixedFlux(faces=mesh.getFacesLeft(), value=0))
from fipy import viewers
viewer = viewers.make(vars=(u))
betaPsi = -N * 2 * Phi / (1 + PhiSq)
A = alpha**2 * c * (1. + c * beta) * betaPsi
D = alpha**2 * (1. + c * beta)**2
dxi = phase.getFaceGrad()._take((1, 0), axis=1) * (-1, 1)
anisotropySource = (A * dxi).getDivergence()
from fipy.terms.transientTerm import TransientTerm
from fipy.terms.explicitDiffusionTerm import ExplicitDiffusionTerm
from fipy.terms.implicitSourceTerm import ImplicitSourceTerm
phaseEq = TransientTerm(tau) == ExplicitDiffusionTerm(D) + \
ImplicitSourceTerm(mVar * ((mVar < 0) - phase)) + \
((mVar > 0.) * mVar * phase + anisotropySource)
from fipy.terms.implicitDiffusionTerm import ImplicitDiffusionTerm
temperatureEq = TransientTerm() == \
ImplicitDiffusionTerm(tempDiffusionCoeff) + \
(phase - phase.getOld()) / timeStepDuration
bench.stop('terms')
phase.updateOld()
temperature.updateOld()
phaseEq.solve(phase, dt=timeStepDuration)
temperatureEq.solve(temperature, dt=timeStepDuration)
steps = 10
bench.start()
for i in range(steps):
phase.updateOld()
temperature.updateOld()
dx = 1.
from fipy.meshes.grid1D import Grid1D
mesh = Grid1D(nx=nx, dx=dx)
from fipy.variables.cellVariable import CellVariable
u = CellVariable(name="solution variable", mesh=mesh, value=0, hasOld=1)
from fipy.terms.implicitDiffusionTerm import ImplicitDiffusionTerm
from fipy.terms.transientTerm import TransientTerm
a = 1
timeStepDuration = 0.01
steps = 1000
eq = ImplicitDiffusionTerm(coeff=(
a, a)) + (u.getOld().getOld() - 2 * u.getOld() + u) / timeStepDuration**2
from fipy.boundaryConditions.fixedValue import FixedValue
from fipy.boundaryConditions.fixedFlux import FixedFlux
from fipy.boundaryConditions.nthOrderBoundaryCondition import NthOrderBoundaryCondition
BCs = (NthOrderBoundaryCondition(faces=mesh.getFacesRight(), value=0, order=2),
NthOrderBoundaryCondition(faces=mesh.getFacesRight(), value=1, order=3),
FixedValue(faces=mesh.getFacesLeft(),
value=0), FixedFlux(faces=mesh.getFacesLeft(), value=0))
from fipy import viewers
viewer = viewers.make(vars=(u))
def hit_continue(Prompt='Hit any key to continue'):
