Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Operator approximation”

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 829 records · Page 46

What Does Reflection from Cloud Sides Tell Us About Vertical Distribution of Cloud Droplet Sizes?

Cloud development, the onset of precipitation and the effect of aerosol on clouds depend on the structure of the cloud profiles of droplet size and phase. Aircraft measurements of cloud profiles are limited in their temporal and spatial extent. Satellites were used to observe cloud tops not cloud profiles with vertical profiles of precipitation-sized droplets anticipated from Cloudsat. The recently proposed CLAIM-3D satellite mission (cloud aerosol interaction mission in 3D) suggests to measure profiles of cloud microphysical properties by retrieving them from the solar and infrared radiation reflected or emitted from cloud sides. Inversion of measurements from the cloud sides requires rigorous understanding of the 3-dimensional (3D) properties of clouds. Here we discuss the reflected sunlight from the cloud sides and top at two wavelengths: one nonabsorbing to solar radiation (0.67 micrometers) and one with liquid water efficient absorption of solar radiation (2.1 micrometers). In contrast to the plane-parallel approximation, a conventional approach to all current operational retrievals, 3D radiative transfer is used for interpreting the observed reflectances. General properties of the radiation reflected from the sides of an isolated cloud are discussed. As a proof of concept, the paper shows a few examples of radiation reflected from cloud fields generated by a simple stochastic cloud model with the prescribed vertically resolved microphysics. To retrieve the information about droplet sizes, we propose to use the probability density function of the droplet size distribution and its first two moments instead of the assumption about fixed values of the droplet effective radius. The retrieval algorithm is based on the Bayesian theorem that combines prior information about cloud structure and microphysics with radiative transfer calculations.

Marshak, Alexander↗

Design Enhancements of the Fourier Kelvin Stellar Interferometer to Enable Detection of Earth Twins

During the last few years, considerable effort has been directed towards very large-scale (> $5 billion) missions to detect and characterize Mars-radius to Earth-radius planets around nearby stars; such as the Terrestrial Planet Finder Interferometer and Darwin missions. However, technological issues such as formation flying and control of systematic noise sources will likely prevent these missions from entering Phase A until at least the end of the next decade. Presently more than 350 planets have been discovered by a variety of techniques, and little is known about the majority of them other than their approximate mass. However, a simplified nulling interferometer operating in the near- to mid-infrared (e.g. approx. 5-15 microns), like the enhanced version of the Fourier Kelvin Stellar Interferometer (FKSI), can characterize the atmospheres of a large sample of the known planets - including Earth twins. Many other scientific problems can be addressed with a system like FKSI, including the studies of debris disks, active galactic nuclei, and low mass companions around nearby stars. We report results of a recent engineering study on an enhanced version of FKSI that includes 1-meter primary mirrors, 20-meter boom length, and an advanced sun shield that will provide a 45-degree FOR and 40K operating temperature for all optics including siderostats.

Barry, Richard K.↗

Small, Fast TES Microcalorimeters with Unprecedented X-ray Spectral Performance

Driven initially by the desire for X-ray microcalorimeter arrays suitable for imaging the dynamic solar corona, we have developed a transition-edge-sensor (TES) microcalorimeter optimization that exhibits a unique combination of high spectral resolving power and a wide X-ray bandpass. These devices have achieved spectral performance of dE approximately 1.3 eV FWHM at 1.5 keV, 1.6 eV at 6 keV, and 2.0 eV at 8 keV, using small TESs (e.g., approximately 35 micron x 35 micron) that operate in a regime in which the superconducting transition is highly current dependent. In order to accommodate high X-ray count rates, the devices sit directly on a solid substrate instead of on membranes, and we use an embedded heatsinking layer to reduce pixel-to-pixel crosstalk. We will present results from devices with a range of TES and absorber sizes, and from device wafers with varied embedded heatsink materials. This contribution will focus on count-rate capabilities, including a discussion of the trade-off between count rate and energy resolution, and the heatsinking design. We will also present preliminary tests of array readout using a code-division multiplexed SQUID readout scheme, which may be necessary to enable large arrays of these fast devices.

Eckart, M. E.↗

High-Speed Imaging of the First Kink Mode Instability in a Magnetoplasmadynamic Thruster

One of the biggest challenges to efficient magnetoplasmadynamic thruster (MPDT) operation is the onset of high-frequency voltage oscillations as the discharge current is increased above a threshold value. The onset regime is closely related to magnetohydrodynamic instabilities known as kink modes. This work documents direct observation of the formation and quasi-steady state behavior of an argon discharge plasma in a MPDT operating at discharge currents of 8 to 10 kA for a pulse length of approximately 4 ms. A high-speed camera images the quasi-steady-state operation of the thruster at 26,143 fps with a frame exposure time of 10 micro s. A 0.9 neutral density filter and 488-nm argon line filter with a 10-nm bandwidth are used on separate trials to capture the time evolution of the discharge plasma. Frame-by-frame analysis of the power flux incident on the CCD sensor shows both the initial discharge plasma formation process and the steady-state behavior of the discharge plasma. Light intensity levels on the order of 4-6 W/m2 indicate radial and azimuthal asymmetries in the concentration of argon plasma in the discharge channel. The plasma concentration exhibits characteristics that suggest the presence of a helical plasma column. This helical behavior has been observed in previous experiments that characterize plasma kink mode instabilities indirectly. Therefore, the direct imaging of these plasma kink modes further supports the link between MPDT onset behavior and the excitation of the magnetohydrodynamic instabilities.

Walker, Jonathan A.↗

Maneuver Design for the Juno Mission: Inner Cruise

The Juno spacecraft launched in August 2011 and, following a successful Earth flyby in October 2013, is on course for a nominal orbit insertion at Jupiter in July 2016. This paper examines the design and execution of deterministic and statistical trajectory correction maneuvers during the first approximately 27 months of post-launch operations that defined the "Inner Cruise" phase of the Juno mission. Topics of emphasis include the two deep space maneuvers, Earth flyby altitude biasing strategy, and the sequence of trajectory correction maneuvers executed in the weeks prior to the successful Earth gravity assist.

Pavlak, Thomas A.↗

Exploring Network-Related Optimization Problems Using Quantum Heuristics

Network-related connectivity optimization problems are underlying a wide range of applications and are also of high computational complexity. We consider studying network optimization problems using two types of quantum heuristics.One is quantum annealing, and the other Quantum Alternating Operator Ansatz, an extension of the Quantum Approximate Optimization Algorithms for gate-model quantum computation, in which a cost-function based unitary and a non-commuting mixing unitary are applied alternately. We present problem mappings for problems of finding the spanning-tree or spanning-graph of a graph that optimizes certain costs, and a variant that further requires the spanning-tree be degree-bounded. With quantum annealing, all constraints are cast into penalty terms in the cost Hamiltonian, and the solution is encoded as the ground state of the Hamiltonian. We provide three mappings to the quadratic unconstrained binary optimization (QUBO) form, compare the resource requirements, and analyze the tradeoffs. For QAOA, we give special focus on the design of mixers based on the constraints presented in the problem, such that the system evolution remains in a subspace of the full Hilbert space where all constraints are satisfied. In the spanning-tree problem, one such hard constraint is that a mixer applied to a spanning-tree needs also be a spanning tree. This involves checking the connectivity of a subgraph, which is a global condition common for most network-related problems. We show how this feature can be efficiently represented in the mixer in a quantum coherent way, based on manipulation of a descendant-matrix and an adjacent matrix. We further develop a mixer for the spanning-graphs based on the spanning-tree mixer.

Wang, Zhihui↗

Design of an Automated Ultrasonic Scanning System for In-Situ Composite Cure Monitoring and Defect Detection

The preliminary design and development of an automated ultrasonic scanning system for in-situ composite cure monitoring and defect detection in the high temperature environment of an oven was completed. This preliminary design is a stepping stone to deployment in the high temperature and high pressure environment of an autoclave, the primary cure method of aerospace grade thermoset composites. Cure monitoring with real-time defect detection during the process could determine when defects form and how they move. In addition, real-time defect detection during cure could assist validating physics-based process models for predicting defects at all stages of the cure cycle. A physics-based process model for predicting porosity and fiber waviness originating during cure is currently under development by the NASA Advanced Composites Project (ACP). For the design, an ultrasonic contact scanner is enclosed in an insulating box that is placed inside an oven during cure. Throughout the cure cycle, the box is nitrogen-cooled to approximately room temperature to maintain a standard operating environment for the scanner. The composite part is mounted on the outside of the box in a vacuum bag on the build/tool plate. The build plate is attached to the bottom surface of the box. The scanner inspects the composite panel through the build plate, tracking the movement of defects introduced during layup and searching for new defects that may form during cure. The focus of this paper is the evaluation and selection of the build plate material and thickness. The selection was based on the required operating temperature of the scanner, the cure temperature of the composite material, thermal conductivity models of the candidate build plates, and a series of ultrasonic attenuation tests. This analysis led to the determination that a 63.5 mm thick build plate of borosilicate glass would be utilized for the system. The borosilicate glass plate was selected as the build plate material due to the low ultrasonic attenuation it demonstrated, its ability to efficiently insulate the scanner while supporting an elevated temperature on the part side of the plate, and the availability of a 63.5 mm thick plate without the need for lamination.

Hudson, Tyler B.↗

Study network-related optimization problems using quantum alternating optimization ansatz

Network-related connectivity optimization problems are underlying a wide range of applications and are also of high computational complexity. We consider studying network optimization problems using two types of quantum heuristics. One is quantum annealing, and the other Quantum Alternating Operator Ansatz, an extension of the Quantum Approximate Optimization Algorithms for gate-model quantum computation, in which a cost-function based unitary and a non-commuting mixing unitary are applied alternately. We present problem mappings for problems of finding the spanning-tree or spanning-graph of a graph that optimizes certain costs, and a variant that further requires the spanning-tree be degree-bounded. With quantum annealing, all constraints are cast into penalty terms in the cost Hamiltonian, and the solution is encoded as the ground state of the Hamiltonian. We provide three mappings to the quadratic unconstrained binary optimization (QUBO) form, compare the resource requirements, and analyze the tradeoffs. For QAOA, we give special focus on the design of mixers based on the constraints presented in the problem, such that the system evolution remains in a subspace of the full Hilbert space where all constraints are satisfied. In the spanning-tree problem, one such hard constraint is that a mixer applied to a spanning-tree needs also be a spanning tree. This involves checking the connectivity of a subgraph, which is a global condition common for most network-related problems. We show how this feature can be efficiently represented in the mixer in a quantum coherent way, based on manipulation of a descendant-matrix and an adjacent matrix. We further develop a mixer for the spanning-graphs based on the spanning-tree mixer.

Zhihui Wang↗

Recent Concept Study for Cryogenic Fluid Management to Support Opposition Class Crewed Missions to Mars

NASA recently completed a mission concept study to evaluate the feasibility and propulsion technology development requirements for reduced travel duration crewed missions to Mars. A high-level goal of the study was to minimize the health impact on the crew caused by the space environment. This was implemented in the study by limiting the crew to a total of approximately-two years of in-space operations and travel time. For the initial mission, the crew would stay about 30 days on the Martian surface. The propulsive demands of such a mission are immense, and the study identified two advanced propulsion options with the potential to meet the mission requirements—both options rely on nuclear fission to provide efficient propulsive energy. One propulsion option was a nuclear electric propulsion (NEP)/Chem Hybrid, with a reactor and energy conversion system powering xenon propellant ion thrusters to provide an efficient, but lower-thrust, push for most of the mission duration. This concept also relied on a liquid oxygen/liquid methane (LO2/LCH4) chemical propulsion stage to provide high thrust for maneuvers while near the Earth and Mars. The second propulsion option was nuclear thermal propulsion (NTP), in which the reactor heats liquid hydrogen (LH2) propellant to expand through a nozzle for thrust at about twice the efficiency of the best chemical propulsion systems. Both vehicle concepts rely on storing large amounts of cryogenic propellant (either LO2/LCH4 or LH2) for multiple years in space without loss, far exceeding state-of-the-art capability. To enable this new capability, the team assumed the use of several advanced cryogenic fluid management (CFM) technologies and analyzed the integrated system performance. This included considering the vehicle-level effects of the size, mass, and power requirements of these CFM elements. Further, the team evaluated the development required to enable such a mission in the mid-2030 s and determined that it was feasible. The paper elaborates on the assumed CFM technologies, provides key analysis results, and illustrates the feasibility of technology development for the proposed solutions to the CFM challenges for each propulsion concept.

cryogenic propellant↗

On-Orbit Performance of the Adiabatic Demagnetization Refrigerator on XRISM

The X-ray Imaging and Spectroscopy Mission (XRISM) observatory was launched on Sep 7, 2023 from Tanegashima Space Center in Japan. Resolve, one its two instruments, performs high-resolution spectroscopy in the soft x-ray band (0.2-13 keV) using a 6x6 microcalorimeter array. The array cooled to 50 mK by a 3-stage ADR, which is linked to both a liquid helium dewar (at <1.2 K) and a Joule-Thomson (JT) cryocooler operating at <4.5 K. While liquid helium is present, two of the ADR stages provide detector cooling at 50 mK and auxiliary cooling to 0.5 K to intercept some parasitic heat loads, while rejecting waste heat to the helium. Once the helium is exhausted, Resolve enters a cryogen-free mode in which all three ADR stages are operated to provide detector cooling at 50 mK and to continuously cool the helium tank to 1.4 K. In this mode, waste heat is rejected to the JT cryocooler. Cryogen-free operation can be sustained as long as the JT and other cryocoolers remain fully operational. At launch, the helium tank contained approximately 35.6 liters of liquid. Within a few days, the helium cooled below 1.15 K, with an estimated 35.0 liters remaining. With an expected time average heat load of 0.68 mW, the helium lifetime was projected to exceed 4 years, but measurements of He volume on orbit suggest a significantly longer lifetime. Details of the ADR’s design and performance are presented.

x-ray astromony↗

A Flexible Power Method for Solving Infinite Dimensional Tensor Eigenvalue Problems

We propose a flexible power method for computing the leftmost, i.e., algebraically smallest, eigenvalue of an infinite dimensional tensor eigenvalue problem, $H x = \lambda x$, where the infinite dimensional symmetric matrix $H$ exhibits a translational invariant structure. We assume the smallest eigenvalue of $H$ is simple and apply a power iteration of $e^{-H}$ with the eigenvector represented in a compact way as a translational invariant infinite Tensor Ring (iTR). Hence, the infinite dimensional eigenvector can be represented by a finite number of iTR cores of finite rank. In order to implement this power iteration, we use a small parameter $t$ so that the infinite matrix-vector operation $e^{-Ht}x$ can efficiently be approximated by the Lie product formula, also known as Suzuki--Trotter splitting, and we employ a low rank approximation through a truncated singular value decomposition on the iTR cores in order to keep the cost of subsequent power iterations bounded. We also use an efficient way for computing the iTR Rayleigh quotient and introduce a finite size iTR residual which is used to monitor the convergence of the Rayleigh quotient and to modify the timestep $t$. In this paper, we discuss 2 different implementations of the flexible power algorithm and illustrate the automatic timestep adaption approach for several numerical examples.

Beeumen, Roel Van↗

Computational Math Problems for a Clean Energy Future

Cutting edge computational mathematics are ubiquitous in renewable energy research. Problems in resilient and reliable electric grid operations, infrastructure planning, wind farm yaw control, and more demand sophisticated and scalable computational tools that enable the transition of renewable energy technologies from proof of concept to deployment into our energy system. The mission of the Computational Science Center at NREL is to lead the lab's efforts to solve energy challenges using high-performance computing (HPC), computational science, applied mathematics, scientific data management, visualization, and informatics. In this poster, we provide a short overview of three areas of computational mathematics research at NREL: wind power scenario generation for stochastic grid operations and infrastructure planning, improved rational function approximations for electromagnetic transients codes, and wind farm yaw control using a combination of the Alternating Direction Method of Multipliers (ADMM) and reinforcement learning (RL). Increasing penetrations of renewable energy into power grids motivate the investigation of new approaches to characterizing uncertainty for five-minute economic dispatch problems. Similarly, as the penetration of distributed energy resources on power grids increases, it becomes important to revisit our methods of modelling transient phenomena, i.e. electromagnetic transients programs. Finally, the combination of ADMM and RL for wind farm yaw control presented here can potentially increase the efficiency of the deployed distributed controllers by orders of magnitude.

ADMM↗

Data-Centric Approach to Capture Non-Polynomial Nonlinear Dynamics

We propose an analytical construction of observable functions in the extended dynamic mode decomposition (EDMD) algorithm. EDMD is a numerical method for approximating the spectral properties of the Koopman operator. The choice of observable functions is fundamental for applying EDMD to nonlinear problems arising in systems and control. Existing methods either start from a set of dictionary functions and look for the subset that best fits the underlying nonlinear dynamics or rely on machine learning algorithms to “learn” observable functions. Conversely, in this paper, we start from the dynamical system model and lift it through the Lie derivatives, rendering it into a polynomial form. This proposed transformation into a polynomial form is exact and provides an adequate set of observable functions. The strength of the proposed approach is its applicability to a broader class of nonlinear dynamical systems, particularly those with nonpolynomial functions and compositions thereof. Moreover, it retains the physical interpretability of the underlying dynamical system and can be readily integrated into existing numerical libraries. We demonstrate the proposed approach with an application to electric power systems. The modeled system consists of a single generator connected to an infinite bus, where nonlinear terms include sine and cosine functions. The results demonstrate the effectiveness of the proposed procedure in off-attractor nonlinear dynamics for estimation and prediction; the observable functions obtained from the proposed construction outperform methods that use dictionary functions comprising monomials or radial basis functions.

extended dynamic mode decomposition↗

Plasma Potential and Langmuir Probe Measurements in the Near-Field Plume of the NASA-457Mv2 Hall Thruster

In order to further the design of future high-power Hall thrusters and provide experimental validation for ongoing modeling efforts, plasma potential and Langmuir probe measurements were performed on the 50-kW NASA-457Mv2. An electrostatic probe array comprised of a near-field Faraday probe, single Langmuir probe, and emissive probe was used to interrogate the near-field plume from approximately 0.1 ? 2.0 mean thruster diameters downstream of the thruster exit plane at the following operating conditions: 300 V, 400 V and 500 V at 30 kW and 500 V at 50 kW. Results have shown that the acceleration zone is limited to within 0.4 mean thruster diameters of the exit plane while the high-temperature region is limited to 0.25 mean thruster diameters from the exit plane at all four operating conditions. Maximum plasma potentials in the near-field at 300 and 400 V were approximately 50 V with respect to cathode potential, while maximum electron temperatures varied from 24 ? 32 eV, depending on operating condition. Isothermal lines at all operating conditions were found to strongly resemble the magnetic field topology in the high-temperature regions. This distribution was found to create regions of high temperature and low density near the magnetic poles, indicating strong, thick sheath formation along these surfaces. The data taken from this study are considered valuable for future design as well as modeling validation.

Shastry, Rohit↗

Plasma Potential and Langmuir Probe Measurements in the Near-field Plume of the NASA-457Mv2 Hall Thruster

In order to further the design of future high-power Hall thrusters and provide experimental validation for ongoing modeling efforts, plasma potential and Langmuir probe measurements were performed on the 50-kW NASA-457Mv2. An electrostatic probe array comprised of a near-field Faraday probe, single Langmuir probe, and emissive probe was used to interrogate the near-field plume from approximately 0.1 - 2.0 mean thruster diameters downstream of the thruster exit plane at the following operating conditions: 300 V, 400 V and 500 V at 30 kW and 500 V at 50 kW. Results have shown that the acceleration zone is limited to within 0.4 mean thruster diameters of the exit plane while the high-temperature region is limited to 0.25 mean thruster diameters from the exit plane at all four operating conditions. Maximum plasma potentials in the near-field at 300 and 400 V were approximately 50 V with respect to cathode potential, while maximum electron temperatures varied from 24 - 32 eV, depending on operating condition. Isothermal lines at all operating conditions were found to strongly resemble the magnetic field topology in the high-temperature regions. This distribution was found to create regions of high temperature and low density near the magnetic poles, indicating strong, thick sheath formation along these surfaces. The data taken from this study are considered valuable for future design as well as modeling validation.

Shastry, Rohit↗

Demonstration of acoustic monitoring for structural health of microreactors: Through use of neural networks and resonant ultrasound spectroscopy

Nuclear microreactors prioritize modularity and portability and are intended to be a cost-effective technology for non-conventional nuclear markets. As such, the development of microreactors into a safe and feasible solution for energy security applications will necessitate the development of non-destructive technologies to monitor the integrity of inaccessible reactors components during operation. This demonstration applies linear and nonlinear acoustic techniques, in combination with machine learning, to detect and classify mechanical changes (stress and damage) in a test article which are broadly representative of potential operating challenges within a functioning microreactor. All necessary data has been collected for this demonstration, with minor experimental issues identified that can be addressed in follow-on work. Motivated by the expected conditions within a functioning microreactor, we have demonstrated our monitoring techniques on a core-block-like test article using an unstructured excitation source that approximates the noisy acoustic environment expected during reactor operation. At all stress states mechanically applied to the test article, a machine learning model using an artificial neural network was able to classify with 100% accuracy whether a 3D laser vibrometry point measurement was made on an intact or artificially defective test article. Further, model predictions about whether the defect interface was rough or smooth were 95% accurate, indicating the ability of acoustic techniques to recover defect characteristics. Resonant ultrasound spectroscopy (RUS) was also applied to the dataset to provide further quantitative insights about material properties. RUS analysis was ultimately hampered by several minor experimental and data issues, limiting results to certain cases for this demonstration. Last, analysis using nonlinear RUS exhibited sensitivity to changing levels of applied stresses for each intact and defective state. As presented in this demonstration, acoustic monitoring exhibits sensitivity to stress changes, which are of concern due to high thermal gradients expected during startup and operation. Further, our techniques distinguish between measurements made on intact and damaged test articles.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

High heat transfer oxidizer heat exchanger design and analysis

The RL10-2B engine, a derivative of the RL10, is capable of multimode thrust operation. This engine operates at two low thrust levels: tank head idle (THI), which is approximately 1 to 2% of full thrust, and pumped idle (PI), which is 10% of full thrust. Operation at THI provides vehicle propellant settling thrust and efficient engine thermal conditioning; PI operation provides vehicle tank pre-pressurization and maneuver thrust for low-g deployment. Stable combustion of the RL10-2B engine during the low thrust operating modes can be accomplished by using a heat exchanger to supply gaseous oxygen to the propellant injector. The oxidizer heat exchanger (OHE) vaporizes the liquid oxygen using hydrogen as the energy source. The design, concept verification testing and analysis for such a heat exchanger is discussed. The design presented uses a high efficiency compact core to vaporize the oxygen, and in the self-contained unit, attenuates any pressure and flow oscillations which result from unstable boiling in the core. This approach is referred to as the high heat transfer design. An alternative approach which prevents unstable boiling of the oxygen by limiting the heat transfer is referred to as the low heat transfer design and is reported in Pratt & Whitney report FR-19135-2.

Kmiec, Thomas D.↗

Two phase capillary pumped heat transfer in the Instrument Thermal Test Bed

An experimental study of the thermal performance of two evaporators installed in the Instrument Thermal Test Bed (ITTB) was conducted. The ITTB was operated as a capillary pumped loop (CPL) with a transport length of approximately 12 meters. Empirical determinations of a general start up procedure, overall heat transfer coefficient, and minimum operating power were accomplished for each evaporator. Additionally, a detailed thermal model was developed for the High Power Spacecraft Thermal Management (HPSTM) evaporator and validated.

Didion, Jeffrey R.↗