proteus 1.9.0
C/C++/Fortran libraries
Loading...
Searching...
No Matches
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
15extern 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
29extern 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);
86extern 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);
96extern 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
110extern 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);
129extern 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);
167extern double smoothedHeaviside(
168 double eps,
169 double phi
170);
171extern double smoothedHeaviside_integral(
172 double eps,
173 double phi
174);
175extern double smoothedDirac(
176 double eps,
177 double phi
178);
179extern 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);
197extern 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);
223extern 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);
254extern 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);
264extern 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);
279extern 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);
294extern 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);
309extern 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);
323extern 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);
345extern 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);
366extern void twophaseLevelSetCoefficientsUpdateVelocity(
367 int nPoints,
368 int nSpace,
369 double v_scale,
370 double *vIn,
371 double *vOut
372);
373extern 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);
387extern 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);
405extern 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);
416extern 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);
426extern 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);
435extern void twophaseSignedDistanceCoefficientsUpdateSignFunction(
436 int nPoints,
437 double eps,
438 double *u_levelSet,
439 double *S
440);
441extern 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);
453extern 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);
466extern 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);
489extern 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);
512extern 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);
535extern 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);
556extern 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);
578extern 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);
598extern 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);
620extern void l2projectScalar(
621 const int nSimplices,
622 const int nPointsPerSimplex,
623 double *dV,
624 double *r
625);
626extern void l2projectVector(
627 const int nSimplices,
628 const int nPointsPerSimplex,
629 const int nSpace,
630 double *dV,
631 double *r
632);
633extern void l2project2Tensor(
634 const int nSimplices,
635 const int nPointsPerSimplex,
636 const int nSpace,
637 double *dV,
638 double *r
639);
640extern 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);
658extern 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);
680extern 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);
698extern 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);
717extern 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);
737extern 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);
756extern 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);
787extern 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);
832extern 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);
857extern 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);
881extern 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);
914extern 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);
949extern 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);
990extern 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);
1039extern 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);
1068extern 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);
1105extern 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);
1140extern 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);
1189extern 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);
1218extern 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);
1255extern 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
1294extern 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
1355extern
1356void 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
1403extern
1404void 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);
1415extern 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);
1426extern 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);
1440extern 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);
1454extern 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);
1468extern 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);
1482extern 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 );
1518extern 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);
1553extern 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);
1580extern 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);
1608extern 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);
1647extern 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);
1685extern 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);
1723extern 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);
1760extern 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);
1769extern 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);
1783extern 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);
1804extern 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
1817extern 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
1826extern 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
1837extern 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
1848extern 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
1862extern 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);
1871extern 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
1910extern 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
1946extern
1947void 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); */
2029extern void groundwaterBiodegradation01EvaluateFC(const int nPoints,
2030 const int nSpace,
2031 const double omega,
2032 const double d_c,
2033 const double d_e,
2034 const double alpha_L,
2035 const double alpha_T,
2036 const double Kox_max,
2037 const double Kox_C,
2038 const double Kox_E,
2039 const double Kox_X,
2040 const double Yield,
2041 const double k_d,
2042 const double *v,
2043 const double *c_c,
2044 const double *c_e,
2045 const double *c_x,
2046 double *m_c,
2047 double *dm_c,
2048 double *m_e,
2049 double *dm_e,
2050 double *m_x,
2051 double *dm_x,
2052 double *f_c,
2053 double *df_c,
2054 double *f_e,
2055 double *df_e,
2056 double *a_c,
2057 double *a_e,
2058 double *r_c,
2059 double *dr_c_dc,
2060 double *dr_c_de,
2061 double *dr_c_dx,
2062 double *r_e,
2063 double *dr_e_dc,
2064 double *dr_e_de,
2065 double *dr_e_dx,
2066 double *r_x,
2067 double *dr_x_dc,
2068 double *dr_x_de,
2069 double *dr_x_dx);
2070extern void groundwaterBryantDawsonIonExEvaluateFC(const int nPoints,
2071 const int nSpace,
2072 const double omega,
2073 const double d_m,
2074 const double d_h,
2075 const double alpha_L,
2076 const double alpha_T,
2077 const double K_m,
2078 const double K_h,
2079 const double K_w,
2080 const double Z_tot,
2081 const double *v,
2082 const double *c_m,
2083 const double *c_h,
2084 double *m_m,
2085 double *dm_m_m,
2086 double *dm_m_h,
2087 double *m_h,
2088 double *dm_h_m,
2089 double *dm_h_h,
2090 double *f_m,
2091 double *df_m,
2092 double *f_h,
2093 double *df_h,
2094 double *a_m,
2095 double *a_h,
2096 double *phi_h,
2097 double *dphi_h,
2098 double *r_m,
2099 double *dr_m_dm,
2100 double *dr_m_dh,
2101 double *r_h,
2102 double *dr_h_dm,
2103 double *dr_h_dh);
2104
2105extern void groundwaterTransportCoefficientsEvaluate_hetMat(const int nSimplex,
2106 const int nPointsPerSimplex,
2107 const int nSpace,
2108 const double d,
2109 const int* materialTypes,
2110 const double *omega_types,
2111 const double *alpha_L_types,
2112 const double *alpha_T_types,
2113 const double *v,
2114 const double *u,
2115 double *m,
2116 double *dm,
2117 double *f,
2118 double *df,
2119 double *a);
2120
2121extern void variablySaturatedGroundwaterTransportCoefficientsEvaluate_hetMat(const int nSimplex,
2122 const int nPointsPerSimplex,
2123 const int nSpace,
2124 const double d,
2125 const int* materialTypes,
2126 const double *theta, /*phase volume fraction*/
2127 const double *alpha_L_types,
2128 const double *alpha_T_types,
2129 const double *v,/*phase darcy velocity*/
2130 const double *u,
2131 double *m,
2132 double *dm,
2133 double *f,
2134 double *df,
2135 double *a);
2136
2137extern void redistanceLevelSetCoefficientsWithWeakPenaltyEvaluate(int nPoints,
2138 int nSpace,
2139 double eps,
2140 double lambda_penalty,
2141 double* u_levelSet,
2142 double* u,
2143 double* grad_u,
2144 double* m,
2145 double* dm,
2146 double* H,
2147 double* dH,
2148 double* r,
2149 double* dr);
2150
2151extern void setWeakDirichletConditionsForLevelSet(int nElements_global,
2152 int nDOF_trial_element,
2153 double epsilon_freeze_factor,
2154 const double *elementDiameter,
2155 const int * u_l2g,
2156 const double *u_dof,
2157 int * freeze_nodes_tmp,
2158 int * weakDirichletConditionFlags);
2159
2160extern void StokesP_3D_Evaluate(const int nPoints,
2161 const double rho,
2162 const double nu,
2163 const double *g,
2164 const double *p,
2165 const double *u,
2166 const double *v,
2167 const double *w,
2168 double *mom_u_acc,
2169 double *dmom_u_acc_u,
2170 double *mom_v_acc,
2171 double *dmom_v_acc_v,
2172 double *mom_w_acc,
2173 double *dmom_w_acc_w,
2174 double *mass_adv,
2175 double *dmass_adv_u,
2176 double *dmass_adv_v,
2177 double *dmass_adv_w,
2178 double *mom_u_adv,
2179 double *dmom_u_adv_p,
2180 double *mom_v_adv,
2181 double *dmom_v_adv_p,
2182 double *mom_w_adv,
2183 double *dmom_w_adv_p,
2184 double *mom_u_diff_ten,
2185 double *mom_v_diff_ten,
2186 double *mom_w_diff_ten,
2187 double *mom_u_source,
2188 double *mom_v_source,
2189 double *mom_w_source);
2190
2191extern void TwophaseNavierStokes_ST_LS_SO_2D_Evaluate_sd(const int nPoints,
2192 const double eps_rho,
2193 const double eps_mu,
2194 const double sigma,
2195 const double rho_0,
2196 const double nu_0,
2197 const double rho_1,
2198 const double nu_1,
2199 const double* g,
2200 const double* phi,
2201 const double* n,
2202 const double* kappa,
2203 const double *p,
2204 const double *grad_p,
2205 const double *u,
2206 const double *v,
2207 double *mom_u_acc,
2208 double *dmom_u_acc_u,
2209 double *mom_v_acc,
2210 double *dmom_v_acc_v,
2211 double *mass_adv,
2212 double *dmass_adv_u,
2213 double *dmass_adv_v,
2214 double *mom_u_adv,
2215 double *dmom_u_adv_u,
2216 double *dmom_u_adv_v,
2217 double *mom_v_adv,
2218 double *dmom_v_adv_u,
2219 double *dmom_v_adv_v,
2220 double *mom_u_diff_ten,
2221 double *mom_v_diff_ten,
2222 double *mom_uv_diff_ten,
2223 double *mom_vu_diff_ten,
2224 double *mom_u_source,
2225 double *mom_v_source,
2226 double *mom_u_ham,
2227 double *dmom_u_ham_grad_p,
2228 double *mom_v_ham,
2229 double *dmom_v_ham_grad_p);
2230
2231extern void ThreephaseNavierStokes_ST_LS_SO_2D_Evaluate(const int nPoints,
2232 const double boundaryPenaltyCoef,
2233 const double volumePenaltyCoef,
2234 const double eps_rho,
2235 const double eps_mu,
2236 const double sigma,
2237 const double rho_0,
2238 const double nu_0,
2239 const double rho_1,
2240 const double nu_1,
2241 const double rho_s,
2242 const double nu_s,
2243 const double* g,
2244 const double* phi,
2245 const double* n,
2246 const double* kappa,
2247 const double* phi_s,
2248 const double* n_s,
2249 const double *p,
2250 const double *grad_p,
2251 const double *u,
2252 const double *v,
2253 double *mom_u_acc,
2254 double *dmom_u_acc_u,
2255 double *mom_v_acc,
2256 double *dmom_v_acc_v,
2257 double *mass_adv,
2258 double *dmass_adv_u,
2259 double *dmass_adv_v,
2260 double *mom_u_adv,
2261 double *dmom_u_adv_u,
2262 double *dmom_u_adv_v,
2263 double *mom_v_adv,
2264 double *dmom_v_adv_u,
2265 double *dmom_v_adv_v,
2266 double *mom_u_diff_ten,
2267 double *mom_v_diff_ten,
2268 double *mom_uv_diff_ten,
2269 double *mom_vu_diff_ten,
2270 double *mom_u_source,
2271 double *dmom_u_source_u,
2272 double *dmom_u_source_v,
2273 double *mom_v_source,
2274 double *dmom_v_source_u,
2275 double *dmom_v_source_v,
2276 double *mom_u_ham,
2277 double *dmom_u_ham_grad_p,
2278 double *mom_v_ham,
2279 double *dmom_v_ham_grad_p);
2280void TwophaseNavierStokes_ST_LS_SO_3D_Evaluate_sd(const int nPoints,
2281 const double eps_rho,
2282 const double eps_mu,
2283 const double sigma,
2284 const double rho_0,
2285 const double nu_0,
2286 const double rho_1,
2287 const double nu_1,
2288 const double* g,
2289 const double* phi,
2290 const double* n,
2291 const double* kappa,
2292 const double *p,
2293 const double *grad_p,
2294 const double *u,
2295 const double *v,
2296 const double *w,
2297 double *mom_u_acc,
2298 double *dmom_u_acc_u,
2299 double *mom_v_acc,
2300 double *dmom_v_acc_v,
2301 double *mom_w_acc,
2302 double *dmom_w_acc_w,
2303 double *mass_adv,
2304 double *dmass_adv_u,
2305 double *dmass_adv_v,
2306 double *dmass_adv_w,
2307 double *mom_u_adv,
2308 double *dmom_u_adv_u,
2309 double *dmom_u_adv_v,
2310 double *dmom_u_adv_w,
2311 double *mom_v_adv,
2312 double *dmom_v_adv_u,
2313 double *dmom_v_adv_v,
2314 double *dmom_v_adv_w,
2315 double *mom_w_adv,
2316 double *dmom_w_adv_u,
2317 double *dmom_w_adv_v,
2318 double *dmom_w_adv_w,
2319 double *mom_u_diff_ten,
2320 double *mom_v_diff_ten,
2321 double *mom_w_diff_ten,
2322 double *mom_uv_diff_ten,
2323 double *mom_uw_diff_ten,
2324 double *mom_vu_diff_ten,
2325 double *mom_vw_diff_ten,
2326 double *mom_wu_diff_ten,
2327 double *mom_wv_diff_ten,
2328 double *mom_u_source,
2329 double *mom_v_source,
2330 double *mom_w_source,
2331 double *mom_u_ham,
2332 double *dmom_u_ham_grad_p,
2333 double *mom_v_ham,
2334 double *dmom_v_ham_grad_p,
2335 double *mom_w_ham,
2336 double *dmom_w_ham_grad_p);
2337extern void conservativeHeadRichardsMualemVanGenuchten_sd_het_linearized_at_saturation(const int nSimplex,
2338 const int nPointsPerSimplex,
2339 const int nSpace,
2340 double linear_break,
2341 const int * rowptr,
2342 const int * colind,
2343 const int* materialTypes,
2344 const double rho,
2345 const double beta,
2346 const double* gravity,
2347 const double* alpha,
2348 const double* n,
2349 const double* thetaR,
2350 const double* thetaSR,
2351 const double* KWs,
2352 double *u,
2353 double *mass,
2354 double *dmass,
2355 double *f,
2356 double *df,
2357 double *a,
2358 double *da,
2359 double* vol_frac);
2360extern void conservativeHeadRichardsMualemVanGenuchten_sd_het(const int nSimplex,
2361 const int nPointsPerSimplex,
2362 const int nSpace,
2363 double pc_eps,
2364 const int * rowptr,
2365 const int * colind,
2366 const int* materialTypes,
2367 const double rho,
2368 const double beta,
2369 const double* gravity,
2370 const double* alpha,
2371 const double* n,
2372 const double* thetaR,
2373 const double* thetaSR,
2374 const double* KWs,
2375 double *u,
2376 double *mass,
2377 double *dmass,
2378 double *f,
2379 double *df,
2380 double *a,
2381 double *da,
2382 double* vol_frac);
2383extern void seepageBrezis(const int nSimplex,
2384 const int nPointsPerSimplex,
2385 const int nSpace,
2386 const int* materialTypes,
2387 const double epsFact,
2388 const double rho,
2389 const double beta,
2390 const double* elementDiameter,
2391 const double* gravity,
2392 const double* alpha,
2393 const double* n,
2394 const double* thetaR,
2395 const double* thetaSR,
2396 const double* KWs,
2397 double *u,
2398 double *mass,
2399 double *dmass,
2400 double *f,
2401 double *df,
2402 double *a,
2403 double *da);
2404extern void MovingMesh_1D_Evaluate(const int nPoints,
2405 const double E0,
2406 const double nu,
2407 const double *g,
2408 const double *det_J,
2409 const double *u,
2410 double *uu_diff_ten,
2411 double *u_force);
2412extern void MovingMesh_2D_Evaluate(const int nPoints,
2413 const double E0,
2414 const double nu,
2415 const double *g,
2416 const double *det_J,
2417 const double *u,
2418 const double *v,
2419 double *uu_diff_ten,double *uv_diff_ten,
2420 double *vu_diff_ten,double *vv_diff_ten,
2421 double *u_force,
2422 double *v_force);
2423extern void MovingMesh_3D_Evaluate(const int nPoints,
2424 const double E0,
2425 const double nu,
2426 const double *g,
2427 const double *det_J,
2428 const double *u,
2429 const double *v,
2430 const double *w,
2431 double *uu_diff_ten,double *uv_diff_ten,double *uw_diff_ten,
2432 double *vu_diff_ten,double *vv_diff_ten,double *vw_diff_ten,
2433 double *wu_diff_ten,double *wv_diff_ten,double *ww_diff_ten,
2434 double *u_force,
2435 double *v_force,
2436 double *w_force);
2437
2438extern void levelSetConservationCoefficientsEvaluate_sd(int nPoints,
2439 double epsHeaviside,
2440 double epsDirac,
2441 double* u_ls,
2442 double* H_vof,
2443 double* u,
2444 double* r,
2445 double* dr);
2446extern void variablySaturatedGroundwaterEnergyTransportCoefficientsEvaluate_hetMat(const int nSimplex,
2447 const int nPointsPerSimplex,
2448 const int nSpace,
2449 const double rho_w,
2450 const double rho_n,
2451 const double specificHeat_w,
2452 const double specificHeat_n,
2453 const int* materialTypes,
2454 const double *theta, /*phase volume fraction*/
2455 const double *thetaS_types,
2456 const double *alpha_L_types,
2457 const double *alpha_T_types,
2458 const double *rho_s_types,
2459 const double *specificHeat_s_types,
2460 const double *lambda_sat_types,
2461 const double *lambda_dry_types,
2462 const double *lambda_aniso_types,
2463 const double *v,/*phase darcy velocity*/
2464 const double *u,
2465 double *m,
2466 double *dm,
2467 double *f,
2468 double *df,
2469 double *a);
2470
2471extern void setSimpleWeakDirichletConditionsForLevelSet(int nElements_global,
2472 int nDOF_trial_element,
2473 double epsilon_freeze_factor,
2474 const double *elementDiameter,
2475 const int * u_l2g,
2476 const double *u_dof,
2477 int * freeze_nodes_tmp,
2478 int * weakDirichletConditionFlags);
2479extern void VolumeAveragedTwophaseNavierStokes_ST_LS_SO_2D_Evaluate(const int nPoints,
2480 const int killNonlinearDrag,
2481 const double eps_rho,
2482 const double eps_mu,
2483 const double sigma,
2484 const double rho_0,
2485 const double nu_0,
2486 const double rho_1,
2487 const double nu_1,
2488 const double* meanGrainSize,
2489 const double* g,
2490 const double* phi,
2491 const double* n,
2492 const double* kappa,
2493 const double *p,
2494 const double *grad_p,
2495 const double *u,
2496 const double *v,
2497 const double *porosity,
2498 double *mom_u_acc,
2499 double *dmom_u_acc_u,
2500 double *mom_v_acc,
2501 double *dmom_v_acc_v,
2502 double *mass_adv,
2503 double *dmass_adv_u,
2504 double *dmass_adv_v,
2505 double *mom_u_adv,
2506 double *dmom_u_adv_u,
2507 double *dmom_u_adv_v,
2508 double *mom_v_adv,
2509 double *dmom_v_adv_u,
2510 double *dmom_v_adv_v,
2511 double *mom_u_diff_ten,
2512 double *mom_v_diff_ten,
2513 double *mom_uv_diff_ten,
2514 double *mom_vu_diff_ten,
2515 double *mom_u_source,
2516 double *mom_v_source,
2517 double *dmom_u_source_u,
2518 double *dmom_u_source_v,
2519 double *dmom_v_source_u,
2520 double *dmom_v_source_v,
2521 double *mom_u_ham,
2522 double *dmom_u_ham_grad_p,
2523 double *mom_v_ham,
2524 double *dmom_v_ham_grad_p);
2525extern void VolumeAveragedTwophaseNavierStokes_ST_LS_SO_3D_Evaluate(const int nPoints,
2526 const int killNonlinearDrag,
2527 const double eps_rho,
2528 const double eps_mu,
2529 const double sigma,
2530 const double rho_0,
2531 const double nu_0,
2532 const double rho_1,
2533 const double nu_1,
2534 const double* meanGrainSize,
2535 const double* g,
2536 const double* phi,
2537 const double* n,
2538 const double* kappa,
2539 const double *p,
2540 const double *grad_p,
2541 const double *u,
2542 const double *v,
2543 const double *w,
2544 const double *porosity,
2545 double *mom_u_acc,
2546 double *dmom_u_acc_u,
2547 double *mom_v_acc,
2548 double *dmom_v_acc_v,
2549 double *mom_w_acc,
2550 double *dmom_w_acc_w,
2551 double *mass_adv,
2552 double *dmass_adv_u,
2553 double *dmass_adv_v,
2554 double *dmass_adv_w,
2555 double *mom_u_adv,
2556 double *dmom_u_adv_u,
2557 double *dmom_u_adv_v,
2558 double *dmom_u_adv_w,
2559 double *mom_v_adv,
2560 double *dmom_v_adv_u,
2561 double *dmom_v_adv_v,
2562 double *dmom_v_adv_w,
2563 double *mom_w_adv,
2564 double *dmom_w_adv_u,
2565 double *dmom_w_adv_v,
2566 double *dmom_w_adv_w,
2567 double *mom_u_diff_ten,
2568 double *mom_v_diff_ten,
2569 double *mom_w_diff_ten,
2570 double *mom_uv_diff_ten,
2571 double *mom_uw_diff_ten,
2572 double *mom_vu_diff_ten,
2573 double *mom_vw_diff_ten,
2574 double *mom_wu_diff_ten,
2575 double *mom_wv_diff_ten,
2576 double *mom_u_source,
2577 double *mom_v_source,
2578 double *mom_w_source,
2579 double *dmom_u_source_u,
2580 double *dmom_u_source_v,
2581 double *dmom_u_source_w,
2582 double *dmom_v_source_u,
2583 double *dmom_v_source_v,
2584 double *dmom_v_source_w,
2585 double *dmom_w_source_u,
2586 double *dmom_w_source_v,
2587 double *dmom_w_source_w,
2588 double *mom_u_ham,
2589 double *dmom_u_ham_grad_p,
2590 double *mom_v_ham,
2591 double *dmom_v_ham_grad_p,
2592 double *mom_w_ham,
2593 double *dmom_w_ham_grad_p);
2594extern void VolumeAveragedTwophaseNavierStokes_ST_LS_SO_3D_Evaluate_sd(const int nPoints,
2595 const int killNonlinearDrag,
2596 const double eps_rho,
2597 const double eps_mu,
2598 const double sigma,
2599 const double rho_0,
2600 const double nu_0,
2601 const double rho_1,
2602 const double nu_1,
2603 const double* meanGrainSize,
2604 const double* g,
2605 const double* phi,
2606 const double* n,
2607 const double* kappa,
2608 const double *p,
2609 const double *grad_p,
2610 const double *u,
2611 const double *v,
2612 const double *w,
2613 const double *porosity,
2614 double *mom_u_acc,
2615 double *dmom_u_acc_u,
2616 double *mom_v_acc,
2617 double *dmom_v_acc_v,
2618 double *mom_w_acc,
2619 double *dmom_w_acc_w,
2620 double *mass_adv,
2621 double *dmass_adv_u,
2622 double *dmass_adv_v,
2623 double *dmass_adv_w,
2624 double *mom_u_adv,
2625 double *dmom_u_adv_u,
2626 double *dmom_u_adv_v,
2627 double *dmom_u_adv_w,
2628 double *mom_v_adv,
2629 double *dmom_v_adv_u,
2630 double *dmom_v_adv_v,
2631 double *dmom_v_adv_w,
2632 double *mom_w_adv,
2633 double *dmom_w_adv_u,
2634 double *dmom_w_adv_v,
2635 double *dmom_w_adv_w,
2636 double *mom_u_diff_ten,
2637 double *mom_v_diff_ten,
2638 double *mom_w_diff_ten,
2639 double *mom_uv_diff_ten,
2640 double *mom_uw_diff_ten,
2641 double *mom_vu_diff_ten,
2642 double *mom_vw_diff_ten,
2643 double *mom_wu_diff_ten,
2644 double *mom_wv_diff_ten,
2645 double *mom_u_source,
2646 double *mom_v_source,
2647 double *mom_w_source,
2648 double *dmom_u_source_u,
2649 double *dmom_u_source_v,
2650 double *dmom_u_source_w,
2651 double *dmom_v_source_u,
2652 double *dmom_v_source_v,
2653 double *dmom_v_source_w,
2654 double *dmom_w_source_u,
2655 double *dmom_w_source_v,
2656 double *dmom_w_source_w,
2657 double *mom_u_ham,
2658 double *dmom_u_ham_grad_p,
2659 double *mom_v_ham,
2660 double *dmom_v_ham_grad_p,
2661 double *mom_w_ham,
2662 double *dmom_w_ham_grad_p);
2663extern void kEpsilon_2D_Evaluate(int nPoints,
2664 int nSpace,
2665 double sigma_k,
2666 double sigma_e,
2667 double c_1,
2668 double c_2,
2669 double c_mu,
2670 double c_e,
2671 double nu,
2672 double *velocity,
2673 double *gradu,
2674 double *gradv,
2675 double *k,
2676 double *epsilon,
2677 double *m_k,
2678 double *dm_k,
2679 double *m_e,
2680 double *dm_e,
2681 double *phi_k,
2682 double *dphi_k,
2683 double *phi_e,
2684 double *dphi_e,
2685 double *f_k,
2686 double *df_k,
2687 double *f_e,
2688 double *df_e,
2689 double *a_k,
2690 double *da_k_dk,
2691 double *da_k_de,
2692 double *a_e,
2693 double *da_e_dk,
2694 double *da_e_de,
2695 double *r_k,
2696 double *dr_k_dk,
2697 double *dr_k_de,
2698 double *r_e,
2699 double *dr_e_dk,
2700 double *dr_e_de);
2701
2702extern void kEpsilon_2D_Evaluate_sd(int nPoints,
2703 int nSpace,
2704 double sigma_k,
2705 double sigma_e,
2706 double c_1,
2707 double c_2,
2708 double c_mu,
2709 double c_e,
2710 double nu,
2711 double *velocity,
2712 double *gradu,
2713 double *gradv,
2714 double *k,
2715 double *epsilon,
2716 double *m_k,
2717 double *dm_k,
2718 double *m_e,
2719 double *dm_e,
2720 double *phi_k,
2721 double *dphi_k,
2722 double *phi_e,
2723 double *dphi_e,
2724 double *f_k,
2725 double *df_k,
2726 double *f_e,
2727 double *df_e,
2728 double *a_k,
2729 double *da_k_dk,
2730 double *da_k_de,
2731 double *a_e,
2732 double *da_e_dk,
2733 double *da_e_de,
2734 double *r_k,
2735 double *dr_k_dk,
2736 double *dr_k_de,
2737 double *r_e,
2738 double *dr_e_dk,
2739 double *dr_e_de);
2740extern void kEpsilon_3D_Evaluate_sd(int nPoints,
2741 int nSpace,
2742 double sigma_k,
2743 double sigma_e,
2744 double c_1,
2745 double c_2,
2746 double c_mu,
2747 double c_e,
2748 double nu,
2749 double *velocity,
2750 double *gradu,
2751 double *gradv,
2752 double *gradw,
2753 double *k,
2754 double *epsilon,
2755 double *m_k,
2756 double *dm_k,
2757 double *m_e,
2758 double *dm_e,
2759 double *phi_k,
2760 double *dphi_k,
2761 double *phi_e,
2762 double *dphi_e,
2763 double *f_k,
2764 double *df_k,
2765 double *f_e,
2766 double *df_e,
2767 double *a_k,
2768 double *da_k_dk,
2769 double *da_k_de,
2770 double *a_e,
2771 double *da_e_dk,
2772 double *da_e_de,
2773 double *r_k,
2774 double *dr_k_dk,
2775 double *dr_k_de,
2776 double *r_e,
2777 double *dr_e_dk,
2778 double *dr_e_de);
2779extern void kEpsilon_3D_Evaluate(int nPoints,
2780 int nSpace,
2781 double sigma_k,
2782 double sigma_e,
2783 double c_1,
2784 double c_2,
2785 double c_mu,
2786 double c_e,
2787 double nu,
2788 double *velocity,
2789 double *gradu,
2790 double *gradv,
2791 double *gradw,
2792 double *k,
2793 double *epsilon,
2794 double *m_k,
2795 double *dm_k,
2796 double *m_e,
2797 double *dm_e,
2798 double *phi_k,
2799 double *dphi_k,
2800 double *phi_e,
2801 double *dphi_e,
2802 double *f_k,
2803 double *df_k,
2804 double *f_e,
2805 double *df_e,
2806 double *a_k,
2807 double *da_k_dk,
2808 double *da_k_de,
2809 double *a_e,
2810 double *da_e_dk,
2811 double *da_e_de,
2812 double *r_k,
2813 double *dr_k_dk,
2814 double *dr_k_de,
2815 double *r_e,
2816 double *dr_e_dk,
2817 double *dr_e_de);
2818extern void kEpsilon_k_2D_Evaluate_sd(int nPoints,
2819 int nSpace,
2820 double sigma_k,
2821 double c_mu,
2822 double nu,
2823 double *velocity,
2824 double *gradu,
2825 double *gradv,
2826 double *k,
2827 double *epsilon,
2828 double *m_k,
2829 double *dm_k,
2830 double *phi_k,
2831 double *dphi_k,
2832 double *f_k,
2833 double *df_k,
2834 double *a_k,
2835 double *da_k_dk,
2836 double *r_k,
2837 double *dr_k_dk);
2838
2839extern void kEpsilon_epsilon_2D_Evaluate_sd(int nPoints,
2840 int nSpace,
2841 double sigma_e,
2842 double c_1,
2843 double c_2,
2844 double c_mu,
2845 double c_e,
2846 double nu,
2847 double *velocity,
2848 double *gradu,
2849 double *gradv,
2850 double *k,
2851 double *epsilon,
2852 double *m_e,
2853 double *dm_e,
2854 double *phi_e,
2855 double *dphi_e,
2856 double *f_e,
2857 double *df_e,
2858 double *a_e,
2859 double *da_e_de,
2860 double *r_e,
2861 double *dr_e_de);
2862
2863extern void kEpsilon_k_3D_Evaluate_sd(int nPoints,
2864 int nSpace,
2865 double sigma_k,
2866 double c_mu,
2867 double nu,
2868 double *velocity,
2869 double *gradu,
2870 double *gradv,
2871 double *gradw,
2872 double *k,
2873 double *epsilon,
2874 double *m_k,
2875 double *dm_k,
2876 double *phi_k,
2877 double *dphi_k,
2878 double *f_k,
2879 double *df_k,
2880 double *a_k,
2881 double *da_k_dk,
2882 double *r_k,
2883 double *dr_k_dk);
2884
2885extern void kEpsilon_epsilon_3D_Evaluate_sd(int nPoints,
2886 int nSpace,
2887 double sigma_e,
2888 double c_1,
2889 double c_2,
2890 double c_mu,
2891 double c_e,
2892 double nu,
2893 double *velocity,
2894 double *gradu,
2895 double *gradv,
2896 double *gradw,
2897 double *k,
2898 double *epsilon,
2899 double *m_e,
2900 double *dm_e,
2901 double *phi_e,
2902 double *dphi_e,
2903 double *f_e,
2904 double *df_e,
2905 double *a_e,
2906 double *da_e_de,
2907 double *r_e,
2908 double *dr_e_de);
2909extern void ReynoldsAveragedNavierStokes_kEpsilon_2D_Update(const int nPoints,
2910 const double nu,
2911 const double c_mu,
2912 const double* k,
2913 const double* grad_k,
2914 const double* epsilon,
2915 double *mom_u_diff_ten,
2916 double *mom_v_diff_ten,
2917 double *mom_uv_diff_ten,
2918 double *mom_vu_diff_ten,
2919 double *mom_u_source,
2920 double *mom_v_source);
2921extern void ReynoldsAveragedNavierStokes_kEpsilon_2D_Update_sd(const int nPoints,
2922 const double rho,
2923 const double nu,
2924 const double c_mu,
2925 const double* k,
2926 const double* grad_k,
2927 const double* epsilon,
2928 double *mom_u_diff_ten,
2929 double *mom_v_diff_ten,
2930 double *mom_uv_diff_ten,
2931 double *mom_vu_diff_ten,
2932 double *mom_u_source,
2933 double *mom_v_source);
2934
2935extern void ReynoldsAveragedNavierStokes_kEpsilon_3D_Update(const int nPoints,
2936 const double nu,
2937 const double c_mu,
2938 const double* k,
2939 const double* grad_k,
2940 const double* epsilon,
2941 double *mom_u_diff_ten,
2942 double *mom_v_diff_ten,
2943 double *mom_w_diff_ten,
2944 double *mom_uv_diff_ten,
2945 double *mom_uw_diff_ten,
2946 double *mom_vu_diff_ten,
2947 double *mom_vw_diff_ten,
2948 double *mom_wu_diff_ten,
2949 double *mom_wv_diff_ten,
2950 double *mom_u_source,
2951 double *mom_v_source,
2952 double *mom_w_source);
2953
2954extern void ReynoldsAveragedNavierStokes_kEpsilon_3D_Update_sd(const int nPoints,
2955 const double nu,
2956 const double c_mu,
2957 const double* k,
2958 const double* grad_k,
2959 const double* epsilon,
2960 double *mom_u_diff_ten,
2961 double *mom_v_diff_ten,
2962 double *mom_w_diff_ten,
2963 double *mom_uv_diff_ten,
2964 double *mom_uw_diff_ten,
2965 double *mom_vu_diff_ten,
2966 double *mom_vw_diff_ten,
2967 double *mom_wu_diff_ten,
2968 double *mom_wv_diff_ten,
2969 double *mom_u_source,
2970 double *mom_v_source,
2971 double *mom_w_source);
2972
2973extern void scriptedSphereMotionSignedDistance(const int nPoints,
2974 const double t,
2975 const int nSpace,
2976 const int nSpheres,
2977 const double * radii,
2978 const double * centers,
2979 const double * x,
2980 double * phi,
2981 double * n);
2982
2983extern void shallowWater_1D_Evaluate(const int nPoints,
2984 const double h_eps,
2985 const double g,
2986 const double bedFrictionCoefficient,
2987 const double bedFrictionPower,
2988 const double eddyViscosity,
2989 const double* x, /*bathymetry elevation in index 1 (y)*/
2990 const double* db_dx,/*db/dx (bed slope) */
2991 const double* h,
2992 const double* hu,
2993 double *H,
2994 double *mass_acc,
2995 double *dmass_acc_dh,
2996 double *mom_acc,
2997 double *dmom_acc_dhu,
2998 double *mass_adv,
2999 double *dmass_adv_dhu,
3000 double *mom_adv,
3001 double *dmom_adv_dh,
3002 double *dmom_adv_dhu,
3003 double *mom_source,
3004 double *dmom_source_dh,
3005 double *dmom_source_dhu,
3006 double *mom_diff);
3007
3008extern void VolumeAveragedTwophaseNavierStokes_ST_LS_SO_2D_Evaluate_sd(const int nPoints,
3009 const int killNonlinearDrag,
3010 const double eps_rho,
3011 const double eps_mu,
3012 const double sigma,
3013 const double rho_0,
3014 const double nu_0,
3015 const double rho_1,
3016 const double nu_1,
3017 const double* meanGrainSize,
3018 const double* g,
3019 const double* phi,
3020 const double* n,
3021 const double* kappa,
3022 const double *p,
3023 const double *grad_p,
3024 const double *u,
3025 const double *v,
3026 const double *porosity,
3027 double *mom_u_acc,
3028 double *dmom_u_acc_u,
3029 double *mom_v_acc,
3030 double *dmom_v_acc_v,
3031 double *mass_adv,
3032 double *dmass_adv_u,
3033 double *dmass_adv_v,
3034 double *mom_u_adv,
3035 double *dmom_u_adv_u,
3036 double *dmom_u_adv_v,
3037 double *mom_v_adv,
3038 double *dmom_v_adv_u,
3039 double *dmom_v_adv_v,
3040 double *mom_u_diff_ten,
3041 double *mom_v_diff_ten,
3042 double *mom_uv_diff_ten,
3043 double *mom_vu_diff_ten,
3044 double *mom_u_source,
3045 double *mom_v_source,
3046 double *dmom_u_source_u,
3047 double *dmom_u_source_v,
3048 double *dmom_v_source_u,
3049 double *dmom_v_source_v,
3050 double *mom_u_ham,
3051 double *dmom_u_ham_grad_p,
3052 double *mom_v_ham,
3053 double *dmom_v_ham_grad_p);
3054extern void shallowWater_2D_Evaluate(const int nPoints,
3055 const double h_eps,
3056 const double g,
3057 const double bedFrictionCoefficient,
3058 const double bedFrictionPower,
3059 const double eddyViscosity,
3060 const double* x, /*bathymetry elevation in index 2 (z)*/
3061 const double* grad_b,/*\nabla b (bed slope) */
3062 const double* h,
3063 const double* hu,
3064 const double* hv,
3065 double *H,
3066 double *mass_acc,
3067 double *dmass_acc_dh,
3068 double *mom_u_acc,
3069 double *dmom_u_acc_dhu,
3070 double *mom_v_acc,
3071 double *dmom_v_acc_dhv,
3072 double *mass_adv,
3073 double *dmass_adv_dhu,
3074 double *dmass_adv_dhv,
3075 double *mom_u_adv,
3076 double *dmom_u_adv_dh,
3077 double *dmom_u_adv_dhu,
3078 double *dmom_u_adv_dhv,
3079 double *mom_v_adv,
3080 double *dmom_v_adv_dh,
3081 double *dmom_v_adv_dhu,
3082 double *dmom_v_adv_dhv,
3083 double *mom_u_diff,
3084 double *mom_v_diff,
3085 double *mom_u_source,
3086 double *dmom_u_source_dh,
3087 double *dmom_u_source_dhu,
3088 double *dmom_u_source_dhv,
3089 double *mom_v_source,
3090 double *dmom_v_source_dh,
3091 double *dmom_v_source_dhu,
3092 double *dmom_v_source_dhv);
3093extern void conservativeHeadRichardsMualemVanGenuchtenHetEvaluateV2withUpwind(const int upwindFlag,
3094 const int computeAverages,
3095 const int nSimplex,
3096 const int nPointsPerSimplex,
3097 const int nSpace,
3098 const int nQuadraturePoints_elementBoundary,
3099 const int* elementBoundaryElementsArray,
3100 const int* quadraturePointToElementBoundary,
3101 const int* materialTypes,
3102 const double rho,
3103 const double beta,
3104 const double* gravity,
3105 const double* alpha,
3106 const double* n_vg,
3107 const double* thetaR,
3108 const double* thetaSR,
3109 const double* KWs,
3110 const double *u,
3111 const double *gradu,
3112 const double *n_global,
3113 const double *dV,
3114 double *mass,
3115 double *dmass,
3116 double *f_avg,
3117 double *df_avg,
3118 double *a_avg,
3119 double *da_avg,
3120 double *f,
3121 double *df,
3122 double *a,
3123 double *da);
3124extern void conservativeHeadRichardsMualemVanGenuchtenHetEvaluateV2withUpwindAndHarm(const int upwindFlag,
3125 const int computeAverages,
3126 const int nSimplex,
3127 const int nPointsPerSimplex,
3128 const int nSpace,
3129 const int nQuadraturePoints_elementBoundary,
3130 const int* elementBoundaryElementsArray,
3131 const int* quadraturePointToElementBoundary,
3132 const int* materialTypes,
3133 const double rho,
3134 const double beta,
3135 const double* gravity,
3136 const double* alpha,
3137 const double* n_vg,
3138 const double* thetaR,
3139 const double* thetaSR,
3140 const double* KWs,
3141 const double *u,
3142 const double *gradu,
3143 const double *n_global,
3144 const double *dV,
3145 double *mass,
3146 double *dmass,
3147 double *f_avg,
3148 double *df_avg,
3149 double *a_avg,
3150 double *da_avg,
3151 double *f,
3152 double *df,
3153 double *a,
3154 double *da);
3155extern void conservativeHeadRichardsMualemVanGenuchtenHetEvaluateV2withUpwindAndHarm_sd(const int upwindFlag,
3156 const int computeAverages,
3157 const int nSimplex,
3158 const int nPointsPerSimplex,
3159 const int nSpace,
3160 const int nQuadraturePoints_elementBoundary,
3161 const int *rowptr,
3162 const int *colind,
3163 const int* elementBoundaryElementsArray,
3164 const int* quadraturePointToElementBoundary,
3165 const int* materialTypes,
3166 const double rho,
3167 const double beta,
3168 const double* gravity,
3169 const double* alpha,
3170 const double* n_vg,
3171 const double* thetaR,
3172 const double* thetaSR,
3173 const double* KWs,
3174 const double *u,
3175 const double *gradu,
3176 const double *n_global,
3177 const double *dV,
3178 double *mass,
3179 double *dmass,
3180 double *f_avg,
3181 double *df_avg,
3182 double *a_avg,
3183 double *da_avg,
3184 double *f,
3185 double *df,
3186 double *a,
3187 double *da);
3188extern void applyContactLineSlip(int nExteriorElementBoundaries_global,
3189 int nQuadraturePoints_elementBoundary,
3190 double eps,
3191 int* isDOFBoundary,
3192 double* phi,
3193 double* advectiveFlux,
3194 double* diffusiveFlux);
3195extern void applyContactLineSlipJacobian(int nExteriorElementBoundaries_global,
3196 int nQuadraturePoints_elementBoundary,
3197 int nDOF_trial_element,
3198 double eps,
3199 int* isDOFBoundary,
3200 double* phi,
3201 double* fluxJacobian);
3202extern void calculateEddyViscosity_Smagorinsky_2D(const int nElements_global,
3203 const int nQuadraturePoints_element,
3204 const double smagorinskyConstant,
3205 const double * h_e,
3206 const double * grad_u,
3207 const double * grad_v,
3208 double * nu_t);
3209extern void calculateEddyViscosity_Smagorinsky_3D(const int nElements_global,
3210 const int nQuadraturePoints_element,
3211 const double smagorinskyConstant,
3212 const double * h_e,
3213 const double * grad_u,
3214 const double * grad_v,
3215 const double * grad_w,
3216 double * nu_t);
3217extern void calculateEddyViscosity_Smagorinsky2P_2D(const int nElements_global,
3218 const int nQuadraturePoints_element,
3219 const double smagorinskyConstant_0,
3220 const double smagorinskyConstant_1,
3221 const double eps,
3222 const double * phi_ls,
3223 const double * h_e,
3224 const double * grad_u,
3225 const double * grad_v,
3226 double * nu_t);
3227extern void calculateEddyViscosity_Smagorinsky2P_3D(const int nElements_global,
3228 const int nQuadraturePoints_element,
3229 const double smagorinskyConstant_0,
3230 const double smagorinskyConstant_1,
3231 const double eps,
3232 const double * phi_ls,
3233 const double * h_e,
3234 const double * grad_u,
3235 const double * grad_v,
3236 const double * grad_w,
3237 double * nu_t);
3238extern void eddyViscosity_2D_Update(const int nPoints,
3239 const double* nu_t,
3240 double *mom_u_diff_ten,
3241 double *mom_v_diff_ten,
3242 double *mom_uv_diff_ten,
3243 double *mom_vu_diff_ten);
3244extern void eddyViscosity_2D_Update_sd(const int nPoints,
3245 const double* nu_t,
3246 double *mom_u_diff_ten,
3247 double *mom_v_diff_ten,
3248 double *mom_uv_diff_ten,
3249 double *mom_vu_diff_ten);
3250extern void eddyViscosity_3D_Update(const int nPoints,
3251 const double* nu_t,
3252 double *mom_u_diff_ten,
3253 double *mom_v_diff_ten,
3254 double *mom_w_diff_ten,
3255 double *mom_uv_diff_ten,
3256 double *mom_uw_diff_ten,
3257 double *mom_vu_diff_ten,
3258 double *mom_vw_diff_ten,
3259 double *mom_wu_diff_ten,
3260 double *mom_wv_diff_ten);
3261extern void eddyViscosity_3D_Update_sd(const int nPoints,
3262 const double* nu_t,
3263 double *mom_u_diff_ten,
3264 double *mom_v_diff_ten,
3265 double *mom_w_diff_ten,
3266 double *mom_uv_diff_ten,
3267 double *mom_uw_diff_ten,
3268 double *mom_vu_diff_ten,
3269 double *mom_vw_diff_ten,
3270 double *mom_wu_diff_ten,
3271 double *mom_wv_diff_ten);
3272extern void calculateWaveFunction3d_ref(//mesh rep
3273 int nElements_global,
3274 int nDOF_element_mesh,
3275 int nQuadraturePoints_element,
3276 const double* mesh_trial_ref,
3277 const double* mesh_dof,
3278 const int* mesh_l2g,
3279 const double* elementDiametersArray,
3280 //
3281 const double* omega_s_x, //source region (rectangular)
3282 const double* omega_s_y,
3283 const double* omega_s_z,
3284 double t,
3285 int waveFlag, //1 secondOrderStokes
3286 //2 solitaryWave
3287 //0 monochromaticWave
3288 double epsFact,
3289 double waveHeight,
3290 double waveCelerity,
3291 double waveFrequency,
3292 double waveNumber,
3293 double waterDepth,
3294 double* source);
3295extern void Mass_2D_Evaluate(const int nPoints,
3296 double rho,
3297 double *p,
3298 double *u,
3299 double *v,
3300 double *mom_p_acc,
3301 double *mom_u_acc,
3302 double *mom_v_acc,
3303 double *dmom_p_acc_p,
3304 double *dmom_u_acc_u,
3305 double *dmom_v_acc_v);
3306extern void Mass_3D_Evaluate(const int nPoints,
3307 double rho,
3308 double *p,
3309 double *u,
3310 double *v,
3311 double *w,
3312 double *mom_p_acc,
3313 double *mom_u_acc,
3314 double *mom_v_acc,
3315 double *mom_w_acc,
3316 double *dmom_p_acc_p,
3317 double *dmom_u_acc_u,
3318 double *dmom_v_acc_v,
3319 double *dmom_w_acc_w);
3320extern void TwoPhaseAdvection_2D_Evaluate(const int nPoints,
3321 const double eps,
3322 const double rho_0,
3323 const double nu_0,
3324 const double rho_1,
3325 const double nu_1,
3326 const double *phi,
3327 const double *p,
3328 const double *u,
3329 const double *v,
3330 double *mass_adv,
3331 double *dmass_adv_p,
3332 double *dmass_adv_u,
3333 double *dmass_adv_v,
3334 double *mom_u_adv,
3335 double *dmom_u_adv_u,
3336 double *dmom_u_adv_v,
3337 double *mom_v_adv,
3338 double *dmom_v_adv_u,
3339 double *dmom_v_adv_v);
3340extern void Laplace_Evaluate2D(const int nPoints,
3341 double *mom_p_diff_ten,
3342 double *mom_u_diff_ten,
3343 double *mom_v_diff_ten);
3344extern void Laplace_Evaluate3D(const int nPoints,
3345 double *mom_p_diff_ten,
3346 double *mom_u_diff_ten,
3347 double *mom_v_diff_ten,
3348 double *mom_w_diff_ten);
3349extern void conservativeHeadRichardsJLeverett(const int nSimplex,
3350 const int nPointsPerSimplex,
3351 const int nSpace,
3352 const int* materialTypes,
3353 const double rho,
3354 const double beta,
3355 const double* gravity,
3356 const double* phi,
3357 const double* psiD,
3358 const double* ns,
3359 const double* nk,
3360 const double* S_wirr,
3361 const double* S_nwr,
3362 const double* kr0,
3363 double *u,
3364 double *mass,
3365 double *dmass,
3366 double *f,
3367 double *df,
3368 double *a,
3369 double *da);
3370extern void conservativeHeadRichardsJLeverettAni(const int nSimplex,
3371 const int nPointsPerSimplex,
3372 const int nSpace,
3373 const int* materialTypes,
3374 const double rho,
3375 const double beta,
3376 const double* gravity,
3377 const double* phi,
3378 const double* psiD,
3379 const double* ns,
3380 const double* nk,
3381 const double* S_wirr,
3382 const double* S_nwr,
3383 const double* kr0x,
3384 const double* kr0y,
3385 const double* kr0z,
3386 double *u,
3387 double *mass,
3388 double *dmass,
3389 double *f,
3390 double *df,
3391 double *a,
3392 double *da);
3393extern void diffusiveWave1DEvaluate(const int nPoints,
3394 const double alpha,
3395 const double gamma,
3396 const double epsilon,
3397 const double* x,
3398 const double* u,
3399 const double* grad_u,
3400 double* m,
3401 double* dm,
3402 double* a,
3403 double* da);
3404extern void diffusiveWave2DEvaluate(const int nd,
3405 const int nPoints,
3406 const double alpha,
3407 const double gamma,
3408 const double epsilon,
3409 const double* x,
3410 const double* u,
3411 const double* grad_u,
3412 double* m,
3413 double* dm,
3414 double* a,
3415 double* da);
3417#endif
nu_0
double nu_0
Definition ErrorResidualMethod.cpp:22
rho_1
double rho_1
Definition ErrorResidualMethod.cpp:22
nu_1
double nu_1
Definition ErrorResidualMethod.cpp:22
rho_0
double rho_0
Definition ErrorResidualMethod.cpp:22
n
Int n
Definition Headers.h:28
r
Double r
Definition Headers.h:83
H
Double H
Definition Headers.h:65
B
Double * B
Definition Headers.h:41
f
Double f
Definition Headers.h:64
u
Double u
Definition Headers.h:89
ns
Int ns
Definition Headers.h:30
v
Double v
Definition Headers.h:95
R
Double R
Definition Headers.h:82
phi
Double phi
Definition Headers.h:76
df
double df(double C, double b, double a, int q, int r)
Definition analyticalSolutions.c:2209
applyContactLineSlip
void applyContactLineSlip(int nExteriorElementBoundaries_global, int nQuadraturePoints_elementBoundary, double eps, int *isDOFBoundary, double *phi, double *advectiveFlux, double *diffusiveFlux)
Definition transportCoefficients.c:11924
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
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
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
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
cLevelSetCoefficientsEvaluate
void cLevelSetCoefficientsEvaluate(int nPoints, int nSpace, double *v, double *u, double *m, double *dm, double *f, double *df)
Definition transportCoefficients.c:1222
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
kEpsilon_epsilon_3D_Evaluate_sd
void kEpsilon_epsilon_3D_Evaluate_sd(int nPoints, int nSpace, double sigma_e, double c_1, double c_2, double c_mu, double c_e, double nu, double *velocity, double *gradu, double *gradv, double *gradw, double *k, double *epsilon, double *m_e, double *dm_e, double *phi_e, double *dphi_e, double *f_e, double *df_e, double *a_e, double *da_e_de, double *r_e, double *dr_e_de)
Definition transportCoefficients.c:10686
eddyViscosity_3D_Update_sd
void eddyViscosity_3D_Update_sd(const int nPoints, const double *nu_t, 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)
Definition transportCoefficients.c:12407
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
levelSetConservationCoefficientsEvaluate_sd
void levelSetConservationCoefficientsEvaluate_sd(int nPoints, double epsHeaviside, double epsDirac, double *u_ls, double *H_vof, double *u, double *r, double *dr)
Definition transportCoefficients.c:8438
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
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
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
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
Mass_3D_Evaluate
void Mass_3D_Evaluate(const int nPoints, double rho, double *p, double *u, double *v, double *w, double *mom_p_acc, double *mom_u_acc, double *mom_v_acc, double *mom_w_acc, double *dmom_p_acc_p, double *dmom_u_acc_u, double *dmom_v_acc_v, double *dmom_w_acc_w)
Definition transportCoefficients.c:12651
calculateWaveFunction3d_ref
void calculateWaveFunction3d_ref(int nElements_global, int nDOF_element_mesh, int nQuadraturePoints_element, const double *mesh_trial_ref, const double *mesh_dof, const int *mesh_l2g, const double *elementDiametersArray, const double *omega_s_x, const double *omega_s_y, const double *omega_s_z, double t, int waveFlag, double epsFact, double waveHeight, double waveCelerity, double waveFrequency, double waveNumber, double waterDepth, double *source)
Definition transportCoefficients.c:12452
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
VolumeAveragedTwophaseNavierStokes_ST_LS_SO_3D_Evaluate_sd
void VolumeAveragedTwophaseNavierStokes_ST_LS_SO_3D_Evaluate_sd(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 *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:9660
TwophaseNavierStokes_ST_LS_SO_2D_Evaluate_sd
void TwophaseNavierStokes_ST_LS_SO_2D_Evaluate_sd(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:2677
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
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
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
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
diffusiveWave1DEvaluate
void diffusiveWave1DEvaluate(const int nPoints, const double alpha, const double gamma, const double epsilon, const double *x, const double *u, const double *grad_u, double *m, double *dm, double *a, double *da)
Definition transportCoefficients.c:11978
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
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
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
shallowWater_2D_Evaluate
void shallowWater_2D_Evaluate(const int nPoints, const double h_eps, const double g, const double bedFrictionCoefficient, const double bedFrictionPower, const double eddyViscosity, const double *x, const double *grad_b, const double *h, const double *hu, const double *hv, double *H, double *mass_acc, double *dmass_acc_dh, double *mom_u_acc, double *dmom_u_acc_dhu, double *mom_v_acc, double *dmom_v_acc_dhv, double *mass_adv, double *dmass_adv_dhu, double *dmass_adv_dhv, double *mom_u_adv, double *dmom_u_adv_dh, double *dmom_u_adv_dhu, double *dmom_u_adv_dhv, double *mom_v_adv, double *dmom_v_adv_dh, double *dmom_v_adv_dhu, double *dmom_v_adv_dhv, double *mom_u_diff, double *mom_v_diff, double *mom_u_source, double *dmom_u_source_dh, double *dmom_u_source_dhu, double *dmom_u_source_dhv, double *mom_v_source, double *dmom_v_source_dh, double *dmom_v_source_dhu, double *dmom_v_source_dhv)
Definition transportCoefficients.c:11180
calculateEddyViscosity_Smagorinsky_3D
void calculateEddyViscosity_Smagorinsky_3D(const int nElements_global, const int nQuadraturePoints_element, const double smagorinskyConstant, const double *h_e, const double *grad_u, const double *grad_v, const double *grad_w, double *nu_t)
Definition transportCoefficients.c:12084
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)
conservativeHeadRichardsJLeverett
void conservativeHeadRichardsJLeverett(const int nSimplex, const int nPointsPerSimplex, const int nSpace, const int *materialTypes, const double rho, const double beta, const double *gravity, const double *phi, const double *psiD, const double *ns, const double *nk, const double *S_wirr, const double *S_nwr, const double *kr0, double *u, double *mass, double *dmass, double *f, double *df, double *a, double *da)
Definition transportCoefficients.c:6143
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
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
groundwaterBryantDawsonIonExEvaluateFC
void groundwaterBryantDawsonIonExEvaluateFC(const int nPoints, const int nSpace, const double omega, const double d_m, const double d_h, const double alpha_L, const double alpha_T, const double K_m, const double K_h, const double K_w, const double Z_tot, const double *v, const double *c_m, const double *c_h, double *m_m, double *dm_m_m, double *dm_m_h, double *m_h, double *dm_h_m, double *dm_h_h, double *f_m, double *df_m, double *f_h, double *df_h, double *a_m, double *a_h, double *phi_h, double *dphi_h, double *r_m, double *dr_m_dm, double *dr_m_dh, double *r_h, double *dr_h_dm, double *dr_h_dh)
Definition transportCoefficients.c:303
seepageBrezis
void seepageBrezis(const int nSimplex, const int nPointsPerSimplex, const int nSpace, const int *materialTypes, const double epsFact, const double rho, const double beta, const double *elementDiameter, 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:6083
eddyViscosity_2D_Update
void eddyViscosity_2D_Update(const int nPoints, const double *nu_t, double *mom_u_diff_ten, double *mom_v_diff_ten, double *mom_uv_diff_ten, double *mom_vu_diff_ten)
Definition transportCoefficients.c:12307
variablySaturatedGroundwaterEnergyTransportCoefficientsEvaluate_hetMat
void variablySaturatedGroundwaterEnergyTransportCoefficientsEvaluate_hetMat(const int nSimplex, const int nPointsPerSimplex, const int nSpace, const double rho_w, const double rho_n, const double specificHeat_w, const double specificHeat_n, const int *materialTypes, const double *theta, const double *thetaS_types, const double *alpha_L_types, const double *alpha_T_types, const double *rho_s_types, const double *specificHeat_s_types, const double *lambda_sat_types, const double *lambda_dry_types, const double *lambda_aniso_types, const double *v, const double *u, double *m, double *dm, double *f, double *df, double *a)
Definition transportCoefficients.c:523
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)
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
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)
setSimpleWeakDirichletConditionsForLevelSet
void setSimpleWeakDirichletConditionsForLevelSet(int nElements_global, int nDOF_trial_element, double epsilon_freeze_factor, const double *elementDiameter, const int *u_l2g, const double *u_dof, int *freeze_nodes_tmp, int *weakDirichletConditionFlags)
Definition transportCoefficients.c:1672
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
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
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
smoothedHeaviside
double smoothedHeaviside(double eps, double phi)
Definition numericalFlux.c:14167
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
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
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
diffusiveWave2DEvaluate
void diffusiveWave2DEvaluate(const int nd, const int nPoints, const double alpha, const double gamma, const double epsilon, const double *x, const double *u, const double *grad_u, double *m, double *dm, double *a, double *da)
Definition transportCoefficients.c:12007
eddyViscosity_2D_Update_sd
void eddyViscosity_2D_Update_sd(const int nPoints, const double *nu_t, double *mom_u_diff_ten, double *mom_v_diff_ten, double *mom_uv_diff_ten, double *mom_vu_diff_ten)
Definition transportCoefficients.c:12340
kEpsilon_3D_Evaluate_sd
void kEpsilon_3D_Evaluate_sd(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 *gradw, 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:10205
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
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
groundwaterBiodegradation01EvaluateFC
void groundwaterBiodegradation01EvaluateFC(const int nPoints, const int nSpace, const double omega, const double d_c, const double d_e, const double alpha_L, const double alpha_T, const double Kox_max, const double Kox_C, const double Kox_E, const double Kox_X, const double Yield, const double k_d, const double *v, const double *c_c, const double *c_e, const double *c_x, double *m_c, double *dm_c, double *m_e, double *dm_e, double *m_x, double *dm_x, double *f_c, double *df_c, double *f_e, double *df_e, double *a_c, double *a_e, double *r_c, double *dr_c_dc, double *dr_c_de, double *dr_c_dx, double *r_e, double *dr_e_dc, double *dr_e_de, double *dr_e_dx, double *r_x, double *dr_x_dc, double *dr_x_de, double *dr_x_dx)
Definition transportCoefficients.c:187
conservativeHeadRichardsMualemVanGenuchten_sd_het_linearized_at_saturation
void conservativeHeadRichardsMualemVanGenuchten_sd_het_linearized_at_saturation(const int nSimplex, const int nPointsPerSimplex, const int nSpace, double linear_break, const int *rowptr, const int *colind, 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, double *vol_frac)
Definition transportCoefficients.c:5847
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
Laplace_Evaluate3D
void Laplace_Evaluate3D(const int nPoints, double *mom_p_diff_ten, double *mom_u_diff_ten, double *mom_v_diff_ten, double *mom_w_diff_ten)
Definition transportCoefficients.c:1822
ThreephaseNavierStokes_ST_LS_SO_2D_Evaluate
void ThreephaseNavierStokes_ST_LS_SO_2D_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, 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 *dmom_u_source_u, double *dmom_u_source_v, double *mom_v_source, 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:3050
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
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
ReynoldsAveragedNavierStokes_kEpsilon_2D_Update
void ReynoldsAveragedNavierStokes_kEpsilon_2D_Update(const int nPoints, const double nu, const double c_mu, const double *k, const double *grad_k, const double *epsilon, 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)
Definition transportCoefficients.c:10808
l2project2Tensor
void l2project2Tensor(const int nSimplices, const int nPointsPerSimplex, const int nSpace, double *dV, double *r)
Definition transportCoefficients.c:5705
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
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
setWeakDirichletConditionsForLevelSet
void setWeakDirichletConditionsForLevelSet(int nElements_global, int nDOF_trial_element, double epsilon_freeze_factor, const double *elementDiameter, const int *u_l2g, const double *u_dof, int *freeze_nodes_tmp, int *weakDirichletConditionFlags)
Definition transportCoefficients.c:1620
kEpsilon_epsilon_2D_Evaluate_sd
void kEpsilon_epsilon_2D_Evaluate_sd(int nPoints, int nSpace, 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_e, double *dm_e, double *phi_e, double *dphi_e, double *f_e, double *df_e, double *a_e, double *da_e_de, double *r_e, double *dr_e_de)
Definition transportCoefficients.c:10526
StokesP_3D_Evaluate
void StokesP_3D_Evaluate(const int nPoints, const double rho, const double nu, const double *g, const double *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_p, double *mom_v_adv, double *dmom_v_adv_p, double *mom_w_adv, double *dmom_w_adv_p, 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)
Definition transportCoefficients.c:2295
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)
Laplace_Evaluate2D
void Laplace_Evaluate2D(const int nPoints, double *mom_p_diff_ten, double *mom_u_diff_ten, double *mom_v_diff_ten)
Definition transportCoefficients.c:1803
conservativeHeadRichardsMualemVanGenuchten_sd_het
void conservativeHeadRichardsMualemVanGenuchten_sd_het(const int nSimplex, const int nPointsPerSimplex, const int nSpace, double pc_eps, const int *rowptr, const int *colind, 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, double *vol_frac)
Definition transportCoefficients.c:5740
conservativeHeadRichardsMualemVanGenuchtenHetEvaluateV2withUpwindAndHarm_sd
void conservativeHeadRichardsMualemVanGenuchtenHetEvaluateV2withUpwindAndHarm_sd(const int upwindFlag, const int computeAverages, const int nSimplex, const int nPointsPerSimplex, const int nSpace, const int nQuadraturePoints_elementBoundary, const int *rowptr, const int *colind, const int *elementBoundaryElementsArray, const int *quadraturePointToElementBoundary, const int *materialTypes, const double rho, const double beta, const double *gravity, const double *alpha, const double *n_vg, const double *thetaR, const double *thetaSR, const double *KWs, const double *u, const double *gradu, const double *n_global, const double *dV, double *mass, double *dmass, double *f_avg, double *df_avg, double *a_avg, double *da_avg, double *f, double *df, double *a, double *da)
Definition transportCoefficients.c:11710
calculateEddyViscosity_Smagorinsky_2D
void calculateEddyViscosity_Smagorinsky_2D(const int nElements_global, const int nQuadraturePoints_element, const double smagorinskyConstant, const double *h_e, const double *grad_u, const double *grad_v, double *nu_t)
Definition transportCoefficients.c:12036
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
twophaseLevelSetCoefficientsUpdateVelocity
void twophaseLevelSetCoefficientsUpdateVelocity(int nPoints, int nSpace, double v_scale, double *vIn, double *vOut)
Definition transportCoefficients.c:1135
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
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
kEpsilon_k_2D_Evaluate_sd
void kEpsilon_k_2D_Evaluate_sd(int nPoints, int nSpace, double sigma_k, double c_mu, double nu, double *velocity, double *gradu, double *gradv, double *k, double *epsilon, double *m_k, double *dm_k, double *phi_k, double *dphi_k, double *f_k, double *df_k, double *a_k, double *da_k_dk, double *r_k, double *dr_k_dk)
Definition transportCoefficients.c:10456
scriptedSphereMotionSignedDistance
void scriptedSphereMotionSignedDistance(const int nPoints, const double t, const int nSpace, const int nSpheres, const double *radii, const double *centers, const double *x, double *phi, double *n)
Definition transportCoefficients.c:11059
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
twophaseSignedDistanceCoefficientsUpdateSignFunction
void twophaseSignedDistanceCoefficientsUpdateSignFunction(int nPoints, double eps, double *u_levelSet, double *S)
Definition transportCoefficients.c:1298
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
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)
conservativeHeadRichardsMualemVanGenuchtenHetEvaluateV2withUpwind
void conservativeHeadRichardsMualemVanGenuchtenHetEvaluateV2withUpwind(const int upwindFlag, const int computeAverages, const int nSimplex, const int nPointsPerSimplex, const int nSpace, const int nQuadraturePoints_elementBoundary, const int *elementBoundaryElementsArray, const int *quadraturePointToElementBoundary, const int *materialTypes, const double rho, const double beta, const double *gravity, const double *alpha, const double *n_vg, const double *thetaR, const double *thetaSR, const double *KWs, const double *u, const double *gradu, const double *n_global, const double *dV, double *mass, double *dmass, double *f_avg, double *df_avg, double *a_avg, double *da_avg, double *f, double *df, double *a, double *da)
Definition transportCoefficients.c:11323
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
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
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
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
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
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
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
conservativeHeadRichardsJLeverettAni
void conservativeHeadRichardsJLeverettAni(const int nSimplex, const int nPointsPerSimplex, const int nSpace, const int *materialTypes, const double rho, const double beta, const double *gravity, const double *phi, const double *psiD, const double *ns, const double *nk, const double *S_wirr, const double *S_nwr, const double *kr0x, const double *kr0y, const double *kr0z, double *u, double *mass, double *dmass, double *f, double *df, double *a, double *da)
Definition transportCoefficients.c:6220
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
l2projectVector
void l2projectVector(const int nSimplices, const int nPointsPerSimplex, const int nSpace, double *dV, double *r)
Definition transportCoefficients.c:5676
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)
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
conservativeHeadRichardsMualemVanGenuchtenHetEvaluateV2withUpwindAndHarm
void conservativeHeadRichardsMualemVanGenuchtenHetEvaluateV2withUpwindAndHarm(const int upwindFlag, const int computeAverages, const int nSimplex, const int nPointsPerSimplex, const int nSpace, const int nQuadraturePoints_elementBoundary, const int *elementBoundaryElementsArray, const int *quadraturePointToElementBoundary, const int *materialTypes, const double rho, const double beta, const double *gravity, const double *alpha, const double *n_vg, const double *thetaR, const double *thetaSR, const double *KWs, const double *u, const double *gradu, const double *n_global, const double *dV, double *mass, double *dmass, double *f_avg, double *df_avg, double *a_avg, double *da_avg, double *f, double *df, double *a, double *da)
Definition transportCoefficients.c:11515
groundwaterTransportCoefficientsEvaluate_hetMat
void groundwaterTransportCoefficientsEvaluate_hetMat(const int nSimplex, const int nPointsPerSimplex, const int nSpace, const double d, const int *materialTypes, const double *omega_types, const double *alpha_L_types, const double *alpha_T_types, const double *v, const double *u, double *m, double *dm, double *f, double *df, double *a)
Definition transportCoefficients.c:413
ReynoldsAveragedNavierStokes_kEpsilon_3D_Update
void ReynoldsAveragedNavierStokes_kEpsilon_3D_Update(const int nPoints, const double nu, const double c_mu, const double *k, const double *grad_k, const double *epsilon, 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)
Definition transportCoefficients.c:10932
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
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
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
ReynoldsAveragedNavierStokes_kEpsilon_2D_Update_sd
void ReynoldsAveragedNavierStokes_kEpsilon_2D_Update_sd(const int nPoints, const double rho, const double nu, const double c_mu, const double *k, const double *grad_k, const double *epsilon, 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)
Definition transportCoefficients.c:10890
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
l2projectScalar
void l2projectScalar(const int nSimplices, const int nPointsPerSimplex, double *dV, double *r)
Definition transportCoefficients.c:5651
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
kEpsilon_k_3D_Evaluate_sd
void kEpsilon_k_3D_Evaluate_sd(int nPoints, int nSpace, double sigma_k, double c_mu, double nu, double *velocity, double *gradu, double *gradv, double *gradw, double *k, double *epsilon, double *m_k, double *dm_k, double *phi_k, double *dphi_k, double *f_k, double *df_k, double *a_k, double *da_k_dk, double *r_k, double *dr_k_dk)
Definition transportCoefficients.c:10606
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
levelSetCurvatureCoefficientsEvaluate
void levelSetCurvatureCoefficientsEvaluate(int nPoints, int nSpace, double *grad_phi, double *u, double *f, double *r, double *dr)
Definition transportCoefficients.c:1270
VolumeAveragedTwophaseNavierStokes_ST_LS_SO_2D_Evaluate_sd
void VolumeAveragedTwophaseNavierStokes_ST_LS_SO_2D_Evaluate_sd(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:9079
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
Mass_2D_Evaluate
void Mass_2D_Evaluate(const int nPoints, double rho, double *p, double *u, double *v, double *mom_p_acc, double *mom_u_acc, double *mom_v_acc, double *dmom_p_acc_p, double *dmom_u_acc_u, double *dmom_v_acc_v)
Definition transportCoefficients.c:12626
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
smoothedHeaviside_integral
double smoothedHeaviside_integral(double eps, double phi)
Definition numericalFlux.c:14181
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
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
kEpsilon_3D_Evaluate
void kEpsilon_3D_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 *gradw, 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:10329
shallowWater_1D_Evaluate
void shallowWater_1D_Evaluate(const int nPoints, const double h_eps, const double g, const double bedFrictionCoefficient, const double bedFrictionPower, const double eddyViscosity, const double *x, const double *db_dx, const double *h, const double *hu, double *H, double *mass_acc, double *dmass_acc_dh, double *mom_acc, double *dmom_acc_dhu, double *mass_adv, double *dmass_adv_dhu, double *mom_adv, double *dmom_adv_dh, double *dmom_adv_dhu, double *mom_source, double *dmom_source_dh, double *dmom_source_dhu, double *mom_diff)
Definition transportCoefficients.c:11107
MovingMesh_3D_Evaluate
void MovingMesh_3D_Evaluate(const int nPoints, const double E0, const double nu, const double *g, const double *det_J, 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:8321
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
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
TwoPhaseAdvection_2D_Evaluate
void TwoPhaseAdvection_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 *phi, const double *p, const double *u, const double *v, double *mass_adv, double *dmass_adv_p, 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)
Definition transportCoefficients.c:12714
TwophaseNavierStokes_ST_LS_SO_3D_Evaluate_sd
void TwophaseNavierStokes_ST_LS_SO_3D_Evaluate_sd(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:3568
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)
ReynoldsAveragedNavierStokes_kEpsilon_3D_Update_sd
void ReynoldsAveragedNavierStokes_kEpsilon_3D_Update_sd(const int nPoints, const double nu, const double c_mu, const double *k, const double *grad_k, const double *epsilon, 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)
Definition transportCoefficients.c:10996
variablySaturatedGroundwaterTransportCoefficientsEvaluate_hetMat
void variablySaturatedGroundwaterTransportCoefficientsEvaluate_hetMat(const int nSimplex, const int nPointsPerSimplex, const int nSpace, const double d, const int *materialTypes, const double *theta, const double *alpha_L_types, const double *alpha_T_types, const double *v, const double *u, double *m, double *dm, double *f, double *df, double *a)
Definition transportCoefficients.c:468
MovingMesh_1D_Evaluate
void MovingMesh_1D_Evaluate(const int nPoints, const double E0, const double nu, const double *g, const double *det_J, const double *u, double *uu_diff_ten, double *u_force)
Definition transportCoefficients.c:8255
redistanceLevelSetCoefficientsWithWeakPenaltyEvaluate
void redistanceLevelSetCoefficientsWithWeakPenaltyEvaluate(int nPoints, int nSpace, double eps, double lambda_penalty, double *u_levelSet, double *u, double *grad_u, double *m, double *dm, double *H, double *dH, double *r, double *dr)
Definition transportCoefficients.c:1443
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
smoothedDirac
double smoothedDirac(double eps, double phi)
Definition numericalFlux.c:14199
kEpsilon_2D_Evaluate_sd
void kEpsilon_2D_Evaluate_sd(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:10090
eddyViscosity_3D_Update
void eddyViscosity_3D_Update(const int nPoints, const double *nu_t, 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)
Definition transportCoefficients.c:12362
VolumeAveragedTwophaseNavierStokes_ST_LS_SO_3D_Evaluate
void VolumeAveragedTwophaseNavierStokes_ST_LS_SO_3D_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 *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:9322
applyContactLineSlipJacobian
void applyContactLineSlipJacobian(int nExteriorElementBoundaries_global, int nQuadraturePoints_elementBoundary, int nDOF_trial_element, double eps, int *isDOFBoundary, double *phi, double *fluxJacobian)
Definition transportCoefficients.c:11948
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
calculateEddyViscosity_Smagorinsky2P_2D
void calculateEddyViscosity_Smagorinsky2P_2D(const int nElements_global, const int nQuadraturePoints_element, const double smagorinskyConstant_0, const double smagorinskyConstant_1, const double eps, const double *phi_ls, const double *h_e, const double *grad_u, const double *grad_v, double *nu_t)
Definition transportCoefficients.c:12165
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
MovingMesh_2D_Evaluate
void MovingMesh_2D_Evaluate(const int nPoints, const double E0, const double nu, const double *g, const double *det_J, 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:8276
calculateEddyViscosity_Smagorinsky2P_3D
void calculateEddyViscosity_Smagorinsky2P_3D(const int nElements_global, const int nQuadraturePoints_element, const double smagorinskyConstant_0, const double smagorinskyConstant_1, const double eps, const double *phi_ls, const double *h_e, const double *grad_u, const double *grad_v, const double *grad_w, double *nu_t)
Definition transportCoefficients.c:12219
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
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
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)
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
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
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)
w
#define w(x)
Definition jf.h:22
pd
#define pd(x)
Definition jf.h:24
Generated on for proteus by doxygen 1.18.0