Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “error indicators”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Technical Development for S-CO 2 Advanced Energy Conversion

This report is divided into four parts. First part of the report describes the methods used to measure and model the flow of supercritical carbon dioxide (S-CO 2 ) through annuli and straight-through labyrinth seals. The effects of shaft eccentricity in small diameter annuli were observed for length-to-hydraulic diameter (L/D) ratios of 6, 12, 143, and 235. Flow rates through tooth-cavity labyrinth seals were measured for inlet pressures of 7.7, 10, and 11 MPa with corresponding inlet densities of 325, 475, and 630 kg/m 3 . Various leakage models were compared to this result to describe their applicability in supercritical carbon dioxide applications. Flow rate measurements were made varying tooth number for labyrinth seals of same total length. Second part of the report describes the computational study performed to understand the leakage through the labyrinth seals using Open source CFD package OpenFOAM. Fluid Property Interpolation Tables (FIT) program was implemented in OpenFOAM to accurately model the properties of CO2 required to solve the governing equations. To predict the flow behavior in the two phase dome Homogeneous Equilibrium Model (HEM) is assumed to be valid. Experimental results for plain orifice (L/D ~ 5) were used to show the capabilities of the FIT model implemented in OpenFOAM. Error analysis indicated that OpenFOAM is capable of predicting experimental data within ±10% error with the majority of data close to ±5% error. Following the validation of computational model, effects of geometrical parameters and operating conditions are isolated from each other and a parametric study was performed in two parts to understand their effects on leakage flow. Third part of the report provides the details of the constructed heat exchanger test facility and presents the experimental results obtained to investigate the effects of buoyancy on heat transfer characteristics of Supercritical carbon dioxide in heating mode. Turbulent flows with Reynolds numbers up to 60,000, at operating pressures of 7.5, 8.1, and 10.2 MPa were tested in a round tube. Local heat transfer coefficients were obtained from measured wall temperatures over a large set of experimental parameters that varied inlet temperature from 20 °C to 55 °C,mass flux from 150 to 350 kg/m 2 s, and a maximum heat flux of 65 KW/m 2 . Horizontal, upward and downward flows were tested to investigate the unusual heat-transfer characteristics to the effect of buoyancy and flow acceleration caused by large variation in density. Final part of this report presents the simplified analysis performed to investigate the possibility of using wet cooling tower option to reject heat from the supercritical carbon dioxide Brayton cycle power convertor for AFR-100 and ABR-1000 plants. A code was developed to estimate the tower dimensions, power and water consumption, and to perform economic analysis. The code developed was verified by comparing the calculations to a vendor quote. The effect of ambient air and water conditions on the sizing and construction of the cooling tower as well as the cooler is studied. Finally, a cost-based optimization technique is used to estimate the optimum water conditions which will improve the plant economics.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Patient motion effects on the quantification of regional myocardial blood flow with dynamic PET imaging

Purpose: Patient motion is a common problem during dynamic positron emission tomography (PET) scans for quantification of myocardial blood flow (MBF). The purpose of this study was to quantify the prevalence of body motion in a clinical setting and evaluate with realistic phantoms the effects of motion on blood flow quantification, including CT attenuation correction (CTAC) artifacts that result from PET–CT misalignment. Methods: A cohort of 236 sequential patients was analyzed for patient motion under resting and peak stress conditions by two independent observers. The presence of motion, affected time-frames, and direction of motion was recorded; discrepancy between observers was resolved by consensus review. Based on these results, patient body motion effects on MBF quantification were characterized using the digital NURBS-based cardiac-torso phantom, with characteristic time activity curves (TACs) assigned to the heart wall (myocardium) and blood regions. Simulated projection data were corrected for attenuation and reconstructed using filtered back-projection. All simulations were performed without noise added, and a single CT image was used for attenuation correction and aligned to the early- or late-frame PET images. Results: In the patient cohort, mild motion of 0.5 ± 0.1 cm occurred in 24% and moderate motion of 1.0 ± 0.3 cm occurred in 38% of patients. Motion in the superior/inferior direction accounted for 45% of all detected motion, with 30% in the superior direction. Anterior/posterior motion was predominant (29%) in the posterior direction. Left/right motion occurred in 24% of cases, with similar proportions in the left and right directions. Computer simulation studies indicated that errors in MBF can approach 500% for scans with severe patient motion (up to 2 cm). The largest errors occurred when the heart wall was shifted left toward the adjacent lung region, resulting in a severe undercorrection for attenuation of the heart wall. Simulations also indicated that the magnitude of MBF errors resulting from motion in the superior/inferior and anterior/posterior directions was similar (up to 250%). Body motion effects were more detrimental for higher resolution PET imaging (2 vs 10 mm full-width at half-maximum), and for motion occurring during the mid-to-late time-frames. Motion correction of the reconstructed dynamic image series resulted in significant reduction in MBF errors, but did not account for the residual PET–CTAC misalignment artifacts. MBF bias was reduced further using global partial-volume correction, and using dynamic alignment of the PET projection data to the CT scan for accurate attenuation correction during image reconstruction. Conclusions: Patient body motion can produce MBF estimation errors up to 500%. To reduce these errors, new motion correction algorithms must be effective in identifying motion in the left/right direction, and in the mid-to-late time-frames, since these conditions produce the largest errors in MBF, particularly for high resolution PET imaging. Ideally, motion correction should be done before or during image reconstruction to eliminate PET-CTAC misalignment artifacts.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Error statistics and scalability of quantum error mitigation formulas

Quantum computing promises advantages over classical computing in many problems. Nevertheless, noise in quantum devices prevents most quantum algorithms from achieving the quantum advantage. Quantum error mitigation provides a variety of protocols to handle such noise using minimal qubit resources. While some of those protocols have been implemented in experiments for a few qubits, it remains unclear whether error mitigation will be effective in quantum circuits with tens to hundreds of qubits. In this paper, we apply statistics principles to quantum error mitigation and analyse the scaling behaviour of its intrinsic error. We find that the error increases linearly O(ϵN) with the gate number N before mitigation and sublinearly O(ϵ'N γ ) after mitigation, where γ ≈ 0.5, ϵ is the error rate of a quantum gate, and ϵ' is a protocol-dependent factor. The $\sqrt{N}$ scaling is a consequence of the law of large numbers, and it indicates that error mitigation can suppress the error by a larger factor in larger circuits. We propose the importance Clifford sampling as a key technique for error mitigation in large circuits to obtain this result.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Chemical composition measurements of the low activity waste (LAW) EPA-Series glasses

In this report, the Savannah River National Laboratory provides chemical analysis results for a series of simulated low activity waste glasses provided by Pacific Northwest National Laboratory as part of an ongoing development task. The measured chemical composition data are reported and compared with the targeted values for each component for each glass. A detailed review showed no indications of errors in the preparation or measurement of the study glasses. All of the measured sums of oxides for the study glasses fell within the interval of 100.2 to 100.8 wt %, indicating recovery of all components. Comparisons of the targeted and measured chemical compositions showed that the measured values for the glasses met the targeted concentrations within 10% for those components present at more than 5 wt %.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2016

Chemical Composition Measurements of LAWA44 Glass Samples

DOE is building the Hanford Tank Waste Treatment and Immobilization Plant (WTP) at the Hanford Site in Washington to remediate 55 million gallons of radioactive waste that is temporarily stored in 177 underground tanks. Both low-activity and high-level wastes will then be vitrified into borosilicate glass using Joule-heated ceramic melters. Efforts are being made to increase the loading of Hanford tank wastes in the glass. One area of work is enhancing waste glass composition/property models and broadening the compositional regions over which those models are applicable. In this report, the Savannah River National Laboratory provides chemical analysis results for several samples of a simulated low-activity waste glass, LAWA44, provided by the Pacific Northwest National Laboratory as part of an ongoing development task. The measured chemical composition data are reported and compared with the targeted values for each component for each glass. A detailed review showed no indications of errors in the preparation or measurement of the study glasses. All of the measured sums of oxides for the study glasses fell within the interval of 97.9 to 102.6 wt %, indicating acceptable recovery of the glass components. Comparisons of the targeted and measured chemical compositions showed that the measured values for the glasses met the targeted concentrations within 10% for those components present at more than 5 wt %. It was noted that the measured B 2 O 3 concentrations are somewhat above the targeted values for the study glasses. No obvious trends were observed with regard to the multiple melting steps used to prepare the study glasses, indicating that any potential effects of volatility were below measurable thresholds.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2014

Estimating uncertainties in statistics computed from direct numerical simulation

Rigorous assessment of uncertainty is crucial to the utility of direct numerical simulation (DNS) results. Uncertainties in the computed statistics arise from two sources: finite statistical sampling and the discretization of the Navier–Stokes equations. Due to the presence of non-trivial sampling error, standard techniques for estimating discretization error (such as Richardson extrapolation) fail or are unreliable. This work provides a systematic and unified approach for estimating these errors. First, a sampling error estimator that accounts for correlation in the input data is developed. Then, this sampling error estimate is used as part of a Bayesian extension of Richardson extrapolation in order to characterize the discretization error. These methods are tested using the Lorenz equations and are shown to perform well. These techniques are then used to investigate the sampling and discretization errors in the DNS of a wall-bounded turbulent flow at Reτ ≈ 180. Both small (Lx/δ × Lz/δ = 4π × 2π) and large (Lx/δ × Lz/δ = 12π × 4π) domain sizes are investigated. For each case, a sequence of meshes was generated by first designing a “nominal” mesh using standard heuristics for wall-bounded simulations. These nominal meshes were then coarsened to generate a sequence of grid resolutions appropriate for the Bayesian Richardson extrapolation method. In addition, the small box case is computationally inexpensive enough to allow simulation on a finer mesh, enabling the results of the extrapolation to be validated in a weak sense. For both cases, it is found that while the sampling uncertainty is large enough to make the order of accuracy difficult to determine, the estimated discretization errors are quite small. This indicates that the commonly used heuristics provide adequate resolution for this class of problems. However, it is also found that, for some quantities, the discretization error is not small relative to sampling error, indicating that the conventional wisdom that sampling error dominates discretization error for this class of simulations needs to be reevaluated.

Oliver, Todd A.↗
DOE OSTI2021

Precision measurement of relative γ-ray intensities from the decay of 61Cu

A discrepancy, well outside reported uncertainties, has been observed between the accepted and measured values of the intensity ratio of the two strongest γ-ray following 61 Cu β + decay. This discrepancy has significant impact since the nat Ni(d,x) 61 Cu reaction has historically been one of only a few IAEA recommendations for use as a deuteron flux monitor and a considerable number of published cross sections measured in ratio to that beam monitor cross section may depend on the choice of either the first or second strongest ray in those calculations. To determine the magnitude of this error most precisely, over a hundred separate measurements of the 283 keV to 656 keV γ-ray emission ratio were collected from seven experiments and a variety of detectors and detection geometries. A weighted average of all these measurements indicates an error in the value listed in the Nuclear Data Sheets of 11% in either the primary or second-highest intensity γ-ray of 61 Cu, potentially introducing an 11% error in 61 Cu production cross section measurements, cross sections using nickel activation as a deuteron beam current monitor, or in dose rates when 61 Cu is used in nuclear medicine. General agreement with the Data Sheets with ten other intensity ratios suggests the most probable error is in the secondary (656 keV) emission, which accordingly should be updated from 10.8% to 9.69%.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Error estimates of finite difference methods for the Dirac equation in the massless and nonrelativistic regime

We present four frequently used finite difference methods and establish the error bounds for the discretization of the Dirac equation in the massless and nonrelativistic regime, involving a small dimensionless parameter 0 < ε &NestedLessLess; 1 inversely proportional to the speed of light. In the massless and nonrelativistic regime, the solution exhibits rapid motion in space and is highly oscillatory in time. Specifically, the wavelength of the propagating waves in time is at O(ε), while in space, it is at O(1) with the wave speed at O(ε -1 ). We adopt one leap-frog, two semi-implicit, and one conservative Crank-Nicolson finite difference methods to numerically discretize the Dirac equation in one dimension and establish rigorously the error estimates which depend explicitly on the time step τ, mesh size h, and the small parameter ε. The error bounds indicate that, to obtain the “correct” numerical solution in the massless and nonrelativistic regime, i.e., 0 < ε &NestedLessLess; 1, all these finite difference methods share the same ε-scalability as time step τ = O(ε 3/2 ) and mesh size h = O(ε 1/2 ). A large number of numerical results are reported to verify the error estimates.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Isotope identity experiments in JET-ILW with H and D L-mode plasmas

NBI-heated L-mode plasmas have been obtained in JET with the Be/W ITER-like wall (JET-ILW) in H and D, with matched profiles of the dimensionless plasma parameters, ρ*, ν*, β and q in the plasma core confinement region and same T i /T e and Z eff . The achieved isotope identity indicates that the confinement scale invariance principle is satisfied in the core confinement region of these plasmas, where the dominant instabilities are Ion Temperature Gradient (ITG) modes. The dimensionless thermal energy confinement time, Ω i τ E ,t h , and the scaled core plasma heat diffusivity, A χ eff /B T , are identical in H and D within error bars, indicating lack of isotope mass dependence of the dimensionless L-mode thermal energy confinement time in JET-ILW. Predictive flux driven simulations with JETTO-TGLF of the H and D identity pair is in very good agreement with experiment for both isotopes: the stiff core heat transport, typical of JET-ILW NBI heated L-modes, overcomes the local gyro-Bohm scaling of gradient-driven TGLF, explaining the lack of isotope mass dependence in the confinement region of these plasmas. The effect of E × B shearing on the predicted heat and particle transport channels is found to be negligible for these low beta and low momentum input plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

A bronchoscopic navigation system using bronchoscope center calibration for accurate registration of electromagnetic tracker and CT volume without markers

Purpose: Various bronchoscopic navigation systems are developed for diagnosis, staging, and treatment of lung and bronchus cancers. To construct electromagnetically navigated bronchoscopy systems, registration of preoperative images and an electromagnetic tracker must be performed. This paper proposes a new marker-free registration method, which uses the centerlines of the bronchial tree and the center of a bronchoscope tip where an electromagnetic sensor is attached, to align preoperative images and electromagnetic tracker systems. Methods: The chest computed tomography (CT) volume (preoperative images) was segmented to extract the bronchial centerlines. An electromagnetic sensor was fixed at the bronchoscope tip surface. A model was designed and printed using a 3D printer to calibrate the relationship between the fixed sensor and the bronchoscope tip center. For each sensor measurement that includes sensor position and orientation information, its corresponding bronchoscope tip center position was calculated. By minimizing the distance between each bronchoscope tip center position and the bronchial centerlines, the spatial alignment of the electromagnetic tracker system and the CT volume was determined. After obtaining the spatial alignment, an electromagnetic navigation bronchoscopy system was established to real-timely track or locate a bronchoscope inside the bronchial tree during bronchoscopic examinations. Results: The electromagnetic navigation bronchoscopy system was validated on a dynamic bronchial phantom that can simulate respiratory motion with a breath rate range of 0–10 min{sup −1}. The fiducial and target registration errors of this navigation system were evaluated. The average fiducial registration error was reduced from 8.7 to 6.6 mm. The average target registration error, which indicates all tracked or navigated bronchoscope position accuracy, was much reduced from 6.8 to 4.5 mm compared to previous registration methods. Conclusions: An electromagnetically navigated bronchoscopy system was constructed with accurate registration of an electromagnetic tracker and the CT volume on the basis of an improved marker-free registration approach that uses the bronchial centerlines and bronchoscope tip center information. The fiducial and target registration errors of our electromagnetic navigation system were about 6.6 and 4.5 mm in dynamic bronchial phantom validation.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Nondestructive Damage Detection of Concrete With Alkali-Silica Reactions Using Coda Wave and Anomaly Detection

An anomaly detection model for early damage detection for concrete structures undergoing alkali-silica reaction (ASR) is presented. It is difficult to detect ASR initiation and early damage without a reference expansion measurement. Coda waves, or the multiply scattered portion of ultrasonic waves, have been found to be indicative of small changes in complex material such as concrete. The relationship between concrete damage and relative velocity change and decorrelation of coda waves has been studied, but a generalized model which detects when damage occurs in a concrete structure is still lacking. The presented method uses features extracted from coda waves to detect early damage in concrete structures. The model uses unsupervised learning and only requires data from undamaged structures for training. During the training process, the reconstruction error of the training data is minimized. When the data collected from damaged concrete structures is used as an input of the model, it returns high reconstruction errors that indicate the occurrence of damage in the structures. The performance of the model is validated using experimental studies and has been shown to generalize across two different ASR specimens.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Evaluation of the Planetary Boundary Layer Height From ERA5 Reanalysis With MOSAiC Observations Over the Arctic Ocean

The planetary boundary layer height (PBLH) is a crucial indicator reflecting the region of the atmosphere characterized by continuous turbulence. Here, we use radiosonde and surface meteorological observations (4–7 times per day, year-round measurements) during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition to derive the PBLH (PBLH MOSAiC ), and further evaluate the PBLH from the ERA5 reanalysis (PBLH ERA5 ). Comparisons between PBLH MOSAiC and PBLH ERA5 from different perspectives reveal that: (a) The overestimation of PBLH ERA5 when the sea ice concentration is >90% is significant with the centered root mean squared error reaching up to 201 m; (b) The difference between the two products is notably pronounced in cold seasons, while it is comparatively diminished in warm seasons; (c) In neutral boundary layers, differences in PBLH ERA5 are larger compared with stable and convective boundary layers. In addition, the analysis of error sources indicates that the bias of PBLH ERA5 is sensitive to the bias of vertical thermal structure and wind speed profiles in ERA5 data sets in all conditions. Finally, we find a Random Forest model effectively reduces the bias of PBLH ERA5 with the index of agreement reaching up to 0.71 in the test data set, while a multiple linear regression demonstrates comparable performance to the Random Forest model.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

MO-D-213-05: Sensitivity of Routine IMRT QA Metrics to Couch and Collimator Rotations

Purpose: To assess the sensitivity of gamma index and other IMRT QA metrics to couch and collimator rotations. Methods: Two brain IMRT plans with couch and/or collimator rotations in one or more of the fields were evaluated using the IBA MatriXX ion chamber array and its associated software (OmniPro-I’mRT). The plans were subjected to routine QA by 1) Creating a composite planar dose in the treatment planning system (TPS) with the couch/collimator rotations and 2) Creating the planar dose after “zeroing” the rotations. Plan deliveries to MatriXX were performed with all rotations set to zero on a Varian 21ex linear accelerator. This in effect created TPS-created planar doses with an induced rotation error. Point dose measurements for the delivered plans were also performed in a solid water phantom. Results: The IMRT QA of the plans with couch and collimator rotations showed clear discrepancies in the planar dose and 2D dose profile overlays. The gamma analysis, however, did pass with the criteria of 3%/3mm (for 95% of the points), albeit with a lower percentage pass rate, when one or two of the fields had a rotation. Similar results were obtained with tighter criteria of 2%/2mm. Other QA metrics such as percentage difference or distance-to-agreement (DTA) histograms produced similar results. The point dose measurements did not obviously indicate the error due to location of dose measurement (on the central axis) and the size of the ion chamber used (0.6 cc). Conclusion: Relying on Gamma analysis, percentage difference, or DTA to determine the passing of an IMRT QA may miss critical errors in the plan delivery due to couch/collimator rotations. A combination of analyses for composite QA plans, or per-beam analysis, would detect these errors.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Correlated $\ n-γ$ angular distributions from the $\ Q$ = 4.4398 MeV 12 C ($\ n, n' γ$) reaction for incident neutron energies from 6.5 MeV to 16.5 MeV

Neutron scattering cross sections and angular distributions represent one of the most glaring sources of uncertainty in calculations of nuclear systems. Even simple nuclei like 12 C show indications of errors in nuclear databases for scattering reactions. Measurements of inelastic neutron scattering have historically measured either the scattered neutrons or the nuclear deexcitation $\ γ$ emission. Only a very small number of experiments attempted correlated measurements of both the neutron and $\ γ$ data simultaneously, even though these $\ n-γ$ correlations could be essential for understanding particle transport in nuclear systems. In this work we describe a measurement of the $\ n, γ$, and correlated $\ n-γ$ angular distributions from the $\ Q$ = 4.4398 MeV 12 C ($\ n, n'γ$) reaction in a single experiment using an EJ-309 liquid scintillator detector array with wide angular coverage, and with a continuous incident neutron energy range from 6.5 to 16.5 MeV. We also provide a thorough covariance description of these results, including normalization of the probability distribution. While the measured n distributions agree well with the relatively large number of available literature measurements, there are comparatively very few measurements of the γ distributions from this reaction. However, our data support the presence of a nonzero α 4 Legendre polynomial component of the γ angular distribution suggested in past measurements, which is currently not incorporated in the ENDF/B-VIII.0 library despite the use of these same literature data for evaluation of the 12 C ($\ n, n'γ$) cross section. The correlated $\ n-γ$ distribution measurements are limited to three measurements at incident neutron energies near 14 MeV. Our results do not generally agree with any of these literature measurements. We observe clear indications of significant changes in the $\ n$ distribution for specific $\ γ$-detection angles and vice versa especially near thresholds for other reaction channels, which shows the potential for significant bias in experiments that, for example, tag on inelastic scattering using a single or small number of $\ γ$ -detection angles and could impact particle transport calculations.

6 ≤ A ≤ 19↗
DOE OSTI2020

Calibration, measurement, and characterization of soil moisture dynamics in a central Amazonian tropical forest

Soil moisture plays a key role in hydrological, biogeochemical, and energy budgets of terrestrial ecosystems. Accurate soil moisture measurements in remote ecosystems such as the Amazon are difficult and limited because of logistical constraints. Time domain reflectometry (TDR) sensors are widely used to monitor soil moisture and require calibration to convert the TDR's dielectric permittivity measurement (K a ) to volumetric water content (θ v ). In this study, our objectives were to develop a field-based calibration of TDR sensors in an old-growth upland forest in the central Amazon, to evaluate the performance of the calibration, and then to apply the calibration to determine the dynamics of soil moisture content within a 14.2 -m-deep vertical soil profile. Depth-specific TDR calibration using local soils in a controlled laboratory setting yielded a novel K a –θ v third-degree polynomial calibration. The sensors were later installed to their specific calibration depth in a 14.2-m pit. The widely used K a –θ v relationship (Topp model) underestimated the site-specific θ v by 22–42%, indicating significant error in the model when applied to these well-structured, clay-rich tropical forest soils. The calibrated wet- and dry-season θ v data showed a variety of depth and temporal variations highlighting the importance of soil textural differentiation, root uptake depths, as well as event to seasonal precipitation effects. Data such as these are greatly needed for improving our understanding of ecohydrological processes within tropical forests and for improving models of these systems in the face of changing environmental conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

MFNets: data efficient all-at-once learning of multifidelity surrogates as directed networks of information sources

We present an approach for constructing a surrogate from ensembles of information sources of varying cost and accuracy. The multifidelity surrogate encodes connections between information sources as a directed acyclic graph, and is trained via gradient-based minimization of a nonlinear least squares objective. While the vast majority of state-of-the-art assumes hierarchical connections between information sources, our approach works with flexibly structured information sources that may not admit a strict hierarchy. The formulation has two advantages: (1) increased data efficiency due to parsimonious multifidelity networks that can be tailored to the application; and (2) no constraints on the training data—we can combine noisy, non-nested evaluations of the information sources. Finally, numerical examples ranging from synthetic to physics-based computational mechanics simulations indicate the error in our approach can be orders-of-magnitude smaller, particularly in the low-data regime, than single-fidelity and hierarchical multifidelity approaches.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Adaptive mesh refinement for steady flows in Nek5000

Adaptive mesh refinement is performed in the framework of the spectral element method augmented by approaches to error estimation and control. The h-refinement technique is used for adapting the mesh, where selected grid elements are split by a quadtree (2D) or octree (3D) structure. Continuity between parent-child elements is enforced by high-order interpolation of the solution across the common faces. Parallel mesh partitioning and grid management respectively, are taken care of by the external libraries ParMETIS and p4est. Two methods are considered for estimating and controlling the error of the solution. The first error estimate is local and based on the spectral properties of the solution on each element. This method gives a local measure of the L-2-norm of the solution over the entire computational domain. The second error estimate uses the dual-weighted residuals method - it is based on and takes into account both the local properties of the solution and the global dependence of the error in the solution via an adjoint problem. The objective of this second approach is to optimize the computation of a given functional of physical interest. The simulations are performed by using the code Nek5000 and three steady-state test cases are studied: a two-dimensional lid-driven cavity at Re = 7, 500, a two-dimensional flow past a cylinder at Re = 40, and a three-dimensional lid-driven cavity at Re = 2,000 with a moving lid tilted by an angle of 30 degrees. Finally, the efficiency of both error estimators is compared in terms of refinement patterns and accuracy on the functional of interest. In the case of the adjoint error estimators, the trend on the error of the functional is shown to be correctly represented up to a multiplicative constant.

42 ENGINEERING↗
DOE OSTI2023

Residual estimation for grid modification in wall-modeled large eddy simulation using unstructured high-order methods

Here, the accuracy and computational cost of a large eddy simulation are highly dependent on the computational grid. Building optimal grids manually from a priori knowledge is not feasible in most practical use cases; instead, solution-adaptive strategies can provide a robust and cost-efficient method to generate a grid with the desired accuracy. We adapt the residual estimation algorithm developed by Toosi and Larsson for Discontinuous Galerkin Spectral Elements Methods (DGSEM) to guide the grid-adaptation process. The core of the method is the computation of the estimated modeling residual using the polynomial basis functions used in DGSEM, and the averaging of the estimated residual over each element. The final method is assessed in multiple channel flow test cases and for the transonic flow over an airfoil, in both cases making use of mortar interfaces between elements with hanging nodes. The method is found to be robust and reliable, and to provide solutions on grids with significantly fewer elements at comparable accuracy compared to when using human-generated grids.

97 MATHEMATICS AND COMPUTING↗