Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “consistent scheme”

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 775 records · Page 43

A Computationally Efficient Algorithm for Computing Convex Hull Prices

Electricity markets worldwide allow participants to bid non-convex production offers. While non-convex offers can more accurately reflect a resource's capabilities, they create challenges for market clearing processes. For example, system operators may execute side payments when a participant’s cost is not covered through energy sale settlements from locational marginal pricing schemes, or when a participant incurs lost opportunity costs to follow the dispatch signal. Convex hull pricing minimizes these and other types of side payments while providing uniform (i.e., locationally and temporally consistent) prices. However, computing convex hull prices involves solving either a large-scale linear program - which in turn requires explicit descriptions of market participants’ convex hulls -or the Lagrangian dual of the corresponding non-convex scheduling problem. Here, we propose a computationally feasible and industrially scalable Benders decomposition approach to computing convex hull prices at least an order of magnitude faster than the current state-of-the-art while leveraging recent advances in convex hull formulations for thermal generating units.

61 RADIATION PROTECTION AND DOSIMETRY↗

Multiresolution Wavelet Based Adaptive Numerical Dissipation Control for Shock-Turbulence Computations

The recently developed essentially fourth-order or higher low dissipative shock-capturing scheme of Yee, Sandham and Djomehri (1999) aimed at minimizing nu- merical dissipations for high speed compressible viscous flows containing shocks, shears and turbulence. To detect non smooth behavior and control the amount of numerical dissipation to be added, Yee et al. employed an artificial compression method (ACM) of Harten (1978) but utilize it in an entirely different context than Harten originally intended. The ACM sensor consists of two tuning parameters and is highly physical problem dependent. To minimize the tuning of parameters and physical problem dependence, new sensors with improved detection properties are proposed. The new sensors are derived from utilizing appropriate non-orthogonal wavelet basis functions and they can be used to completely switch to the extra numerical dissipation outside shock layers. The non-dissipative spatial base scheme of arbitrarily high order of accuracy can be maintained without compromising its stability at all parts of the domain where the solution is smooth. Two types of redundant non-orthogonal wavelet basis functions are considered. One is the B-spline wavelet (Mallat & Zhong 1992) used by Gerritsen and Olsson (1996) in an adaptive mesh refinement method, to determine regions where re nement should be done. The other is the modification of the multiresolution method of Harten (1995) by converting it to a new, redundant, non-orthogonal wavelet. The wavelet sensor is then obtained by computing the estimated Lipschitz exponent of a chosen physical quantity (or vector) to be sensed on a chosen wavelet basis function. Both wavelet sensors can be viewed as dual purpose adaptive methods leading to dynamic numerical dissipation control and improved grid adaptation indicators. Consequently, they are useful not only for shock-turbulence computations but also for computational aeroacoustics and numerical combustion. In addition, these sensors are scheme independent and can be stand alone options for numerical algorithm other than the Yee et al. scheme.

Sjoegreen, B.↗

Conceptual Design of a Thermal–to-14-MeV Neutron Conversion Device for Use in the Advanced Test Reactor

Here, this paper details the conceptual design of a thermal-to-14-MeV neutron converter consisting of a mixture of lithium and deuterium in a blanket material. Such a device operates a two-step reaction, first generating tritons via thermal neutron absorption in the tritium breeding material, and in the second step, high-energy neutrons are produced either via deuterium-tritium fusion reaction or with tritium reacting with lithium. A thermal-to-14-MeV neutron converter significantly hardens the neutron spectrum by virtually removing thermal neutrons and adding a high-energy 14 MeV component to the neutron spectrum. While similar concepts have been previously proposed and tested in other reactors, the unique characteristics of the Advanced Test Reactor (ATR), namely its important thermal flux (up to $10^{15}$ n.cm$^{-2}$.s$^{-1}$) makes it markedly attractive for obtaining a very large fast neutron flux, usable for irradiation studies under neutron flux conditions prototypical of fusion reactors. The paper provides a description of a new computational scheme developed for handling the coupled neutron-triton transport mechanism using the Geant4 toolkit. Resulting neutron spectra and high-energy neutron yields are summarized for different irradiation positions and potential neutron breeder materials. Maximum predicted thermal-to-14-MeV neutron yields are on the order of $2.10^{-4}$, which is consistent with previous studies found in the literature. Thus, when placed inside the ATR, such a neutron converter will be providing the largest high-energy neutron source available for activation and irradiation studies of materials foreseen for use in fusion reactors. Future steps will involve qualifying the computational scheme using the ATR critical facility using activation foil measurements.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Transported PDF Modeling of Compressible Turbulent Reactive Flows by using the Eulerian Monte Carlo Fields Method

Although the transported probability density function (PDF) method has been developed for decades, its application has been mainly focused on the low-Mach number flow problems. This work extends the transported PDF method to compressible flow problems. The Eulerian Monte Carlo fields (EMCF) solution method is employed to solve the transported PDF equation for compressible flow problems. A pseudo stagnation enthalpy is introduced and its stochastic partial differential equation is derived to ensure total energy conservation numerically. A new mixing model called interaction by partial exchange with mean (IPEM) is introduced to expand the available choices of mixing models for the EMCF method. The consistency of the EMCF method is examined for solving the transported PDF equation. Numerical implementation details are discussed, such as the density coupling between the compressible flow solver and the EMCF solver, discretization schemes for the mixing terms and the stochastic terms. The implemented compressible flow solver coupled with the EMCF solver is verified and validated in a series of test cases with increasing level of complexity, ranging from a statistically one-dimensional turbulent mixing layer to a self-excited resonance model rocket combustor. It is observed that in general with the increase of compressibility, there is an increase in the sensitivity of the modeling results to the different models and algorithms. This makes it necessary to develop a thorough understanding of the model sensitivity in order to develop a robust and accurate simulation solver for highly compressible turbulent reactive flows. The thermo-acoustic instability inside the model rocket combustor case is captured reasonably, which demonstrates the overall capability of the developed compressible turbulent combustion solver based on the transported PDF method.

42 ENGINEERING↗

Generation of entangled-photons by a quantum dot cascade source in polarized cavities: Using cavity resonances to boost signals and preserve the entanglements

Motivated by recent advances in the development of single photon emitters for quantum information sciences, here we design and formulate a quantum cascade model that describes cascade emission by a quantum dot (QD) in a cavity structure while preserving entanglement that stores information needed for single photon emission. The theoretical approach is based on a photonic structure that consists of two orthogonal cavities in which resonance with either the first or second of the two emitted photons is possible, leading to amplification and rerouting of the entangled light. The cavity–QD scheme uses a four-level cascade emitter that involves three levels for each polarization, leading to two spatially entangled photons for each polarization. By solving the Schrodinger equation, we identify the characteristic properties of the system, which can be used in conjunction with optimization techniques to achieve the “best” design relative to a set of prioritized criteria or constraints in our optical system. The theoretical investigations include an analysis of emission spectra in addition to the joint spectral density profile, and the results demonstrate the ability of the cavities to act as frequency filters for the photons that make up the entanglements and to modify entanglement properties. The results provide new opportunities for the experimental design and engineering of on-demand single photon sources.

Chemistry↗

Perturbatively including inhomogeneities in axion inflation

Axion inflation, i.e. an axion-like inflaton coupled to an Abelian gauge field through a Chern-Simons interaction, comes with a rich and testable phenomenology. This is particularly true in the strong backreaction regime, where the gauge field production heavily impacts the axion dynamics. Lattice simulations have recently demonstrated the importance of accounting for inhomogeneities of the axion field in this regime. We propose a perturbative scheme to account for these inhomogeneities while maintaining high computational efficiency. Our goal is to accurately capture deviations from the homogeneous axion field approximation within the perturbative regime as well as self-consistently determine the onset of the nonperturbative regime.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: cosmological implications from multitracer BAO analysis with galaxies and voids

We construct cosmic void catalogues with the dive void finder upon SDSS BOSS DR12 and eBOSS DR16 galaxy samples with BAO reconstruction, and perform a joint BAO analysis using different types of galaxies and the corresponding voids. The BAO peak is evident for the galaxy–galaxy, galaxy–void, and void–void correlation functions of all data sets, including the ones cross-correlating LRG and ELG samples. Two multitracer BAO fitting schemes are tested, one combining the galaxy and void correlation functions with a weight applied to voids, and the other using a single BAO dilation parameter for all clustering measurements. Both methods produce consistent results with mock catalogues, and on average ~10 per cent improvements of the BAO statistical uncertainties are observed for all samples, compared to the results from galaxies alone. By combining the clustering of galaxies and voids, the uncertainties of BAO measurements from the SDSS data are reduced by 5–15 per cent, yielding 0.9 per cent, 0.8 per cent, 1.1 per cent, 2.3 per cent, and 2.9 per cent constraints on the distance $D_{_{\rm V}}(z)$, at effective redshifts 0.38, 0.51, 0.70, 0.77, and 0.85, respectively. When combined with BAO measurements from SDSS MGS, QSO, and Lyα samples, as well as the BBN results, we obtain $H_0 = 67.58 \pm 0.91\, {\rm km}\, {\rm s}^{-1}\, {\rm Mpc}^{-1}$, Ω m = 0.290 ± 0.015, and $\Omega _\Lambda h^2 = 0.3241 \pm 0.0079$ in the flat-ΛCDM framework, where the 1σ uncertainties are around 6 per cent, 6 per cent, and 17 per cent smaller respectively, compared to constraints from the corresponding anisotropic BAO measurements without voids and LRG–ELG cross-correlations.

79 ASTRONOMY AND ASTROPHYSICS↗

Generalized parton distributions from symbolic regression

AI/ML informed Symbolic Regression is the next stage of scientific modeling. We utilize a highly customizable symbolic regression package “PySR” to model the x and t dependence of the flavor isovector combination Hu-d(x, t, ) at = 0. These PySR models were trained on GPD results provided by both Lattice QCD and phenomenological sources GGL, GK, and VGG. We demonstrate, for the first time, the consistency and systematic convergence of Symbolic Regression by quantifying the disparate models through their Taylor expansion coefficients. In addition to PySR penalizing models with higher complexity and mean-squared error, we implement schemes that test specific physics hypotheses, including force-factorized x and t dependence and Regge behavior in PySR GPDs. We show that PySR can identify factorizing GPD sources based on their response to the Force-Factorized model. Here, knowing the precise behavior of the GPDs, and their uncertainties in a wide range in x and t, crucially impacts our ability to concretely and quantitatively predict hadronic spatial distributions and their derived quantities.

Dotson, Andrew [New Mexico State University, Las C↗

Coupled rotor and fuselage equations of motion

The governing equations of motion of a helicopter rotor coupled to a rigid body fuselage are derived. A consistent formulation is used to derive nonlinear periodic coefficient equations of motion which are used to study coupled rotor/fuselage dynamics in forward flight. Rotor/fuselage coupling is documented and the importance of an ordering scheme in deriving nonlinear equations of motion is reviewed. The nature of the final equations and the use of multiblade coordinates are discussed.

Warmbrodt, W.↗

Stress and Fracture Analyses Under Elastic-plastic and Creep Conditions: Some Basic Developments and Computational Approaches

A new hybrid-stress finite element algorith, suitable for analyses of large quasi-static deformations of inelastic solids, is presented. Principal variables in the formulation are the nominal stress-rate and spin. A such, a consistent reformulation of the constitutive equation is necessary, and is discussed. The finite element equations give rise to an initial value problem. Time integration has been accomplished by Euler and Runge-Kutta schemes and the superior accuracy of the higher order schemes is noted. In the course of integration of stress in time, it has been demonstrated that classical schemes such as Euler's and Runge-Kutta may lead to strong frame-dependence. As a remedy, modified integration schemes are proposed and the potential of the new schemes for suppressing frame dependence of numerically integrated stress is demonstrated. The topic of the development of valid creep fracture criteria is also addressed.

Reed, K. W.↗

Combination probe for hi-frequency unsteady aerodynamic measurements

A combination probe for time-resolved measurements of unsteady compressible flows in transonic wind tunnels is described. The probe measures stagnation (total) temperature and pressure, static pressure, and flow angles in two planes. From these, the fluctuating mass flux, Mach number, and velocity as well as their components in three directions can be deduced. The combination probe consists of a dual hot-wire aspirating temperature and pressure probe mounted piggyback with a high-frequency angle probe. The angle probe has four surface-mounted silicon pressure sensors. A scheme is described for retrieving from the four pressure signals the stagnation and static pressures, Mach number, and flow angles in two planes. The calibrations forming the base for this procedure, obtained from steady-state tests, are given. Typical data obtained in the Karman vortex street shed from a cylinder and at the exit of a Mach 0.4 air jet are presented.

Ng, W. F.↗

A combination probe for high-frequency unsteady aerodynamic measurements in transonic wind tunnels

A combination probe for time-resolved measurements for unsteady compressible flows is described. The probe measures stagnation (total) temperature and pressure, static pressure, and flow angles in two planes. From these, the fluctuating mass flux, Mach number, and velocity as well as their components in three directions can be deduced. The combination probe consists of a dual hot-wire aspirating temperature and pressure probe mounted piggyback with a high-frequency angle probe. The angle probe has four surface-mounted silicon pressure sensors. A scheme is described for retrieving from the four pressure signals the stagnation and static pressures, Mach number, and flow angles in two planes. The calibrations forming the basis for this procedure, obtained from steady-state tests, are given. Typical data obtained in the Karman vortex street shed from a cylinder and at the exit of a Mach-0.4 air jet are presented.

Ng, Wing F.↗

Force control of a multi-arm robot system

A force-compensated control method for multiple manipulators is presented that allows coordinated manipulation of a jointly grasped object. In this scheme, each arm independently carries out the motions required to realize the desired motion of a prescribed point on the manipulated object. The approach has been implemented and demonstrated on a laboratory system consisting of two industrial, computer-controller manipulators.

Alberts, Thomas E.↗

PARAS program: Phased array radio astronomy from space

An orbiting radio telescope is proposed which, when operated in a Very Long Baseline Interferometry (VLBLI) scheme, would allow higher (than currently available) angular resolution and dynamic range in the maps, and the ability of observing rapidly changing astronomical sources. Using a passive phases array technology, the proposed design consists of 656 hexagonal modules forming a 150 meter diameter dish. Each observatory module is largely autonomous, having its own photovoltaic power supply and low-noise receiver and processor for phase shifting. The signals received by the modules are channeled via fiber optics to the central control computer in the central bus module. After processing and multiplexing, the data is transmitted to telemetry stations on the ground. The truss frame supporting each observatory pane is a hybrid structure consisting of a bottom graphite/epoxy tubular triangle and rigidized inflatable Kevlar tubes connecting the top observatory panel and bottom triangle. Attitude control and stationkeeping functions are performed by a system of momentum wheels in the bus and four propulsion modules located at the compass points on the periphery of the observatory dish. Each propulsion module has four monopropellant thrusters and six hydrazine arcjets, the latter supported by a nuclear reactor. The total mass of the spacecraft is 22,060 kg.

Jakubowski, Antoni K.↗

On controlling nonlinear dissipation in high order filter methods for ideal and non-ideal MHD

The newly developed adaptive numerical dissipation control in spatially high order filter schemes for the compressible Euler and Navier-Stokes equations has been recently extended to the ideal and non-ideal magnetohydrodynamics (MHD) equations. These filter schemes are applicable to complex unsteady MHD high-speed shock/shear/turbulence problems. They also provide a natural and efficient way for the minimization of Div(B) numerical error. The adaptive numerical dissipation mechanism consists of automatic detection of different flow features as distinct sensors to signal the appropriate type and amount of numerical dissipation/filter where needed and leave the rest of the region free from numerical dissipation contamination. The numerical dissipation considered consists of high order linear dissipation for the suppression of high frequency oscillation and the nonlinear dissipative portion of high-resolution shock-capturing methods for discontinuity capturing. The applicable nonlinear dissipative portion of high-resolution shock-capturing methods is very general. The objective of this paper is to investigate the performance of three commonly used types of nonlinear numerical dissipation for both the ideal and non-ideal MHD.

Yee, H. C.↗

Adjoint-Based Algorithms for Adaptation and Design Optimizations on Unstructured Grids

Schemes based on discrete adjoint algorithms present several exciting opportunities for significantly advancing the current state of the art in computational fluid dynamics. Such methods provide an extremely efficient means for obtaining discretely consistent sensitivity information for hundreds of design variables, opening the door to rigorous, automated design optimization of complex aerospace configuration using the Navier-Stokes equation. Moreover, the discrete adjoint formulation provides a mathematically rigorous foundation for mesh adaptation and systematic reduction of spatial discretization error. Error estimates are also an inherent by-product of an adjoint-based approach, valuable information that is virtually non-existent in today's large-scale CFD simulations. An overview of the adjoint-based algorithm work at NASA Langley Research Center is presented, with examples demonstrating the potential impact on complex computational problems related to design optimization as well as mesh adaptation.

Nielsen, Eric J.↗

Multiple Embedded Processors for Fault-Tolerant Computing

A fault-tolerant computer architecture has been conceived in an effort to reduce vulnerability to single-event upsets (spurious bit flips caused by impingement of energetic ionizing particles or photons). As in some prior fault-tolerant architectures, the redundancy needed for fault tolerance is obtained by use of multiple processors in one computer. Unlike prior architectures, the multiple processors are embedded in a single field-programmable gate array (FPGA). What makes this new approach practical is the recent commercial availability of FPGAs that are capable of having multiple embedded processors. A working prototype (see figure) consists of two embedded IBM PowerPC 405 processor cores and a comparator built on a Xilinx Virtex-II Pro FPGA. This relatively simple instantiation of the architecture implements an error-detection scheme. A planned future version, incorporating four processors and two comparators, would correct some errors in addition to detecting them.

Bolotin, Gary↗

The Baryon Mapping Experiment (BMX), a 21cm intensity mapping pathfinder

The Baryon Mapping eXperiment (BMX) is an interferometric array designed as a pathfinder for a future post- reionization 21 cm intensity mapping survey. It consists of four 4-meter parabolic reflectors each having offset pyramidal horn feed, quad-ridge orthomode transducer, temperature-stabilized RF amplification and filtering, and pulsed noise injection diode. An undersampling readout scheme uses 8-bit digitizers running at 1.1 Gsamples/sec to provide access to signals from 1.1 - 1.55 GHz (third Nyquist zone), corresponding to HI emission from sources at redshift 0 < z < 0.3. An FX correlator is implemented in GPU and generates 28 GB/day of time-ordered visibility data. About 7,000 hours of data were collected from Jan. 2019 - May 2020, and we will present results on system performance including sensitivity, beam mapping studies, observations of bright celestial targets, and system electronics upgrades. BMX is a pathfinder for the proposed PUMA intensity mapping survey in the 2030s.

79 ASTRONOMY AND ASTROPHYSICS↗