Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “error field”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,225 records · Page 68

A theoretical investigation of the sound radiation fields associated with a Bellmouth inlet

Analytical results are obtained by numerical integration of a cylindrically symmetric integral representation of the external solutions of the Helmholtz equation. The accuracy of this method is checked by comparisons of computed results with 'exact' solutions generated by the point source method. In all cases, the average error for the amplitude and phase of the points calculated in the field is found to be less than ten percent. Theoretical studies which model experiments run for the NASA Langely Bellmouth inlet configuration are presented and comparisons are made with the experimental results. In all cases very good agreement is obtained between the experimental and theoretically calculated values.

Meyer, W. L.↗

Intercomparison of Deep Learning Model Architectures for Atmospheric River Prediction

With a rapid surge in the application of machine learning (ML) for a diverse range of tasks in climate science, the present study addresses a challenge for climate scientists when selecting the optimal ML or deep learning (DL) architecture for a given application. In particular, a DL intercomparison study was performed with a focus on forecasting the position of atmospheric rivers (ARs) on short-range time scales (up to 5-day lead times). AR predictions from multiple DL architectures, including various types of convolutional autoencoders and a vision transformer (ViT), were compared against ECMWF ERA5 reanalysis and hindcasts from a global climate model. DL models with similar trainable parameters were trained on ERA5 reanalysis data and AR positions derived from a thresholding algorithm to ensure a fair comparison among the DL models. Each model’s performance and accuracy in forecasting AR location and key input fields within a 5-day window were assessed using metrics of root-mean-square error, anomaly correlation, and mean intersection over union. The ViT architecture outperformed other autoencoder models in most of the metrics. Incorporating additional meteorological fields only yielded slight improvements in forecasting certain fields at longer lead times. The results also suggest that a smaller number of input time steps or smaller number of autoregressive steps can achieve better prediction skills, while also improving the overall computational efficiency. This research offers valuable insights into the strengths and weaknesses of different DL techniques for AR forecasting, hopefully guiding the development of improved models for forecasting this phenomenon.

54 ENVIRONMENTAL SCIENCES↗

Confidence set inference with a prior quadratic bound

In the uniqueness part of a geophysical inverse problem, the observer wants to predict all likely values of P unknown numerical properties z = (z sub 1,...,z sub p) of the earth from measurement of D other numerical properties y(0)=(y sub 1(0),...,y sub D(0)) knowledge of the statistical distribution of the random errors in y(0). The data space Y containing y(0) is D-dimensional, so when the model space X is infinite-dimensional the linear uniqueness problem usually is insoluble without prior information about the correct earth model x. If that information is a quadratic bound on x (e.g., energy or dissipation rate), Bayesian inference (BI) and stochastic inversion (SI) inject spurious structure into x, implied by neither the data nor the quadratic bound. Confidence set inference (CSI) provides an alternative inversion technique free of this objection. CSI is illustrated in the problem of estimating the geomagnetic field B at the core-mantle boundary (CMB) from components of B measured on or above the earth's surface. Neither the heat flow nor the energy bound is strong enough to permit estimation of B(r) at single points on the CMB, but the heat flow bound permits estimation of uniform averages of B(r) over discs on the CMB, and both bounds permit weighted disc-averages with continous weighting kernels. Both bounds also permit estimation of low-degree Gauss coefficients at the CMB. The heat flow bound resolves them up to degree 8 if the crustal field at satellite altitudes must be treated as a systematic error, but can resolve to degree 11 under the most favorable statistical treatment of the crust. These two limits produce circles of confusion on the CMB with diameters of 25 deg and 19 deg respectively.

Backus, George E.↗

Principles of error detection and error correction codes

Report is reviewed which considers theoretical basis of groups, rings, fields, and vector spaces, and their relationship to algebraic coding theory. Report serves as summary for engineers and scientists involved in data handling and processing systems.

Perlman, M.↗

Measurement of static pressure on aircraft

Existing data on the errors involved in the measurement of static pressure by means of static-pressure tubes and fuselage vents are presented. The errors associated with the various design features of static-pressure tubes are discussed for the condition of zero angle of attack and for the case where the tube is inclined to flow. Errors which result from variations in the configuration of static-pressure vents are also presented. Errors due to the position of a static-pressure tube in the flow field of the airplane are given for locations ahead of the fuselage nose, ahead of the wing tip, and ahead of the vertical tail fin. The errors of static-pressure vents on the fuselage of an airplane are also presented. Various methods of calibrating static-pressure installations in flight are briefly discussed.

Gracey, William↗

The inverse gravimetric problem in gravity modelling

One of the main purposes of geodesy is to determine the gravity field of the Earth in the space outside its physical surface. This purpose can be pursued without any particular knowledge of the internal density even if the exact shape of the physical surface of the Earth is not known, though this seems to entangle the two domains, as it was in the old Stoke's theory before the appearance of Molodensky's approach. Nevertheless, even when large, dense and homogeneous data sets are available, it was always recognized that subtracting from the gravity field the effect of the outer layer of the masses (topographic effect) yields a much smoother field. This is obviously more important when a sparse data set is bad so that any smoothing of the gravity field helps in interpolating between the data without raising the modeling error, this approach is generally followed because it has become very cheap in terms of computing time since the appearance of spectral techniques. The mathematical description of the Inverse Gravimetric Problem (IGP) is dominated mainly by two principles, which in loose terms can be formulated as follows: the knowledge of the external gravity field determines mainly the lateral variations of the density; and the deeper the density anomaly giving rise to a gravity anomaly, the more improperly posed is the problem of recovering the former from the latter. The statistical relation between rho and n (and its inverse) is also investigated in its general form, proving that degree cross-covariances have to be introduced to describe the behavior of rho. The problem of the simultaneous estimate of a spherical anomalous potential and of the external, topographic masses is addressed criticizing the choice of the mixed collection approach.

Sanso, F.↗

Optical and X-ray observations of the PSR 1913 + 16 field

We have obtained medium-resolution spectra in the yellow of seven stars in the error box of the binary radio pulsar PSR 1913 + 16, using the image-tube scanner at the Lick 3-m telescope. None of the objects is spectroscopically abnormal and thus a strong candidate for membership in the system. The data are used to establish the extinction in the field, which is modest for this low galactic latitude. For three of the stars, limits of the order V exceeding 23 mag on optical pulsations with the frequency of the radio pulsar are derived. X-ray observations of the field have been obtained in the 2.5-7.5 keV band, using the UCL experiment aboard OAO Copernicus. No source is detected, with a 2 sigma upper limit on the time-averaged flux of 1.0 times 10 to the -10th power ergs per sq cm per sec.

Davidsen, A.↗

Achieving Accuracy Requirements for Forest Biomass Mapping: A Spaceborne Data Fusion Method for Estimating Forest Biomass and Lidar Sampling Error

The synergistic use of active and passive remote sensing (i.e., data fusion) demonstrates the ability of spaceborne light detection and ranging (LiDAR), synthetic aperture radar (SAR) and multispectral imagery for achieving the accuracy requirements of a global forest biomass mapping mission (+/-20 Mg/ha or 20%, the greater of the two, for at least 80% of grid cells). A data fusion approach also provides a means to extend 3D information from discrete spaceborne LiDAR measurements of forest structure across scales much larger than that of the LiDAR footprint. For estimating biomass, these measurements mix a number of errors including those associated with LiDAR footprint sampling over regional-global extents. A general framework for mapping above ground live forest biomass density (AGB) with a data fusion approach is presented and verified using data from NASA field campaigns near Howland, ME, USA, to assess AGB and LiDAR sampling errors across a regionally representative landscape. We combined SAR and Landsat-derived optical (passive optical) image data to identify contiguous areas (>0.5 ha) that are relatively homogenous in remote sensing metrics (forest patches). We used this image-derived data with simulated spaceborne LiDAR derived from orbit and cloud cover simulations and airborne data from NASA's Laser Vegetation Imaging Sensor (LVIS) to compute AGB and estimate LiDAR sampling error for forest patches and 100 m, 250 m, 500 m, and 1 km grid cells. At both the patch and grid scales, we evaluated differences in AGB estimation and sampling error from the combined use of LiDAR with both SAR and passive optical and with either SAR or passive optical alone. First, this data fusion approach demonstrates that incorporating forest patches into the AGB mapping framework can provide sub-grid forest information for coarser grid-level AGB reporting. Second, a data fusion approach for estimating AGB using simulated spaceborne LiDAR with SAR and passive optical image combinations reduced forest AGB sampling errors 12%-38% from those where LiDAR is used with SAR or passive optical alone. In absolute terms, sampling errors were reduced from 14-40 Mg/ha to 11-28 Mg/ha across all grid scales and prediction methods, where minimum sampling errors were 11, 15, 18, and 22 Mg/ha for 1 km, 500 m, 250 m, and 100 m grid scales, respectively. Third, spaceborne global scale accuracy requirements were achieved whereby at least 80% of the grid cells at 100 m, 250 m, 500 m, and 1 km grid levels met AGB accuracy requirements using a combination of passive optical and SAR along with machine learning methods to predict vegetation structure metrics for forested areas without LiDAR samples. Finally, using either passive optical or SAR, accuracy requirements were met at the 500 m and 250 m grid level, respectively..

LiDAR↗

Brighter-fatter Effect in Near-infrared Detectors—III. Fourier-domain Treatment of Flat Field Correlations and Application to WFIRST

Weak gravitational lensing has emerged as a leading probe of the growth of cosmic structure. However, the shear signal is very small and accurate measurement depends critically on our ability to understand how non-ideal instrumental effects affect astronomical images. The Wide-Field Infrared Survey Telescope (WFIRST) will fly a focal plane containing 18 Teledyne H4RG-10 near-infrared detector arrays, which present different instrument calibration challenges from previous weak lensing observations. Previous work [Paper I: Hirata & Choi, PASP, 132, 014501 (2020); and Paper II: Choi & Hirata, PASP, 132, 014502 (2020)] has shown that correlation functions of flat field images, including cross-correlations between different time slices that are enabled by the non-destructive read capability of the infrared detectors, are effective tools for disentangling linear and nonlinear inter-pixel capacitance (IPC) and the brighter-fatter effect (BFE). Here we present a Fourier-domain treatment of the flat field correlations, which allows us to expand the previous formalism to all orders in IPC, BFE, and classical nonlinearity. We show that biases in simulated flat field analyses in Paper I are greatly reduced through the use of this formalism. We then apply this updated formalism to flat field data from three WFIRST flight candidate detectors, and explore the robustness to variations in the analysis. We find that the BFE is present in all three detectors, and that its contribution to the flat field correlations dominates over the nonlinear IPC, in accordance with the results from Paper II on a development detector. The magnitude of the BFE is such that the effective area of a pixel is increased by (3.54 ± 0.03) × 10 -7 for every electron deposited in a neighboring pixel (sensor chip assembly [SCA] 20829, statistical error, not IPC-deconvolved). We compare IPC maps from flat field autocorrelation measurements to those obtained from the single pixel reset method and find a median difference of 0.113% for SCA 20829. After further diagnosis of this difference, we ascribe it largely to an additional source of cross-talk, the vertical trailing pixel effect, and recommend further work to develop a model for this effect. These results represent a significant step toward calibration of the non-ideal effects in WFIRST detectors.

79 ASTRONOMY AND ASTROPHYSICS↗

The CCSDS Next Generation Space Data Link Protocol (NGSLP)

The CCSDS space link protocols i.e., Telemetry (TM), Telecommand (TC), Advanced Orbiting Systems (AOS) were developed in the early growth period of the space program. They were designed to meet the needs of the early missions, be compatible with the available technology and focused on the specific link environments. Digital technology was in its infancy and spacecraft power and mass issues enforced severe constraints on flight implementations. Therefore the Telecommand protocol was designed around a simple Bose, Hocquenghem, Chaudhuri (BCH) code that provided little coding gain and limited error detection but was relatively simple to decode on board. The infusion of the concatenated Convolutional and Reed-Solomon codes5 for telemetry was a major milestone and transformed telemetry applications by providing them the ability to more efficiently utilize the telemetry link and its ability to deliver user data. The ability to significantly lower the error rates on the telemetry links enabled the use of packet telemetry and data compression. The infusion of the high performance codes for telemetry was enabled by the advent of digital processing, but it was limited to earth based systems supporting telemetry. The latest CCSDS space link protocol, Proximity-1 was developed in early 2000 to meet the needs of short-range, bi-directional, fixed or mobile radio links characterized by short time delays, moderate but not weak signals, and short independent sessions. Proximity-1 has been successfully deployed on both NASA and ESA missions at Mars and is planned to be utilized by all Mars missions in development. A new age has arisen, one that now provides the means to perform advanced digital processing in spacecraft systems enabling the use of improved transponders, digital correlators, and high performance forward error correcting codes for all communications links. Flight transponders utilizing digital technology have emerged and can efficiently provide the means to make the next leap in performance for space link communications. Field Programmable Gate Arrays (FPGAs) provide the capability to incorporate high performance forward error correcting codes implemented within software transponders providing improved performance in data transfer, ranging, link security, and time correlation. Given these synergistic technological breakthroughs, the time has come to take advantage of them in applying them to both on going (e.g., command, telemetry) and emerging (e.g., space link security, optical communication) space link applications. However one of the constraining factors within the Data Link Layer in realizing these performance gains is the lack of a generic transfer frame format and common supporting services amongst the existing CCSDS link layer protocols. Currently each of the four CCSDS link layer protocols (TM, TC, AOS, and Proximity-1) have unique formats and services which prohibits their reuse across the totality of all space link applications of CCSDS member space agencies. For example, Mars missions. These missions implement their proximity data link layer using the Proximity-1 frame format and the services it supports but is still required to support the direct from Earth (TC) protocols and the Direct To Earth (AOS/TM) protocols. The prime purpose of this paper, is to describe a new general purpose CCSDS Data Link layer protocol, the NGSLP that will provide the required services along with a common transfer frame format for all the CCSDS space links (ground to/from space and space to space links) targeted for emerging missions after a CCSDS agency-wide coordinated date. This paper will also describe related options that can be included for the Coding and Synchronization sub-layer of the Data Link layer to extend the capacities of the link and additionally provide an independence of the transfer frame sub-layer from the coding sublayer. This feature will provide missions the option of running either the currently performed synchronous coding and transfer frame data link or an asynchronous coding/frame data link, in which the transfer frame length is independent of the block size of the code. The benefits from the elimination of this constraint (frame synchronized to the code block) will simplify the interface between the transponder and the data handling equipment and reduce implementation costs and complexities. The benefits include: inclusion of encoders/decoders into transmitters and receivers without regard to data link protocols, providing the ability to insert latency sensitive messages into the link to support launch, landing/docking, telerobotics. and Variable Coded Modulation (VCM). In addition the ability to transfer different sized frames can provide a backup for delivering stored anomaly engineering data simultaneously with real time data, or relaying of frames from various sources onto a trunk line for delivery to Earth.

Kazz, Greg J.↗

Studying electroweak few-body observables in chiral effective field theory

The use of nuclei to study electroweak probes is becoming increasingly relevant experimentally. The success of dark matter and neutrino experiments strongly depends on the ability to control nuclear effects in order to extract the fundamental parameters associated with external probes. Therefore, reliable theoretical calculations of nuclear structure and reactions, with well-controlled errors, are crucial for the success of experimental efforts. Currently, chiral effective field theory ($\chi$EFT) coupled with {\it ab-initio} methods represents one of the best approaches that fulfills these requirements. To use this approach as a tool for studying fundamental physics, it is essential to validate it against experimental data for which the calculations are well under control, such as the elastic scattering of electrons on nuclei. In this proceeding, I will present recent developments in the fitting of electromagnetic currents derived using $\chi$EFT and the calculation of electromagnetic form factors of light nuclei. The results of these calculations demonstrate the strength of the theory in describing the interaction of nuclei with electromagnetic probes over a broad range of momentum transfers and highlight the robustness of $\chi$EFT for analyzing future experimental data aimed at extracting fundamental parameters.

Gnech, Alex [Old Dominion Univ., Norfolk, VA (Unit↗

An LES-based airborne Doppler lidar simulator and its application to wind profiling in inhomogeneous flow conditions

Wind profiling by Doppler lidar is common practice and highly useful in a wide range of applications. Airborne Doppler lidar can provide additional insights relative to ground-based systems by allowing for spatially distributed and targeted measurements. Providing a link between theory and measurement, a first large eddy simulation (LES)-based airborne Doppler lidar simulator (ADLS) has been developed. Simulated measurements are conducted based on LES wind fields, considering the coordinate and geometric transformations applicable to real-world measurements. The ADLS provides added value as the input truth used to create the measurements is known exactly, which is nearly impossible in real-world situations. Thus, valuable insight can be gained into measurement system characteristics as well as retrieval strategies. As an example application, airborne Doppler lidar wind profiling is investigated using the ADLS. For commonly used airborne velocity azimuth display (AVAD) techniques, flow homogeneity is assumed throughout the retrieval volume, a condition which is violated in turbulent boundary layer flow. Assuming an ideal measurement system, the ADLS allows to isolate and evaluate the error in wind profiling which occurs due to the violation of the flow homogeneity assumption. Overall, the ADLS demonstrates that wind profiling is possible in turbulent wind field conditions with reasonable errors (root mean squared error of 0.36 m s –1 for wind speed when using a commonly used system setup and retrieval strategy for the conditions investigated). Nevertheless, flow inhomogeneity, e.g., due to boundary layer turbulence, can cause an important contribution to wind profiling error and is non-negligible. Results suggest that airborne Doppler lidar wind profiling at low wind speeds (<5 m s –1 ) can be biased, if conducted in regions of inhomogeneous flow conditions.

17 WIND ENERGY↗

Correction of Dual-PRF Doppler Velocity Outliers in the Presence of Aliasing

Abstract In Doppler weather radars, the presence of unfolding errors or outliers is a well-known quality issue for radial velocity fields estimated using the dual–pulse repetition frequency (PRF) technique. Postprocessing methods have been developed to correct dual-PRF outliers, but these need prior application of a dealiasing algorithm for an adequate correction. This paper presents an alternative procedure based on circular statistics that corrects dual-PRF errors in the presence of extended Nyquist aliasing. The correction potential of the proposed method is quantitatively tested by means of velocity field simulations and is exemplified in the application to real cases, including severe storm events. The comparison with two other existing correction methods indicates an improved performance in the correction of clustered outliers. The technique proposed is well suited for real-time applications requiring high-quality Doppler radar velocity fields, such as wind shear and mesocyclone detection algorithms, or assimilation in numerical weather prediction models.

58 GEOSCIENCES↗

The effects of sampling location and turbulence on discharge estimates in short converging turbine intakes

Standards provide recommendations for best practices when installing current meters to measure fluid flow in closed conduits. A central guideline requires the velocity distribution to be regular and the flow steady. Because of the nature of the short converging intakes typical of low-head hydroturbines, these assumptions may be invalid if current meters are intended to be used to estimate discharge. Usual concerns are (1) the effects of the number of devices, (2) the sampling location and (3) the high turbulence caused by blockage from submersible traveling screens usually deployed for safe downstream fish passage. These three effects were examined in the present study by using 3D simulated flow fields in both steady-state and transient modes. In the process of describing an application at an existing hydroturbine intake at Ice Harbor Dam, the present work outlined the methods involved, which combined computational fluid dynamics, laboratory measurements in physical models of the hydroturbine, and current meter performance evaluations in experimental settings. The main conclusions in this specific application were that a steady-state flow field sufficed to determine the adequate number of meters and their location, and that both the transverse velocity and turbulence intensity had a small impact on estimate errors. However, while it may not be possible to extrapolate these findings to other field conditions and measuring devices, the study laid out a path to conduct similar assessments in other applications.

computational fluid dynamics (CFD)↗

Time-resolved spectroscopy of the Mercury 6 3P1 state

The time-resolved fluorescence was observed from the Hg 6 3P1 state under the influence of the earth's magnetic field and with applied fields of up to 14 G. Modulation of the fluorescence decay signal was observed as a function of both time and space and can be interpreted in terms of a classical precession of the excited atom about the magnetic field or as quantum beats resulting from interference between coherently populated Zeeman sublevels. This modulation was studied for each of the five resolvable components of the hyperfine structure separately. The fluorescence from the even isotopes was determined to be almost completely modulated while the fluorescence from the odd isotopes was only partially modulated. The frequency of modulation of the fluorescence from the mercury-202 isotope was observed as a function of the applied magnetic field and a value for the Lande factor of 1.46 + or - 0.03 was obtained. This is within experimental error of the accepted value of 1.486. In addition, the frequency of modulation as a function of applied magnetic field was determined for each of the three resolvable components with more than one contributing isotopic hyperfine line. An investigation of the effect of radiation trapping on the degree modulation was also made.

Halstead, J. A.↗

Isobaric Reconstruction of the Baryonic Acoustic Oscillation

In this Letter, we report a significant recovery of the linear baryonic acoustic oscillation (BAO) signature by applying the isobaric reconstruction algorithm to the nonlinear matter density field. Assuming only the longitudinal component of the displacement being cosmologically relevant, this algorithm iteratively solves the coordinate transform between the Lagrangian and Eulerian frames without requiring any specific knowledge of the dynamics. For dark matter field, it produces the nonlinear displacement potential with very high fidelity. The reconstruction error at the pixel level is within a few percent and is caused only by the emergence of the transverse component after the shell-crossing. As it circumvents the strongest nonlinearity of the density evolution, the reconstructed field is well described by linear theory and immune from the bulk-flow smearing of the BAO signature. Therefore, this algorithm could significantly improve the measurement accuracy of the sound horizon scale s. For a perfect large-scale structure survey at redshift zero without Poisson or instrumental noise, the fractional error Δs/s is reduced by a factor of ∼2.7, very close to the ideal limit with the linear power spectrum and Gaussian covariance matrix.

79 ASTRONOMY AND ASTROPHYSICS↗

WE-AB-BRB-10: Filmless QA of CyberKnife MLC-Collimated and Iris-Collimated Fields

Purpose: Current methods of CK field shape QA is based on the use of radiochromic film. Though accurate results can be attained, these methods are prone to error, time consuming, and expensive. The techniques described herein perform similar QA using the FOIL Detector (Field, Output, and Image Localization). A key feature of this in-house QA solution, and central to this study, is an aSi flat-panel detector which provides the user with the means to perform accurate, immediate, and quantitative field analysis. Methods: The FOIL detector is automatically aligned in the CK beam using fiducial markers implanted within the detector case. Once the system is aligned, a treatment plan is delivered which irradiates the flat-panel imager using the field being tested. The current study tests each of the clinically-used fields shaped using the Iris variable-aperture collimation system using a plan which takes 6 minutes to deliver. The user is immediately provided with field diameter and beam profile, as well as a comparison to baseline values. Additionally, the detector is used to acquire and analyze leaf positions of the InCise multi-leaf collimation system. Results: Using a 6-minute plan consisting of 11 beams of 25MU-per-beam, the FOIL detector provided the user with a quantitative analysis of all clinically-used field shapes. The FOIL detector was also able to clearly resolve field edge junctions in a picket fence test, including slight over-travel of individual leaves as well as inter-leaf leakage. Conclusion: The FOIL system provided comparable field diameter and profile data when compared to methods using film; providing results much faster and with 5% of the MU used for film. When used with the MLC system, the FOIL detector provided the means for immediate quantification of the performance of the system through analysis of leaf positions in a picket fence test field. Author is the President/Owner of Spectrum Medical Physics, LLC, a company which maintains contracts with Siemens Healthcare and Standard Imaging, Inc.

60 APPLIED LIFE SCIENCES↗

Discretization Error Estimation and Control for Farfield Acoustic Signatures

We investigate the utility of adjoint-based error estimates for sonic boom farfield simulations governed by solutions of the augmented Burgers’ equation. Solution of this nonlinear system uses operator splitting with a second-order finite volume discretization in space and second-order Runge-Kutta time marching, while the absorption and molecular relaxation are solved using second-order central differencing. The discretization error in selected ground sonic boom cost functionals is estimated using the method of adjoint-weighted residuals. Key elements of the implementation process are emphasized with details provided on the practical aspects as appliedto the sonic boom farfield propagation. We establish the accuracy of the adjoint solutions usingcomplex step and finite difference approaches, and examine the accuracy of the error estimates using analytical N-wave solutions. We then apply it to a pressure waveform corresponding to the X-59 research aircraft. The investigations demonstrate that the method of adjoint-weighted residuals accurately predicts the level of discretization error present in sonic boom farfield simulations while offering insight into which features of the near field signal are the primary drivers of ground noise metrics. The numerical results indicate that at sampling frequencies as low as50kHz, discretization error in the propagation is under 0.01 dB[A] for realistically complex examples.

CST↗