Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “convergence analysis”

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,207 records · Page 67

Effects of threshold resummation for large-x PDF in large momentum effective theory

Parton distribution functions (PDFs) at large x are challenging to extract from experimental data, yet they are essential for understanding hadron structure and searching for new physics beyond the Standard Model. Within the framework of the large momentum P z expansion of lattice quasi-PDFs, we investigate large x PDFs, where the matching coefficient is factorized into the hard kernel, related to the active quark momentum xP z , and the threshold soft function, associated with the spectator momentum (1 − x)P z . The renormalization group equation of the soft function enables the resummation of the threshold double logarithms α k ln 2k (1 − x), which is crucial for a reliable and controllable calculation of large x PDFs. Our analysis with pion valence PDFs indicates that perturbative matching breaks down when the spectator momentum (1 − x)P z approaches Λ QCD , but remains valid when both xP z and (1 − x)P z are much larger than Λ QCD . Additionally, we incorporate leading renormalon resummation within the threshold framework, demonstrating good perturbative convergence in the region where both spectator and active quark momenta are perturbative scales.

factorization↗

Coexisting kagome and heavy fermion flat bands in YbCr 6 Ge 6

Flat bands, electronic states with nearly dispersionless energy-momentum structure, provide fertile ground for unconventional quantum phases. Recent observations of flat bands at the Fermi level in kagome metals open the possibility of unifying topology and correlation-driven heavy-fermion physics. Here we show that topology and heavy-fermion correlations coexist in the layered kagome metal YbCr 6 Ge 6 . At high temperatures, an intrinsic kagome flat band—arising from frustrated hopping on the kagome lattice—dominates the Fermi level. Upon cooling, localized Yb 4f-states hybridize with the topological kagome flat bands, transforming this state into momentum-independent Kondo resonance states across the entire Brillouin zone. Topological analysis of the hybridization gaps reveals filling-tunable weak and strong topological Kondo-insulating regimes, and identifies a topological Dirac–Kondo semimetal. Taken together, these results identify YbCr 6 Ge 6 as a prototype of a topological heavy-fermion system and a platform where geometric frustration, strong correlations, and topology converge, with broad implications for correlated quantum matter.

36 MATERIALS SCIENCE↗

Assessing MP2 frozen natural orbitals in relativistic correlated electronic structure calculations

The high computational scaling with the basis set size and the number of correlated electrons is a bottleneck limiting applications of coupled cluster algorithms, in particular for calculations based on two- or four-component relativistic Hamiltonians, which often employ uncontracted basis sets. This problem may be alleviated by replacing canonical Hartree–Fock virtual orbitals by natural orbitals (NOs). Here, in this paper, we describe the implementation of a module for generating NOs for correlated wavefunctions and, in particular, second order Møller–Plesset perturbation frozen natural orbitals (MP2FNOs) as a component of our novel implementation of relativistic coupled cluster theory for massively parallel architectures [Pototschnig et al. J. Chem. Theory Comput. 17, 5509, (2021)]. Our implementation can manipulate complex or quaternion density matrices, thus allowing for the generation of both Kramers-restricted and Kramers-unrestricted MP2FNOs. Furthermore, NOs are re-expressed in the parent atomic orbital (AO) basis, allowing for generating coupled cluster singles and doubles NOs in the AO basis for further analysis. By investigating the truncation errors of MP2FNOs for both the correlation energy and molecular properties—electric field gradients at the nuclei, electric dipole and quadrupole moments for hydrogen halides HX (X = F–Ts), and parity-violating energy differences for H 2 Z 2 (Z = O–Se)—we find MP2FNOs accelerate the convergence of the correlation energy in a roughly uniform manner across the Periodic Table. It is possible to obtain reliable estimates for both energies and the molecular properties considered with virtual molecular orbital spaces truncated to about half the size of the full spaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dark Energy Survey Year 3 results: Imprints of cosmic voids and superclusters in the Planck CMB lensing map

ABSTRACT The CMB lensing signal from cosmic voids and superclusters probes the growth of structure in the low-redshift cosmic web. In this analysis, we cross-correlated the Planck CMB lensing map with voids detected in the Dark Energy Survey Year 3 (Y3) data set (∼5000 deg2), expanding on previous measurements that used Y1 catalogues (∼1300 deg2). Given the increased statistical power compared to Y1 data, we report a 6.6σ detection of negative CMB convergence (κ) imprints using approximately 3600 voids detected from a redMaGiC luminous red galaxy sample. However, the measured signal is lower than expected from the MICE N-body simulation that is based on the ΛCDM model (parameters Ωm = 0.25, σ8 = 0.8), and the discrepancy is associated mostly with the void centre region. Considering the full void lensing profile, we fit an amplitude $A_{\kappa }=\kappa _{{\rm DES}}/\kappa _{{\rm MICE}}$ to a simulation-based template with fixed shape and found a moderate 2σ deviation in the signal with Aκ ≈ 0.79 ± 0.12. We also examined the WebSky simulation that is based on a Planck 2018 ΛCDM cosmology, but the results were even less consistent given the slightly higher matter density fluctuations than in MICE. We then identified superclusters in the DES and the MICE catalogues, and detected their imprints at the 8.4σ level; again with a lower-than-expected Aκ = 0.84 ± 0.10 amplitude. The combination of voids and superclusters yields a 10.3σ detection with an Aκ = 0.82 ± 0.08 constraint on the CMB lensing amplitude, thus the overall signal is 2.3σ weaker than expected from MICE.

79 ASTRONOMY AND ASTROPHYSICS↗

A Phenomenological Model for Cavitation

A phenomenological model of cavitation is presented, based on observations that both large relative negative pressures and large negative time derivatives of pressure are required for cavitation onset. We simulated two cavitation experiments to generate cavitation scaling parameters for relative pressure drop and rate of pressure drop. Our results show the model, while simple, is effective at reproducing results from laboratory experiments of cavitation. The parameters were then used in conjunction with a human surrogate computational model to predict, at any position within the head, the probability of intracranial cavitation caused by exposure to a blast event. The results suggest that the magnitude of blast overpressure observed in field data is sufficient to cause intracranial cavitation. Our analysis indicates that the helmeted head, when compared to the unhelmeted head configuration, results in a decrease but not elimination of cavitation exposure. When density functions of cavitation probability versus cumulative brain volume are combined with an injury severity model, the results show helmet efficacy at low and moderate risk levels. However, the convergence of unhelmeted and helmeted probability density functions at high-to-excessive risk thresholds indicates the helmet offers diminishing protection at elevated exposure levels, relative to the unhelmeted baseline. Future investigation and collaboration with neuroscience subject matter experts are needed to contextualize the current computational results. While the present work contributes specific and quantified predictions of intracranial cavitation location and severity, more research is required to apply our results to clinical settings with population-based brain injury subjects and controls. The relationship between our intracranial cavitation predictions with their anticipated clinical sequelae remains a topic in need of exploration.

60 APPLIED LIFE SCIENCES↗

Blockchain for Fault-Tolerant Grid Operations

Radial topology and vast geographic coverage make distribution systems prone to widespread power outages upon the failure of a single (or multiple) upstream component. Fault-handling algorithms depend heavily on correct estimations of the system’s state to effectively isolate the affected area and reduce the number of affected customers while maintaining operational safety. The work described here leverages the core features of distributed, consensus-based decision-making processes and the immutability of blockchain, and demonstrates their value in improving fault-tolerant grid operations. In this work, blockchain was used to create a trusted data-sharing platform that enables independent actors to reconstruct the system state; this enables distributed resources to make intelligent decisions with limited knowledge. Although the process requires data sharing, its algorithms have been designed to limit the amount of private information that is exchanged, which helps preserve business-sensitive data and maintain customer privacy. In addition, by reducing the information that must be shared, the communication requirements are also reduced; (however, an in-depth analysis of the communication requirements is beyond the scope of this project). The proposed use cases are intended to represent a foundational basis for third parties to develop functional solutions that can eventually be deployed in the field. To further provide guidance, the envisioned use cases have incorporated design requirements that consider the blockchain characteristics and a need to limit information from surrounding resources, which preserve the assumption and the possibility that such resources could belong to different entities. This report presents a detailed design of the three use cases with the tools needed to enable the analysis being tested. The implemented gross error detection method can detect mismatches when the error exceeds 3.8 times the sensor’s rated accuracy. Detection of the circuit breaker state successfully identified the correct states across all simulation tests. A distribution-system power-flow solution in the simulator OpenDSS generally possesses a convergency tolerance of 0.01% on the voltage magnitude. The evaluation of possible reconnection using voltage magnitude—preserving the data ownership—has a voltage magnitude difference smaller than 0.001% from the OpenDSS result. The results preserving data ownership have a difference within the expected power flow tolerance with full knowledge of the system, which surpasses expectations.

24 POWER TRANSMISSION AND DISTRIBUTION↗

The reduction of remote sensing data by visual means

Issues likely to be of concern to educators called upon to teach courses involving the reduction (interpretation) of remotely sensed data by visual means are considered. Topics covered include: (1) information requirements of those using remotely-sensed data; (2) educational concepts involved in teaching students how to generate the desired information from a visual analysis of the data; (3) principles and techniques specific to the photointerpretation process; (4) concepts involved in the making of photographic measurements, as dictated by the geometry of remote sensing imagery; (5) the nature of the various kinds of mapping, plotting, and photointerpretation equipment; and (6) some special considerations with respect to the convergence of evidence and other principles involved in the interpretation of photographs. A recommended procedure for determining the usefulness of any given type of aerial or space photography in relation to the inventory of natural resources is included.

Colwell, R. N.↗

Multigrid applications to three-dimensional elliptic coordinate generation

The multigrid algorithm was applied to solve the coupled set of elliptic quasilinear partial differential equations associated with three-dimensional coordinate generation. The results indicate that the multigrid scheme is more than twice as fast as conventional relaxation schemes on moderate-size grids. Convergence factors of order 0.90 per work unit were achieved on 36,000-point grids. The paper covers the form of transformation, develops the set of generation equations, and gives details on the multigrid approach used. Included are a development of the full-approximation storage scheme, details of the smoothing-rate analysis, and a section devoted to rational programming techniques applicable to the multigrid algorithm.

Thames, F. C.↗

Boundary layer development as a function of chamber pressure in the NASA Lewis 1030:1 area ratio rocket nozzle

Through the use of theoretical predictions of fluid properties and experimental heat transfer and thrust measurements, the zones of laminar, transitional, and turbulent boundary layer flow were defined for the NASA Lewis 1030:1 area ratio rocket nozzle. Tests were performed on the nozzle at chamber pressures from 350 to 100 psia. For these conditions, the throat diameter Reynolds numbers varied from 300,000 to 1 million. The propellants used were gaseous hydrogen and gaseous oxygen. Thrust measurements and nozzle outer wall temperature measurements were taken during the 3-sec test runs. Comparison of experimental heat transfer and thrust data with the corresponding predictions from the Two-Dimensional Kinetics (TDK) nozzle analysis program indicated laminar flow in the nozzle at a throat diameter Reynolds number of 320,000 or chamber pressure of 360 psia. Comparison of experimental and predicted heat transfer data indicated transitional flow up to and including a chamber pressure of 1000 psia. Predicted values of the axisymmetric acceleration parameter within the convergent and divergent nozzle were consistent with the above results. Based upon an extrapolation of the heat transfer data and predicted distributions of the axisymmetric acceleration parameter, transitional flow was predicted up to a throat diameter Reynolds number of 220,000 or 2600-psia chamber pressure. Above 2600-psia chamber pressure, fully developed turbulent flow was predicted.

Smith, Tamara A.↗

Dual adaptive control: Design principles and applications

The design of an actively adaptive dual controller based on an approximation of the stochastic dynamic programming equation for a multi-step horizon is presented. A dual controller that can enhance identification of the system while controlling it at the same time is derived for multi-dimensional problems. This dual controller uses sensitivity functions of the expected future cost with respect to the parameter uncertainties. A passively adaptive cautious controller and the actively adaptive dual controller are examined. In many instances, the cautious controller is seen to turn off while the latter avoids the turn-off of the control and the slow convergence of the parameter estimates, characteristic of the cautious controller. The algorithms have been applied to a multi-variable static model which represents a simplified linear version of the relationship between the vibration output and the higher harmonic control input for a helicopter. Monte Carlo comparisons based on parametric and nonparametric statistical analysis indicate the superiority of the dual controller over the baseline controller.

Mookerjee, Purusottam↗

Boundary layer development as a function of chamber pressure in the NASA Lewis 1030:1 area ratio rocket nozzle

Through the use of theoretical predictions of fluid properties and experimental heat transfer and thrust measurements, the zones of laminar, transitional, and turbulent boundary layer flow were defined for the NASA Lewis 1039:1 area ratio rocket nozzle. Tests were performed on the nozzle at chamber pressures from 350 to 100 psia. For these conditions, the throat diameter Reynolds numbers varied from 300,000 to 1 million. The propellants used were gaseous hydrogen and gaseous oxygen. Thrust measurements and nozzle outer wall temperature measurements were taken during the 3-sec test runs. Comparison of experimental heat transfer and thrust data with the corresponding predictions from the Two-Dimensional Kinetics (TDK) nozzle analysis program indicated laminar flow in the nozzle at a throat diameter Reynolds number of 320,000 or chamber pressure of 360 psia. Comparison of experimental and predicted heat transfer data indicated transitional flow up to and including a chamber pressure of 1000 psia. Predicted values of the axisymmetric acceleration parameter within the convergent and divergent nozzle were consistent with the above results. Based upon an extrapolation of the heat transfer data and predicted distributions of the axisymmetric acceleration parameter, transitional flow was predicted up to a throat diameter Reynolds number of 220,000 or 2600-psia chamber pressure. Above 2600-psia chamber pressure, fully developed turbulent flow was predicted.

Smith, Tamara A.↗

Blocking reduction of Landsat Thematic Mapper JPEG browse images using optimal PSNR estimated spectra adaptive postfiltering

Two representative sample images of Band 4 of the Landsat Thematic Mapper are compressed with the JPEG algorithm at 8:1, 16:1 and 24:1 Compression Ratios for experimental browsing purposes. We then apply the Optimal PSNR Estimated Spectra Adaptive Postfiltering (ESAP) algorithm to reduce the DCT blocking distortion. ESAP reduces the blocking distortion while preserving most of the image's edge information by adaptively postfiltering the decoded image using the block's spectral information already obtainable from each block's DCT coefficients. The algorithm iteratively applied a one dimensional log-sigmoid weighting function to the separable interpolated local block estimated spectra of the decoded image until it converges to the optimal PSNR with respect to the original using a 2-D steepest ascent search. Convergence is obtained in a few iterations for integer parameters. The optimal logsig parameters are transmitted to the decoder as a negligible byte of overhead data. A unique maxima is guaranteed due to the 2-D asymptotic exponential overshoot shape of the surface generated by the algorithm. ESAP is based on a DFT analysis of the DCT basis functions. It is implemented with pixel-by-pixel spatially adaptive separable FIR postfilters. PSNR objective improvements between 0.4 to 0.8 dB are shown together with their corresponding optimal PSNR adaptive postfiltered images.

Linares, Irving↗

A Discrete Constraint for Entropy Conservation and Sound Waves in Cloud-Resolving Modeling

Ideal cloud-resolving models contain little-accumulative errors. When their domain is so large that synoptic large-scale circulations are accommodated, they can be used for the simulation of the interaction between convective clouds and the large-scale circulations. This paper sets up a framework for the models, using moist entropy as a prognostic variable and employing conservative numerical schemes. The models possess no accumulative errors of thermodynamic variables when they comply with a discrete constraint on entropy conservation and sound waves. Alternatively speaking, the discrete constraint is related to the correct representation of the large-scale convergence and advection of moist entropy. Since air density is involved in entropy conservation and sound waves, the challenge is how to compute sound waves efficiently under the constraint. To address the challenge, a compensation method is introduced on the basis of a reference isothermal atmosphere whose governing equations are solved analytically. Stability analysis and numerical experiments show that the method allows the models to integrate efficiently with a large time step.

Zeng, Xi-Ping↗

Fluctuations in Conjunction Miss Distance Projections as Time Approaches Time of Closest Approach

A responsibility of the Trajectory Operations Officer is to ensure that the International Space Station (ISS) avoids colliding with debris. United States Space Command (USSPACECOM) tracks and catalogs a portion of the debris in Earth orbit, but only objects with a perigee less than 600 km and a radar cross section (RCS) greater than 10 cm-objects that, in fact, represent only a small fraction of the objects in Earth orbit. To accommodate for this, the ISS uses shielding to protect against collisions with smaller objects. This study provides a better understanding of how quickly, and to what degree, USSPACECOM projections tend to converge to the final, true miss distance. The information included is formulated to better predict the behavior of miss distance data during real-time operations. It was determined that the driving components, in order of impact on miss distance fluctuations, are energy dissipation rate (EDR), RCS, and inclination. Data used in this analysis, calculations made, and conclusions drawn are stored in Microsoft Excel log sheets. A separate log sheet, created for each conjunction, contains information such as predicted miss distances, apogee and perigee of debris orbit, EDR, RCS, inclination, tracks and observations, statistical data, and other evaluation/orbital parameters.

Christian, John A., III↗

Flight Results of the NF-15B Intelligent Flight Control System (IFCS) Aircraft with Adaptation to a Longitudinally Destabilized Plant

Adaptive flight control systems have the potential to be resilient to extreme changes in airplane behavior. Extreme changes could be a result of a system failure or of damage to the airplane. The goal for the adaptive system is to provide an increase in survivability in the event that these extreme changes occur. A direct adaptive neural-network-based flight control system was developed for the National Aeronautics and Space Administration NF-15B Intelligent Flight Control System airplane. The adaptive element was incorporated into a dynamic inversion controller with explicit reference model-following. As a test the system was subjected to an abrupt change in plant stability simulating a destabilizing failure. Flight evaluations were performed with and without neural network adaptation. The results of these flight tests are presented. Comparison with simulation predictions and analysis of the performance of the adaptation system are discussed. The performance of the adaptation system is assessed in terms of its ability to stabilize the vehicle and reestablish good onboard reference model-following. Flight evaluation with the simulated destabilizing failure and adaptation engaged showed improvement in the vehicle stability margins. The convergent properties of this initial system warrant additional improvement since continued maneuvering caused continued adaptation change. Compared to the non-adaptive system the adaptive system provided better closed-loop behavior with improved matching of the onboard reference model. A detailed discussion of the flight results is presented.

Bosworth, John T.↗

Air Pollution Over the Ganges Basin and Northwest Bay of Bengal in the Early Postmonsoon Season Based on NASA MERRAero Data

The MERRA Aerosol Reanalysis (MERRAero) has been recently developed at NASA's Global Modeling Assimilation Office. This reanalysis is based on a version of the Goddard Earth Observing System-5 (GEOS-5) model radiatively coupled with Goddard Chemistry, Aerosol, Radiation, and Transport aerosols, and it includes assimilation of bias-corrected aerosol optical thickness (AOT) from the Moderate Resolution Imaging Spectroradiometer (MODIS) sensor on both Terra and Aqua satellites. In October over the period 2002-2009, MERRAero showed that AOT was lower over the east of the Ganges basin than over the northwest of the Ganges basin: this was despite the fact that the east of the Ganges basin should have produced higher anthropogenic aerosol emissions because of higher population density, increased industrial output, and transportation. This is evidence that higher aerosol emissions do not always correspond to higher AOT over the areas where the effects of meteorological factors on AOT dominate those of aerosol emissions. MODIS AOT assimilation was essential for correcting modeled AOT mainly over the northwest of the Ganges basin, where AOT increments were maximal. Over the east of the Ganges basin and northwest Bay of Bengal (BoB), AOT increments were low and MODIS AOT assimilation did not contribute significantly to modeled AOT. Our analysis showed that increasing AOT trends over northwest BoB (exceeding those over the east of the Ganges basin) were reproduced by GEOS-5, not because of MODIS AOT assimilation butmainly because of the model capability of reproducing meteorological factors contributing to AOT trends. Moreover, vertically integrated aerosol mass flux was sensitive to wind convergence causing aerosol accumulation over northwest BoB.

AOT↗

A Study of Parallel Scalability and Dynamic Workload Balancing in GlennICE

The Glenn Icing Computational Environment (GlennICE) is a computational tool designed to calculate ice growth on complex three-dimensional geometries. It utilizes user-supplied computational fluid dynamics solutions for the geometry of interest. Key developments include advancements in convergence of collection efficiency, trajectory optimization, and refinement methodology. These improvements have significantly enhanced GlennICE’s efficiency for practical engineering applications. A recent study focused on benchmarking GlennICE’s scalability in a parallel environment using static scheduling. Findings indicated a potential twofold increase in efficiency through workload balance enhancements. This paper presents an analysis of the solver’s new workload balancing improvements, incorporating shared memory and dynamic scheduling routines. Results demonstrate a highly efficient and consistent algorithm across high-performance computing clusters.

Computational Icing↗

A Study of Parallel Scalability and Dynamic Workload Balancing in GlennICE

The Glenn Icing Computational Environment (GlennICE) is a computational tool designed to calculate ice growth on complex three-dimensional geometries. It utilizes user-supplied computational fluid dynamics solutions for the geometry of interest. Key developments include advancements in convergence of collection efficiency, trajectory optimization, and refinement methodology. These improvements have significantly enhanced GlennICE’s efficiency for practical engineering applications. A recent study focused on benchmarking GlennICE’s scalability in a parallel environment using static scheduling. Findings indicated a potential twofold increase in efficiency through workload balance enhancements. This paper presents an analysis of the solver’s new workload balancing improvements, incorporating shared memory and dynamic scheduling routines. Results demonstrate a highly efficient and consistent algorithm across high-performance computing clusters.

Computational Icing↗