proteus  1.8.1
C/C++/Fortran libraries
transportCoefficients.h
Go to the documentation of this file.
1 #ifndef TRANSPORTCOEFFICIENTS_H
2 #define TRANSPORTCOEFFICIENTS_H
3 
4 #include <math.h>
5 #include <assert.h>
6 #include <stdio.h>
7 #include <memory.h>
8 #include <stdlib.h>
9 
15 extern void groundwaterTransportCoefficientsEvaluate(const int nPoints,
16  const int nSpace,
17  const double omega,
18  const double d,
19  const double alpha_L,
20  const double alpha_T,
21  const double *v,
22  const double *u,
23  double *m,
24  double *dm,
25  double *f,
26  double *df,
27  double *a);
28 
29 extern void TwophaseNavierStokes_ST_LS_SO_3D_Evaluate(const int nPoints,
30  const double eps_rho,
31  const double eps_mu,
32  const double sigma,
33  const double rho_0,
34  const double nu_0,
35  const double rho_1,
36  const double nu_1,
37  const double* g,
38  const double* phi,
39  const double* n,
40  const double* kappa,
41  const double *p,
42  const double *grad_p,
43  const double *u,
44  const double *v,
45  const double *w,
46  double *mom_u_acc,
47  double *dmom_u_acc_u,
48  double *mom_v_acc,
49  double *dmom_v_acc_v,
50  double *mom_w_acc,
51  double *dmom_w_acc_w,
52  double *mass_adv,
53  double *dmass_adv_u,
54  double *dmass_adv_v,
55  double *dmass_adv_w,
56  double *mom_u_adv,
57  double *dmom_u_adv_u,
58  double *dmom_u_adv_v,
59  double *dmom_u_adv_w,
60  double *mom_v_adv,
61  double *dmom_v_adv_u,
62  double *dmom_v_adv_v,
63  double *dmom_v_adv_w,
64  double *mom_w_adv,
65  double *dmom_w_adv_u,
66  double *dmom_w_adv_v,
67  double *dmom_w_adv_w,
68  double *mom_u_diff_ten,
69  double *mom_v_diff_ten,
70  double *mom_w_diff_ten,
71  double *mom_uv_diff_ten,
72  double *mom_uw_diff_ten,
73  double *mom_vu_diff_ten,
74  double *mom_vw_diff_ten,
75  double *mom_wu_diff_ten,
76  double *mom_wv_diff_ten,
77  double *mom_u_source,
78  double *mom_v_source,
79  double *mom_w_source,
80  double *mom_u_ham,
81  double *dmom_u_ham_grad_p,
82  double *mom_v_ham,
83  double *dmom_v_ham_grad_p,
84  double *mom_w_ham,
85  double *dmom_w_ham_grad_p);
86 extern void VOFCoefficientsEvaluate(int nPoints,
87  int nSpace,
88  double eps,
89  double* v,
90  double* phi,
91  double* u,
92  double* m,
93  double* dm,
94  double* f,
95  double* df);
96 extern void redistanceLevelSetSandFCoefficientsEvaluate(int nSimplex,
97  int nPointsPerSimplex,
98  int nSpace,
99  double eps,
100  double* u_levelSet,
101  double* dV,
102  double* u,
103  double* grad_u,
104  double* m,
105  double* dm,
106  double* H,
107  double* dH,
108  double* r);
109 
110 extern void conservativeHeadRichardsMualemVanGenuchtenHetEvaluateV2(const int nSimplex,
111  const int nPointsPerSimplex,
112  const int nSpace,
113  const int* materialTypes,
114  const double rho,
115  const double beta,
116  const double* gravity,
117  const double* alpha,
118  const double* n,
119  const double* thetaR,
120  const double* thetaSR,
121  const double* KWs,
122  double *u,
123  double *mass,
124  double *dmass,
125  double *f,
126  double *df,
127  double *a,
128  double *da);
129 extern void kEpsilon_2D_Evaluate(int nPoints,
130  int nSpace,
131  double sigma_k,
132  double sigma_e,
133  double c_1,
134  double c_2,
135  double c_mu,
136  double c_e,
137  double nu,
138  double *velocity,
139  double *gradu,
140  double *gradv,
141  double *k,
142  double *epsilon,
143  double *m_k,
144  double *dm_k,
145  double *m_e,
146  double *dm_e,
147  double *phi_k,
148  double *dphi_k,
149  double *phi_e,
150  double *dphi_e,
151  double *f_k,
152  double *df_k,
153  double *f_e,
154  double *df_e,
155  double *a_k,
156  double *da_k_dk,
157  double *da_k_de,
158  double *a_e,
159  double *da_e_dk,
160  double *da_e_de,
161  double *r_k,
162  double *dr_k_dk,
163  double *dr_k_de,
164  double *r_e,
165  double *dr_e_dk,
166  double *dr_e_de);
167 extern double smoothedHeaviside(
168  double eps,
169  double phi
170 );
171 extern double smoothedHeaviside_integral(
172  double eps,
173  double phi
174 );
175 extern double smoothedDirac(
176  double eps,
177  double phi
178 );
179 extern void linearADR_ConstantCoefficientsEvaluate(
180  const int nPoints,
181  const int nSpace,
182  const double M,
183  const double *A,
184  const double *B,
185  const double C,
186  const double t,
187  const double *x,
188  const double *u,
189  double *m,
190  double *dm,
191  double *f,
192  double *df,
193  double *a,
194  double *r,
195  double *dr
196 );
197 extern void nonlinearADR_pqrstEvaluate(
198  const int nPoints,
199  const int nSpace,
200  const double M,
201  const double *A,
202  const double *B,
203  const double C,
204  const double p_pow,
205  const double q_pow,
206  const double r_pow,
207  const double s_pow,
208  const double t_pow,
209  const double t,
210  const double *x,
211  const double *u,
212  double *m,
213  double *dm,
214  double *f,
215  double *df,
216  double *a,
217  double *da,
218  double *phi,
219  double *dphi,
220  double *r,
221  double *dr
222 );
223 extern void nonlinearADR_pqrstDualEvaluate(
224  const int nPoints,
225  const int nSpace,
226  const double M,
227  const double *A,
228  const double *B,
229  const double C,
230  const double p1,
231  const double q1,
232  const double r1,
233  const double s1,
234  const double t1,
235  const double p2,
236  const double q2,
237  const double r2,
238  const double s2,
239  const double t2,
240  const double t,
241  const double *x,
242  const double *u,
243  double *m,
244  double *dm,
245  double *f,
246  double *df,
247  double *a,
248  double *da,
249  double *phi,
250  double *dphi,
251  double *r,
252  double *dr
253 );
254 extern void unitSquareRotationEvaluate(
255  const int nPoints,
256  const int nSpace,
257  const double *x,
258  const double *u,
259  double *m,
260  double *dm,
261  double *f,
262  double *df
263 );
264 extern void rotatingPulseVelEvaluate(
265  const int nPoints,
266  const int nSpace,
267  const double self_a,
268  const double *x,
269  const double *u,
270  double *m,
271  double *dm,
272  double *f,
273  double *df,
274  double *a,
275  double *da,
276  double *phi,
277  double *dphi
278 );
279 extern void disRotatingPulseVelEvaluate(
280  const int nPoints,
281  const int nSpace,
282  const double self_a,
283  const double *x,
284  const double *u,
285  double *m,
286  double *dm,
287  double *f,
288  double *df,
289  double *a,
290  double *da,
291  double *phi,
292  double *dphi
293 );
294 extern void disVelEvaluate(
295  const int nPoints,
296  const int nSpace,
297  const double self_a,
298  const double *x,
299  const double *u,
300  double *m,
301  double *dm,
302  double *f,
303  double *df,
304  double *a,
305  double *da,
306  double *phi,
307  double *dphi
308 );
309 extern void burgersDiagonalVelEvaluate(
310  const int nPoints,
311  const int nSpace,
312  const double self_a,
313  const double * self_v,
314  const double *u,
315  double *m,
316  double *dm,
317  double *f,
318  double *df,
319  double *a,
320  double *phi,
321  double *dphi
322 );
323 extern void twophasePotentialFlowEvaluate(
324  int nPoints,
325  int nSpace,
326  double *M,
327  double *A,
328  double *B,
329  double *Bcon,
330  double *C,
331  double t,
332  double *x,
333  double *u,
334  double *m,
335  double *dm,
336  double *f,
337  double *df,
338  double *a,
339  double *da,
340  double *phi,
341  double *dphi,
342  double *r,
343  double *dr
344 );
345 extern void twophasePotentialFlowUpdateFreeSurface(
346  int nPoints,
347  int nSpace,
348  double eps,
349  double *u_levelSet,
350  double M1,
351  double M2,
352  double *M,
353  double *A1,
354  double *A2,
355  double *A,
356  double *B1,
357  double *B2,
358  double *B,
359  double *Bcon1,
360  double *Bcon2,
361  double *Bcon,
362  double C1,
363  double C2,
364  double *C
365 );
366 extern void twophaseLevelSetCoefficientsUpdateVelocity(
367  int nPoints,
368  int nSpace,
369  double v_scale,
370  double *vIn,
371  double *vOut
372 );
373 extern void twophaseLevelSetCoefficientsEvaluate(
374  int nPoints,
375  int nSpace,
376  double *B,
377  double t,
378  double *x,
379  double *u,
380  double *grad_u,
381  double *m,
382  double *dm,
383  double *h,
384  double *dh,
385  double *rh
386 );
387 extern void twophaseLevelSetCoefficientsEvaluateCI(
388  int nPoints,
389  int nSpace,
390  double *B,
391  double t,
392  double *x,
393  double *u,
394  double *m,
395  double *dm,
396  double *f,
397  double *df,
398  double *a,
399  double *da,
400  double *phi,
401  double *dphi,
402  double *r,
403  double *dr
404 );
405 extern void ncLevelSetCoefficientsEvaluate(
406  int nPoints,
407  int nSpace,
408  double *v,
409  double *u,
410  double *grad_u,
411  double *m,
412  double *dm,
413  double *H,
414  double *dH
415 );
416 extern void cLevelSetCoefficientsEvaluate(
417  int nPoints,
418  int nSpace,
419  double *v,
420  double *u,
421  double *m,
422  double *dm,
423  double *f,
424  double *df
425 );
426 extern void levelSetCurvatureCoefficientsEvaluate(
427  int nPoints,
428  int nSpace,
429  double *grad_phi,
430  double *u,
431  double *f,
432  double *r,
433  double *dr
434 );
435 extern void twophaseSignedDistanceCoefficientsUpdateSignFunction(
436  int nPoints,
437  double eps,
438  double *u_levelSet,
439  double *S
440 );
441 extern void twophaseSignedDistanceCoefficientsEvaluate(
442  int nPoints,
443  int nSpace,
444  double *S,
445  double *u,
446  double *grad_u,
447  double *m,
448  double *dm,
449  double *h,
450  double *dh,
451  double *rh
452 );
453 extern void redistanceLevelSetCoefficientsEvaluate(
454  int nPoints,
455  int nSpace,
456  double eps,
457  double *u_levelSet,
458  double *u,
459  double *grad_u,
460  double *m,
461  double *dm,
462  double *H,
463  double *dH,
464  double *r
465 );
466 extern void darcySharpInterfaceFlowEvaluate(
467  int nPoints,
468  int nSpace,
469  double Km,
470  double rhoM,
471  double Kp,
472  double rhoP,
473  double eps,
474  double *gravity_u,
475  double *u,
476  double *gradu,
477  double *u_levelSet,
478  double *phi_pot,
479  double *a,
480  double *f,
481  double *r,
482  double *m,
483  double *dphi_pot,
484  double *da,
485  double *df,
486  double *dr,
487  double *dm
488 );
489 extern void darcySharpInterfaceFlowImEvaluate(
490  int nPoints,
491  int nSpace,
492  double Km,
493  double rhoM,
494  double Kp,
495  double rhoP,
496  double eps,
497  double *gravity_u,
498  double *u,
499  double *gradu,
500  double *u_levelSet,
501  double *phi_pot,
502  double *a,
503  double *f,
504  double *r,
505  double *m,
506  double *dphi_pot,
507  double *da,
508  double *df,
509  double *dr,
510  double *dm
511 );
512 extern void conservativeHeadRichardsMualemVanGenuchtenHomEvaluate(
513  const int nPoints,
514  const int nSpace,
515  const double rho,
516  const double beta,
517  const double *gravity,
518  const double *x,
519  const double alpha,
520  const double n,
521  const double m,
522  const double thetaR,
523  const double thetaSR,
524  const double KWs,
525  double *u,
526  double *mass,
527  double *dmass,
528  double *f,
529  double *df,
530  double *a,
531  double *da,
532  double *phi,
533  double *dphi
534 );
535 extern void conservativeHeadRichardsL2projMualemVanGenuchtenHomEvaluate(
536  const int nSimplices,
537  const int nPointsPerSimplex,
538  const int nSpace,
539  const double rho,
540  const double *gravity,
541  const double alpha,
542  const double n,
543  const double m,
544  const double thetaR,
545  const double thetaSR,
546  const double KWs,
547  double *dV,
548  double *u,
549  double *mass,
550  double *dmass,
551  double *f,
552  double *df,
553  double *a,
554  double *da
555 );
556 extern void conservativeHeadRichardsL2projBndMualemVanGenuchtenHomEvaluate(
557  const int nElements,
558  const int nElementBoundaries_element,
559  const int nPointsPerElementBoundary,
560  const int nSpace,
561  const double rho,
562  const double *gravity,
563  const double alpha,
564  const double n,
565  const double m,
566  const double thetaR,
567  const double thetaSR,
568  const double KWs,
569  double *dV,
570  double *u,
571  double *mass,
572  double *dmass,
573  double *f,
574  double *df,
575  double *a,
576  double *da
577 );
578 extern void conservativeHeadRichardsL2projMualemVanGenuchtenHetEvaluate(
579  const int nSimplices,
580  const int nPointsPerSimplex,
581  const int nSpace,
582  const double rho,
583  const double *gravity,
584  const double *alpha,
585  const double *n,
586  const double *thetaR,
587  const double *thetaSR,
588  const double *KWs,
589  double *dV,
590  double *u,
591  double *mass,
592  double *dmass,
593  double *f,
594  double *df,
595  double *a,
596  double *da
597 );
598 extern void conservativeTotalHeadRichardsMualemVanGenuchtenHomEvaluate(
599  const int nPoints,
600  const int nSpace,
601  const double rho,
602  const double *gravity,
603  const double *x,
604  const double alpha,
605  const double n,
606  const double m,
607  const double thetaR,
608  const double thetaSR,
609  const double KWs,
610  double *u,
611  double *mass,
612  double *dmass,
613  double *f,
614  double *df,
615  double *a,
616  double *da,
617  double *phi,
618  double *dphi
619 );
620 extern void l2projectScalar(
621  const int nSimplices,
622  const int nPointsPerSimplex,
623  double *dV,
624  double *r
625 );
626 extern void l2projectVector(
627  const int nSimplices,
628  const int nPointsPerSimplex,
629  const int nSpace,
630  double *dV,
631  double *r
632 );
633 extern void l2project2Tensor(
634  const int nSimplices,
635  const int nPointsPerSimplex,
636  const int nSpace,
637  double *dV,
638  double *r
639 );
640 extern void conservativeHeadRichardsMualemVanGenuchtenHetEvaluate(
641  const int nPoints,
642  const int nSpace,
643  const double rho,
644  const double *gravity,
645  const double *alpha,
646  const double *n,
647  const double *thetaR,
648  const double *thetaSR,
649  const double *KWs,
650  double *u,
651  double *mass,
652  double *dmass,
653  double *f,
654  double *df,
655  double *a,
656  double *da
657 );
658 extern void conservativeSatRichardsMualemVanGenuchtenHomEvaluate(
659  const int nPoints,
660  const int nSpace,
661  const double rho,
662  const double *gravity,
663  const double *x,
664  const double alpha,
665  const double n,
666  const double m,
667  const double thetaR,
668  const double thetaSR,
669  const double KWs,
670  double *u,
671  double *mass,
672  double *dmass,
673  double *f,
674  double *df,
675  double *a,
676  double *da,
677  double *phi,
678  double *dphi
679 );
680 extern void conservativeHeadRichardsBrooksCoreyBurdineHetEvaluate(
681  const int nPoints,
682  const int nSpace,
683  const double rho,
684  const double *gravity,
685  const double *lambda,
686  const double *pd,
687  const double *thetaR,
688  const double *thetaS,
689  const double *KWs,
690  double *u,
691  double *mass,
692  double *dmass,
693  double *f,
694  double *df,
695  double *a,
696  double *da
697 );
698 extern void conservativeHeadRichardsBrooksCoreyBurdineHomEvaluate(
699  const int nPoints,
700  const int nSpace,
701  const double rho,
702  const double beta,
703  const double *gravity,
704  const double lambda,
705  const double pd,
706  const double thetaR,
707  const double thetaSR,
708  const double KWs,
709  double *u,
710  double *mass,
711  double *dmass,
712  double *f,
713  double *df,
714  double *a,
715  double *da
716 );
717 extern void conservativeSatRichardsBrooksCoreyBurdineHomEvaluate(
718  const int nPoints,
719  const int nSpace,
720  const double rho,
721  const double *gravity,
722  const double lambda,
723  const double pd,
724  const double thetaR,
725  const double thetaSR,
726  const double KWs,
727  double *u,
728  double *mass,
729  double *dmass,
730  double *f,
731  double *df,
732  double *a,
733  double *da,
734  double *phi,
735  double *dphi
736 );
737 extern void conservativeHeadRichardsBCBfromMVGHomEvaluate(
738  const int nPoints,
739  const int nSpace,
740  const double rho,
741  const double *gravity,
742  const double alpha,
743  const double n,
744  const double m,
745  const double thetaR,
746  const double thetaSR,
747  const double KWs,
748  double *u,
749  double *mass,
750  double *dmass,
751  double *f,
752  double *df,
753  double *a,
754  double *da
755 );
756 extern void NavierStokes_2D_Evaluate(
757  const int nPoints,
758  const double rho,
759  const double nu,
760  const double *g,
761  const double *p,
762  const double *grad_p,
763  const double *u,
764  const double *v,
765  double *mom_u_acc,
766  double *dmom_u_acc_u,
767  double *mom_v_acc,
768  double *dmom_v_acc_v,
769  double *mass_adv,
770  double *dmass_adv_u,
771  double *dmass_adv_v,
772  double *mom_u_adv,
773  double *dmom_u_adv_u,
774  double *dmom_u_adv_v,
775  double *mom_v_adv,
776  double *dmom_v_adv_u,
777  double *dmom_v_adv_v,
778  double *mom_u_diff_ten,
779  double *mom_v_diff_ten,
780  double *mom_u_source,
781  double *mom_v_source,
782  double *mom_u_ham,
783  double *dmom_u_ham_grad_p,
784  double *mom_v_ham,
785  double *dmom_v_ham_grad_p
786 );
787 extern void NavierStokes_3D_Evaluate(
788  const int nPoints,
789  const double rho,
790  const double nu,
791  const double *g,
792  const double *p,
793  const double *grad_p,
794  const double *u,
795  const double *v,
796  const double *w,
797  double *mom_u_acc,
798  double *dmom_u_acc_u,
799  double *mom_v_acc,
800  double *dmom_v_acc_v,
801  double *mom_w_acc,
802  double *dmom_w_acc_w,
803  double *mass_adv,
804  double *dmass_adv_u,
805  double *dmass_adv_v,
806  double *dmass_adv_w,
807  double *mom_u_adv,
808  double *dmom_u_adv_u,
809  double *dmom_u_adv_v,
810  double *dmom_u_adv_w,
811  double *mom_v_adv,
812  double *dmom_v_adv_u,
813  double *dmom_v_adv_v,
814  double *dmom_v_adv_w,
815  double *mom_w_adv,
816  double *dmom_w_adv_u,
817  double *dmom_w_adv_v,
818  double *dmom_w_adv_w,
819  double *mom_u_diff_ten,
820  double *mom_v_diff_ten,
821  double *mom_w_diff_ten,
822  double *mom_u_source,
823  double *mom_v_source,
824  double *mom_w_source,
825  double *mom_u_ham,
826  double *dmom_u_ham_grad_p,
827  double *mom_v_ham,
828  double *dmom_v_ham_grad_p,
829  double *mom_w_ham,
830  double *dmom_w_ham_grad_p
831 );
832 extern void Stokes_2D_Evaluate(
833  const int nPoints,
834  const double rho,
835  const double nu,
836  const double *g,
837  const double *p,
838  const double *grad_p,
839  const double *u,
840  const double *v,
841  double *mom_u_acc,
842  double *dmom_u_acc_u,
843  double *mom_v_acc,
844  double *dmom_v_acc_v,
845  double *mass_adv,
846  double *dmass_adv_u,
847  double *dmass_adv_v,
848  double *mom_u_diff_ten,
849  double *mom_v_diff_ten,
850  double *mom_u_source,
851  double *mom_v_source,
852  double *mom_u_ham,
853  double *dmom_u_ham_grad_p,
854  double *mom_v_ham,
855  double *dmom_v_ham_grad_p
856 );
857 extern void StokesP_2D_Evaluate(
858  const int nPoints,
859  const double rho,
860  const double nu,
861  const double *g,
862  const double *p,
863  const double *u,
864  const double *v,
865  double *mom_u_acc,
866  double *dmom_u_acc_u,
867  double *mom_v_acc,
868  double *dmom_v_acc_v,
869  double *mass_adv,
870  double *dmass_adv_u,
871  double *dmass_adv_v,
872  double *mom_u_adv,
873  double *dmom_u_adv_p,
874  double *mom_v_adv,
875  double *dmom_v_adv_p,
876  double *mom_u_diff_ten,
877  double *mom_v_diff_ten,
878  double *mom_u_source,
879  double *mom_v_source
880 );
881 extern void Stokes_3D_Evaluate(
882  const int nPoints,
883  const double rho,
884  const double nu,
885  const double *g,
886  const double *p,
887  const double *grad_p,
888  const double *u,
889  const double *v,
890  const double *w,
891  double *mom_u_acc,
892  double *dmom_u_acc_u,
893  double *mom_v_acc,
894  double *dmom_v_acc_v,
895  double *mom_w_acc,
896  double *dmom_w_acc_w,
897  double *mass_adv,
898  double *dmass_adv_u,
899  double *dmass_adv_v,
900  double *dmass_adv_w,
901  double *mom_u_diff_ten,
902  double *mom_v_diff_ten,
903  double *mom_w_diff_ten,
904  double *mom_u_source,
905  double *mom_v_source,
906  double *mom_w_source,
907  double *mom_u_ham,
908  double *dmom_u_ham_grad_p,
909  double *mom_v_ham,
910  double *dmom_v_ham_grad_p,
911  double *mom_w_ham,
912  double *dmom_w_ham_grad_p
913 );
914 extern void TwophaseNavierStokes_LS_SO_2D_Evaluate(
915  const int nPoints,
916  const double eps,
917  const double rho_0,
918  const double nu_0,
919  const double rho_1,
920  const double nu_1,
921  const double *g,
922  const double *phi,
923  const double *p,
924  const double *grad_p,
925  const double *u,
926  const double *v,
927  double *mom_u_acc,
928  double *dmom_u_acc_u,
929  double *mom_v_acc,
930  double *dmom_v_acc_v,
931  double *mass_adv,
932  double *dmass_adv_u,
933  double *dmass_adv_v,
934  double *mom_u_adv,
935  double *dmom_u_adv_u,
936  double *dmom_u_adv_v,
937  double *mom_v_adv,
938  double *dmom_v_adv_u,
939  double *dmom_v_adv_v,
940  double *mom_u_diff_ten,
941  double *mom_v_diff_ten,
942  double *mom_u_source,
943  double *mom_v_source,
944  double *mom_u_ham,
945  double *dmom_u_ham_grad_p,
946  double *mom_v_ham,
947  double *dmom_v_ham_grad_p
948 );
949 extern void TwophaseNavierStokes_ST_LS_SO_2D_Evaluate(
950  const int nPoints,
951  const double eps_rho,
952  const double eps_mu,
953  const double sigma,
954  const double rho_0,
955  const double nu_0,
956  const double rho_1,
957  const double nu_1,
958  const double *g,
959  const double *phi,
960  const double *n,
961  const double *kappa,
962  const double *p,
963  const double *grad_p,
964  const double *u,
965  const double *v,
966  double *mom_u_acc,
967  double *dmom_u_acc_u,
968  double *mom_v_acc,
969  double *dmom_v_acc_v,
970  double *mass_adv,
971  double *dmass_adv_u,
972  double *dmass_adv_v,
973  double *mom_u_adv,
974  double *dmom_u_adv_u,
975  double *dmom_u_adv_v,
976  double *mom_v_adv,
977  double *dmom_v_adv_u,
978  double *dmom_v_adv_v,
979  double *mom_u_diff_ten,
980  double *mom_v_diff_ten,
981  double *mom_uv_diff_ten,
982  double *mom_vu_diff_ten,
983  double *mom_u_source,
984  double *mom_v_source,
985  double *mom_u_ham,
986  double *dmom_u_ham_grad_p,
987  double *mom_v_ham,
988  double *dmom_v_ham_grad_p
989 );
990 extern void TwophaseNavierStokes_LS_SO_3D_Evaluate(
991  const int nPoints,
992  const double eps,
993  const double rho_0,
994  const double nu_0,
995  const double rho_1,
996  const double nu_1,
997  const double *g,
998  const double *phi,
999  const double *p,
1000  const double *grad_p,
1001  const double *u,
1002  const double *v,
1003  const double *w,
1004  double *mom_u_acc,
1005  double *dmom_u_acc_u,
1006  double *mom_v_acc,
1007  double *dmom_v_acc_v,
1008  double *mom_w_acc,
1009  double *dmom_w_acc_w,
1010  double *mass_adv,
1011  double *dmass_adv_u,
1012  double *dmass_adv_v,
1013  double *dmass_adv_w,
1014  double *mom_u_adv,
1015  double *dmom_u_adv_u,
1016  double *dmom_u_adv_v,
1017  double *dmom_u_adv_w,
1018  double *mom_v_adv,
1019  double *dmom_v_adv_u,
1020  double *dmom_v_adv_v,
1021  double *dmom_v_adv_w,
1022  double *mom_w_adv,
1023  double *dmom_w_adv_u,
1024  double *dmom_w_adv_v,
1025  double *dmom_w_adv_w,
1026  double *mom_u_diff_ten,
1027  double *mom_v_diff_ten,
1028  double *mom_w_diff_ten,
1029  double *mom_u_source,
1030  double *mom_v_source,
1031  double *mom_w_source,
1032  double *mom_u_ham,
1033  double *dmom_u_ham_grad_p,
1034  double *mom_v_ham,
1035  double *dmom_v_ham_grad_p,
1036  double *mom_w_ham,
1037  double *dmom_w_ham_grad_p
1038 );
1039 extern void TwophaseStokes_LS_SO_2D_Evaluate(
1040  const int nPoints,
1041  const double eps,
1042  const double rho_0,
1043  const double nu_0,
1044  const double rho_1,
1045  const double nu_1,
1046  const double *g,
1047  const double *phi,
1048  const double *p,
1049  const double *grad_p,
1050  const double *u,
1051  const double *v,
1052  double *mom_u_acc,
1053  double *dmom_u_acc_u,
1054  double *mom_v_acc,
1055  double *dmom_v_acc_v,
1056  double *mass_adv,
1057  double *dmass_adv_u,
1058  double *dmass_adv_v,
1059  double *mom_u_diff_ten,
1060  double *mom_v_diff_ten,
1061  double *mom_u_source,
1062  double *mom_v_source,
1063  double *mom_u_ham,
1064  double *dmom_u_ham_grad_p,
1065  double *mom_v_ham,
1066  double *dmom_v_ham_grad_p
1067 );
1068 extern void TwophaseStokes_LS_SO_3D_Evaluate(
1069  const int nPoints,
1070  const double eps,
1071  const double rho_0,
1072  const double nu_0,
1073  const double rho_1,
1074  const double nu_1,
1075  const double *g,
1076  const double *phi,
1077  const double *p,
1078  const double *grad_p,
1079  const double *u,
1080  const double *v,
1081  const double *w,
1082  double *mom_u_acc,
1083  double *dmom_u_acc_u,
1084  double *mom_v_acc,
1085  double *dmom_v_acc_v,
1086  double *mom_w_acc,
1087  double *dmom_w_acc_w,
1088  double *mass_adv,
1089  double *dmass_adv_u,
1090  double *dmass_adv_v,
1091  double *dmass_adv_w,
1092  double *mom_u_diff_ten,
1093  double *mom_v_diff_ten,
1094  double *mom_w_diff_ten,
1095  double *mom_u_source,
1096  double *mom_v_source,
1097  double *mom_w_source,
1098  double *mom_u_ham,
1099  double *dmom_u_ham_grad_p,
1100  double *mom_v_ham,
1101  double *dmom_v_ham_grad_p,
1102  double *mom_w_ham,
1103  double *dmom_w_ham_grad_p
1104 );
1105 extern void TwophaseNavierStokes_VOF_SO_2D_Evaluate(
1106  const int nPoints,
1107  const double eps,
1108  const double rho_0,
1109  const double nu_0,
1110  const double rho_1,
1111  const double nu_1,
1112  const double *g,
1113  const double *vof,
1114  const double *p,
1115  const double *grad_p,
1116  const double *u,
1117  const double *v,
1118  double *mom_u_acc,
1119  double *dmom_u_acc_u,
1120  double *mom_v_acc,
1121  double *dmom_v_acc_v,
1122  double *mass_adv,
1123  double *dmass_adv_u,
1124  double *dmass_adv_v,
1125  double *mom_u_adv,
1126  double *dmom_u_adv_u,
1127  double *dmom_u_adv_v,
1128  double *mom_v_adv,
1129  double *dmom_v_adv_u,
1130  double *dmom_v_adv_v,
1131  double *mom_u_diff_ten,
1132  double *mom_v_diff_ten,
1133  double *mom_u_source,
1134  double *mom_v_source,
1135  double *mom_u_ham,
1136  double *dmom_u_ham_grad_p,
1137  double *mom_v_ham,
1138  double *dmom_v_ham_grad_p
1139 );
1140 extern void TwophaseNavierStokes_VOF_SO_3D_Evaluate(
1141  const int nPoints,
1142  const double eps,
1143  const double rho_0,
1144  const double nu_0,
1145  const double rho_1,
1146  const double nu_1,
1147  const double *g,
1148  const double *vof,
1149  const double *p,
1150  const double *grad_p,
1151  const double *u,
1152  const double *v,
1153  const double *w,
1154  double *mom_u_acc,
1155  double *dmom_u_acc_u,
1156  double *mom_v_acc,
1157  double *dmom_v_acc_v,
1158  double *mom_w_acc,
1159  double *dmom_w_acc_w,
1160  double *mass_adv,
1161  double *dmass_adv_u,
1162  double *dmass_adv_v,
1163  double *dmass_adv_w,
1164  double *mom_u_adv,
1165  double *dmom_u_adv_u,
1166  double *dmom_u_adv_v,
1167  double *dmom_u_adv_w,
1168  double *mom_v_adv,
1169  double *dmom_v_adv_u,
1170  double *dmom_v_adv_v,
1171  double *dmom_v_adv_w,
1172  double *mom_w_adv,
1173  double *dmom_w_adv_u,
1174  double *dmom_w_adv_v,
1175  double *dmom_w_adv_w,
1176  double *mom_u_diff_ten,
1177  double *mom_v_diff_ten,
1178  double *mom_w_diff_ten,
1179  double *mom_u_source,
1180  double *mom_v_source,
1181  double *mom_w_source,
1182  double *mom_u_ham,
1183  double *dmom_u_ham_grad_p,
1184  double *mom_v_ham,
1185  double *dmom_v_ham_grad_p,
1186  double *mom_w_ham,
1187  double *dmom_w_ham_grad_p
1188 );
1189 extern void TwophaseStokes_VOF_SO_2D_Evaluate(
1190  const int nPoints,
1191  const double eps,
1192  const double rho_0,
1193  const double nu_0,
1194  const double rho_1,
1195  const double nu_1,
1196  const double *g,
1197  const double *vof,
1198  const double *p,
1199  const double *grad_p,
1200  const double *u,
1201  const double *v,
1202  double *mom_u_acc,
1203  double *dmom_u_acc_u,
1204  double *mom_v_acc,
1205  double *dmom_v_acc_v,
1206  double *mass_adv,
1207  double *dmass_adv_u,
1208  double *dmass_adv_v,
1209  double *mom_u_diff_ten,
1210  double *mom_v_diff_ten,
1211  double *mom_u_source,
1212  double *mom_v_source,
1213  double *mom_u_ham,
1214  double *dmom_u_ham_grad_p,
1215  double *mom_v_ham,
1216  double *dmom_v_ham_grad_p
1217 );
1218 extern void TwophaseStokes_VOF_SO_3D_Evaluate(
1219  const int nPoints,
1220  const double eps,
1221  const double rho_0,
1222  const double nu_0,
1223  const double rho_1,
1224  const double nu_1,
1225  const double *g,
1226  const double *vof,
1227  const double *p,
1228  const double *grad_p,
1229  const double *u,
1230  const double *v,
1231  const double *w,
1232  double *mom_u_acc,
1233  double *dmom_u_acc_u,
1234  double *mom_v_acc,
1235  double *dmom_v_acc_v,
1236  double *mom_w_acc,
1237  double *dmom_w_acc_w,
1238  double *mass_adv,
1239  double *dmass_adv_u,
1240  double *dmass_adv_v,
1241  double *dmass_adv_w,
1242  double *mom_u_diff_ten,
1243  double *mom_v_diff_ten,
1244  double *mom_w_diff_ten,
1245  double *mom_u_source,
1246  double *mom_v_source,
1247  double *mom_w_source,
1248  double *mom_u_ham,
1249  double *dmom_u_ham_grad_p,
1250  double *mom_v_ham,
1251  double *dmom_v_ham_grad_p,
1252  double *mom_w_ham,
1253  double *dmom_w_ham_grad_p
1254 );
1255 extern void VolumeAveragedNavierStokesFullDevStress_2D_Evaluate(const int nPoints,
1256  const double rho,
1257  const double mu,
1258  const double *meanGrainSize,
1259  const double *g,
1260  const double *p,
1261  const double *grad_p,
1262  const double *u,
1263  const double *v,
1264  const double *porosity,
1265  double *mom_u_acc,
1266  double *dmom_u_acc_u,
1267  double *mom_v_acc,
1268  double *dmom_v_acc_v,
1269  double *mass_adv,
1270  double *dmass_adv_u,
1271  double *dmass_adv_v,
1272  double *mom_u_adv,
1273  double *dmom_u_adv_u,
1274  double *dmom_u_adv_v,
1275  double *mom_v_adv,
1276  double *dmom_v_adv_u,
1277  double *dmom_v_adv_v,
1278  double *mom_u_diff_ten,
1279  double *mom_v_diff_ten,
1280  double *mom_uv_diff_ten,
1281  double *mom_vu_diff_ten,
1282  double *mom_u_source,
1283  double *mom_v_source,
1284  double *dmom_u_source_u,
1285  double *dmom_u_source_v,
1286  double *dmom_v_source_u,
1287  double *dmom_v_source_v,
1288  double *mom_u_ham,
1289  double *dmom_u_ham_grad_p,
1290  double *mom_v_ham,
1291  double *dmom_v_ham_grad_p);
1292 
1293 
1294 extern void VolumeAveragedNavierStokesFullDevStress_3D_Evaluate(const int nPoints,
1295  const double rho,
1296  const double mu,
1297  const double *meanGrainSize,
1298  const double *g,
1299  const double *p,
1300  const double *grad_p,
1301  const double *u,
1302  const double *v,
1303  const double *w,
1304  const double *porosity,
1305  double *mom_u_acc,
1306  double *dmom_u_acc_u,
1307  double *mom_v_acc,
1308  double *dmom_v_acc_v,
1309  double *mom_w_acc,
1310  double *dmom_w_acc_w,
1311  double *mass_adv,
1312  double *dmass_adv_u,
1313  double *dmass_adv_v,
1314  double *dmass_adv_w,
1315  double *mom_u_adv,
1316  double *dmom_u_adv_u,
1317  double *dmom_u_adv_v,
1318  double *dmom_u_adv_w,
1319  double *mom_v_adv,
1320  double *dmom_v_adv_u,
1321  double *dmom_v_adv_v,
1322  double *dmom_v_adv_w,
1323  double *mom_w_adv,
1324  double *dmom_w_adv_u,
1325  double *dmom_w_adv_v,
1326  double *dmom_w_adv_w,
1327  double *mom_u_diff_ten,
1328  double *mom_v_diff_ten,
1329  double *mom_w_diff_ten,
1330  double *mom_uv_diff_ten,
1331  double *mom_uw_diff_ten,
1332  double *mom_vu_diff_ten,
1333  double *mom_vw_diff_ten,
1334  double *mom_wu_diff_ten,
1335  double *mom_wv_diff_ten,
1336  double *mom_u_source,
1337  double *mom_v_source,
1338  double *mom_w_source,
1339  double *dmom_u_source_u,
1340  double *dmom_u_source_v,
1341  double *dmom_u_source_w,
1342  double *dmom_v_source_u,
1343  double *dmom_v_source_v,
1344  double *dmom_v_source_w,
1345  double *dmom_w_source_u,
1346  double *dmom_w_source_v,
1347  double *dmom_w_source_w,
1348  double *mom_u_ham,
1349  double *dmom_u_ham_grad_p,
1350  double *mom_v_ham,
1351  double *dmom_v_ham_grad_p,
1352  double *mom_w_ham,
1353  double *dmom_w_ham_grad_p);
1354 
1355 extern
1356 void VolumeAveragedTwophaseNavierStokes_ST_LS_SO_2D_Evaluate(const int nPoints,
1357  const int killNonlinearDrag,
1358  const double eps_rho,
1359  const double eps_mu,
1360  const double sigma,
1361  const double rho_0,
1362  const double nu_0,
1363  const double rho_1,
1364  const double nu_1,
1365  const double* meanGrainSize,
1366  const double* g,
1367  const double* phi,
1368  const double* n,
1369  const double* kappa,
1370  const double *p,
1371  const double *grad_p,
1372  const double *u,
1373  const double *v,
1374  const double *porosity,
1375  double *mom_u_acc,
1376  double *dmom_u_acc_u,
1377  double *mom_v_acc,
1378  double *dmom_v_acc_v,
1379  double *mass_adv,
1380  double *dmass_adv_u,
1381  double *dmass_adv_v,
1382  double *mom_u_adv,
1383  double *dmom_u_adv_u,
1384  double *dmom_u_adv_v,
1385  double *mom_v_adv,
1386  double *dmom_v_adv_u,
1387  double *dmom_v_adv_v,
1388  double *mom_u_diff_ten,
1389  double *mom_v_diff_ten,
1390  double *mom_uv_diff_ten,
1391  double *mom_vu_diff_ten,
1392  double *mom_u_source,
1393  double *mom_v_source,
1394  double *dmom_u_source_u,
1395  double *dmom_u_source_v,
1396  double *dmom_v_source_u,
1397  double *dmom_v_source_v,
1398  double *mom_u_ham,
1399  double *dmom_u_ham_grad_p,
1400  double *mom_v_ham,
1401  double *dmom_v_ham_grad_p);
1402 
1403 extern
1404 void VolumeAveragedVOFCoefficientsEvaluate(int nPoints,
1405  int nSpace,
1406  double eps,
1407  double* v,
1408  double* phi,
1409  double* porosity,
1410  double* u,
1411  double* m,
1412  double* dm,
1413  double* f,
1414  double* df);
1415 extern void unitSquareVortexEvaluate(
1416  const int nPoints,
1417  const int nSpace,
1418  double t,
1419  const double *x,
1420  const double *u,
1421  double *m,
1422  double *dm,
1423  double *f,
1424  double *df
1425 );
1426 extern void constantVelocityLevelSetEvaluate(
1427  const int nPoints,
1428  const int nSpace,
1429  const double *b,
1430  const double *x,
1431  const double *u,
1432  const double *gradu,
1433  double *m,
1434  double *dm,
1435  double *f,
1436  double *df,
1437  double *H,
1438  double *dH
1439 );
1440 extern void constantNormalVelocityLevelSetEvaluate(
1441  const int nPoints,
1442  const int nSpace,
1443  double b,
1444  const double *x,
1445  const double *u,
1446  const double *gradu,
1447  double *m,
1448  double *dm,
1449  double *f,
1450  double *df,
1451  double *H,
1452  double *dH
1453 );
1454 extern void unitSquareVortexLevelSetEvaluate(
1455  const int nPoints,
1456  const int nSpace,
1457  double t,
1458  const double *x,
1459  const double *u,
1460  const double *gradu,
1461  double *m,
1462  double *dm,
1463  double *f,
1464  double *df,
1465  double *H,
1466  double *dH
1467 );
1468 extern void unitSquareRotationLevelSetEvaluate(
1469  const int nPoints,
1470  const int nSpace,
1471  double t,
1472  const double *x,
1473  const double *u,
1474  const double *gradu,
1475  double *m,
1476  double *dm,
1477  double *f,
1478  double *df,
1479  double *H,
1480  double *dH
1481 );
1482 extern void FractionalFlowPhaseForm_saturationEvaluate(const int nPoints,
1483  const int nSpace,
1484  const int nc,
1485  const int pskModelFlag,
1486  const double Kbar,
1487  const double rhon,
1488  const double rhow,
1489  const double *g,
1490  const double g_norm,
1491  const double alpha,
1492  const double bc_lambda,
1493  const double bc_pd,
1494  const double mvg_n,
1495  const double mvg_m,
1496  const double omega,
1497  const double mun,
1498  const double muw,
1499  const double sw_min,
1500  const double sw_max,
1501  const double M,
1502  const double R,
1503  const double Temp,
1504  const double p_o,
1505  const double b,
1506  double *u,
1507  double *m,
1508  double *dm,
1509  double *phi,
1510  double *dphi,
1511  double *f,
1512  double *df,
1513  double *a,
1514  double *da,
1515  double *q_t,
1516  double *psiw
1517  );
1518 extern void FractionalFlowPhaseForm_potentialEvaluate(const int nPoints,
1519  const int nSpace,
1520  const int nc,
1521  const int pskModelFlag,
1522  const double Kbar,
1523  const double rhon,
1524  const double rhow,
1525  const double *g,
1526  const double g_norm,
1527  const double alpha,
1528  const double bc_lambda,
1529  const double bc_pd,
1530  const double mvg_n,
1531  const double mvg_m,
1532  const double omega,
1533  const double mun,
1534  const double muw,
1535  const double sw_min,
1536  const double sw_max,
1537  const double M,
1538  const double R,
1539  const double Temp,
1540  const double p_o,
1541  const double b,
1542  double *u,
1543  double *m,
1544  double *dm,
1545  double *phi,
1546  double *dphi,
1547  double *f,
1548  double *df,
1549  double *a,
1550  double *da,
1551  double *s_w,
1552  double *grad_psic);
1553 extern void FractionalFlowPhaseForm_saturationHetEvaluate(const int nPoints,
1554  const int nSpace,
1555  const int nc,
1556  const int pskModelFlag,
1557  const double *Kbar,
1558  const double rhon,
1559  const double rhow,
1560  const double *g,
1561  const double *alpha,
1562  const double *bc_lambda,
1563  const double *bc_pd,
1564  const double *mvg_m,
1565  const double *thetaS,
1566  const double *thetaR,
1567  const double mun,
1568  const double muw,
1569  const double b,
1570  double *u,
1571  double *m,
1572  double *dm,
1573  double *phi,
1574  double *dphi,
1575  double *f,
1576  double *df,
1577  double *a,
1578  double *da,
1579  double *q_t);
1580 extern void FractionalFlowPhaseForm_potentialHetEvaluate(const int nPoints,
1581  const int nSpace,
1582  const int nc,
1583  const int pskModelFlag,
1584  const double *Kbar,
1585  const double rhon,
1586  const double rhow,
1587  const double *g,
1588  const double *alpha,
1589  const double *bc_lambda,
1590  const double *bc_pd,
1591  const double *mvg_m,
1592  const double *thetaS,
1593  const double *thetaR,
1594  const double mun,
1595  const double muw,
1596  const double b,
1597  double *u,
1598  double *m,
1599  double *dm,
1600  double *phi,
1601  double *dphi,
1602  double *f,
1603  double *df,
1604  double *a,
1605  double *da,
1606  double *s_w,
1607  double *grad_psic);
1608 extern void TwophaseDarcyFC_Evaluate(
1609  const int nPoints,
1610  const int nSpace,
1611  const int pskModelFlag,
1612  const double Kbar,
1613  const double rhon,
1614  const double rhow,
1615  const double *g,
1616  const double *x,
1617  const double alpha,
1618  const double bc_lambda,
1619  const double bc_pd,
1620  const double mvg_n,
1621  const double mvg_m,
1622  const double omega,
1623  const double omega_r,
1624  const double mun,
1625  const double muw,
1626  const double b,
1627  double *sw,
1628  double *psiw,
1629  double *mw,
1630  double *dmw,
1631  double *mn,
1632  double *dmn,
1633  double *phi_psiw,
1634  double *dphi_psiw_dpsiw,
1635  double *phi_psin,
1636  double *dphi_psin_dpsiw,
1637  double *dphi_psin_dsw,
1638  double *fw,
1639  double *dfw,
1640  double *fn,
1641  double *dfn,
1642  double *aw,
1643  double *daw,
1644  double *an,
1645  double *dan
1646 );
1647 extern void TwophaseFFDarcyFC_Evaluate(
1648  const int nPoints,
1649  const int nSpace,
1650  const int pskModelFlag,
1651  const double Kbar,
1652  const double rhon,
1653  const double rhow,
1654  const double *g,
1655  const double *x,
1656  const double alpha,
1657  const double bc_lambda,
1658  const double bc_pd,
1659  const double mvg_n,
1660  const double mvg_m,
1661  const double omega,
1662  const double omega_r,
1663  const double mun,
1664  const double muw,
1665  const double b,
1666  double *sw,
1667  double *psiw,
1668  double *mw,
1669  double *dmw_dsw,
1670  double *mm,
1671  double *dmm_dsw,
1672  double *phi_psic,
1673  double *dphi_psic_dsw,
1674  double *phi_psiw,
1675  double *dphi_psiw_dpsiw,
1676  double *fm,
1677  double *dfm_dsw,
1678  double *aw_psiw,
1679  double *daw_psiw_dsw,
1680  double *am_psiw,
1681  double *dam_psiw_dsw,
1682  double *am_psic,
1683  double *dam_psic_dsw
1684 );
1685 extern void TwophaseDarcyFCHet_Evaluate(
1686  const int nPoints,
1687  const int nSpace,
1688  const int pskModelFlag,
1689  const double *Kbar,
1690  const double rhon,
1691  const double rhow,
1692  const double *g,
1693  const double *x,
1694  const double *alpha,
1695  const double *bc_lambda,
1696  const double *bc_pd,
1697  const double *mvg_m,
1698  const double *omega,
1699  const double *omega_r,
1700  const double mun,
1701  const double muw,
1702  const double b,
1703  double *sw,
1704  double *psiw,
1705  double *mw,
1706  double *dmw,
1707  double *mn,
1708  double *dmn,
1709  double *phi_psiw,
1710  double *dphi_psiw_dpsiw,
1711  double *phi_psin,
1712  double *dphi_psin_dpsiw,
1713  double *dphi_psin_dsw,
1714  double *fw,
1715  double *dfw,
1716  double *fn,
1717  double *dfn,
1718  double *aw,
1719  double *daw,
1720  double *an,
1721  double *dan
1722 );
1723 extern void TwophaseFFDarcyFCHet_Evaluate(
1724  const int nPoints,
1725  const int nSpace,
1726  const int pskModelFlag,
1727  const double *Kbar,
1728  const double rhon,
1729  const double rhow,
1730  const double *g,
1731  const double *x,
1732  const double *alpha,
1733  const double *bc_lambda,
1734  const double *bc_pd,
1735  const double *mvg_m,
1736  const double *omega,
1737  const double *omega_r,
1738  const double mun,
1739  const double muw,
1740  const double b,
1741  double *sw,
1742  double *psiw,
1743  double *mw,
1744  double *dmw_dsw,
1745  double *mm,
1746  double *dmm_dsw,
1747  double *phi_psic,
1748  double *dphi_psic_dsw,
1749  double *phi_psiw,
1750  double *dphi_psiw_dpsiw,
1751  double *fm,
1752  double *dfm_dsw,
1753  double *aw_psiw,
1754  double *daw_psiw_dsw,
1755  double *am_psiw,
1756  double *dam_psiw_dsw,
1757  double *am_psic,
1758  double *dam_psic_dsw
1759 );
1760 extern void LinearElasticity_1D_Evaluate(
1761  const int nPoints,
1762  const double E,
1763  const double nu,
1764  const double *g,
1765  const double *u,
1766  double *uu_diff_ten,
1767  double *u_force
1768 );
1769 extern void LinearElasticity_2D_Evaluate(
1770  const int nPoints,
1771  const double E,
1772  const double nu,
1773  const double *g,
1774  const double *u,
1775  const double *v,
1776  double *uu_diff_ten,
1777  double *uv_diff_ten,
1778  double *vu_diff_ten,
1779  double *vv_diff_ten,
1780  double *u_force,
1781  double *v_force
1782 );
1783 extern void LinearElasticity_3D_Evaluate(
1784  const int nPoints,
1785  const double E,
1786  const double nu,
1787  const double *g,
1788  const double *u,
1789  const double *v,
1790  const double *w,
1791  double *uu_diff_ten,
1792  double *uv_diff_ten,
1793  double *uw_diff_ten,
1794  double *vu_diff_ten,
1795  double *vv_diff_ten,
1796  double *vw_diff_ten,
1797  double *wu_diff_ten,
1798  double *wv_diff_ten,
1799  double *ww_diff_ten,
1800  double *u_force,
1801  double *v_force,
1802  double *w_force
1803 );
1804 extern void levelSetConservationCoefficientsEvaluate(int nPoints,
1805  int nSpace,
1806  double epsHeaviside,
1807  double epsDirac,
1808  double epsDiffusion,
1809  double* u_ls,
1810  double* H_vof,
1811  double* u,
1812  double* r,
1813  double* dr,
1814  double* a);
1815 
1816 
1817 extern void unitCubeRotationEvaluate(const int nPoints,
1818  const int nSpace,
1819  const double *x,
1820  const double *u,
1821  double *m,
1822  double *dm,
1823  double *f,
1824  double *df);
1825 
1826 extern void evaluateBuckleyLeverettLiuExample(int nPoints,
1827  int nSpace,
1828  const double * x,
1829  const double * u,
1830  double * m,
1831  double * dm,
1832  double * f,
1833  double * df,
1834  double * a);
1835 
1836 
1837 extern void eikonalEquationEvaluate(int nPoints,
1838  int nSpace,
1839  double rhs,
1840  double* u,
1841  double* grad_u,
1842  double* m,
1843  double* dm,
1844  double* H,
1845  double* dH,
1846  double* r);
1847 
1848 extern void burgersDiagonalVelHJEvaluate(const int nPoints,
1849  const int nSpace,
1850  const double self_a,
1851  const double *self_v,
1852  const double *u,
1853  const double *grad_u,
1854  double *m,
1855  double *dm,
1856  double *H,
1857  double *dH,
1858  double *a,
1859  double *phi,
1860  double *dphi);
1861 
1862 extern void HJBurgersEvaluate(const int nPoints,
1863  const int nSpace,
1864  const double offset,
1865  const double *u,
1866  const double *gradu,
1867  double *m,
1868  double *dm,
1869  double *H,
1870  double *dH);
1871 extern void kEpsilon_2D_Evaluate(int nPoints,
1872  int nSpace,
1873  double sigma_k,
1874  double sigma_e,
1875  double c_1,
1876  double c_2,
1877  double c_mu,
1878  double c_e,
1879  double nu,
1880  double *velocity,
1881  double *gradu,
1882  double *gradv,
1883  double *k,
1884  double *epsilon,
1885  double *m_k,
1886  double *dm_k,
1887  double *m_e,
1888  double *dm_e,
1889  double *phi_k,
1890  double *dphi_k,
1891  double *phi_e,
1892  double *dphi_e,
1893  double *f_k,
1894  double *df_k,
1895  double *f_e,
1896  double *df_e,
1897  double *a_k,
1898  double *da_k_dk,
1899  double *da_k_de,
1900  double *a_e,
1901  double *da_e_dk,
1902  double *da_e_de,
1903  double *r_k,
1904  double *dr_k_dk,
1905  double *dr_k_de,
1906  double *r_e,
1907  double *dr_e_dk,
1908  double *dr_e_de);
1909 
1910 extern void ReynoldsAveragedNavierStokes_kEpsilon_2D_Evaluate(const int nPoints,
1911  const double rho,
1912  const double nu,
1913  const double c_mu,
1914  const double* k,
1915  const double* grad_k,
1916  const double* epsilon,
1917  const double* g,
1918  const double *p,
1919  const double *grad_p,
1920  const double *u,
1921  const double *v,
1922  double *mom_u_acc,
1923  double *dmom_u_acc_u,
1924  double *mom_v_acc,
1925  double *dmom_v_acc_v,
1926  double *mass_adv,
1927  double *dmass_adv_u,
1928  double *dmass_adv_v,
1929  double *mom_u_adv,
1930  double *dmom_u_adv_u,
1931  double *dmom_u_adv_v,
1932  double *mom_v_adv,
1933  double *dmom_v_adv_u,
1934  double *dmom_v_adv_v,
1935  double *mom_u_diff_ten,
1936  double *mom_v_diff_ten,
1937  double *mom_uv_diff_ten,
1938  double *mom_vu_diff_ten,
1939  double *mom_u_source,
1940  double *mom_v_source,
1941  double *mom_u_ham,
1942  double *dmom_u_ham_grad_p,
1943  double *mom_v_ham,
1944  double *dmom_v_ham_grad_p);
1945 
1946 extern
1947 void ThreephaseNavierStokes_ST_LS_SO_3D_Evaluate(const int nPoints,
1948  const double boundaryPenaltyCoef,
1949  const double volumePenaltyCoef,
1950  const double eps_rho,
1951  const double eps_mu,
1952  const double sigma,
1953  const double rho_0,
1954  const double nu_0,
1955  const double rho_1,
1956  const double nu_1,
1957  const double rho_s,
1958  const double nu_s,
1959  const double* g,
1960  const double* phi,
1961  const double* n,
1962  const double* kappa,
1963  const double* phi_s,
1964  const double* n_s,
1965  const double *p,
1966  const double *grad_p,
1967  const double *u,
1968  const double *v,
1969  const double *w,
1970  double *mom_u_acc,
1971  double *dmom_u_acc_u,
1972  double *mom_v_acc,
1973  double *dmom_v_acc_v,
1974  double *mom_w_acc,
1975  double *dmom_w_acc_w,
1976  double *mass_adv,
1977  double *dmass_adv_u,
1978  double *dmass_adv_v,
1979  double *dmass_adv_w,
1980  double *mom_u_adv,
1981  double *dmom_u_adv_u,
1982  double *dmom_u_adv_v,
1983  double *dmom_u_adv_w,
1984  double *mom_v_adv,
1985  double *dmom_v_adv_u,
1986  double *dmom_v_adv_v,
1987  double *dmom_v_adv_w,
1988  double *mom_w_adv,
1989  double *dmom_w_adv_u,
1990  double *dmom_w_adv_v,
1991  double *dmom_w_adv_w,
1992  double *mom_u_diff_ten,
1993  double *mom_v_diff_ten,
1994  double *mom_w_diff_ten,
1995  double *mom_uv_diff_ten,
1996  double *mom_uw_diff_ten,
1997  double *mom_vu_diff_ten,
1998  double *mom_vw_diff_ten,
1999  double *mom_wu_diff_ten,
2000  double *mom_wv_diff_ten,
2001  double *mom_u_source,
2002  double *dmom_u_source_u,
2003  double *dmom_u_source_v,
2004  double *dmom_u_source_w,
2005  double *mom_v_source,
2006  double *dmom_v_source_u,
2007  double *dmom_v_source_v,
2008  double *dmom_v_source_w,
2009  double *mom_w_source,
2010  double *dmom_w_source_u,
2011  double *dmom_w_source_v,
2012  double *dmom_w_source_w,
2013  double *mom_u_ham,
2014  double *dmom_u_ham_grad_p,
2015  double *mom_v_ham,
2016  double *dmom_v_ham_grad_p,
2017  double *mom_w_ham,
2018  double *dmom_w_ham_grad_p);
2019 
2020 /* int findInterval(const double* vertices, int nv, double x, int* ival, double tol); */
2021 
2022 /* void piecewiseLinearTableLookup(double x, */
2023 /* int nv, */
2024 /* int* start, */
2025 /* double* y, */
2026 /* double* dy, */
2027 /* const double* xv, */
2028 /* const double* yv); */
2029 
2031 #endif
TwophaseDarcyFCHet_Evaluate
void TwophaseDarcyFCHet_Evaluate(const int nPoints, const int nSpace, const int pskModelFlag, const double *Kbar, const double rhon, const double rhow, const double *g, const double *x, const double *alpha, const double *bc_lambda, const double *bc_pd, const double *mvg_m, const double *omega, const double *omega_r, const double mun, const double muw, const double b, double *sw, double *psiw, double *mw, double *dmw, double *mn, double *dmn, double *phi_psiw, double *dphi_psiw_dpsiw, double *phi_psin, double *dphi_psin_dpsiw, double *dphi_psin_dsw, double *fw, double *dfw, double *fn, double *dfn, double *aw, double *daw, double *an, double *dan)
ReynoldsAveragedNavierStokes_kEpsilon_2D_Evaluate
void ReynoldsAveragedNavierStokes_kEpsilon_2D_Evaluate(const int nPoints, const double rho, const double nu, const double c_mu, const double *k, const double *grad_k, const double *epsilon, const double *g, const double *p, const double *grad_p, const double *u, const double *v, double *mom_u_acc, double *dmom_u_acc_u, double *mom_v_acc, double *dmom_v_acc_v, double *mass_adv, double *dmass_adv_u, double *dmass_adv_v, double *mom_u_adv, double *dmom_u_adv_u, double *dmom_u_adv_v, double *mom_v_adv, double *dmom_v_adv_u, double *dmom_v_adv_v, double *mom_u_diff_ten, double *mom_v_diff_ten, double *mom_uv_diff_ten, double *mom_vu_diff_ten, double *mom_u_source, double *mom_v_source, double *mom_u_ham, double *dmom_u_ham_grad_p, double *mom_v_ham, double *dmom_v_ham_grad_p)
Stokes_3D_Evaluate
void Stokes_3D_Evaluate(const int nPoints, const double rho, const double nu, const double *g, const double *p, const double *grad_p, const double *u, const double *v, const double *w, double *mom_u_acc, double *dmom_u_acc_u, double *mom_v_acc, double *dmom_v_acc_v, double *mom_w_acc, double *dmom_w_acc_w, double *mass_adv, double *dmass_adv_u, double *dmass_adv_v, double *dmass_adv_w, double *mom_u_diff_ten, double *mom_v_diff_ten, double *mom_w_diff_ten, double *mom_u_source, double *mom_v_source, double *mom_w_source, double *mom_u_ham, double *dmom_u_ham_grad_p, double *mom_v_ham, double *dmom_v_ham_grad_p, double *mom_w_ham, double *dmom_w_ham_grad_p)
Definition: transportCoefficients.c:2206
smoothedHeaviside
double smoothedHeaviside(double eps, double phi)
Definition: numericalFlux.c:14167
TwophaseNavierStokes_LS_SO_2D_Evaluate
void TwophaseNavierStokes_LS_SO_2D_Evaluate(const int nPoints, const double eps, const double rho_0, const double nu_0, const double rho_1, const double nu_1, const double *g, const double *phi, const double *p, const double *grad_p, const double *u, const double *v, double *mom_u_acc, double *dmom_u_acc_u, double *mom_v_acc, double *dmom_v_acc_v, double *mass_adv, double *dmass_adv_u, double *dmass_adv_v, double *mom_u_adv, double *dmom_u_adv_u, double *dmom_u_adv_v, double *mom_v_adv, double *dmom_v_adv_u, double *dmom_v_adv_v, double *mom_u_diff_ten, double *mom_v_diff_ten, double *mom_u_source, double *mom_v_source, double *mom_u_ham, double *dmom_u_ham_grad_p, double *mom_v_ham, double *dmom_v_ham_grad_p)
Definition: transportCoefficients.c:2383
conservativeHeadRichardsMualemVanGenuchtenHetEvaluate
void conservativeHeadRichardsMualemVanGenuchtenHetEvaluate(const int nPoints, const int nSpace, const double rho, const double *gravity, const double *alpha, const double *n, const double *thetaR, const double *thetaSR, const double *KWs, double *u, double *mass, double *dmass, double *f, double *df, double *a, double *da)
Definition: transportCoefficients.c:6337
redistanceLevelSetSandFCoefficientsEvaluate
void redistanceLevelSetSandFCoefficientsEvaluate(int nSimplex, int nPointsPerSimplex, int nSpace, double eps, double *u_levelSet, double *dV, double *u, double *grad_u, double *m, double *dm, double *H, double *dH, double *r)
Definition: transportCoefficients.c:1512
kEpsilon_2D_Evaluate
void kEpsilon_2D_Evaluate(int nPoints, int nSpace, double sigma_k, double sigma_e, double c_1, double c_2, double c_mu, double c_e, double nu, double *velocity, double *gradu, double *gradv, double *k, double *epsilon, double *m_k, double *dm_k, double *m_e, double *dm_e, double *phi_k, double *dphi_k, double *phi_e, double *dphi_e, double *f_k, double *df_k, double *f_e, double *df_e, double *a_k, double *da_k_dk, double *da_k_de, double *a_e, double *da_e_dk, double *da_e_de, double *r_k, double *dr_k_dk, double *dr_k_de, double *r_e, double *dr_e_dk, double *dr_e_de)
Definition: transportCoefficients.c:9975
w
#define w(x)
Definition: jf.h:22
rotatingPulseVelEvaluate
void rotatingPulseVelEvaluate(const int nPoints, const int nSpace, const double self_a, const double *x, const double *u, double *m, double *dm, double *f, double *df, double *a, double *da, double *phi, double *dphi)
Definition: transportCoefficients.c:818
B
Double * B
Definition: Headers.h:41
l2project2Tensor
void l2project2Tensor(const int nSimplices, const int nPointsPerSimplex, const int nSpace, double *dV, double *r)
Definition: transportCoefficients.c:5705
twophaseLevelSetCoefficientsEvaluate
void twophaseLevelSetCoefficientsEvaluate(int nPoints, int nSpace, double *B, double t, double *x, double *u, double *grad_u, double *m, double *dm, double *h, double *dh, double *rh)
Definition: transportCoefficients.c:1147
l2projectVector
void l2projectVector(const int nSimplices, const int nPointsPerSimplex, const int nSpace, double *dV, double *r)
Definition: transportCoefficients.c:5676
conservativeSatRichardsMualemVanGenuchtenHomEvaluate
void conservativeSatRichardsMualemVanGenuchtenHomEvaluate(const int nPoints, const int nSpace, const double rho, const double *gravity, const double *x, const double alpha, const double n, const double m, const double thetaR, const double thetaSR, const double KWs, double *u, double *mass, double *dmass, double *f, double *df, double *a, double *da, double *phi, double *dphi)
Definition: transportCoefficients.c:6418
TwophaseStokes_LS_SO_3D_Evaluate
void TwophaseStokes_LS_SO_3D_Evaluate(const int nPoints, const double eps, const double rho_0, const double nu_0, const double rho_1, const double nu_1, const double *g, const double *phi, const double *p, const double *grad_p, const double *u, const double *v, const double *w, double *mom_u_acc, double *dmom_u_acc_u, double *mom_v_acc, double *dmom_v_acc_v, double *mom_w_acc, double *dmom_w_acc_w, double *mass_adv, double *dmass_adv_u, double *dmass_adv_v, double *dmass_adv_w, double *mom_u_diff_ten, double *mom_v_diff_ten, double *mom_w_diff_ten, double *mom_u_source, double *mom_v_source, double *mom_w_source, double *mom_u_ham, double *dmom_u_ham_grad_p, double *mom_v_ham, double *dmom_v_ham_grad_p, double *mom_w_ham, double *dmom_w_ham_grad_p)
Definition: transportCoefficients.c:4203
NavierStokes_2D_Evaluate
void NavierStokes_2D_Evaluate(const int nPoints, const double rho, const double nu, const double *g, const double *p, const double *grad_p, const double *u, const double *v, double *mom_u_acc, double *dmom_u_acc_u, double *mom_v_acc, double *dmom_v_acc_v, double *mass_adv, double *dmass_adv_u, double *dmass_adv_v, double *mom_u_adv, double *dmom_u_adv_u, double *dmom_u_adv_v, double *mom_v_adv, double *dmom_v_adv_u, double *dmom_v_adv_v, double *mom_u_diff_ten, double *mom_v_diff_ten, double *mom_u_source, double *mom_v_source, double *mom_u_ham, double *dmom_u_ham_grad_p, double *mom_v_ham, double *dmom_v_ham_grad_p)
Definition: transportCoefficients.c:1850
TwophaseNavierStokes_LS_SO_3D_Evaluate
void TwophaseNavierStokes_LS_SO_3D_Evaluate(const int nPoints, const double eps, const double rho_0, const double nu_0, const double rho_1, const double nu_1, const double *g, const double *phi, const double *p, const double *grad_p, const double *u, const double *v, const double *w, double *mom_u_acc, double *dmom_u_acc_u, double *mom_v_acc, double *dmom_v_acc_v, double *mom_w_acc, double *dmom_w_acc_w, double *mass_adv, double *dmass_adv_u, double *dmass_adv_v, double *dmass_adv_w, double *mom_u_adv, double *dmom_u_adv_u, double *dmom_u_adv_v, double *dmom_u_adv_w, double *mom_v_adv, double *dmom_v_adv_u, double *dmom_v_adv_v, double *dmom_v_adv_w, double *mom_w_adv, double *dmom_w_adv_u, double *dmom_w_adv_v, double *dmom_w_adv_w, double *mom_u_diff_ten, double *mom_v_diff_ten, double *mom_w_diff_ten, double *mom_u_source, double *mom_v_source, double *mom_w_source, double *mom_u_ham, double *dmom_u_ham_grad_p, double *mom_v_ham, double *dmom_v_ham_grad_p, double *mom_w_ham, double *dmom_w_ham_grad_p)
Definition: transportCoefficients.c:3981
TwophaseFFDarcyFCHet_Evaluate
void TwophaseFFDarcyFCHet_Evaluate(const int nPoints, const int nSpace, const int pskModelFlag, const double *Kbar, const double rhon, const double rhow, const double *g, const double *x, const double *alpha, const double *bc_lambda, const double *bc_pd, const double *mvg_m, const double *omega, const double *omega_r, const double mun, const double muw, const double b, double *sw, double *psiw, double *mw, double *dmw_dsw, double *mm, double *dmm_dsw, double *phi_psic, double *dphi_psic_dsw, double *phi_psiw, double *dphi_psiw_dpsiw, double *fm, double *dfm_dsw, double *aw_psiw, double *daw_psiw_dsw, double *am_psiw, double *dam_psiw_dsw, double *am_psic, double *dam_psic_dsw)
StokesP_2D_Evaluate
void StokesP_2D_Evaluate(const int nPoints, const double rho, const double nu, const double *g, const double *p, const double *u, const double *v, double *mom_u_acc, double *dmom_u_acc_u, double *mom_v_acc, double *dmom_v_acc_v, double *mass_adv, double *dmass_adv_u, double *dmass_adv_v, double *mom_u_adv, double *dmom_u_adv_p, double *mom_v_adv, double *dmom_v_adv_p, double *mom_u_diff_ten, double *mom_v_diff_ten, double *mom_u_source, double *mom_v_source)
Definition: transportCoefficients.c:2143
TwophaseNavierStokes_ST_LS_SO_3D_Evaluate
void TwophaseNavierStokes_ST_LS_SO_3D_Evaluate(const int nPoints, const double eps_rho, const double eps_mu, const double sigma, const double rho_0, const double nu_0, const double rho_1, const double nu_1, const double *g, const double *phi, const double *n, const double *kappa, const double *p, const double *grad_p, const double *u, const double *v, const double *w, double *mom_u_acc, double *dmom_u_acc_u, double *mom_v_acc, double *dmom_v_acc_v, double *mom_w_acc, double *dmom_w_acc_w, double *mass_adv, double *dmass_adv_u, double *dmass_adv_v, double *dmass_adv_w, double *mom_u_adv, double *dmom_u_adv_u, double *dmom_u_adv_v, double *dmom_u_adv_w, double *mom_v_adv, double *dmom_v_adv_u, double *dmom_v_adv_v, double *dmom_v_adv_w, double *mom_w_adv, double *dmom_w_adv_u, double *dmom_w_adv_v, double *dmom_w_adv_w, double *mom_u_diff_ten, double *mom_v_diff_ten, double *mom_w_diff_ten, double *mom_uv_diff_ten, double *mom_uw_diff_ten, double *mom_vu_diff_ten, double *mom_vw_diff_ten, double *mom_wu_diff_ten, double *mom_wv_diff_ten, double *mom_u_source, double *mom_v_source, double *mom_w_source, double *mom_u_ham, double *dmom_u_ham_grad_p, double *mom_v_ham, double *dmom_v_ham_grad_p, double *mom_w_ham, double *dmom_w_ham_grad_p)
Definition: transportCoefficients.c:3279
f
Double f
Definition: Headers.h:64
unitSquareVortexEvaluate
void unitSquareVortexEvaluate(const int nPoints, const int nSpace, double t, const double *x, const double *u, double *m, double *dm, double *f, double *df)
Definition: transportCoefficients.c:4717
darcySharpInterfaceFlowImEvaluate
void darcySharpInterfaceFlowImEvaluate(int nPoints, int nSpace, double Km, double rhoM, double Kp, double rhoP, double eps, double *gravity_u, double *u, double *gradu, double *u_levelSet, double *phi_pot, double *a, double *f, double *r, double *m, double *dphi_pot, double *da, double *df, double *dr, double *dm)
Definition: transportCoefficients.c:1750
linearADR_ConstantCoefficientsEvaluate
void linearADR_ConstantCoefficientsEvaluate(const int nPoints, const int nSpace, const double M, const double *A, const double *B, const double C, const double t, const double *x, const double *u, double *m, double *dm, double *f, double *df, double *a, double *r, double *dr)
Definition: transportCoefficients.c:101
nonlinearADR_pqrstEvaluate
void nonlinearADR_pqrstEvaluate(const int nPoints, const int nSpace, const double M, const double *A, const double *B, const double C, const double p_pow, const double q_pow, const double r_pow, const double s_pow, const double t_pow, const double t, const double *x, const double *u, double *m, double *dm, double *f, double *df, double *a, double *da, double *phi, double *dphi, double *r, double *dr)
Definition: transportCoefficients.c:597
Stokes_2D_Evaluate
void Stokes_2D_Evaluate(const int nPoints, const double rho, const double nu, const double *g, const double *p, const double *grad_p, const double *u, const double *v, double *mom_u_acc, double *dmom_u_acc_u, double *mom_v_acc, double *dmom_v_acc_v, double *mass_adv, double *dmass_adv_u, double *dmass_adv_v, double *mom_u_diff_ten, double *mom_v_diff_ten, double *mom_u_source, double *mom_v_source, double *mom_u_ham, double *dmom_u_ham_grad_p, double *mom_v_ham, double *dmom_v_ham_grad_p)
Definition: transportCoefficients.c:2079
conservativeHeadRichardsL2projMualemVanGenuchtenHomEvaluate
void conservativeHeadRichardsL2projMualemVanGenuchtenHomEvaluate(const int nSimplices, const int nPointsPerSimplex, const int nSpace, const double rho, const double *gravity, const double alpha, const double n, const double m, const double thetaR, const double thetaSR, const double KWs, double *dV, double *u, double *mass, double *dmass, double *f, double *df, double *a, double *da)
Definition: transportCoefficients.c:5095
n
Int n
Definition: Headers.h:28
NavierStokes_3D_Evaluate
void NavierStokes_3D_Evaluate(const int nPoints, const double rho, const double nu, const double *g, const double *p, const double *grad_p, const double *u, const double *v, const double *w, double *mom_u_acc, double *dmom_u_acc_u, double *mom_v_acc, double *dmom_v_acc_v, double *mom_w_acc, double *dmom_w_acc_w, double *mass_adv, double *dmass_adv_u, double *dmass_adv_v, double *dmass_adv_w, double *mom_u_adv, double *dmom_u_adv_u, double *dmom_u_adv_v, double *dmom_u_adv_w, double *mom_v_adv, double *dmom_v_adv_u, double *dmom_v_adv_v, double *dmom_v_adv_w, double *mom_w_adv, double *dmom_w_adv_u, double *dmom_w_adv_v, double *dmom_w_adv_w, double *mom_u_diff_ten, double *mom_v_diff_ten, double *mom_w_diff_ten, double *mom_u_source, double *mom_v_source, double *mom_w_source, double *mom_u_ham, double *dmom_u_ham_grad_p, double *mom_v_ham, double *dmom_v_ham_grad_p, double *mom_w_ham, double *dmom_w_ham_grad_p)
Definition: transportCoefficients.c:1938
redistanceLevelSetCoefficientsEvaluate
void redistanceLevelSetCoefficientsEvaluate(int nPoints, int nSpace, double eps, double *u_levelSet, double *u, double *grad_u, double *m, double *dm, double *H, double *dH, double *r)
Definition: transportCoefficients.c:1384
df
double df(double C, double b, double a, int q, int r)
Definition: analyticalSolutions.c:2209
phi
Double phi
Definition: Headers.h:76
TwophaseDarcyFC_Evaluate
void TwophaseDarcyFC_Evaluate(const int nPoints, const int nSpace, const int pskModelFlag, const double Kbar, const double rhon, const double rhow, const double *g, const double *x, const double alpha, const double bc_lambda, const double bc_pd, const double mvg_n, const double mvg_m, const double omega, const double omega_r, const double mun, const double muw, const double b, double *sw, double *psiw, double *mw, double *dmw, double *mn, double *dmn, double *phi_psiw, double *dphi_psiw_dpsiw, double *phi_psin, double *dphi_psin_dpsiw, double *dphi_psin_dsw, double *fw, double *dfw, double *fn, double *dfn, double *aw, double *daw, double *an, double *dan)
constantVelocityLevelSetEvaluate
void constantVelocityLevelSetEvaluate(const int nPoints, const int nSpace, const double *b, const double *x, const double *u, const double *gradu, double *m, double *dm, double *f, double *df, double *H, double *dH)
Definition: transportCoefficients.c:4745
LinearElasticity_3D_Evaluate
void LinearElasticity_3D_Evaluate(const int nPoints, const double E, const double nu, const double *g, const double *u, const double *v, const double *w, double *uu_diff_ten, double *uv_diff_ten, double *uw_diff_ten, double *vu_diff_ten, double *vv_diff_ten, double *vw_diff_ten, double *wu_diff_ten, double *wv_diff_ten, double *ww_diff_ten, double *u_force, double *v_force, double *w_force)
Definition: transportCoefficients.c:8186
TwophaseFFDarcyFC_Evaluate
void TwophaseFFDarcyFC_Evaluate(const int nPoints, const int nSpace, const int pskModelFlag, const double Kbar, const double rhon, const double rhow, const double *g, const double *x, const double alpha, const double bc_lambda, const double bc_pd, const double mvg_n, const double mvg_m, const double omega, const double omega_r, const double mun, const double muw, const double b, double *sw, double *psiw, double *mw, double *dmw_dsw, double *mm, double *dmm_dsw, double *phi_psic, double *dphi_psic_dsw, double *phi_psiw, double *dphi_psiw_dpsiw, double *fm, double *dfm_dsw, double *aw_psiw, double *daw_psiw_dsw, double *am_psiw, double *dam_psiw_dsw, double *am_psic, double *dam_psic_dsw)
cLevelSetCoefficientsEvaluate
void cLevelSetCoefficientsEvaluate(int nPoints, int nSpace, double *v, double *u, double *m, double *dm, double *f, double *df)
Definition: transportCoefficients.c:1222
twophasePotentialFlowEvaluate
void twophasePotentialFlowEvaluate(int nPoints, int nSpace, double *M, double *A, double *B, double *Bcon, double *C, double t, double *x, double *u, double *m, double *dm, double *f, double *df, double *a, double *da, double *phi, double *dphi, double *r, double *dr)
Definition: transportCoefficients.c:1027
VolumeAveragedTwophaseNavierStokes_ST_LS_SO_2D_Evaluate
void VolumeAveragedTwophaseNavierStokes_ST_LS_SO_2D_Evaluate(const int nPoints, const int killNonlinearDrag, const double eps_rho, const double eps_mu, const double sigma, const double rho_0, const double nu_0, const double rho_1, const double nu_1, const double *meanGrainSize, const double *g, const double *phi, const double *n, const double *kappa, const double *p, const double *grad_p, const double *u, const double *v, const double *porosity, double *mom_u_acc, double *dmom_u_acc_u, double *mom_v_acc, double *dmom_v_acc_v, double *mass_adv, double *dmass_adv_u, double *dmass_adv_v, double *mom_u_adv, double *dmom_u_adv_u, double *dmom_u_adv_v, double *mom_v_adv, double *dmom_v_adv_u, double *dmom_v_adv_v, double *mom_u_diff_ten, double *mom_v_diff_ten, double *mom_uv_diff_ten, double *mom_vu_diff_ten, double *mom_u_source, double *mom_v_source, double *dmom_u_source_u, double *dmom_u_source_v, double *dmom_v_source_u, double *dmom_v_source_v, double *mom_u_ham, double *dmom_u_ham_grad_p, double *mom_v_ham, double *dmom_v_ham_grad_p)
Definition: transportCoefficients.c:8832
R
Double R
Definition: Headers.h:82
darcySharpInterfaceFlowEvaluate
void darcySharpInterfaceFlowEvaluate(int nPoints, int nSpace, double Km, double rhoM, double Kp, double rhoP, double eps, double *gravity_u, double *u, double *gradu, double *u_levelSet, double *phi_pot, double *a, double *f, double *r, double *m, double *dphi_pot, double *da, double *df, double *dr, double *dm)
Definition: transportCoefficients.c:1698
H
Double H
Definition: Headers.h:65
conservativeTotalHeadRichardsMualemVanGenuchtenHomEvaluate
void conservativeTotalHeadRichardsMualemVanGenuchtenHomEvaluate(const int nPoints, const int nSpace, const double rho, const double *gravity, const double *x, const double alpha, const double n, const double m, const double thetaR, const double thetaSR, const double KWs, double *u, double *mass, double *dmass, double *f, double *df, double *a, double *da, double *phi, double *dphi)
Definition: transportCoefficients.c:5550
HJBurgersEvaluate
void HJBurgersEvaluate(const int nPoints, const int nSpace, const double offset, const double *u, const double *gradu, double *m, double *dm, double *H, double *dH)
Definition: transportCoefficients.c:4885
nu_0
double nu_0
Definition: ErrorResidualMethod.cpp:22
levelSetConservationCoefficientsEvaluate
void levelSetConservationCoefficientsEvaluate(int nPoints, int nSpace, double epsHeaviside, double epsDirac, double epsDiffusion, double *u_ls, double *H_vof, double *u, double *r, double *dr, double *a)
Definition: transportCoefficients.c:8393
conservativeHeadRichardsMualemVanGenuchtenHetEvaluateV2
void conservativeHeadRichardsMualemVanGenuchtenHetEvaluateV2(const int nSimplex, const int nPointsPerSimplex, const int nSpace, const int *materialTypes, const double rho, const double beta, const double *gravity, const double *alpha, const double *n, const double *thetaR, const double *thetaSR, const double *KWs, double *u, double *mass, double *dmass, double *f, double *df, double *a, double *da)
Definition: transportCoefficients.c:5990
v
Double v
Definition: Headers.h:95
conservativeHeadRichardsBrooksCoreyBurdineHetEvaluate
void conservativeHeadRichardsBrooksCoreyBurdineHetEvaluate(const int nPoints, const int nSpace, const double rho, const double *gravity, const double *lambda, const double *pd, const double *thetaR, const double *thetaS, const double *KWs, double *u, double *mass, double *dmass, double *f, double *df, double *a, double *da)
Definition: transportCoefficients.c:6648
nu_1
double nu_1
Definition: ErrorResidualMethod.cpp:22
FractionalFlowPhaseForm_potentialHetEvaluate
void FractionalFlowPhaseForm_potentialHetEvaluate(const int nPoints, const int nSpace, const int nc, const int pskModelFlag, const double *Kbar, const double rhon, const double rhow, const double *g, const double *alpha, const double *bc_lambda, const double *bc_pd, const double *mvg_m, const double *thetaS, const double *thetaR, const double mun, const double muw, const double b, double *u, double *m, double *dm, double *phi, double *dphi, double *f, double *df, double *a, double *da, double *s_w, double *grad_psic)
TwophaseStokes_VOF_SO_3D_Evaluate
void TwophaseStokes_VOF_SO_3D_Evaluate(const int nPoints, const double eps, const double rho_0, const double nu_0, const double rho_1, const double nu_1, const double *g, const double *vof, const double *p, const double *grad_p, const double *u, const double *v, const double *w, double *mom_u_acc, double *dmom_u_acc_u, double *mom_v_acc, double *dmom_v_acc_v, double *mom_w_acc, double *dmom_w_acc_w, double *mass_adv, double *dmass_adv_u, double *dmass_adv_v, double *dmass_adv_w, double *mom_u_diff_ten, double *mom_v_diff_ten, double *mom_w_diff_ten, double *mom_u_source, double *mom_v_source, double *mom_w_source, double *mom_u_ham, double *dmom_u_ham_grad_p, double *mom_v_ham, double *dmom_v_ham_grad_p, double *mom_w_ham, double *dmom_w_ham_grad_p)
Definition: transportCoefficients.c:4620
twophaseLevelSetCoefficientsUpdateVelocity
void twophaseLevelSetCoefficientsUpdateVelocity(int nPoints, int nSpace, double v_scale, double *vIn, double *vOut)
Definition: transportCoefficients.c:1135
conservativeHeadRichardsL2projMualemVanGenuchtenHetEvaluate
void conservativeHeadRichardsL2projMualemVanGenuchtenHetEvaluate(const int nSimplices, const int nPointsPerSimplex, const int nSpace, const double rho, const double *gravity, const double *alpha, const double *n, const double *thetaR, const double *thetaSR, const double *KWs, double *dV, double *u, double *mass, double *dmass, double *f, double *df, double *a, double *da)
Definition: transportCoefficients.c:5412
unitSquareRotationEvaluate
void unitSquareRotationEvaluate(const int nPoints, const int nSpace, const double *x, const double *u, double *m, double *dm, double *f, double *df)
Definition: transportCoefficients.c:767
twophaseSignedDistanceCoefficientsUpdateSignFunction
void twophaseSignedDistanceCoefficientsUpdateSignFunction(int nPoints, double eps, double *u_levelSet, double *S)
Definition: transportCoefficients.c:1298
conservativeHeadRichardsBCBfromMVGHomEvaluate
void conservativeHeadRichardsBCBfromMVGHomEvaluate(const int nPoints, const int nSpace, const double rho, const double *gravity, const double alpha, const double n, const double m, const double thetaR, const double thetaSR, const double KWs, double *u, double *mass, double *dmass, double *f, double *df, double *a, double *da)
Definition: transportCoefficients.c:6864
TwophaseNavierStokes_VOF_SO_2D_Evaluate
void TwophaseNavierStokes_VOF_SO_2D_Evaluate(const int nPoints, const double eps, const double rho_0, const double nu_0, const double rho_1, const double nu_1, const double *g, const double *vof, const double *p, const double *grad_p, const double *u, const double *v, double *mom_u_acc, double *dmom_u_acc_u, double *mom_v_acc, double *dmom_v_acc_v, double *mass_adv, double *dmass_adv_u, double *dmass_adv_v, double *mom_u_adv, double *dmom_u_adv_u, double *dmom_u_adv_v, double *mom_v_adv, double *dmom_v_adv_u, double *dmom_v_adv_v, double *mom_u_diff_ten, double *mom_v_diff_ten, double *mom_u_source, double *mom_v_source, double *mom_u_ham, double *dmom_u_ham_grad_p, double *mom_v_ham, double *dmom_v_ham_grad_p)
Definition: transportCoefficients.c:4300
burgersDiagonalVelEvaluate
void burgersDiagonalVelEvaluate(const int nPoints, const int nSpace, const double self_a, const double *self_v, const double *u, double *m, double *dm, double *f, double *df, double *a, double *phi, double *dphi)
Definition: transportCoefficients.c:953
u
Double u
Definition: Headers.h:89
smoothedDirac
double smoothedDirac(double eps, double phi)
Definition: numericalFlux.c:14199
LinearElasticity_1D_Evaluate
void LinearElasticity_1D_Evaluate(const int nPoints, const double E, const double nu, const double *g, const double *u, double *uu_diff_ten, double *u_force)
Definition: transportCoefficients.c:8124
ThreephaseNavierStokes_ST_LS_SO_3D_Evaluate
void ThreephaseNavierStokes_ST_LS_SO_3D_Evaluate(const int nPoints, const double boundaryPenaltyCoef, const double volumePenaltyCoef, const double eps_rho, const double eps_mu, const double sigma, const double rho_0, const double nu_0, const double rho_1, const double nu_1, const double rho_s, const double nu_s, const double *g, const double *phi, const double *n, const double *kappa, const double *phi_s, const double *n_s, const double *p, const double *grad_p, const double *u, const double *v, const double *w, double *mom_u_acc, double *dmom_u_acc_u, double *mom_v_acc, double *dmom_v_acc_v, double *mom_w_acc, double *dmom_w_acc_w, double *mass_adv, double *dmass_adv_u, double *dmass_adv_v, double *dmass_adv_w, double *mom_u_adv, double *dmom_u_adv_u, double *dmom_u_adv_v, double *dmom_u_adv_w, double *mom_v_adv, double *dmom_v_adv_u, double *dmom_v_adv_v, double *dmom_v_adv_w, double *mom_w_adv, double *dmom_w_adv_u, double *dmom_w_adv_v, double *dmom_w_adv_w, double *mom_u_diff_ten, double *mom_v_diff_ten, double *mom_w_diff_ten, double *mom_uv_diff_ten, double *mom_uw_diff_ten, double *mom_vu_diff_ten, double *mom_vw_diff_ten, double *mom_wu_diff_ten, double *mom_wv_diff_ten, double *mom_u_source, double *dmom_u_source_u, double *dmom_u_source_v, double *dmom_u_source_w, double *mom_v_source, double *dmom_v_source_u, double *dmom_v_source_v, double *dmom_v_source_w, double *mom_w_source, double *dmom_w_source_u, double *dmom_w_source_v, double *dmom_w_source_w, double *mom_u_ham, double *dmom_u_ham_grad_p, double *mom_v_ham, double *dmom_v_ham_grad_p, double *mom_w_ham, double *dmom_w_ham_grad_p)
Definition: transportCoefficients.c:3749
ncLevelSetCoefficientsEvaluate
void ncLevelSetCoefficientsEvaluate(int nPoints, int nSpace, double *v, double *u, double *grad_u, double *m, double *dm, double *H, double *dH)
Definition: transportCoefficients.c:1198
TwophaseStokes_LS_SO_2D_Evaluate
void TwophaseStokes_LS_SO_2D_Evaluate(const int nPoints, const double eps, const double rho_0, const double nu_0, const double rho_1, const double nu_1, const double *g, const double *phi, const double *p, const double *grad_p, const double *u, const double *v, double *mom_u_acc, double *dmom_u_acc_u, double *mom_v_acc, double *dmom_v_acc_v, double *mass_adv, double *dmass_adv_u, double *dmass_adv_v, double *mom_u_diff_ten, double *mom_v_diff_ten, double *mom_u_source, double *mom_v_source, double *mom_u_ham, double *dmom_u_ham_grad_p, double *mom_v_ham, double *dmom_v_ham_grad_p)
Definition: transportCoefficients.c:4129
rho_1
double rho_1
Definition: ErrorResidualMethod.cpp:22
FractionalFlowPhaseForm_saturationHetEvaluate
void FractionalFlowPhaseForm_saturationHetEvaluate(const int nPoints, const int nSpace, const int nc, const int pskModelFlag, const double *Kbar, const double rhon, const double rhow, const double *g, const double *alpha, const double *bc_lambda, const double *bc_pd, const double *mvg_m, const double *thetaS, const double *thetaR, const double mun, const double muw, const double b, double *u, double *m, double *dm, double *phi, double *dphi, double *f, double *df, double *a, double *da, double *q_t)
TwophaseNavierStokes_ST_LS_SO_2D_Evaluate
void TwophaseNavierStokes_ST_LS_SO_2D_Evaluate(const int nPoints, const double eps_rho, const double eps_mu, const double sigma, const double rho_0, const double nu_0, const double rho_1, const double nu_1, const double *g, const double *phi, const double *n, const double *kappa, const double *p, const double *grad_p, const double *u, const double *v, double *mom_u_acc, double *dmom_u_acc_u, double *mom_v_acc, double *dmom_v_acc_v, double *mass_adv, double *dmass_adv_u, double *dmass_adv_v, double *mom_u_adv, double *dmom_u_adv_u, double *dmom_u_adv_v, double *mom_v_adv, double *dmom_v_adv_u, double *dmom_v_adv_v, double *mom_u_diff_ten, double *mom_v_diff_ten, double *mom_uv_diff_ten, double *mom_vu_diff_ten, double *mom_u_source, double *mom_v_source, double *mom_u_ham, double *dmom_u_ham_grad_p, double *mom_v_ham, double *dmom_v_ham_grad_p)
Definition: transportCoefficients.c:2481
VolumeAveragedNavierStokesFullDevStress_3D_Evaluate
void VolumeAveragedNavierStokesFullDevStress_3D_Evaluate(const int nPoints, const double rho, const double mu, const double *meanGrainSize, const double *g, const double *p, const double *grad_p, const double *u, const double *v, const double *w, const double *porosity, double *mom_u_acc, double *dmom_u_acc_u, double *mom_v_acc, double *dmom_v_acc_v, double *mom_w_acc, double *dmom_w_acc_w, double *mass_adv, double *dmass_adv_u, double *dmass_adv_v, double *dmass_adv_w, double *mom_u_adv, double *dmom_u_adv_u, double *dmom_u_adv_v, double *dmom_u_adv_w, double *mom_v_adv, double *dmom_v_adv_u, double *dmom_v_adv_v, double *dmom_v_adv_w, double *mom_w_adv, double *dmom_w_adv_u, double *dmom_w_adv_v, double *dmom_w_adv_w, double *mom_u_diff_ten, double *mom_v_diff_ten, double *mom_w_diff_ten, double *mom_uv_diff_ten, double *mom_uw_diff_ten, double *mom_vu_diff_ten, double *mom_vw_diff_ten, double *mom_wu_diff_ten, double *mom_wv_diff_ten, double *mom_u_source, double *mom_v_source, double *mom_w_source, double *dmom_u_source_u, double *dmom_u_source_v, double *dmom_u_source_w, double *dmom_v_source_u, double *dmom_v_source_v, double *dmom_v_source_w, double *dmom_w_source_u, double *dmom_w_source_v, double *dmom_w_source_w, double *mom_u_ham, double *dmom_u_ham_grad_p, double *mom_v_ham, double *dmom_v_ham_grad_p, double *mom_w_ham, double *dmom_w_ham_grad_p)
Definition: transportCoefficients.c:8630
LinearElasticity_2D_Evaluate
void LinearElasticity_2D_Evaluate(const int nPoints, const double E, const double nu, const double *g, const double *u, const double *v, double *uu_diff_ten, double *uv_diff_ten, double *vu_diff_ten, double *vv_diff_ten, double *u_force, double *v_force)
Definition: transportCoefficients.c:8143
conservativeHeadRichardsMualemVanGenuchtenHomEvaluate
void conservativeHeadRichardsMualemVanGenuchtenHomEvaluate(const int nPoints, const int nSpace, const double rho, const double beta, const double *gravity, const double *x, const double alpha, const double n, const double m, const double thetaR, const double thetaSR, const double KWs, double *u, double *mass, double *dmass, double *f, double *df, double *a, double *da, double *phi, double *dphi)
Definition: transportCoefficients.c:5007
unitCubeRotationEvaluate
void unitCubeRotationEvaluate(const int nPoints, const int nSpace, const double *x, const double *u, double *m, double *dm, double *f, double *df)
Definition: transportCoefficients.c:791
evaluateBuckleyLeverettLiuExample
void evaluateBuckleyLeverettLiuExample(int nPoints, int nSpace, const double *x, const double *u, double *m, double *dm, double *f, double *df, double *a)
Definition: transportCoefficients.c:8455
VolumeAveragedNavierStokesFullDevStress_2D_Evaluate
void VolumeAveragedNavierStokesFullDevStress_2D_Evaluate(const int nPoints, const double rho, const double mu, const double *meanGrainSize, const double *g, const double *p, const double *grad_p, const double *u, const double *v, const double *porosity, double *mom_u_acc, double *dmom_u_acc_u, double *mom_v_acc, double *dmom_v_acc_v, double *mass_adv, double *dmass_adv_u, double *dmass_adv_v, double *mom_u_adv, double *dmom_u_adv_u, double *dmom_u_adv_v, double *mom_v_adv, double *dmom_v_adv_u, double *dmom_v_adv_v, double *mom_u_diff_ten, double *mom_v_diff_ten, double *mom_uv_diff_ten, double *mom_vu_diff_ten, double *mom_u_source, double *mom_v_source, double *dmom_u_source_u, double *dmom_u_source_v, double *dmom_v_source_u, double *dmom_v_source_v, double *mom_u_ham, double *dmom_u_ham_grad_p, double *mom_v_ham, double *dmom_v_ham_grad_p)
Definition: transportCoefficients.c:8506
VOFCoefficientsEvaluate
void VOFCoefficientsEvaluate(int nPoints, int nSpace, double eps, double *v, double *phi, double *u, double *m, double *dm, double *f, double *df)
Definition: transportCoefficients.c:1244
conservativeHeadRichardsL2projBndMualemVanGenuchtenHomEvaluate
void conservativeHeadRichardsL2projBndMualemVanGenuchtenHomEvaluate(const int nElements, const int nElementBoundaries_element, const int nPointsPerElementBoundary, const int nSpace, const double rho, const double *gravity, const double alpha, const double n, const double m, const double thetaR, const double thetaSR, const double KWs, double *dV, double *u, double *mass, double *dmass, double *f, double *df, double *a, double *da)
Definition: transportCoefficients.c:5243
VolumeAveragedVOFCoefficientsEvaluate
void VolumeAveragedVOFCoefficientsEvaluate(int nPoints, int nSpace, double eps, double *v, double *phi, double *porosity, double *u, double *m, double *dm, double *f, double *df)
Definition: transportCoefficients.c:9890
pd
#define pd(x)
Definition: jf.h:24
eikonalEquationEvaluate
void eikonalEquationEvaluate(int nPoints, int nSpace, double rhs, double *u, double *grad_u, double *m, double *dm, double *H, double *dH, double *r)
Definition: transportCoefficients.c:1354
conservativeSatRichardsBrooksCoreyBurdineHomEvaluate
void conservativeSatRichardsBrooksCoreyBurdineHomEvaluate(const int nPoints, const int nSpace, const double rho, const double *gravity, const double lambda, const double pd, const double thetaR, const double thetaSR, const double KWs, double *u, double *mass, double *dmass, double *f, double *df, double *a, double *da, double *phi, double *dphi)
Definition: transportCoefficients.c:6798
conservativeHeadRichardsBrooksCoreyBurdineHomEvaluate
void conservativeHeadRichardsBrooksCoreyBurdineHomEvaluate(const int nPoints, const int nSpace, const double rho, const double beta, const double *gravity, const double lambda, const double pd, const double thetaR, const double thetaSR, const double KWs, double *u, double *mass, double *dmass, double *f, double *df, double *a, double *da)
Definition: transportCoefficients.c:6706
twophasePotentialFlowUpdateFreeSurface
void twophasePotentialFlowUpdateFreeSurface(int nPoints, int nSpace, double eps, double *u_levelSet, double M1, double M2, double *M, double *A1, double *A2, double *A, double *B1, double *B2, double *B, double *Bcon1, double *Bcon2, double *Bcon, double C1, double C2, double *C)
Definition: transportCoefficients.c:1073
r
Double r
Definition: Headers.h:83
groundwaterTransportCoefficientsEvaluate
void groundwaterTransportCoefficientsEvaluate(const int nPoints, const int nSpace, const double omega, const double d, const double alpha_L, const double alpha_T, const double *v, const double *u, double *m, double *dm, double *f, double *df, double *a)
Definition: transportCoefficients.c:140
FractionalFlowPhaseForm_potentialEvaluate
void FractionalFlowPhaseForm_potentialEvaluate(const int nPoints, const int nSpace, const int nc, const int pskModelFlag, const double Kbar, const double rhon, const double rhow, const double *g, const double g_norm, const double alpha, const double bc_lambda, const double bc_pd, const double mvg_n, const double mvg_m, const double omega, const double mun, const double muw, const double sw_min, const double sw_max, const double M, const double R, const double Temp, const double p_o, const double b, double *u, double *m, double *dm, double *phi, double *dphi, double *f, double *df, double *a, double *da, double *s_w, double *grad_psic)
burgersDiagonalVelHJEvaluate
void burgersDiagonalVelHJEvaluate(const int nPoints, const int nSpace, const double self_a, const double *self_v, const double *u, const double *grad_u, double *m, double *dm, double *H, double *dH, double *a, double *phi, double *dphi)
Definition: transportCoefficients.c:989
nonlinearADR_pqrstDualEvaluate
void nonlinearADR_pqrstDualEvaluate(const int nPoints, const int nSpace, const double M, const double *A, const double *B, const double C, const double p1, const double q1, const double r1, const double s1, const double t1, const double p2, const double q2, const double r2, const double s2, const double t2, const double t, const double *x, const double *u, double *m, double *dm, double *f, double *df, double *a, double *da, double *phi, double *dphi, double *r, double *dr)
Definition: transportCoefficients.c:664
l2projectScalar
void l2projectScalar(const int nSimplices, const int nPointsPerSimplex, double *dV, double *r)
Definition: transportCoefficients.c:5651
smoothedHeaviside_integral
double smoothedHeaviside_integral(double eps, double phi)
Definition: numericalFlux.c:14181
FractionalFlowPhaseForm_saturationEvaluate
void FractionalFlowPhaseForm_saturationEvaluate(const int nPoints, const int nSpace, const int nc, const int pskModelFlag, const double Kbar, const double rhon, const double rhow, const double *g, const double g_norm, const double alpha, const double bc_lambda, const double bc_pd, const double mvg_n, const double mvg_m, const double omega, const double mun, const double muw, const double sw_min, const double sw_max, const double M, const double R, const double Temp, const double p_o, const double b, double *u, double *m, double *dm, double *phi, double *dphi, double *f, double *df, double *a, double *da, double *q_t, double *psiw)
disVelEvaluate
void disVelEvaluate(const int nPoints, const int nSpace, const double self_a, const double *x, const double *u, double *m, double *dm, double *f, double *df, double *a, double *da, double *phi, double *dphi)
Definition: transportCoefficients.c:908
TwophaseStokes_VOF_SO_2D_Evaluate
void TwophaseStokes_VOF_SO_2D_Evaluate(const int nPoints, const double eps, const double rho_0, const double nu_0, const double rho_1, const double nu_1, const double *g, const double *vof, const double *p, const double *grad_p, const double *u, const double *v, double *mom_u_acc, double *dmom_u_acc_u, double *mom_v_acc, double *dmom_v_acc_v, double *mass_adv, double *dmass_adv_u, double *dmass_adv_v, double *mom_u_diff_ten, double *mom_v_diff_ten, double *mom_u_source, double *mom_v_source, double *mom_u_ham, double *dmom_u_ham_grad_p, double *mom_v_ham, double *dmom_v_ham_grad_p)
Definition: transportCoefficients.c:4546
rho_0
double rho_0
Definition: ErrorResidualMethod.cpp:22
disRotatingPulseVelEvaluate
void disRotatingPulseVelEvaluate(const int nPoints, const int nSpace, const double self_a, const double *x, const double *u, double *m, double *dm, double *f, double *df, double *a, double *da, double *phi, double *dphi)
Definition: transportCoefficients.c:859
constantNormalVelocityLevelSetEvaluate
void constantNormalVelocityLevelSetEvaluate(const int nPoints, const int nSpace, double b, const double *x, const double *u, const double *gradu, double *m, double *dm, double *f, double *df, double *H, double *dH)
Definition: transportCoefficients.c:4776
TwophaseNavierStokes_VOF_SO_3D_Evaluate
void TwophaseNavierStokes_VOF_SO_3D_Evaluate(const int nPoints, const double eps, const double rho_0, const double nu_0, const double rho_1, const double nu_1, const double *g, const double *vof, const double *p, const double *grad_p, const double *u, const double *v, const double *w, double *mom_u_acc, double *dmom_u_acc_u, double *mom_v_acc, double *dmom_v_acc_v, double *mom_w_acc, double *dmom_w_acc_w, double *mass_adv, double *dmass_adv_u, double *dmass_adv_v, double *dmass_adv_w, double *mom_u_adv, double *dmom_u_adv_u, double *dmom_u_adv_v, double *dmom_u_adv_w, double *mom_v_adv, double *dmom_v_adv_u, double *dmom_v_adv_v, double *dmom_v_adv_w, double *mom_w_adv, double *dmom_w_adv_u, double *dmom_w_adv_v, double *dmom_w_adv_w, double *mom_u_diff_ten, double *mom_v_diff_ten, double *mom_w_diff_ten, double *mom_u_source, double *mom_v_source, double *mom_w_source, double *mom_u_ham, double *dmom_u_ham_grad_p, double *mom_v_ham, double *dmom_v_ham_grad_p, double *mom_w_ham, double *dmom_w_ham_grad_p)
Definition: transportCoefficients.c:4398
unitSquareRotationLevelSetEvaluate
void unitSquareRotationLevelSetEvaluate(const int nPoints, const int nSpace, double t, const double *x, const double *u, const double *gradu, double *m, double *dm, double *f, double *df, double *H, double *dH)
Definition: transportCoefficients.c:4848
twophaseSignedDistanceCoefficientsEvaluate
void twophaseSignedDistanceCoefficientsEvaluate(int nPoints, int nSpace, double *S, double *u, double *grad_u, double *m, double *dm, double *h, double *dh, double *rh)
Definition: transportCoefficients.c:1325
levelSetCurvatureCoefficientsEvaluate
void levelSetCurvatureCoefficientsEvaluate(int nPoints, int nSpace, double *grad_phi, double *u, double *f, double *r, double *dr)
Definition: transportCoefficients.c:1270
twophaseLevelSetCoefficientsEvaluateCI
void twophaseLevelSetCoefficientsEvaluateCI(int nPoints, int nSpace, double *B, double t, double *x, double *u, double *m, double *dm, double *f, double *df, double *a, double *da, double *phi, double *dphi, double *r, double *dr)
Definition: transportCoefficients.c:1173
unitSquareVortexLevelSetEvaluate
void unitSquareVortexLevelSetEvaluate(const int nPoints, const int nSpace, double t, const double *x, const double *u, const double *gradu, double *m, double *dm, double *f, double *df, double *H, double *dH)
Definition: transportCoefficients.c:4812
Generated on Fri Jul 1 2022 11:17:46 for proteus by doxygen 1.8.20