You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
Copy file name to clipboardExpand all lines: docs/documentation/case.md
+95-35Lines changed: 95 additions & 35 deletions
Display the source diff
Display the rich diff
Original file line number
Diff line number
Diff line change
@@ -522,6 +522,9 @@ To restart the simulation from $k$-th time step, see [Restarting Cases](running.
522
522
|`probe_wrt`| Logical | Write the flow chosen probes data files for each time step |
523
523
|`num_probes`| Integer | Number of probes |
524
524
|`probe(i)%[x,y,z]`| Real | Coordinates of probe $i$ |
525
+
|`output_partial_domain`| Logical | Output part of the domain |
526
+
|`[x,y,z]_output%beg`| Real | Beginning of the output domain in the [x,y,z]-direction |
527
+
|`[x,y,z]_output%end`| Real | End of the output domain in the [x,y,z]-direction |
525
528
526
529
The table lists formatted database output parameters. The parameters define variables that are outputted from simulation and file types and formats of data as well as options for post-processing.
527
530
@@ -549,6 +552,11 @@ If `file_per_process` is true, then pre_process, simulation, and post_process mu
549
552
550
553
-`probe_wrt` activates the output of state variables at coordinates specified by `probe(i)%[x;y,z]`.
551
554
555
+
-`output_partial_domain` activates the output of part of the domain specified by `[x,y,z]_output%beg` and `[x,y,z]_output%end`.
556
+
This is useful for large domains where only a portion of the domain is of interest.
557
+
It is not supported when `precision = 1` and `format = 1`.
558
+
It also cannot be enabled with `flux_wrt`, `heat_ratio_wrt`, `pres_inf_wrt`, `c_wrt`, `omega_wrt`, `ib`, `schlieren_wrt`, or `qm_wrt`.
559
+
552
560
### 8. Acoustic Source {#acoustic-source}
553
561
554
562
| Parameter | Type | Description |
@@ -630,56 +638,73 @@ Details of the transducer acoustic source model can be found in [Maeda and Colon
630
638
631
639
-`%%bb_lowest_freq` specifies the lower frequency bound of the broadband acoustic wave. The upper frequency bound will be calculated as `%%bb_lowest_freq + %%bb_num_freq * %%bb_bandwidth`. The wave is no longer broadband below the lower bound and above the upper bound.
|`mu_l0` * | Real | Liquid viscosity (only specify in liquid phase) |
649
+
|`ss` * | Real | Surface tension (only specify in liquid phase) |
650
+
|`pv` * | Real | Vapor pressure (only specify in liquid phase) |
651
+
|`gamma_v` † | Real | Specific heat ratio |
652
+
|`M_v` † | Real | Molecular weight |
653
+
|`mu_v` † | Real | Viscosity |
654
+
|`k_v` † | Real | Thermal conductivity |
655
+
|`cp_v` † | Real | Specific heat capacity |
656
+
657
+
These options work only for gas-liquid two-component flows.
658
+
Component indexes are required to be 1 for liquid and 2 for gas.
659
+
660
+
-\* These parameters should be prepended with patch index $1$ that is filled with liquid: `fluid_pp(1)%`.
661
+
- † These parameters should be prepended with patch indexes filled with liquid and gas: `fluid_pp(1)%` and `fluid_pp(2)%`.
662
+
663
+
This table lists the sub-grid bubble model parameters, which can be utilized in both the ensemble-averaged and volume-averaged bubble models.
664
+
665
+
-`bubbles_euler` activates the ensemble-averaged bubble model.
666
+
667
+
-`bubbles_lagrange` activates the volume-averaged bubble model.
668
+
669
+
-`bubble_model` specified a model for spherical bubble dynamics by an integer of 1 and 2.
670
+
`bubble_model = 1`, `2`, and `3` correspond to the Gilmore, Keller-Miksis, and Rayleigh-Plesset models.
671
+
672
+
-`mu_l0`, `ss`, and `pv`, `gamma_v`, `M_v`, `mu_v`, `k_v`, and `cp_v` specify simulation parameters for the non-polytropic gas compression model.
673
+
`mu_l0`, `ss`, and `pv` correspond to the liquid viscosity, surface tension, and vapor pressure, respectively.
674
+
`gamma_v`, `M_v`, `mu_v`, `k_v`, and `cp_v` specify the specific heat ratio, molecular weight, viscosity, thermal conductivity and specific heat capacity of a chosen component (`cp_v` only for ``bubbles_lagrange = 'T'``).
675
+
Implementation of the parameters into the model follow [Ando (2010)](references.md#Ando10).
|`polydisperse`| Logical | Polydispersity in equilibrium bubble radius R0|
685
+
|`polydisperse`| Logical | Polydispersity in equilibrium bubble radius R0|
643
686
|`nb`| Integer | Number of bins: [1] Monodisperse; [$>1$] Polydisperse |
644
687
|`poly_sigma`| Real | Standard deviation for probability density function of polydisperse bubble populations |
645
688
|`R0_type`| Integer | Quadrature rule for probability density function of polydisperse bubble populations |
646
689
|`Ca`| Real | Cavitation number |
647
690
|`Web`| Real | Weber number |
648
691
|`Re_inv`| Real | Inverse Reynolds number |
649
-
|`mu_l0` * | Real | Liquid viscosity (only specify in liquid phase) |
650
-
|`ss` * | Real | Surface tension (only specify in liquid phase) |
651
-
|`pv` * | Real | Vapor pressure (only specify in liquid phase) |
652
-
|`gamma_v` † | Real | Specific heat ratio |
653
-
|`M_v` † | Real | Molecular weight |
654
-
|`mu_v` † | Real | Viscosity |
655
-
|`k_v` † | Real | Thermal conductivity |
656
-
|`qbmm`| Logical | Quadrature by method of moments|
657
-
|`dist_type`| Integer | Joint probability density function for bubble radius and velocity (only for ``qbmm = 'T'``)|
658
-
|`sigR`| Real | Standard deviation for the probability density function of bubble radius (only for ``qbmm = 'T'``) |
659
-
|`sigV`| Real | Standard deviation for the probability density function of bubble velocity (only for ``qbmm = 'T'``) |
692
+
|`qbmm`| Logical | Quadrature by method of moments |
693
+
|`dist_type`| Integer | Joint probability density function for bubble radius and velocity (only for ``qbmm = 'T'``) |
694
+
|`sigR`| Real | Standard deviation for the probability density function of bubble radius (only for ``qbmm = 'T'``) (EE)|
695
+
|`sigV`| Real | Standard deviation for the probability density function of bubble velocity (only for ``qbmm = 'T'``) (EE)|
660
696
|`rhoRV`| Real | Correlation coefficient for the joint probability density function of bubble radius and velocity (only for ``qbmm = 'T'``) |
661
697
662
-
These options work only for gas-liquid two-component flows.
663
-
Component indexes are required to be 1 for liquid and 2 for gas.
664
-
665
-
-\* These parameters should be prepended with patch index $1$ that is filled with liquid: `fluid_pp(1)%`.
666
-
- † These parameters should be prepended with patch indexes filled with liquid and gas: `fluid_pp(1)%` and `fluid_pp(2)%`.
667
-
668
698
This table lists the ensemble-averaged bubble model parameters.
669
699
670
-
-`bubbles` activates the ensemble-averaged bubble model.
671
-
672
-
-`bubble_model` specified a model for spherical bubble dynamics by an integer of 1 and 2.
673
-
`bubble_model = 1`, `2`, and `3` correspond to the Gilmore, Keller-Miksis, and Rayleigh-Plesset models.
674
-
675
700
-`polytropic` activates polytropic gas compression in the bubble.
676
701
When ``polytropic = 'F'``, the gas compression is modeled as non-polytropic due to heat and mass transfer across the bubble wall with constant heat and mass transfer coefficients based on ([Preston et al., 2007](references.md#Preston07)).
677
702
678
-
-`polydisperse` activates polydispersity in the bubble model through a probability density function (PDF) of the equilibrium bubble radius.
679
-
680
703
-`thermal` specifies a model for heat transfer across the bubble interface by an integer from 1 through 3.
681
704
`thermal = 1`, `2`, and `3` correspond to no heat transfer (adiabatic gas compression), isothermal heat transfer, and heat transfer with a constant heat transfer coefficient based on [Preston et al., 2007](references.md#Preston07), respectively.
682
705
706
+
-`polydisperse` activates polydispersity in the bubble model through a probability density function (PDF) of the equilibrium bubble radius.
707
+
683
708
-`R0ref` specifies the reference bubble radius.
684
709
685
710
-`nb` specifies the number of discrete bins that define the probability density function (PDF) of the equilibrium bubble radius.
@@ -691,11 +716,6 @@ When ``polytropic = 'F'``, the gas compression is modeled as non-polytropic due
691
716
692
717
-`Ca`, `Web`, and `Re_inv` respectively specify the Cavitation number, Weber number, and the inverse Reynolds number that characterize the offset of the gas pressure from the vapor pressure, surface tension, and liquid viscosity when the polytropic gas compression model is used.
693
718
694
-
-`mu_l0`, `ss`, and `pv`, `gamma_v`, `M_v`, `mu_v`, and `k_v` specify simulation parameters for the non-polytropic gas compression model.
695
-
`mu_l0`, `ss`, and `pv` correspond to the liquid viscosity, surface tension, and vapor pressure, respectively.
696
-
`gamma_v`, `M_v`, `mu_v`, and `k_v` specify the specific heat ratio, molecular weight, viscosity, and thermal conductivity of a chosen component.
697
-
Implementation of the parameters into the model follow [Ando (2010)](references.md#Ando10).
698
-
699
719
-`qbmm` activates quadrature by method of moments, which assumes a PDF for bubble radius and velocity.
700
720
701
721
-`dist_type` specifies the initial joint PDF of initial bubble radius and bubble velocity required in qbmm. `dist_type = 1` and `2` correspond to binormal and lognormal-normal distributions respectively.
@@ -706,6 +726,46 @@ Implementation of the parameters into the model follow [Ando (2010)](references.
706
726
707
727
-`rhoRV` specifies the correlation coefficient of the joint PDF of bubble radius and bubble velocity required in the QBMM feature.
708
728
729
+
#### 9.2 Volume-Averaged Bubble Model
730
+
731
+
| Parameter | Type | Description |
732
+
| ---: | :---: | :--- |
733
+
|`bubbles_lagrange`| Logical | Lagrangian subgrid bubble model switch |
734
+
|`nBubs_glb`| Integer | Global number of bubbles |
|`heatTransfer_model`| Logical | Activates the interface heat transfer model |
740
+
|`massTransfer_model`| Logical | Activates the interface mass transfer model |
741
+
|`write_bubbles`| Logical | Write files to track the bubble evolution each time step |
742
+
|`write_bubbles_stats`| Logical | Write the maximum and minimum radius of each bubble |
743
+
|`epsilonb`| Real | Standard deviation scaling for the gaussian function |
744
+
|`charwidth`| Real | Domain virtual depth (z direction, for 2D simulations) |
745
+
|`valmaxvoid`| Real | Maximum void fraction permitted |
746
+
|`c0`| Real | Reference speed |
747
+
|`rho0`| Real | Reference density |
748
+
|`T0`| Real | Reference temperature |
749
+
|`x0`| Real | Reference length |
750
+
|`Thost`| Real | Temperature of the surrounding liquid (host) |
751
+
|`diffcoefvap`| Real | Vapor diffusivity in the gas |
752
+
|`rkck_adap_dt`| Logical | Activates the adaptive rkck time stepping algorithm |
753
+
|`rkck_tolerance`| Real | Admissible error truncation tolerance in the rkck stepper |
754
+
755
+
-`nBubs_glb` Total number of bubbles. Their initial conditions need to be specified in the ./input/lag_bubbles.dat file. See the example cases for additional information.
756
+
757
+
-`solver_approach` Specifies the Euler-Lagrange coupling method: [1] enables a one-way coupling approach, where the bubbles do not influence the Eulerian field. [2] activates the two-way coupling approach based on [Maeda and Colonius (2018)](references.md#Maeda18), where the effect of the bubbles is added in the Eulerian field as source terms.
758
+
759
+
-`cluster_type` Specifies method to find p_inf (pressure that drives the bubble dynamics): [1] activates the bilinear interpolation of the pressure field, while [2] enables the bubble dynamic closure based on [Maeda and Colonius (2018)](references.md#Maeda18), the full model is obtained when `pressure_corrector` is true.
760
+
761
+
-`smooth_type` Specifies the smoothening method of projecting the lagrangian bubbles in the Eulerian field: [1] activates the gaussian kernel function described in [Maeda and Colonius (2018)](references.md#Maeda18), while [2] activates the delta kernel function where the effect of the bubble is only seen in the specific bubble location cell.
762
+
763
+
-`heatTransfer_model` Activates the heat transfer model at the bubble's interface based on ([Preston et al., 2007](references.md#Preston07)).
764
+
765
+
-`massTransfer_model` Activates the mass transfer model at the bubble's interface based on ([Preston et al., 2007](references.md#Preston07)).
766
+
767
+
-`rkck_adap_dt` Activates the adaptive 4th/5th order Runge—Kutta–Cash–Karp (RKCK) time-stepping algorithm (requires `time_stepper ==4`). A maximum error between the 4th and 5th order Runge-Kutta-Cash-Karp solutions for the same time step size is calculated. If the error is smaller than a tolerance (`rkck_tolerance`), then the algorithm employs the 5th order solution, while if not, both eulerian/lagrangian variables are re-calculated with a smaller time step size.
Copy file name to clipboardExpand all lines: docs/documentation/references.md
+2Lines changed: 2 additions & 0 deletions
Display the source diff
Display the rich diff
Original file line number
Diff line number
Diff line change
@@ -30,6 +30,8 @@
30
30
31
31
- <aid="Maeda17">Maeda, K. and Colonius, T. (2017). A source term approach for generation of one-way acoustic waves in the euler and navier–stokes equations. Wave Motion, 75:36–49.</a>
32
32
33
+
- <aid="Maeda18">Maeda, K. and Colonius, T. (2018). Eulerian–lagrangian method for simulation of cloud cavitation. Journal of computational physics, 371:994–1017.</a>
34
+
33
35
- <aid="Meng16">Meng, J. C. C. (2016). Numerical simulations of droplet aerobreakup. PhD thesis, California Institute of Technology.</a>
34
36
35
37
- <aid="Pirozzoli13">Pirozzoli, S., and Colonius, T. (2013). Generalized characteristic relaxation boundary conditions for unsteady compressible flow simulations. Journal of Computational Physics, 248:109-126.</a>
0 commit comments