SEARCH · Engineering Papers

Results for “successive linear 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4

All sources NASA NTRS DOE OSTI

NASA NTRS1974

Changes in the performance characteristics of a GaAs near infrared light emitting diode when exposed to various current and thermal stresses

The changes that occurred in the optical and electrical characteristics of a near infrared, GaAs light emitting diode, when operated under various levels and combinations of current and thermal stresses are discussed. A total of forty parts were operated for two thousand hours under eight different sets of dc current and ambient temperature conditions. Degradation in the radiant optical power of these devices was thirty-four percent when operated at their rated current and an ambient temperature of 298K (25 C). Derating the current and/or the thermal stress reduced the degradation of this parameter in approximately a linear manner. All degraded devices behaved similarly, exhibiting rapid nonlinear degradation followed by a gradual linear degradation and finally a period of stable operation. An attempt was made to correlate initial device condition to degradation during stress testing, but met with little success.

Thomas, E. F., Jr.↗

NASA NTRS1994

Sensitivity derivatives for three dimensional supersonic Euler code using incremental iterative strategy

In a recent work, an incremental strategy was proposed to iteratively solve the very large systems of linear equations that are required to obtain quasianalytical sensitivity derivatives from advanced computational fluid dynamics (CFD) codes. The technique was sucessfully demonstrated for two large two-dimensional problems: a subsonic and a transonic airfoil. The principal feature of this incremental iterative stategy is that it allows the use of the identical approximate coefficient matrix operator and algorithm to solve the nonlinear flow and the linear sensitivity equations; at convergence, the accuracy of the sensitivity derivatives is not compromised. This feature allows a comparatively straightforward extension of the methodology to three-dimensional problems; this extension is successfully demonstrated in the present study for a space-marching solution of the three-dimensional Euler equations over a Mach 2.4 blended wing-body configuration.

Korivi, Vamshi Mohan↗

DOE OSTI2021

Success of digital adiabatic simulation with large Trotter step

The simulation of adiabatic evolution has deep connections with adiabatic quantum computation, the quantum approximate optimization algorithm, and adiabatic state preparation. Here we address the error analysis problem in quantum simulation of adiabatic process using Trotter formulas. Here we show that with additional conditions, the circuit depth can be linear in simulation time T. The improvement comes from the observation that the fidelity error here can't be estimated by the norm distance between evolution operators. This phenomenon is termed the robustness of discretization in digital adiabatic simulation. It can be explained in three steps, from analytical and numerical evidence: (1) The fidelity error should be estimated by applying adiabatic theorem on the effective Hamiltonian instead. (2) Because of the specialty of Riemann-Lebesgue lemma, most adiabatic process is naturally robust against discretization. (3) As the Trotter step gets larger, the spectral gap of effective Hamiltonian tends to close, which results in the failure of digital adiabatic simulation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Striking luminescence phenomena of carbon dots and their applications as a double ratiometric fluorescence probes for H2S detection

Here, we report some new observations on the luminescence of carbon dots (CDs) that are passivated with polyethylene glycol, and explored them, for the first time, as highly sensitive H2S detection probe. The as-prepared CDs have an ultraviolet (UV) emission at 350 nm and a green emission at 540 nm when excited at 270 nm, but they have only a green emission at 523 nm when excited at 365 nm. As Na2S is added to the CDs, the UV emission decreases gradually, while the green emission increases slightly, and interestingly a new emission appears at 455 nm is increased linearly in intensity with the increase of Na2S concentration. The blue emission at 455 nm is from the intrinsic core of CDs, and the green emissions at 540 or 523 nm are from their surface states, while the UV emission at 350 nm is from the anchored PEG. Both the emission intensity ratios of F-450/F-350 and F-455/F-523 show an excellent linear relationship with the Na2S concentration in the range of 0-800 mu M, and the detection limit estimated from both of them is approximately 7.0 nM, indicating our method is the most sensitive H2S CDs based detection method reported thus far. Furthermore, the detection was tested successfully for H2S imaging in live cells. We propose that the interactions of H2S with the CDs and the anchored ligands, as well as the energy transfer between the ligands and the CDs, are responsible for the luminescence responses for H2S detection. Our studies enrich the world of carbon dots with plenty of potential applications. (C) 2020 Elsevier Ltd. All rights reserved.

Bioimaging↗

NASA NTRS1989

Projection methods for the numerical solution of Markov chain models

Projection methods for computing stationary probability distributions for Markov chain models are presented. A general projection method is a method which seeks an approximation from a subspace of small dimension to the original problem. Thus, the original matrix problem of size N is approximated by one of dimension m, typically much smaller than N. A particularly successful class of methods based on this principle is that of Krylov subspace methods which utilize subspaces of the form span(v,av,...,A(exp m-1)v). These methods are effective in solving linear systems and eigenvalue problems (Lanczos, Arnoldi,...) as well as nonlinear equations. They can be combined with more traditional iterative methods such as successive overrelaxation, symmetric successive overrelaxation, or with incomplete factorization methods to enhance convergence.

Saad, Youcef↗

DOE OSTI2025

McCormick envelopes in mixed-integer PDE-constrained optimization

McCormick envelopes are a standard tool for deriving convex relaxations of optimization problems that involve polynomial terms. Such McCormick relaxations provide lower bounds, for example, in branch-and-bound procedures for mixed-integer nonlinear programs but have not gained much attention in PDE-constrained optimization so far. This lack of attention may be due to the distributed nature of such problems, which on the one hand leads to infinitely many linear constraints (generally state constraints that may be difficult to handle) in addition to the state equation for a pointwise formulation of the McCormick envelopes and renders bound-tightening procedures that successively improve the resulting convex relaxations computationally intractable. We analyze McCormick envelopes for a model problem class that is governed by a semilinear PDE involving a bilinearity and integrality constraints. We approximate the nonlinearity and in turn the McCormick envelopes by averaging the involved terms over the cells of a partition of the computational domain on which the PDE is defined. This yields convex relaxations that underestimate the original problem up to an a priori error estimate that depends on the mesh size of the discretization. These approximate McCormick relaxations can be improved by means of an optimization-based bound-tightening procedure. We show that their minimizers converge to minimizers to a limit problem with a pointwise formulation of the McCormick envelopes when driving the mesh size to zero. We provide a computational example, for which we certify all of our imposed assumptions. The results point to both the potential of the methodology and the gaps in the research that need to be closed. Our methodology provides a framework first for obtaining pointwise underestimators for nonconvexities and second for approximating them with finitely many linear inequalities in an infinite-dimensional setting.

Approximations and Expansions↗

NASA NTRS1975

Cosmic ray propagation in interplanetary space

The validity of the test-particle picture, the approximation of static fields, and the spatial-diffusion approximation are discussed in a general way before specific technical assumptions are introduced. It is argued that the spatial-diffusion equation for the intensity per unit energy has a much wider range of applicability than the kinetic (Fokker-Planck) equation it is derived from. This gives strong weight to the phenomenological propagation theory. The general success (and possible failure at small energies) of the phenomenological theory for the modulation of galactic cosmic rays and solar events is described. Apparent effects such as the 'free boundary' are given disproportionate weight since they establish the connection with the detailed plasma physics of the solar wind. Greatest attention is paid to the pitch-angle diffusion theory. A general theory is presented which removes the well-known secularities of the quasi-linear approximation. The possible breakdown of any pitch-angle diffusion theory at very small energies is perhaps connected with the observed 'turn up' of the spectrum at low energies. A first attempt to derive the spatial dependence of the diffusion coefficient in the solar cavity, using such a divergence free scattering theory, is described and compared with recent observations out to 5 AU.

Voelk, H. J.↗

NASA NTRS2024

COLUTA: Custom 8-Channel 15-bit 40-MSPS ADC for the ATLAS Liquid Argon Calorimeter Readout

The stringent radiation-tolerance, quantization requirements, and the need for seamless integration of the on-detector readout electronics chain require the design of a full-custom analog-to-digital converter (ADC) for the 182,468 channels in the ATLAS Liquid Argon (LAr) calorimeter at the Large Hadron Collider (LHC) at CERN, as part of the High-Luminosity LHC (HL-LHC) upgrade. Each of the 8 channels in the prototype 65 nm CMOS ADC with 15-bit resolution and >68 dB signal-to-noise-and-distortion ratio (SNDR), or equivalently >11 effective number of bits (ENOB), consists of a Multiplying-DAC (MDAC) and a successive-approximation (SAR) ADC. A digital data processing unit (DDPU) calibrates and serially transmits the quantized data. Initial lab characterization of 18 chips shows a minimum SNDR of 69.5 dB at full-scale at about 5 MHz and a maximum power consumption of 1.17 W; differential non-linearity (DNL) measurements show no missing codes. Additional validation of the ADC is ongoing before approximately 70k chips can be mass-produced for the upgrade.

Rui Xu↗

DOE OSTI2021

Modeling Contact Angle vs. Temperature for the Quartz-Water-Decane System

Innovative approaches are needed to improve the efficiency of oil recovery technologies to meet the growing demands of fossil-fuel based energy consumption. Enhanced oil recovery (EOR) methods such as low-salinity waterflooding and chemically tuned waterflooding aim to optimize the reservoir’s wetting properties, detaching oil globules from rock surfaces and allowing easier oil flow through pore throats. This wetting behavior is commonly quantified by contact angle measurements of the rock-oil-brine interface, which have been thoroughly investigated and theorized for many systems at ambient temperatures and pressures. However, few studies exist for extending contact angle theories away from ambient conditions. In this paper, we model the contact angles of a quartz-water-decane system at elevated temperatures using the surface tension component (STC) approach. Temperature-dependent van der Waals [Lifshitz-van der Waals (LW)] interactions and hydrogen-bonding (acid-base) interactions were calculated and are incorporated into the model for the quartz-water-decane interface. Additionally, the Hough and White procedure was used to create temperature-dependent dielectric functions of quartz, water, and normal decane for calculations of Hamaker coefficients. Hamaker coefficients calculated this way are highly linear with temperature and agree well with Israelachvili’s approximation. The acid-base interactions likely contribute the most to system wettability changes. Resulting contact angles of the quartz-water-decane system shift from water-wet (16°) to slightly water-wet (57.4°) as temperature increases. The model was also successfully verified for the quartz-air-water system. Our results can be used in future studies to determine optimal injected water compositions for specific rock-oil-brine and other systems with consideration of reservoir temperature.

02 PETROLEUM↗

NASA NTRS1981

Recursive algorithms for two-dimensional smoothing using bicubic hermite polynomial

It is noted that in the past, smoothing splines originated from approximation theory have been successfully applied to data filtering and image smoothing problems. Even though the nonrecursive technique of smoothing splines gives an optimal solution, the amount of computation increases rapidly with the size of the two-dimensional data. A derivation is presented here of quarter-plane filtering algorithms that provide smoothed estimates of function values and their derivatives by fitting two-dimensional smoothing splines in a recursive manner. The derivation procedure sheds light on specific problems encountered in two-dimensional filtering problems. What is more, the amount of computation for this recursive processor increases only linearly with the size of the two-dimensional data. Because of certain approximations introduced in its derivation, this recursive processor becomes suboptimal.

Kim, C. S.↗

NASA NTRS1984

Development and testing of stable, invariant, isoparametric curvilinear 2- and 3-D hybrid-stress elements

Linear and quadratic Serendipity hybrid-stress elements are examined in respect of stability, coordinate invariance, and optimality. A formulation based upon symmetry group theory successfully addresses these issues in undistorted geometries and is fully detailed for plane elements. The resulting least-order stable invariant stress polynomials can be applied as astute approximations in distorted cases through a variety of tensor components and variational principles. A distortion sensitivity study for two- and three-dimensional elements provides favorable numerical comparisons with the assumed displacement method.

Punch, E. F.↗

NASA NTRS2002

Integrated Stirling Convertor and Hall Thruster Test Conducted

An important aspect of implementing Stirling Radioisotope Generators on future NASA missions is the integration of the generator and controller with potential spacecraft loads. Some recent studies have indicated that the combination of Stirling Radioisotope Generators and electric propulsion devices offer significant trip time and payload fraction benefits for deep space missions. A test was devised to begin to understand the interactions between Stirling generators and electric thrusters. An electrically heated RG- 350 (350-W output) Stirling convertor, designed and built by Stirling Technology Company of Kennewick, Washington, under a NASA Small Business Innovation Research agreement, was coupled to a 300-W SPT-50 Hall-effect thruster built for NASA by the Moscow Aviation Institute (RIAME). The RG-350 and the SPT-50 shown, were installed in adjacent vacuum chamber ports at NASA Glenn Research Center's Electric Propulsion Laboratory, Vacuum Facility 8. The Stirling electrical controller interfaced directly with the Hall thruster power-processing unit, both of which were located outside of the vacuum chamber. The power-processing unit accepted the 48 Vdc output from the Stirling controller and distributed the power to all the loads of the SPT-50, including the magnets, keeper, heater, and discharge. On February 28, 2001, the Glenn test team successfully operated the Hall-effect thruster with the Stirling convertor. This is the world's first known test of a dynamic power source with electric propulsion. The RG-350 successfully managed the transition from the purely resistive load bank within the Stirling controller to the highly capacitive power-processing unit load. At the time of the demonstration, the Stirling convertor was operating at a hot temperature of 530 C and a cold temperature of -6 C. The linear alternator was producing approximately 250 W at 109 Vac, while the power-processing unit was drawing 175 W at 48 Vdc. The majority of power was delivered to the Hall thruster discharge circuit operating at 115 Vdc and 0.9 A. Testing planned for late 2001 will examine the possibility of directly driving the Hall thruster discharge circuit using rectified and filtered output from the Stirling alternator.

Mason, Lee S.↗

DOE OSTI2025

Deep Koopman operators for causal discovery

Causal discovery aims to identify cause-effect mechanisms for better scientific understanding, explainable decision-making, and more accurate modeling. Standard statistical frameworks, such as Granger causality, lack the ability to quantify causal relationships in nonlinear dynamics due to the presence of complex feedback mechanisms, timescale mixing, and nonstationarity. Thus, applying these methods to study causal dynamics in real-world systems, such as the Earth, is a major challenge. Addressing this shortcoming, we leverage deep learning and a Koopman operator-theoretic formalism to present a class of causal discovery algorithms. Kausal uses deep Koopman operator methods to approximate nonlinear dynamics in a linearized vector space in which traditional causal inference methods such as Granger causality can be more easily applied. Our idealized experiments demonstrate Kausal’s superior ability in discovering and characterizing causal signals compared to existing deep learning and non-deep learning state-of-the-art approaches. Finally, the successful identification of major El Niño and La Niña events in observations showcases Kausal’s skill to handle real-world applications.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1994

Delay Line Detectors for the UVCS and Sumer Instruments on the SOHO Satellite

Microchannel plate based detectors with cross delay line image readout have been rapidly implemented for the SUMER and UVCS instruments aboard the Solar Orbiting Heliospheric Observatory (SOHO) mission to be launched in July 1995. In October 1993 a fast track program to build and characterize detectors and detector control electronics was initiated. We present the detector system design for the SOHO UVCS and SUMER detector programs, and results from the detector test program. Two deliverable detectors have been built at this point, a demonstration model for UVCS, and the flight Ly alpha detector for UVCS, both of which are to be delivered in the next few weeks. Test results have also been obtained with one other demonstration detector system. The detector format is 26mm x 9mm, with 1024 x 360 digitized pixels, using a low resistance Z stack of microchannel plates (MCP's) and a multilayer cross delay line anode (XDL). This configuration provides gains of approximately 2 x 10(exp 7) with good pulse height distributions (less than 50% FWHM) under uniform flood illumination, and background levels typical for this configuration (approximately 0.6 event cm (exp -2)sec(exp -1)). Local counting rates up to about 400 events/pixel/sec have been achieved with no degradation of the MCP gain. The detector and event encoding electronics achieves about 25 millimeter FVHM with good linearity (plus or minus approximately 1 pixel) and is stable to high global counting rates (greater than 4 x 10(exp 5) events sec(exp -1)). Flat field images are dominated by MCP fixed pattern noise and are stable, but the MCP multifiber modulation usually expected is uncharacteristically absent. The detector and electronics have also successfully passed both thermal vacuum and vibration tests.

Seigmund, O. H. W.↗

NASA NTRS2000

Nondestructive Determination of Bond Strength

Although many nondestructive techniques have been applied to detect disbonds in adhesive joints, no absolutely reliable nondestructive method has been developed to detect poor adhesion and evaluate the strength of bonded joints prior to the present work which used nonlinear ultrasonic methods to investigate adhesive bond cure conditions. Previously, a variety of linear and nonlinear ultrasonic methods with water coupling had been used to study aluminum-adhesive-aluminum laminates, prepared under different adhesive curing conditions, for possible bond strength determination. Therefore, in the course of this research effort, a variety of finite-amplitude experimental methods which could possibly differentiate various cure conditions were investigated, including normal and oblique incidence approaches based on nonlinear harmonic generation as well as several non-collinear two-wave interaction approaches. Test samples were mechanically scanned in various ways with respect to the focus of a transmitting transducer operated at several variable excitation frequencies and excitation levels. Even when powerful sample-related resonances were exploited by means of a frequency scanning approach, it was very difficult to isolate the nonlinear characteristics of adhesive bonds. However, a multi-frequency multi-power approach was quite successful and reliable. Ultrasonic tone burst signals at increasing power levels, over a wide frequency range, were transmitted through each bond specimen to determine its excitation dependent nonlinear harmonic resonance behavior. Relative amplitude changes were observed particularly in the higher harmonic spectral data and analyzed using a local displacement and strain analysis in the linear approximation. Two analysis approaches of the excitation-dependent data at specific resonances were found to be quite promising. One of these approaches may represent a very robust algorithm for classifying an adhesive bond as being properly cured or not. Another approach, in addition to differentiation between various cure conditions, may even provide information with respect to the bond strength. Several technical papers were published during the course of this research and a summary is presented in the Ph.D. dissertation of Tobias P. Berndt, a graduate student financially supported by this NASA Grant.

Source record↗

NASA NTRS1994

Optimizing the Zeldovich approximation

We have recently learned that the Zeldovich approximation can be successfully used for a far wider range of gravitational instability scenarios than formerly proposed; we study here how to extend this range. In previous work (Coles, Melott and Shandarin 1993, hereafter CMS) we studied the accuracy of several analytic approximations to gravitational clustering in the mildly nonlinear regime. We found that what we called the 'truncated Zeldovich approximation' (TZA) was better than any other (except in one case the ordinary Zeldovich approximation) over a wide range from linear to mildly nonlinear (sigma approximately 3) regimes. TZA was specified by setting Fourier amplitudes equal to zero for all wavenumbers greater than k(sub nl), where k(sub nl) marks the transition to the nonlinear regime. Here, we study the cross correlation of generalized TZA with a group of n-body simulations for three shapes of window function: sharp k-truncation (as in CMS), a tophat in coordinate space, or a Gaussian. We also study the variation in the crosscorrelation as a function of initial truncation scale within each type. We find that k-truncation, which was so much better than other things tried in CMS, is the worst of these three window shapes. We find that a Gaussian window e(exp(-k(exp 2)/2k(exp 2, sub G))) applied to the initial Fourier amplitudes is the best choice. It produces a greatly improved crosscorrelation in those cases which most needed improvement, e.g. those with more small-scale power in the initial conditions. The optimum choice of kG for the Gaussian window is (a somewhat spectrum-dependent) 1 to 1.5 times k(sub nl). Although all three windows produce similar power spectra and density distribution functions after application of the Zeldovich approximation, the agreement of the phases of the Fourier components with the n-body simulation is better for the Gaussian window. We therefore ascribe the success of the best-choice Gaussian window to its superior treatment of phases in the nonlinear regime. We also report on the accuracy of particle positions and velocities produced by TZA.

Melott, Adrian L.↗

DOE OSTI2022

Application of machine learning and artificial intelligence to extend EFIT equilibrium reconstruction

Recent progress in the application of machine learning (ML)/artificial intelligence (AI) algorithms to improve the Equilibrium Fitting (EFIT) code equilibrium reconstruction for fusion data analysis applications is presented. A device-independent portable core equilibrium solver capable of computing or reconstructing equilibrium for different tokamaks has been created to facilitate adaptation of ML/AI algorithms. A large EFIT database comprising of DIII-D magnetic, motional Stark effect, and kinetic reconstruction data has been generated for developments of EFIT model-order-reduction (MOR) surrogate models to reconstruct approximate equilibrium solutions. Furthermore, a neural-network MOR surrogate model has been successfully trained and tested using the magnetically reconstructed datasets with encouraging results. Other progress includes developments of a Gaussian process Bayesian framework that can adapt its many hyperparameters to improve processing of experimental input data and a 3D perturbed equilibrium database from toroidal full magnetohydrodynamic linear response modeling using the Magnetohydrodynamic Resistive Spectrum - Feedback (MARS-F) code for developments of 3D-MOR surrogate models.

Gaussian process↗

NASA NTRS1985

Air fluorescence detection of large air showers below the horizon

In the interest of exploring the cosmic ray spectrum at energies greater than 10 to the 18th power eV, where flux rates at the Earth's surface drop below 100 yr(-1) km(-2) sr(-1), cosmic ray physicists have been forced to construct ever larger detectors in order to collect useful amounts of data in reasonable lengths of time. At present, the ultimate example of this trend is the Fly's Eye system in Utah, which uses the atmosphere around an array of skyward-looking photomultiplier tubes. The air acts as a scintillator to give detecting areas as large as 5000 square kilometers sr (for highest energy events). This experiment has revealed structure (and a possible cutoff) in the ultra-high energy region above 10 o the 19th power eV. The success of the Fly's Eye experiment provides impetus for continuing the development of larger detectors to make accessible even higher energies. However, due to the rapidly falling flux, a tenfold increase in observable energy would call for a hundredfold increase in the detecting area. But, the cost of expanding the Fly's Eye detecting area will approximately scale linearly with area. It is for these reasons that the authors have proposed a new approach to using the atmosphere as a scintillator; one which will require fewer photomultipliers, less hardware (thus being less extensive), yet will provide position and shower size information.

Halverson, P.↗