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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
NASA NTRS1997

The Climate Signal in Regional Moisture Fluxes: A Comparison of Three Global Data Assimilation Products

This study assesses the quality of estimates of climate variability in moisture flux and convergence from three assimilated data sets: two are reanalysis products generated at the Goddard Data Assimilation Office (DAO) and the National Centers for Environmental Prediction/National Centers for Atmospheric Research (NCEPJNCAR), and the third consists of the operational analyses generated at the European Center for Medium Range Forecasts (ECMWF). The regions under study (the United States Great Plains, the Indian monsoon region, and Argentina east of the Andes) are characterized by frequent low level jets (LLJs) and other interannual low level wind variations tied to the large-scale flow. While the emphasis is on the reanalysis products, the comparison with the operational product is provided to help assess the improvements gained from a fixed analysis system. All three analyses capture the main moisture flux anomalies associated with selected extreme climate (drought and flood) events during the period 1985-93. The correspondence is strongest over the Great Plains and weakest over the Indian monsoon region reflecting differences in the observational coverage. For the reanalysis products, the uncertainties in the lower tropospheric winds is by far the dominant source of the discrepancies in the moisture flux anomalies in the middle latitude regions. Only in the Indian Monsoon region, where interannual variability in the low level winds is comparatively small, does the moisture bias play a substantial role. In contrast, the comparisons with the operational product show differences in moisture which are comparable torhe differences in the wind in all three regions. Compared with the fluxes, the anomalous moisture convergences show substantially larger differences among the three products. The best agreement occurs over the Great Plains region where all three products show vertically-integrated moisture convergence during the floods and divergence during the drought with differences in magnitude of about 25%. The reanalysis products, in particular, show good agreement in depicting the different roles of the mean flow and transients during the flood and drought periods. Differences between the three products in the other two regions exceed 100% reflecting differences in the low level jets and the large scale circulation patterns. The operational product tends to have locally larger amplitude convergence fields which average out in area-mean budgets: this appears to be at least in part due to errors in the surface pressure fields and aliasing from the higher resolution of the original ECMWF fields. On average, the reanalysis products show higher coherence with each other than with the operational product in the estimates of interannual variability. This result is less clear in the Indian monsoon region where differences in the input observations appears to be an important factor. The agreement in the anomalous convergence patterns is, however, still rather poor even over relatively data dense regions such as the United States Great Plains. These differences are attributed to deficiencies in the assimilating GCM's representations of the planetary boundary layer and orography, and a global observing system incapable of resolving the highly confined low level winds associated with the climate anomalies.

Min, Wei↗
DOE OSTI2015

SU-C-BRD-06: Sensitivity Study of An Automated System to Acquire and Analyze EPID Exit Dose Images

Purpose: The dosimetric consequences of errors in patient setup or beam delivery and anatomical changes are not readily known. A new product, PerFRACTION (Sun Nuclear Corporation), is designed to detect these errors by comparing the EPID exit dose image from each field of each fraction to those from baseline fraction images. This work investigates the sensitivity of PerFRACTION to detect the deviation of induced errors in a variety of realistic scenarios. Methods: Eight plans were created mimicking potential delivery or setup errors. The plans consisted of a nominal field and the field with an induced error. These were used to irradiate the EPID simulating multiple fractions with and without the error. Integrated EPID images were acquired in clinical mode and saved in ARIA. PerFRACTION automatically pulls the images into its database and performs the user defined comparison. In some cases, images were manually pushed to PerFRACTION. We varied the distance-to-agreement or dose tolerance until PerFRACTION showed failing pixels in the affected region and recorded the values. We induced errors of 1mm and greater in jaw, MLC, and couch position, 2 degree collimation rotation (patient yaw), and 0.5% to 1.5% in machine output. Both static and arc fields with the rails in or out were also acquired and compared. Results: PerFRACTION detected position errors of the jaws, MLC, and couch with an accuracy of better than 0.5 mm, and 0.2 degrees for collimator and gantry error. PerFRACTION detected a machine output error within 0.2% and detected the change in rail position. Conclusion: A new automated system for monitoring daily treatments for machine or patient variations from the first fraction using integrated EPID images was found to be sensitive enough to detect small positional, angular, and dosimetric errors within 0.5mm, 0.2 degrees, and 0.2%, respectively. Sun Nuclear Corporation has provided a software license for the product described.

60 APPLIED LIFE SCIENCES↗
NASA NTRS2006

Terrestrial Water Mass Load Changes from Gravity Recovery and Climate Experiment (GRACE)

Recent studies show that data from the Gravity Recovery and Climate Experiment (GRACE) is promising for basin- to global-scale water cycle research. This study provides varied assessments of errors associated with GRACE water storage estimates. Thirteen monthly GRACE gravity solutions from August 2002 to December 2004 are examined, along with synthesized GRACE gravity fields for the same period that incorporate simulated errors. The synthetic GRACE fields are calculated using numerical climate models and GRACE internal error estimates. We consider the influence of measurement noise, spatial leakage error, and atmospheric and ocean dealiasing (AOD) model error as the major contributors to the error budget. Leakage error arises from the limited range of GRACE spherical harmonics not corrupted by noise. AOD model error is due to imperfect correction for atmosphere and ocean mass redistribution applied during GRACE processing. Four methods of forming water storage estimates from GRACE spherical harmonics (four different basin filters) are applied to both GRACE and synthetic data. Two basin filters use Gaussian smoothing, and the other two are dynamic basin filters which use knowledge of geographical locations where water storage variations are expected. Global maps of measurement noise, leakage error, and AOD model errors are estimated for each basin filter. Dynamic basin filters yield the smallest errors and highest signal-to-noise ratio. Within 12 selected basins, GRACE and synthetic data show similar amplitudes of water storage change. Using 53 river basins, covering most of Earth's land surface excluding Antarctica and Greenland, we document how error changes with basin size, latitude, and shape. Leakage error is most affected by basin size and latitude, and AOD model error is most dependent on basin latitude.

Seo, K.-W.↗
DOE OSTI2021

A multiresolution adaptive wavelet method for nonlinear partial differential equations

We report the multiscale complexity of modern problems in computational science and engineering can prohibit the use of traditional numerical methods in multi-dimensional simulations. Therefore, novel algorithms are required in these situations to solve partial differential equations (PDEs) with features evolving on a wide range of spatial and temporal scales. To meet these challenges, we present a multiresolution wavelet algorithm to solve PDEs with significant data compression and explicit error control. We discretize in space by projecting fields and spatial derivative operators onto wavelet basis functions. We provide error estimates for the wavelet representation of fields and their derivatives. Then, our estimates are used to construct a sparse multiresolution discretization which guarantees the prescribed accuracy. Additionally, we embed a predictor-corrector procedure within the temporal integration to dynamically adapt the computational grid and maintain the accuracy of the solution of the PDE as it evolves. We present examples to highlight the accuracy and adaptivity of our approach.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

A new field solver for modeling of relativistic particle-laser interactions using the particle-in-cell algorithm

A customized finite-difference field solver for the particle-in-cell (PIC) algorithm that provides higher fidelity for wave-particle interactions in intense electromagnetic waves is presented. In many problems of interest, particles with relativistic energies interact with intense electromagnetic fields that have phase velocities near the speed of light. Numerical errors can arise due to (1) dispersion errors in the phase velocity of the wave, (2) the staggering in time between the electric and magnetic fields and between particle velocity and position and (3) errors in the time derivative in the momentum advance. Errors of the first two kinds are analyzed in detail. It is shown that by using field solvers with different -space operators in Faraday’s and Ampere’s law, the dispersion errors and magnetic field time-staggering errors in the particle pusher can be simultaneously removed for electromagnetic waves moving primarily in a specific direction. Here, the new algorithm was implemented into Osiris by using customized higher-order finite-difference operators. Schemes using the proposed solver in combination with different particle pushers are compared through PIC simulation. It is shown that the use of the new algorithm, together with an analytic particle pusher (assuming constant fields over a time step), can lead to accurate modeling of the motion of a single electron in an intense laser field with normalized vector potentials, eA / mc 2 , exceeding for typical cell sizes and time steps.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Propagation of Noise Uncertainty Through Virtual Strain Gauge Formulations for 2D Digital Image Correlation

The effect of displacement uncertainty is examined on 2-dimensional strain, calculated using linear surfaces fitted to the displacement field. A classical engineering error propagation method is used to calculate uncertainty in Green-Lagrangian strain calculations. The derived uncertainty is compared to a Monte Carlo simulation and discrepancies under 2% are seen between these two methods. The effect of virtual strain gauge size, displacement uncertainty, and boundaries on the region of interest on the strain uncertainty are considered. Here, an exponential decay relationship is observed between strain uncertainty and virtual strain gauge size, while a linear relationship is seen between strain and displacement uncertainty. For boundaries in the region of interest, strain uncertainty is affected by the reduced number of points available to perform the regression.

42 ENGINEERING↗
DOE OSTI2014

Asymptotic safety on the lattice: The nonlinear O(N) sigma model

We study the non-perturbative renormalization group flow of the nonlinear O(N) sigma model in two and three spacetime dimensions using a scheme that combines an effective local Hybrid Monte Carlo update routine, blockspin transformations and a Monte Carlo demon method. In two dimensions our results verify perturbative renormalizability. In three dimensions, we determine the flow diagram of the theory for various N and different truncations and find a non-trivial fixed point, which indicates non-perturbative renormalizability. It is related to the well-studied phase transition of the O(N) universality class and characterizes the continuum physics of the model. We compare the obtained renormalization group flows with recent investigations by means of the Functional Renormalization Group. - Highlights: • We compute the global flow diagram from lattice simulations. • Truncation of the effective action may lead to artificial fixed points. • Our optimization scheme effectively reduces truncation errors. • The Wilson–Fisher fixed point is identified as a non-trivial UV fixed point. • Asymptotic safety is realized for the 3D nonlinear sigma model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Computation of stationary 3D halo currents in fusion devices with accuracy control

This paper addresses the calculation of the resistive distribution of halo currents in three-dimensional structures of large magnetic confinement fusion machines. A Neumann electrokinetic problem is solved on a geometry so complicated that complementarity is used to monitor the discretization error. An irrotational electric field is obtained by a geometric formulation based on the electric scalar potential, whereas three geometric formulations are compared to obtain a solenoidal current density: a formulation based on the electric vector potential and two geometric formulations inspired from mixed and mixed-hybrid Finite Elements. The electric vector potential formulation is usually considered impractical since an enormous computing power is wasted by the topological pre-processing it requires. To solve this challenging problem, we present novel algorithms based on lazy cohomology generators that enable to save orders of magnitude computational time with respect to all other state-of-the-art solutions proposed in literature. Believing that our results are useful in other fields of scientific computing, the proposed algorithm is presented as a detailed pseudocode in such a way that it can be easily implemented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Comparison of Machine Learning Approaches for Prediction of the Equivalent Alkane Carbon Number for Microemulsions Based on Molecular Properties

The chemical properties of oils are vital in the design of microemulsion systems. The hydrophilic–lipophilic difference equation used to predict microemulsions’ phase behavior expresses the oils’ physiochemical properties as the equivalent alkane carbon number (EACN). The experimental determination of EACN requires knowledge of the temperature dependence of the microemulsion system and the effects of different surfactant concentrations. Thus, the experimental determination is time-intensive and tedious, requiring days to months for proper separations. Furthermore, the experiments require high purity of chemicals because microemulsions are sensitive to impurities. Our work focuses on the quick and reliable predictions of the EACN with machine learning (ML) models. Due to the immaturity of ML chemical predictions, we compare three graph neural networks (GNNs) and a gradient-boosted tree algorithm, known as XGBoost. The GNNs use the molecular structures represented as simplified molecular-input line-entry system (SMILES) codes for the initial input, which allows us to assess whether geometry optimization is necessary for reliable results. The XGBoost model also begins with the SMILES representations of the molecules but uses molecular descriptors instead of geometry optimizations. As a result, the best model tested (crystal graph convolutional neural network with Merck molecular force field-94) has an error of 1.15 EACN units of the true EACN for unknown data with the errors skewed toward zero and an R² score of 0.9

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

Investigating Tropical Versus Extratropical Influences on the Southern Hemisphere Tropical Edge in the Unified Model

Abstract Since the late 1970s, observations have shown a widening of the tropical Hadley cell (HC) circulation. State‐of‐the‐art climate models reproduce the general trend along with a projected continuous expansion. Discrepancies in expansion rates of observation‐ and model‐based studies have been attributed to differences in applied methods, natural variability and model shortcomings. Furthermore, the driving influence of tropical or extratropical processes on these changes is not well understood. All of this highlights the dynamical mechanisms and the region of origin controlling the tropical width are still insufficiently understood. Here we examine the influence of systematic model biases of the atmosphere‐only Unified Model (UM) onto the simulation of the Southern Hemisphere (SH) tropical edge. We utilize nudged experiments with prescribed sea surface temperatures, where potential temperature and horizontal winds are relaxed back to ERA‐Interim reanalysis for a 20‐year period in selected regions. Correcting model biases in the tropics and extratropics separately allows us to dissect the dominant remote impacts of present model errors onto the SH tropical edge simulation. The experiments are applied to established tropical width metrics ranging from near‐surface to upper‐level metrics capturing the poleward flank of the HC. We find both regions work remotely to reduce errors in the UM fields and location of the tropical edge. Surprisingly, correcting the extratropical biases, south of 45°S, more consistently improves the tropical width across the metrics and seasons than nudging the tropics (10°N–10°S). These findings demonstrate the substantial role of extratropical influences in locating the SH tropical edge.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Tolerance analyses of a quadrupole magnet for advanced photon source upgrade

Given physics requirements, the mechanical fabrication and assembly tolerances for storage ring magnets can be calculated using analytical methods [1, 2]. However, this method is not easy for complicated magnet designs [1]. In this paper, a novel method is proposed to determine fabrication and assembly tolerances consistent with physics requirements, through a combination of magnetic and mechanical tolerance analyses. In this study, finite element analysis using OPERA is conducted to estimate the effect of fabrication and assembly errors on the magnetic field of a quadrupole magnet and to determine the allowable tolerances to achieve the specified magnetic performances. Based on the study, allowable fabrication and assembly tolerances for the quadrupole assembly are specified for the mechanical design of the quadrupole magnet. Next, to achieve the required assembly level tolerances, mechanical tolerance stackup analyses using a 3D tolerance analysis package are carried out to determine the part and subassembly level fabrication tolerances. This method can be used to determine the tolerances for design of other individual magnets and of magnet strings.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

A unified thermionic and thermionic-field emission (TE–TFE) model for ideal Schottky reverse-bias leakage current

We present a unified thermionic emission (TE) and thermionic-field emission (TFE) model for the ideal reverse-bias leakage current in Schottky junctions. The unified TE–TFE analytical model advances upon previous analytical TFE models by Murphy–Good and Padovani–Stratton, which are the two most widely adopted models by the community, in two major aspects: (i) the applicability of the TFE expression therein is extended to near-zero surface electric fields by an error-function correction, allowing for the calculation of the total current by a nontrivial sum of TE (over-the-barrier current) and TFE (below-the-barrier current) contributions; therefore, an accurate description of the TE-to-TFE transition region is captured analytically for the first time; (ii) image-force lowering is considered with much-simpler correction terms. Comparisons with the reference numerical model show that the unified TE–TFE model has excellent accuracy, as well as a 10 000× reduction in computation time. The unified model is further tested against experimental data from Schottky barrier diodes based on Si, 4H-SiC, GaN, and β-Ga 2 O 3 , revealing accurate extractions of barrier heights and correct descriptions of the ideal reverse leakage characteristics. With the extended applicable range, improved accuracy, and computational efficiency, the unified TE–TFE model is highly valuable for the design and analysis of devices based on Schottky junctions, as well as for potential integration in technology computer-aided design (TCAD) tools.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Uncovering multiple populations with washington photometry. I. The globular cluster NGC 1851

The analysis of multiple populations (MPs) in globular clusters (GCs) has become a forefront area of research in astronomy. Multiple red giant branches (RGBs), subgiant branches (SGBs), and even main sequences (MSs) have now been observed photometrically in many GCs, while broad abundance distributions of certain elements have been detected spectroscopically in most, if not all, GCs. UV photometry has been crucial in discovering and analyzing these MPs, but the Johnson U and the Stromgren and Sloan u filters that have generally been used are relatively inefficient and very sensitive to reddening and atmospheric extinction. In contrast, the Washington C filter is much broader and redder than these competing UV filters, making it far more efficient at detecting MPs and much less sensitive to reddening and extinction. Here, we investigate the use of the Washington system to uncover MPs using only a 1 m telescope. Our analysis of the well-studied GC NGC 1851 finds that the C filter is both very efficient and effective at detecting its previously discovered MPs in the RGB and SGB. Remarkably, we have also detected an intrinsically broad MS best characterized by two distinct but heavily overlapping populations that cannot be explained by binaries, field stars, or photometric errors. The MS distribution is in very good agreement with that seen on the RGB, with ∼30% of the stars belonging to the second population. There is also evidence for two sequences in the red horizontal branch, but this appears to be unrelated to the MPs in this cluster. Neither of these latter phenomena have been observed previously in this cluster. The redder MS stars are also more centrally concentrated than the blue MS. This is the first time MPs in an MS have been discovered from the ground, and using only a 1 m telescope. The Washington system thus proves to be a very powerful tool for investigating MPs, and holds particular promise for extragalactic objects where photons are limited.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Linearised Fokker–Planck collision model for gyrokinetic simulations

We introduce a gyrokinetic, linearised Fokker–Planck collision model that satisfies conservation laws and is accurate at arbitrary collisionalities. The differential test-particle component of the operator is exact; the integral field-particle component is approximated using a spherical harmonic and a modified Laguerre polynomial expansion developed by Hirshman and Sigmar (1976 Phys. Fluids 19 1532). The numerical methods of the implementation in the δf-gyrokinetic code stella (Barnes et al 2019 J. Comput. Phys. 391 365–80) are discussed, and conservation properties of the operator are demonstrated. The collision model is then benchmarked against the collision model of the gyrokinetic solver GS2 in the limiting cases of a reduced test-particle collision operator and energy- and momentum-conserving operator. The accuracy of the full collision model is investigated by solving the parallel Spitzer-Härm problem for the transport coefficients. It is shown that retaining collisional energy flux and higher-order terms in the field-particle operator reduces errors in the transport coefficients from 10%–25% for a simple momentum- and energy-conserving model to under 1%.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Cross-correlating 21cm intensity maps with Lyman Break Galaxies in the post-reionization era

We investigate the cross-correlation between the spatial distribution of Lyman Break Galaxies (LBGs) and the 21cm intensity mapping signal at z∼[3–5]. At these redshifts, galactic feedback is supposed to only marginally affect the matter power spectrum, and the neutral hydrogen distribution is independently constrained by quasar spectra. Using a high resolution N-body simulation, populated with neutral hydrogen a posteriori, we forecast for the expected LBG-21cm cross-spectrum and its error for a 21cm field observed by the Square Kilometre Array (SKA1-LOW and SKA1-MID), combined with a spectroscopic LBG survey with the same volume. The cross power can be detected with a signal-to-noise ratio (SNR) up to ∼10 times higher (and down to ∼ 4 times smaller scales) than the 21cm auto-spectrum for this set-up, with the SNR depending only very weakly on redshift and the LBG population. We also show that while both the 21cm auto- and LBG-21cm cross-spectra can be reliably recovered after the cleaning of smooth-spectrum foreground contamination, only the cross-power is robust to problematic non-smooth foregrounds like polarized synchrotron emission.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Quantum Simulation of the First-Quantized Pauli-Fierz Hamiltonian

We provide an explicit recursive divide-and-conquer approach for simulating quantum dynamics and derive a discrete first-quantized nonrelativistic QED Hamiltonian based on the many-particle Pauli-Fierz Hamiltonian. We apply this recursive divide-and-conquer algorithm to this Hamiltonian and compare it to a concrete simulation algorithm that uses qubitization. Our divide-and-conquer algorithm, using lowest-order Trotterization, scales for fixed grid spacing as O ~ ( Λ N 2 η 2 t 2 / ϵ ) for grid size N , η particles, simulation time t , field cutoff Λ , and error ϵ . Our qubitization algorithm scales as O ~ ( N ( η + N ) ( η + Λ 2 ) t log ( 1 / ϵ ) ) . This shows that even a naive partitioning and low-order splitting formula can yield, through our divide-and-conquer formalism, superior scaling to qubitization for large Λ . We compare the relative costs of these two algorithms on systems that are relevant for applications such as the spontaneous emission of photons and the photoionization of electrons. We observe that for different parameter regimes, one method can be favored over the other. Finally, we give new algorithmic and circuit-level techniques for gate optimization, including a new way of implementing a group of multicontrolled- X gates that can be used for better analysis of circuit cost. Published by the American Physical Society 2024

Mukhopadhyay, Priyanka (ORCID:0000000164639100)↗
DOE OSTI2024

Ambient-Frequency-Data Based System-Level Inertia Estimation Using Physical Equation and its Practice on Hawaii Islands

Here, in this paper, a practical ambient-frequency-data-based inertia estimation method using a physical equation is proposed and validated by real measurement data from Hawaii island grids. With high renewable penetration, accurate inertia estimation is important and urgent. Ambient frequency oscillation always exists in power grids, so the proposed method has advantages of real-time inertia estimation and no need for additional disturbances. This paper first developed the physical equation for inertia estimation to offer a clear mechanism for easy implementation in practice. To apply the proposed method in actual grids, a practical method to extract the ambient frequency oscillation from frequency measurement is further proposed. The inertia estimation using the physical equation is validated by KIUC simulation data and HECO field data, in which error rates are around 2% and 8%, respectively. Practical inertia estimation is challenging due to the large amounts of resources contributing to the power grid's effective inertia, but the method provided in this paper can offer a novel way for practical inertia estimation, which can help renewable penetration to boost carbon-free grid.

14 SOLAR ENERGY↗