CAA_OUTPUT
Specifies output parameters for Computational Aero-Acoustics (CAA) post-processing.
Type
AcuSolve Command
Syntax
CAA_OUTPUT("name") {parameters...}
Qualifier
User-given name.
Parameters
- coordinates or coord (array) [no default]
- An array of sample points. Each row contains a unique identification number and the x, y, and z coordinates of a point in the global xyz coordinate system.
- directions or dir (array) [no default]
- An array of outward normal direction vectors corresponding to the coordinates array. Each row contains an identification number number and the x, y, and z components of a vector in the global xyz coordinate system. Required by some of the output quantities.
- gap (real) >=0 [=0]
- Tolerance for sample points falling outside of all elements.
- output_coordinates (boolean) [=off]
- Flag specifying whether to output the coordinates.
- output_velocity (boolean) [=off]
- Flag specifying whether to output the velocity.
- output_pressure (boolean) [=off]
- Flag specifying whether to output the pressure.
- output_momentum_stress (boolean) [=off]
- Flag specifying whether to output the momentum stress.
- output_total_stress (boolean) [=off]
- Flag specifying whether to output the total stress.
- output_lighthill_stress (boolean) [=off]
- Flag specifying whether to output the Lighthill stress.
- output_reduced_lighthill_stress (boolean) [=off]
- Flag specifying whether to output the reduced Lighthill stress.
- output_div_momentum_stress (boolean) [=off]
- Flag specifying whether to output the divergence of the momentum stress.
- output_div_total_stress (boolean) [=off]
- Flag specifying whether to output the divergence of the total stress.
- output_div_lighthill_stress (boolean) [=off]
- Flag specifying whether to output the divergence of the Lighthill stress.
- output_div_reduced_lighthill_stress (boolean) [=off]
- Flag specifying whether to output the divergence of the reduced Lighthill stress.
- output_normal_momentum_flux (boolean) [=off]
- Flag specifying whether to output the normal component of the momentum flux.
- output_normal_momentum_flux_rate (boolean) [=off]
- Flag specifying whether to output the normal component of the momentum flux rate.
- output_normal_div_total_stress (boolean) [=off]
- Flag specifying whether to output the normal component of the divergence of the total stress.
- output_fwh_monopole (boolean) [=off]
- Flag specifying whether to output the monopole term of the Ffowcs Williams Hawkings acoustics equations.
- output_fwh_dipole (boolean) [=off]
- Flag specifying whether to output the dipole term of the Ffowcs Williams Hawkings acoustics equations.
- output_ei_na (boolean) [=off]
- Flag specifying whether to output the element-integrated volume.
- output_ei_na_velocity (boolean) [=off]
- Flag specifying whether to output the element-integrated velocity.
- output_ei_na_pressure (boolean) [=off]
- Flag specifying whether to output the element-integrated pressure.
- output_ei_na_momentum_stress (boolean) [=off]
- Flag specifying whether to output the element-integrated momentum stress.
- output_ei_na_total_stress (boolean) [=off]
- Flag specifying whether to output the element-integrated total stress.
- output_ei_na_lighthill_stress (boolean) [=off]
- Flag specifying whether to output the element-integrated Lighthill stress.
- output_ei_na_reduced_lighthill_stress (boolean) [=off]
- Flag specifying whether to output the element-integrated reduced Lighthill stress.
- output_ei_na_div_momentum_stress (boolean) [=off]
- Flag specifying whether to output the element-integrated divergence of the momentum stress.
- output_ei_na_div_total_stress (boolean) [=off]
- Flag specifying whether to output the element-integrated divergence of the total stress.
- output_ei_na_div_lighthill_stress (boolean) [=off]
- Flag specifying whether to output the element-integrated divergence of the Lighthill stress.
- output_ei_na_div_reduced_lighthill_stress (boolean) [=off]
- Flag specifying whether to output the element-integrated divergence of the reduced Lighthill stress.
- output_ei_gna_div_momentum_stress (boolean) [=off]
- Flag specifying whether to output the element-integrated shape function gradient contracted with the divergence of the momentum stress.
- output_ei_gna_div_total_stress (boolean) [=off]
- Flag specifying whether to output the element-integrated shape function gradient contracted with the divergence of the total stress.
- output_ei_gna_div_lighthill_stress (boolean) [=off]
- Flag specifying whether to output the element-integrated shape function gradient contracted with the divergence of the Lighthill stress.
- output_ei_gna_div_reduced_lighthill_stress (boolean) [=off]
- Flag specifying whether to output the element-integrated shape function gradient contracted with the divergence of the reduced Lighthill stress.
- output_si_na (boolean) [=off]
- Flag specifying whether to output the surface-integrated area.
- output_si_na_normal (boolean) [=off]
- Flag specifying whether to output the surface-integrated normal.
- output_si_na_pressure (boolean) [=off]
- Flag specifying whether to output the surface-integrated pressure.
- output_si_na_normal_pressure (boolean) [=off]
- Flag specifying whether to output the surface-integrated normal-pressure.
- output_si_na_normal_momentum_flux (boolean) [=off]
- Flag specifying whether to output the surface-integrated normal component of the momentum flux.
- output_si_na_normal_momentum_flux_rate (boolean) [=off]
- Flag specifying whether to output the surface-integrated normal component of the momentum flux rate.
- output_si_na_normal_div_total_stress (boolean) [=off]
- Flag specifying whether to output the surface-integrated normal component of the divergence of the total stress.
- output_si_na_fwh_monopole (boolean) [=off]
- Flag specifying whether to output the surface-integrated monopole term of the Ffowcs Williams Hawkings acoustics equations.
- output_si_na_fwh_dipole (boolean) [=off]
- Flag specifying whether to output the surface-integrated dipole term of the Ffowcs Williams Hawkings acoustics equations.
- reference_density or ref_dens (real) [=1.225]
- Reference density. Required by some of the output quantities.
- reference_pressure or ref_pres (real) [=0]
- Reference pressure. Required by some of the output quantities.
- reference_sound_speed or sound (real) [=330]
- Reference sound speed. Required by some of the output quantities.
- output_frequency or out_freq (integer) >=0 [=1000]
- Time step frequency at which to output CAA data. If zero, this option is ignored.
- output_time_interval or out_intv (real) >=0 [=0]
- Time frequency at which to output CAA data. If zero, this option is ignored.
- num_saved_states or save (integer) >=0 [=0]
- Number of CAA data files to keep. If zero, all are kept.
Description
CAA_OUTPUT( "sample points" ) {
coordinates = { 1, 0., 0., 0. ;
2, 1., 0., 0. ;
11, 0., 1., 0. ;
12, 1., 1., 0. ; }
directions = { 1, 0., 0., 1. ;
2, 0., 0., 1. ;
11, 0., 0., 1. ;
12, 0., 0., 1. ; }
gap = 0
output_coordinates = on
output_velocity = on
output_pressure = off
output_momentum_stress = off
output_total_stress = off
output_lighthill_stress = off
output_reduced_lighthill_stress = off
output_div_momentum_stress = off
output_div_total_stress = off
output_div_lighthill_stress = off
output_div_reduced_lighthill_stress = off
output_normal_momentum_flux = off
output_normal_momentum_flux_rate = off
output_normal_div_total_stress = off
output_fwh_monopole = off
output_fwh_dipole = off
output_ei_na = off
output_ei_na_velocity = off
output_ei_na_pressure = off
output_ei_na_momentum_stress = off
output_ei_na_total_stress = off
output_ei_na_lighthill_stress = off
output_ei_na_reduced_lighthill_stress = off
output_ei_na_div_momentum_stress = off
output_ei_na_div_total_stress = off
output_ei_na_div_lighthill_stress = off
output_ei_na_div_reduced_lighthill_stress = off
output_ei_gna_div_momentum_stress = off
output_ei_gna_div_total_stress = off
output_ei_gna_div_lighthill_stress = off
output_ei_gna_div_reduced_lighthill_stress = off
output_si_na = off
output_si_na_normal = off
output_si_na_normal_pressure = off
output_si_na_normal_momentum_flux = off
output_si_na_normal_momentum_flux_rate = off
output_si_na_normal_div_total_stress = off
output_si_na_fwh_monopole = off
output_si_na_fwh_dipole = off
reference_density = 1.225
reference_pressure = 0
reference_sound_speed = 330
output_frequency = 1
output_time_interval = 0
num_saved_states = 0
}
A set of sample points is specified by the coordinates parameter, which is a four-column array corresponding to an identification number and the x, y, and z coordinates of each sample point, one row per point. The identification number may be any arbitrary but unique (within the set) integer. Negative numbers are acceptable. The coordinates are in the global xyz coordinate system. This is the same format as the coordinates parameter of the COORDINATE command, but the sample points do not have to be a subset of the nodal coordinates nor related to them in any other way. The coordinates parameter is mandatory. Within AcuPrep, an element is found for each sample point and interpolation shape functions are computed. If a sample point falls outside all fluid elements, then the closest element face is found, up to a distance specified by gap. An error is issued if no element is found.
Some of the output fields require a direction vector. If needed, the directions parameter must be given as a four column array: identification number, and the x, y, and z components of the direction. The identification numbers must be the same set as given by coordinates. The parameters reference_density, reference_pressure and reference_sound_speed provide the required data for certain fields. Any of these extra parameters that a particular field requires are noted in the table below.
1 0. 0. 0.
2 1. 0. 0.
11 0. 1. 0.
12 1. 1. 0.
1 0. 0. 1.
2 0. 0. 1.
11 0. 0. 1.
12 0. 0. 1.
CAA_OUTPUT( "sample points" ) {
coordinates = Read( "channel.caa.crd" )
directions = Read( "channel.caa.dir" )
...
}
Since the volume of data written on disk by CAA_OUTPUT could be very large, this command provides the option of selecting which data are to be written on disk via the parameters output_coordinates through output_si_na_fwh_dipole. The default is to not output any fields. The output quantities and their required data are defined in the following table:
Name | DOFs | Definition | Data |
---|---|---|---|
coordinates | 3 | x | |
velocity | 3 | u | |
pressure | 1 | p | |
momentum_stress | 6 | ![]() |
|
total_stress | 6 | ![]() |
|
lighthill_stress | 6 | ![]() |
![]() |
reduced_lighthill_stress | 6 | ![]() |
|
div_momentum_stress | 3 | ![]() |
|
div_total_stress | 3 | ![]() |
|
div_lighthill_stress | 3 | ![]() |
![]() |
div_reduced_lighthill_stress | 3 | ![]() |
|
normal_momentum_flux | 1 | ![]() |
![]() |
normal_momentum_flux_rate | 1 | ![]() |
![]() |
normal_div_total_stress | 1 | ![]() |
![]() |
fwh_monopole | 1 | ![]() |
![]() |
fwh_dipole | 1 | ![]() |
![]() |
ei_na | 1 | ![]() |
|
ei_na_velocity | 3 | ![]() |
|
ei_na_pressure | 1 | ![]() |
|
ei_na_momentum_stress | 6 | ![]() |
|
ei_na_total_stress | 6 | ![]() |
|
ei_na_lighthill_stress | 6 | ![]() |
![]() |
ei_na_reduced_lighthill_stress | 6 | ![]() |
|
ei_na_div_momentum_stress | 3 | ![]() |
|
ei_na_div_total_stress | 3 | ![]() |
|
ei_na_div_lighthill_stress | 3 | ![]() |
![]() |
ei_na_div_reduced_lighthill_stress | 3 | ![]() |
|
ei_gna_div_momentum_stress | 1 | ![]() |
|
ei_gna_div_total_stress | 1 | ![]() |
|
ei_gna_div_lighthill_stress | 1 | ![]() |
![]() |
ei_gna_div_reduced_lighthill_stress | 1 | ![]() |
|
si_na | 1 | ![]() |
|
si_na_normal | 3 | ![]() |
|
si_na_pressure | 1 | ![]() |
|
si_na_normal_pressure | 3 | ![]() |
|
si_na_normal_momentum_flux | 1 | ![]() |
|
si_na_normal_momentum_flux_rate | 1 | ![]() |
![]() |
si_na_normal_div_total_stress | 1 | ![]() |
![]() |
si_na_fwh_monopole | 1 | ![]() |
![]() |
si_na_fwh_dipole | 1 | ![]() |
![]() |
where x is the coordinate vector;
u is the velocity;
v is the mesh velocity;
n is the outward normal direction, given by
directions; p is the pressure; is the reference pressure,
given by reference_pressure; ρ is the density;
is the reference density,
given by reference_density; c is the reference speed
of sound, given by reference_sound_speed;
d is the Kronecker delta (identity matrix);
is the viscous stress tensor
(including turbulence effects); and Na is the element shape
function. The stress-type output quantities with 6 DOFs are ordered as xx, yy, zz, xy, yz, zx.
Due to the mesh velocity term the fwh_monopole variables do not contain the
usual time derivative.
The variables starting with "ei_" are element-integrated. The elements that contribute to the integration are specified by the CAA_ELEMENT_OUTPUT commands that reference the qualifier of the CAA_OUTPUT command.
CAA_OUTPUT( "sample points" ) {
...
output_ei_na = on
output_ei_na_lighthill_stress = on
output_si_na = on
output_si_na_fwh_monopole = on
...
}
CAA_ELEMENT_OUTPUT( "sample elements" ) {
caa_output = "sample points"
shape = eight_node_brick
elements = Read( "sample.elm" )
}
CAA_SURFACE_OUTPUT( "sample surfaces" ) {
caa_output = "sample points"
shape = four_node_quad
partial_surfaces = Read( "sample.srf" )
}
The parameter output_frequency determines how often the CAA data will be output, in number of steps. output_time_interval does the same, except in time. If either of these parameters is non-zero, the CAA data will be output at the end of the run. If both are zero, no data are written to disk.
Run times may not coincide with output_time_interval. In these cases the data are output at every time step which passes through a multiple of output_time_interval.
The num_saved_states parameter indicates the number of CAA data files to save. Once the (num_saved_states + 1)th file is written, the first one is removed. When running in parallel, once all files are written (in parallel), then the old ones are removed (in parallel).
Once the CAA data have been written to disk, they can be translated to other formats using the AcuTrans program; see the AcuSolve Programs Reference Manual for details. Access to the CAA results are provided also via the C-callable API library, libadb.a (Acusim Data Base), and its python/perl/shell bindings/access programs. See variables "nOaas" through "oaaValues" (15 variables total) in doc/adb .txt.