Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reduced basis method”

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 433 records · Page 24

PDG: A Composite Method Based on the Resolution of the Identity

The Gaussian-3 (G3) composite approach for thermochemical properties is revisited in light of the enhanced computational efficiency and reduced memory costs by applying the resolution-of-the-identity (RI) approximation for two-electron repulsion integrals (ERIs) to the computationally demanding component methods in the G3 model: the energy and gradient computations via the second-order Møller–Plesset perturbation theory (MP2) and the energy computations using the coupled-cluster singles–doubles method augmented with noniterative triples corrections [CCSD(T)]. Efficient implementation of the RI-based methods is achieved by employing a hybrid distributed/shared memory model based on MPI and OpenMP. The new variant of the G3 composite approach based on the RI approximation is termed the RI-G3 scheme, or alternatively the PDG method. The accuracy of the new RI-G3/PDG scheme is compared to the “standard” G3 composite approach that employs the memory-expensive four-center ERIs in the MP2 and CCSD(T) calculations. Taking the computation of the heats of formation of the closed-shell molecules in the G3/99 test set as a test case, it is demonstrated that the RI approximation introduces negligible changes to the mean absolute errors relative to the standard G3 model (less than 0.1 kcal/mol), while the standard deviations remain unaltered. Furthermore, the efficiency and memory requirements for the RI-MP2 and RI-CCSD(T) methods are compared to the standard MP2 and CCSD(T) approaches, respectively. The hybrid MPI/OpenMP-based RI-MP2 energy plus gradient computation is found to attain a 7.5× speedup over the standard MP2 calculations. For the most demanding CCSD(T) calculations, the application of the RI approximation is found to nearly halve the memory demand, confer about a 4–5× speedup for the CCSD iterations, and reduce the computational time for the compute-intensive triples correction step by several hours.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A hybrid surrogate modeling framework for the Digital Twin of a Fluoride-salt-cooled High-temperature Reactor (FHR)

While nuclear energy is a non-greenhouse-gas emitting energy source, expensive operational costs due to the high-level of safety requirements decreases their competitiveness in the sustainable energy market. Advanced reactor concepts paired with Digital Twins aim to increase the commercialization gains of nuclear energy by reducing operational costs, increasing reactor reliability and enhancing power generation. To support Digital Twin tasks such as real-time autonomous control, proactive maintenance monitoring or optimizing power demand operations, a fast and accurate virtual representation of the Nuclear Power Plant (NPP) is required. The computational cost of high-fidelity, physics-based models are unsuitable for real-time analysis or scalability. Here, in this work, a hybrid surrogate modeling framework is developed fora Fluoride-salt-cooled High-temperature Reactor (FHR) that leverages physics-inspired models for key reactor components and uses data-driven methods for rapid system state space prediction. The Xenon reactivity feedback model is integrated to inform the surrogate model about the reactor core and the homologous pump theory model is the basis for representing pump degradation. Using a detailed, two dimensional thermal hydraulics model to generate data on the FHR, we train a network of Vectorized Autoregressive Moving-Average with eXogenous input (VARMAX) models to predict the remaining state values. The result is a surrogate model that provides a detailed reactor state representation of 41 system states and a pump degradation analysis. The framework is applied to Load Follows profiles, yielding high accuracy and a speedup that is more than 4000x faster compared to the higher- fidelity thermal hydraulics model, enabling real-time operational intelligence and applications in long horizon predictions. While the surrogate model framework is demonstrated for the particular case of FHR, the hybrid physical/data-driven modeling approach including the network of surrogates and the underlying modularity has the potential to be applied to other physical asset systems.

Digital Twins↗

Application of automated iterative target detection for standoff hyperspectral imaging

The utility of hyperspectral imaging (HSI) has been well established for a wide array of applications but has generated a need for automated screening of high volumes of large HSI cubes. We report two important automated algorithms for more efficient standoff processing: atmospheric correction and target detection. The atmospheric correction method is based on a fast asymmetric least squares approach that is applied on a pixel-by-pixel basis. Here, the correction can be applied to entire images without manually identifying regions of interest and utilizes only in-scene information, no ancillary modeling of the atmosphere is required. An iterative target detection approach is also introduced which demonstrates faster speeds relative to moving window approaches. The target detection algorithm classifies each pixel as true target detections, near target detections, clutter, and no-calls. The algorithms were tested on forty images of twenty-two solid mineral targets placed at a 14-meter standoff distance allowing general observations on expected detection performance for a variety of minerals. In addition to identifying anomalous pixels, the inclusion of “no-calls” reduced the number of false detections significantly.

47 OTHER INSTRUMENTATION↗

Application of a low order panel method to complex three-dimensional internal flow problems

An evaluation of the ability of a low order panel method to predict complex three-dimensional internal flow fields was made. The computer code VSAERO was used as a basis for the evaluation. Guidelines for modeling internal flow geometries were determined and the effects of varying the boundary conditions and the use of numerical approximations on the solutions accuracy were studied. Several test cases were run and the results were compared with theoretical or experimental results. Modeling an internal flow geometry as a closed box with normal velocities specified on an inlet and exit face provided accurate results and gave the user control over the boundary conditions. The values of the boundary conditions greatly influenced the amount of leakage an internal flow geometry suffered and could be adjusted to eliminate leakage. The use of the far-field approximation to reduce computation time influenced the accuracy of a solution and was coupled with the values of the boundary conditions needed to eliminate leakage. The error induced in the influence coefficients by using the far-field approximation was found to be dependent on the type of influence coefficient, the far-field radius, and the aspect ratio of the panels.

Ashby, D. L.↗

Method and apparatus for white-light dispersed-fringe interferometric measurement of corneal topography

An novel interferometric apparatus and method for measuring the topography of aspheric surfaces, without requiring any form of scanning or phase shifting. The apparatus and method of the present invention utilize a white-light interferometer, such as a white-light Twyman-Green interferometer, combined with a means for dispersing a polychromatic interference pattern, using a fiber-optic bundle and a disperser such as a prism for determining the monochromatic spectral intensities of the polychromatic interference pattern which intensities uniquely define the optical path differences or OPD between the surface under test and a reference surface such as a reference sphere. Consequently, the present invention comprises a snapshot approach to measuring aspheric surface topographies such as the human cornea, thereby obviating vibration sensitive scanning which would otherwise reduce the accuracy of the measurement. The invention utilizes a polychromatic interference pattern in the pupil image plane, which is dispersed on a point-wise basis, by using a special area-to-line fiber-optic manifold, onto a CCD or other type detector comprising a plurality of columns of pixels. Each such column is dedicated to a single point of the fringe pattern for enabling determination of the spectral content of the pattern. The auto-correlation of the dispersed spectrum of the fringe pattern is uniquely characteristic of a particular optical path difference between the surface under test and a reference surface.

Hochberg, Eric B.↗

Chemical application of diffusion quantum Monte Carlo

The diffusion quantum Monte Carlo (QMC) method gives a stochastic solution to the Schroedinger equation. This approach is receiving increasing attention in chemical applications as a result of its high accuracy. However, reducing statistical uncertainty remains a priority because chemical effects are often obtained as small differences of large numbers. As an example, the single-triplet splitting of the energy of the methylene molecule CH sub 2 is given. The QMC algorithm was implemented on the CYBER 205, first as a direct transcription of the algorithm running on the VAX 11/780, and second by explicitly writing vector code for all loops longer than a crossover length C. The speed of the codes relative to one another as a function of C, and relative to the VAX, are discussed. The computational time dependence obtained versus the number of basis functions is discussed and this is compared with that obtained from traditional quantum chemistry codes and that obtained from traditional computer architectures.

Reynolds, P. J.↗

Enhancement of the no-vent fill process

This paper presents an analytical and experimental evaluation of an enhanced techniques for no-vent fill. The method entails injecting liquid through the top of the receiver vessel, thereby increasing surface area and agitation of the ullage/liquid interface. Both of these factors promote condensation induced ullage collapse, and reduce compressive impedance to the incoming liquid. The enhanced process was analyzed by modifying the surface area algorithm of an existing tank thermodynamic code to model a downward-pointing, conical jet impringing on a steadily rising liquid surface. Transient pressure and temperature measurements from several tests with Freon-114 were input into the revised model to calculate condensation rate as a function of fill level. By expressing these rates in dimensionless form (i.e., in terms of Stanton number and Prandtl number), an empirical correlation similar to the submerged jet model of Brown and Sonin (1989) was derived. This provided a basis for developing an expression which relates top fill to bottom fill performance.

Vaughan, David A.↗

High Speed Civil Transport (HSCT) Isolated Nacelle Transonic Boattail Drag Study and Results Using Computational Fluid Dynamics (CFD)

Nozzle boattail drag is significant for the High Speed Civil Transport (HSCT) and can be as high as 25% of the overall propulsion system thrust at transonic conditions. Thus, nozzle boattail drag has the potential to create a thrust-drag pinch and can reduce HSCT aircraft aerodynamic efficiencies at transonic operating conditions. In order to accurately predict HSCT performance, it is imperative that nozzle boattail drag be accurately predicted. Previous methods to predict HSCT nozzle boattail drag were suspect in the transonic regime. In addition, previous prediction methods were unable to account for complex nozzle geometry and were not flexible enough for engine cycle trade studies. A computational fluid dynamics (CFD) effort was conducted by NASA and McDonnell Douglas to evaluate the magnitude and characteristics of HSCT nozzle boattail drag at transonic conditions. A team of engineers used various CFD codes and provided consistent, accurate boattail drag coefficient predictions for a family of HSCT nozzle configurations. The CFD results were incorporated into a nozzle drag database that encompassed the entire HSCT flight regime and provided the basis for an accurate and flexible prediction methodology.

Midea, Anthony C.↗

High Speed Civil Transport (HSCT) Isolated Nacelle Transonic Boattail Drag Study and Results Using Computational Fluid Dynamics (CFD)

Nozzle boattail drag is significant for the High Speed Civil Transport (HSCT) and can be as high as 25 percent of the overall propulsion system thrust at transonic conditions. Thus, nozzle boattail drag has the potential to create a thrust drag pinch and can reduce HSCT aircraft aerodynamic efficiencies at transonic operating conditions. In order to accurately predict HSCT performance, it is imperative that nozzle boattail drag be accurately predicted. Previous methods to predict HSCT nozzle boattail drag were suspect in the transonic regime. In addition, previous prediction methods were unable to account for complex nozzle geometry and were not flexible enough for engine cycle trade studies. A computational fluid dynamics (CFD) effort was conducted by NASA and McDonnell Douglas to evaluate the magnitude and characteristics of HSCT nozzle boattail drag at transonic conditions. A team of engineers used various CFD codes and provided consistent, accurate boattail drag coefficient predictions for a family of HSCT nozzle configurations. The CFD results were incorporated into a nozzle drag database that encompassed the entire HSCT flight regime and provided the basis for an accurate and flexible prediction methodology.

Midea, Anthony C.↗

Investigation of the NO reduction by CO reaction over oxidized and reduced NiO x /CeO 2 catalysts

CeO 2 -supported NiO x catalysts have been widely studied in various catalytic reactions including NO reduction by CO. This work is mainly focused on investigation of the impact of catalyst synthesis conditions (e.g., oxidation and reduction) on the physicochemical properties of NiO x /CeO 2 catalysts and the catalytic response for the NO reduction by CO reaction. The oxide NiO x /CeO 2 sample was prepared by an incipient wetness impregnation (IWI) method and reduced under hydrogen reduction treatment at high temperatures (500 and 700 °C). The physicochemical properties of the synthesized samples were characterized by BET analysis, Raman spectroscopy, XRD, XPS, EELS and high-resolution transmission electron microscopy (HR-TEM). The results showed that higher reduction temperature led to the decrease in specific surface area (SSA), fewer oxygen vacancy/defect site, larger crystallite size of the CeO 2 support, and formation of metallic Ni on the surface. The oxidized NiO x /CeO 2 catalyst showed the highest catalytic activity, indicating that the presence of oxygen vacancy/defect sites, Ni 2+ oxidation state, and smaller crystallite size are believed to enhance the catalytic activity. In situ DRIFTS confirmed the generation of several intermediate species, such as nitrate, carbonate, and N 2 O. Finally, on the basis of in situ DRIFTS and activity results, the possible reaction mechanism of NO reduction by CO over NiO x /CeO 2 was proposed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

GPU acceleration of rank-reduced coupled-cluster singles and doubles

Here, we have developed a graphical processing unit (GPU) accelerated implementation of our recently introduced rank-reduced coupled-cluster singles and doubles (RR-CCSD) method. RR-CCSD introduces a low-rank approximation of the doubles amplitudes. This is combined with a low-rank approximation of the electron repulsion integrals via Cholesky decomposition. The result of these two low-rank approximations is the replacement of the usual fourth-order CCSD tensors with products of second- and third-order tensors. In our implementation, only a single fourth-order tensor must be constructed as an intermediate during the solution of the amplitude equations. Owing in large part to the compression of the doubles amplitudes, the GPU-accelerated implementation shows excellent parallel efficiency (95% on eight GPUs). Our implementation can solve the RR-CCSD equations for up to 400 electrons and 1550 basis functions—roughly 50% larger than the largest canonical CCSD computations that have been performed on any hardware. In addition to increased scalability, the RR-CCSD computations are faster than the corresponding CCSD computations for all but the smallest molecules. We test the accuracy of RR-CCSD for a variety of chemical systems including up to 1000 basis functions and determine that accuracy to better than 0.1% error in the correlation energy can be achieved with roughly 95% compression of the ov space for the largest systems considered. We also demonstrate that conformational energies can be predicted to be within 0.1 kcal mol -1 with efficient compression applied to the wavefunction. Finally, we find that low-rank approximations of the CCSD doubles amplitudes used in the similarity transformation of the Hamiltonian prior to a conventional equation-of-motion CCSD computation will not introduce significant errors (on the order of a few hundredths of an electronvolt) into the resulting excitation energies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

AN OVERVIEW OF PLASMA PROCESSING OF SRF CAVITIES AT JLAB

Plasma processing is a common technique where the free oxygen produced in a low-pressure RF plasma breaks down and removes hydrocarbons from surfaces. This increases the work function and reduces the secondary emission coefficient of the treated surfaces. Jefferson Lab has an ongoing R&D program in plasma processing. The experimental program investigated processing using argon/oxygen and helium/oxygen gas mixtures. The initial focus of the effort was processing C100 cavities by injecting RF power into the HOM coupler ports. We also developed the methods for establishing a plasma C75 cavities where the RF power is injected via the fundamental power-coupler. As part of the process development we processed, three C100 cryomodules in our off-line cryomodule test facility. In May 2023 we processed four C100 cryomodules in-situ in the CEBAF accelerator with the cryomodules returning to an operational status in Sept. 2023. The cumulative improvement in field emission free operation, as measured on a cavity by cavity basis, was 59 MeV or 24%. We recently started processing 7 cryomodules in the CEBAF accelerator in August 2024. Methods systems and results from processing cryomodules and individual cavities in the vertical test will be presented. Current status and future plans will also be presented.

Ganey, Tiffany↗

Analytical evaluation of a surface integral expressing the coupling between interior and exterior volumes in a FE-IE approach

Recently an approach which combines the finite element technique and an integral equation to determine the fields scattered by inhomogeneous bodies of complicated shape has been proposed. Basically, a mathematical surface which encloses the scatterers is introduced, thus dividing the space into an interior and an exterior volume, in which the finite element technique and an integral equation for EM scattering, respectively, are applied. The integral equation is set up for the tangential components of the fields at the surface, while the interior volume the unknowns are the total fields. Continuity of the tangential fields at the boundary, as required by Maxwell's equations, is imposed, thus coupling the two methods to obtain a consistent solution. The coupling term is expressed by a surface integral formed by the dot product of a FE basis function and an IE testing function, or viceversa. By choosing the boundary to be a surface of revolution and by making a convenient selection of IE basis (testing) functions, it is possible to evaluate the integrals analytically on surfaces such as curved triangles, curved quadrilaterals and curved pentagons. We will illustrate the salient steps involved in setting up and carrying out these integrals and discuss what class of basis (testing) functions and analytic surfaces of revolution they are applicable to. Analytic calculations offer the advantage of better accuracy than purely numerical ones, and, when combined with them, often shed light on issues of numerical convergence and limiting values. Furthermore, they may reduce computation time and storage requirements.

Zuffada, C.↗

Efficient core-excited state orbital perspective on calculating x-ray absorption transitions in determinant framework

X-ray absorption spectroscopy (XAS) is an explicit probe of the unoccupied electronic structure of materials and an invaluable tool for fingerprinting various electronic properties and phenomena. Computational methods capable of simulating and analyzing such spectra are therefore in high demand for complementing experimental results and for extracting valuable insights therefrom. In particular, a recently proposed first-principles approach titled many-body XAS (MBXAS), which approximates the final (initial) state as a Slater determinant constructed from Kohn-Sham (KS) orbitals optimized in the absence (presence) of the relevant core electron has shown promising prospects in evaluating the transition amplitudes. In this article, we show that the MBXAS approach can be rederived using a transition operator expressed entirely in the basis of core-excited state KS orbitals and that this reformulation offers substantial practical and conceptual advantages. In addition to circumventing previous issues of convergence with respect to the number of unoccupied ground-state orbitals, the aforementioned representation reduces the computational expense by rendering the calculation of such orbitals unnecessary altogether. Furthermore, the reformulated approach also provides a direct pathway for comparing the many-body approximation with the so-called single-particle treatment and indicates the relative importance in observed XAS intensity of the relaxation of the valence occupied subspace induced by the core excitation. Finally, using the core-excited state basis, we define auxiliary orbitals for x-ray absorption and demonstrate their utility in explaining the spectral intensity by contrasting them with single-particle approximations to the excited state.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

A modified Monte Carlo model for the ionospheric heating rates.

A Monte Carlo method is adopted as a basis for the derivation of the photoelectron-heat input into the ionospheric plasma. Since a great number of Monte Carlo runs are required normally for the computation of the heating rates, this approach is modified in an attempt to minimize the computation time. The heat-input distributions are computed for arbitrarily small source elements that are spaced apart at distances corresponding to the photoelectron dissipation range. By means of a nonlinear interpolation procedure their individual heating-rate distributions are utilized to produce synthetic ones that fill the gaps between the Monte Carlo generated distributions. By varying these gaps and the corresponding number of Monte Carlo runs the accuracy of the results is tested to verify the validity of this procedure. It is concluded that this model can reduce the computation time by as much as an order of magnitude, thus improving the feasibility of including Monte Carlo calculations in self-consistent ionosphere models.

Mayr, H. G.↗

Novel Engineering and Fabrication Techniques Tested in Low-Noise- Research Fan Blades

A major source of fan noise in commercial turbofan engines is the interaction of the wake from the fan blades with the stationary vanes (stators) directly behind them. The Trailing Edge Blowing (TEB) project team at the NASA Glenn Research Center designed and fabricated new fan blades to study the effects of fan trailing edge blowing as a potential noise-reduction concept. The intent is to fill the rotor wake by supplying air to the rotor blade trailing edge at the proper conditions to minimize the wake deficit, and thus generate less noise. The TEB hardware is designed for the Active Noise Control Fan (ANCF) test rig in Glenn's Aeroacoustic Propulsion Laboratory. For this test, the air is fed from an external supply through the shaft of the rig. It is distributed to the base of each blade through an impeller, where it is forced into a plenum at the core of each blade. In actual engine configuration, air would most likely be bled from the compressor, but only at times when noise is an issue, such as takeoffs and landings. Glenn researchers designed and manufactured the blades in-house, using new techniques and concepts. The skins, which were designed for maximum strength in the directions of highest stress, were molded from multiple layers of carbon fiber. Considerable use was made of rapid prototyping techniques, such as laser sintering. The core was sintered from a lightweight polymer, and the retainer was CNC-machined (computer numerical control machined) from aluminum. All the components were joined with a cold-cure aerospace adhesive. These techniques and processes reduced the overall cost and allowed the new concept to be studied much sooner than would be possible using traditional fabrication methods. Since this test rig did not support the use of blade-monitoring techniques such as strain gauges, extensive bench testing was required to qualify the design. The blades were examined using a variety of methods including holography, pull tests (cyclic and failure), shake tests, rap tests, and nondestructive inspection. Acoustic testing of the ANCF fan using TEB has been ongoing since January of 2001. The fan has completed about 100 hr of testing with no structural, vibrational, or fatigue problems. Far-field acoustic measurements, in-duct mode measurements, precise hot wire surveys, and detailed performance measurements are providing data for evaluating the concept. The far-field noise data show that tone noise was reduced significantly with the initial ANCF TEB fan design. In addition, a significant reduction in unsteady stator loading has been measured, indicating the potential for stator broadband noise reduction. The acoustic benefits will be assessed and important design parameters identified to improve the ability to fully exploit any benefit provided by this technique. On the basis of the success of trailing edge blowing, Glenn plans to continue this research with a higher speed, higher pressure ratio fan operating in an acoustic wind tunnel to simulate flight conditions.

Cunningham, Cameron C.↗

Synchrotron Photoionization Study of the Diisopropyl Ether Oxidation

Abstract Scientific evidence has shown oxygenates help to reduce dangerous pollutants arising from burning fossil fuel in the automotive sector. For this reason, their use as additives has spread widely. The aim of this work consists in providing a comprehensive identification of the main primary oxidation products of diisopropyl ether (DIPE), one of the most promising among etheric oxygenates. The Cl‐initiated oxidation of DIPE is examinated by using a vacuum ultraviolet (VUV) synchrotron radiation at the Advanced Light Source (ALS) of the Lawrence Berkeley National Laboratory (LBNL). Products are identified on the basis of their mass‐to‐charge ratio, shape of photoionization spectra, adiabatic ionization energies, and chemical kinetic profiles, at three different temperatures (298, 550, and 650 K). Acetone, propanal, propene, and isopropyl acetate have been identified as major reaction products. Acetone is the main primary product. Theoretical calculations using the composite CBS‐QB3 method provided useful tools to validate the postulated reaction mechanisms leading to experimentally observed species. The formation of other species is also discussed.

Giustini, Andrea↗

Composition-Nanoarchitecture-Performance Analysis of High Energy Density Electrodeposited Silicon for Lithium-Ion Battery Anodes

Electrodeposition of silicon (Si) was previously demonstrated as a promising method for fabricating 3D-structured lithium-ion battery anodes. However, the relationship between the electrochemical performance and chemical composition of the relatively impure electrodeposited silicon is not well understood. Here, we report the electrodeposition of a Si-dominant active material (EDEP-Si) onto 3D-structured nickel (Ni) scaffolds and systematically compare the electrochemical properties, elemental composition, atomistic Si coordination, and molecular structure of EDEP-Si with high-purity amorphous Si grown via static chemical vapor deposition. Despite the considerable amount of carbon (9–11 at %) and oxygen (42–44 at %) present in EDEP-Si, the cycling stability and high reversible specific capacity are remarkably similar to those of CVD-Si on a silicon basis (~2400 mA h/g-Si after 100 cycles). The primary difference is that EDEP-Si exhibits reduced cycling efficiency over the first 10–20 cycles. Reactions between carbon and, more importantly, oxygen in EDEP-Si with lithium are likely responsible for the reduced early cycle performance and lower capacity of the total deposit. Finally, our observations suggest ultrapure Si is not necessary for high electrochemical access to reversible charge storage, although limiting the presence of incorporated impurity species would improve energy density and first cycle efficiency.

25 ENERGY STORAGE↗