Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical approximations”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,675 records · Page 93
All sourcesNASA NTRSDOE OSTI
NASA NTRS1998

Detection of Non-Symmetrical Damage in Smart Plate-Like Structures

A two-dimensional model for in-plane vibrations of a cantilever plate with a non-symmetrical damage is used in the context of defect identification in materials with piezoelectric ceramic patches bonded to their surface. These patches can act both as actuators and sensors in a self-analyzing fashion, which is a characteristic of smart materials. A Galerkin method is used to approximate the dynamic response of these structures. The natural frequency shifts due to the damage are estimated numerically and compared to experimental data obtained from tests on cantilever aluminum plate-like structures damaged at different locations with defects of different depths. The damage location and extent are determined by an enhanced least square identification method. Efficacy of the frequency shift based algorithms is demonstrated using experimental data.

Blanks, H. T.↗
NASA NTRS1999

Evidence for Methane Segregation at the Surface of Pluto

In May 1995, a set of spectrophotometric curves of the system Pluto-Charon was recorded with the UKIRT telescope equipped with the spectrometer CGS4. As for the previous observations, the spectra cover a part of the near infrared range, between 1.4 and 2.55 micrometers, but with a higher resolution of approximately 700. In both the 1992 and 1995 data, the existence of solid methane is confirmed by numerous absorption bands, and the carbon monoxide and the nitrogen ices are identified by their respective signatures at 2.35 and 2.15 um. The solid nitrogen seems to be the principal icy component and forms a matrix in which the CH4 and CO molecules are diluted. However a spectroscopic analysis of the 1995 observations indicates that pure methane may coexist with its diluted phase in N2. In order to derive the horizontal and vertical distribution of these different species and to obtain some quantitative information about their characteristics, we have modeled the spectrum of May 15 that corresponds to the maximum of Pluto's visible light curve. This was achieved by means of a radiative transfer algorithm dealing with compact and stratified media. Among the various representations we have tested to describe the surface of Pluto, only a geographical mixture of three distinct units explains all the significant structures of the analyzed spectrum. The first unit is a thin granular layer of pure CH4 covering a compact polycrystalline substratum of N2-CH4-CO, which are in a molecular mixture (concentrations of and CO of the order of 0.45%, 0.1-0.2% respectively). It covers about 70% of the observed area and corresponds to volatile deposits that are sublimating under solar illumination. The second unit is either (a) a single thick layer of pure granular methane or (b) a unit similar to the first unit but with the two components inverted (i.e. with CH4 forming a substratum and the N2-CH4-CO mixture a superficial layer of fine grains). Covering 20% of the surface, it represents some old surfaces that have been sublimated for a long time, and eventually recovered later by very small amounts of fresh deposits of the molecular mixture N2-CH4-CO. Finally, the third unit may result from the condensation of very fine grains of nearly pure N2. It covers the remainder of the surface (about 10%). All these results allow a better understanding of the processes of deposition, metamorphism, sublimation and transport affecting the different ices detected on Pluto during its climatic cycles.

Doute, S.↗
NASA NTRS2000

Radiative Transfer and Satellite Remote Sensing of Cirrus Clouds Using FIRE-2-IFO Data

Under the support of the NASA grant, we have developed a new geometric-optics model (GOM2) for the calculation of the single-scattering and polarization properties for arbitrarily oriented hexagonal ice crystals. From comparisons with the results computed by the finite difference time domain (FDTD) method, we show that the novel geometric-optics can be applied to the computation of the extinction cross section and single-scattering albedo for ice crystals with size parameters along the minimum dimension as small as approximately 6. We demonstrate that the present model converges to the conventional ray tracing method for large size parameters and produces single-scattering results close to those computed by the FDTD method for size parameters along the minimum dimension smaller than approximately 20. We demonstrate that neither the conventional geometric optics method nor the Lorenz-Mie theory can be used to approximate the scattering, absorption, and polarization features for hexagonal ice crystals with size parameters from approximately 5 to 20. On the satellite remote sensing algorithm development and validation, we have developed a numerical scheme to identify multilayer cirrus cloud systems using AVHRR data. We have applied this scheme to the satellite data collected over the FIRE-2-IFO area during nine overpasses within seven observation dates. Determination of the threshold values used in the detection scheme are based on statistical analyses of these satellite data.

Source record↗
NASA NTRS2002

Quantitative X-ray - UV Line and Continuum Spectroscopy with Application to AGN: State-Specific Hydrogenic Recombination Cooling Coefficients for a Wide Range of Conditions

Recombination cooling, in which a free electron emits light while being captured to an ion, is an important cooling process in photoionized clouds that are optically thick or have low metallicity. State specific rather than total recombination cooling rates are needed since the hydrogen atom tends to become optically thick in high-density regimes such as Active Galactic Nuclei. This paper builds upon previous work to derive the cooling rate over the full temperature range where the process can be a significant contributor in a photoionized plasma. We exploit the fact that the recombination and cooling rates are given by intrinsically similar formulae to express the cooling rate in terms of the closely related radiative recombination rate. We give an especially simple but accurate approximation that works for any high hydrogenic level and can be conveniently employed in large-scale numerical simulations.

LaMothe, J.↗
NASA NTRS2004

Verification of Java Programs using Symbolic Execution and Invariant Generation

Software verification is recognized as an important and difficult problem. We present a norel framework, based on symbolic execution, for the automated verification of software. The framework uses annotations in the form of method specifications an3 loop invariants. We present a novel iterative technique that uses invariant strengthening and approximation for discovering these loop invariants automatically. The technique handles different types of data (e.g. boolean and numeric constraints, dynamically allocated structures and arrays) and it allows for checking universally quantified formulas. Our framework is built on top of the Java PathFinder model checking toolset and it was used for the verification of several non-trivial Java programs.

Pasareanu, Corina↗
NASA NTRS2004

Calculation of the Full Scattering Amplitude without Partial Wave Decomposition: Inclusion of Exchange - 2

The development of a practical method of accurately calculating the full scattering amplitude, without making a partial wave decomposition is continued. The method is developed in the context of electron-hydrogen scattering, and here exchange is dealt with by considering e-H scattering in the static exchange approximation. The Schroedinger equation in this approximation can be simplified to a set of coupled integro-differential equations. The equations are solved numerically for the full scattering wave function. The scattering amplitude can most accurately be calculated from an integral expression for the amplitude; that integral can be formally simplified, and then evaluated using the numerically determined wave function. The results are essentially identical to converged partial wave results.

Shertzer, Janine↗
DOE OSTI2023

Exponential acceleration of macroscopic quantum tunneling in a Floquet Ising model

The exponential suppression of macroscopic quantum tunneling (MQT) in the number of elements to be reconfigured is an essential element of broken symmetry phases. This suppression is also a core bottleneck in quantum algorithms, such as traversing an energy landscape in optimization, and adiabatic state preparation more generally. In this work, we demonstrate exponential acceleration of MQT through Floquet engineering with the application of a uniform, high frequency transverse drive field. Using the ferromagnetic phase of the transverse field Ising model in one and two dimensions as a prototypical example, we identify three phenomenological regimes as a function of drive strength. For weak drives, the system exhibits exponentially decaying tunneling rates but robust magnetic order; in the crossover regime at intermediate drive strength, we find polynomial decay of tunnelling alongside vanishing magnetic order; and at very strong drive strengths both the Rabi frequency and time-averaged magnetic order are approximately constant with increasing system size. We support these claims with extensive full wavefunction and tensor network numerical simulations, and theoretical analysis. An experimental test of these results presents a technologically important and novel scientific question accessible on NISQ-era quantum computers.

Grattan, George↗
NASA NTRS2002

Transport of Space Environment Electrons: A Simplified Rapid-Analysis Computational Procedure

A computational procedure for describing transport of electrons in condensed media has been formulated for application to effects and exposures from spectral distributions typical of electrons trapped in planetary magnetic fields. The procedure is based on earlier parameterizations established from numerous electron beam experiments. New parameterizations have been derived that logically extend the domain of application to low molecular weight (high hydrogen content) materials and higher energies (approximately 50 MeV). The production and transport of high energy photons (bremsstrahlung) generated in the electron transport processes have also been modeled using tabulated values of photon production cross sections. A primary purpose for developing the procedure has been to provide a means for rapidly performing numerous repetitive calculations essential for electron radiation exposure assessments for complex space structures. Several favorable comparisons have been made with previous calculations for typical space environment spectra, which have indicated that accuracy has not been substantially compromised at the expense of computational speed.

Nealy, John E.↗
NASA NTRS1988

Periodic response of nonlinear systems

A procedure is developed to determine approximate periodic solutions of autonomous and non-autonomous systems. The trignometric collocation method (TCM) is formalized to allow for the analysis of relatively small order systems directly in physical coordinates. The TCM is extended to large order systems by utilizing modal analysis in a component mode synthesis strategy. The procedure was coded and verified by several check cases. Numerical results for two small order mechanical systems and one large order rotor dynamic system are presented. The method allows for the possibility of approximating periodic responses for large order forced and self-excited nonlinear systems.

Nataraj, C.↗
DOE OSTI2023

How good are learning-based control v.s. model-based control for load shifting? Investigations on a single zone building energy system

Both model predictive control (MPC) and deep reinforcement learning control (DRL) have been presented as a way to approximate the true optimality of a dynamic programming problem, and these two have shown significant operational cost saving potentials for building energy systems. Furthermore, there is still a lack of in-depth quantitative studies on their approximation levels to the true optimality, especially in the building energy domain. To fill in the gap, this paper provides a numerical framework that enables the evaluation of the optimality levels of different controllers for building energy systems. This framework is then used to comprehensively compare the optimal control performance of both MPC and DRL controllers with given computation budgets for a single zone fan coil unit system. Note the optimality is estimated based on a user-specific selection of trade-off weights among energy costs, thermal comfort and control slew rates. Compared with the best optimality we can find through expensive optimization simulations, the best DRL agent can maximally approximate the optimality by 96.54%, which outperforms the best MPC whose optimality level is 90.11%. However, due to the stochasticity, the DRL agent is only expected to approximate the optimality by 90.42%, which is almost equivalent to the best MPC. Except for Proximal Policy Optimization (PPO), all DRL agents can have a better approximation to the optimality than the best MPC, and are expected to have better approximation than the MPC with a prediction horizon of 32 steps (15 min per step). In terms of reducing energy cost and thermal discomfort, MPC can outperform the rule-based control (RBC) by 18.47%–25.44%. DRL can be expected to outperform RBC by 18.95%–25.65% ,and the best DRL control policy can outperform RBC by 20.29%–29.72%. Although the comparison of the optimality level is performed in a perfect setting, e.g., MPC assumes perfect models, and DRL assumes a perfect offline training process and online deployment process, this can shed insight on their capabilities of approximating to the original dynamic programming problem.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS2016

The XXL Survey I. Scientific Motivations - Xmm-Newton Observing Plan - Follow-up Observations and Simulation Programme

The quest for the cosmological parameters that describe our universe continues to motivate the scientific community to undertake very large survey initiatives across the electromagnetic spectrum. Over the past two decades, the Chandra and XMM-Newton observatories have supported numerous studies of X-ray-selected clusters of galaxies, active galactic nuclei (AGNs), and the X-ray background. The present paper is the first in a series reporting results of the XXL-XMM survey; it comes at a time when the Planck mission results are being finalized. Aims. We present the XXL Survey, the largest XMM programme totaling some 6.9 Ms to date and involving an international consortium of roughly 100 members. The XXL Survey covers two extragalactic areas of 25 deg2 each at a point-source sensitivity of approx. 5 x 10(exp 15) erg/s/sq cm in the [0.5-2] keV band (completeness limit). The surveys main goals are to provide constraints on the dark energy equation of state from the space-time-distribution of clusters of galaxies and to serve as a pathfinder for future, wide-area X-ray missions. We review science objectives, including cluster studies, AGN evolution, and large-scale structure, that are being conducted with the support of approximately 30 follow-up programs. Methods. We describe the 542 XMM observations along with the associated multi- and numerical simulation programmes. We give a detailed account of the X-ray processing steps and describe innovative tools being developed for the cosmological analysis. Results. The paper provides a thorough evaluation of the X-ray data, including quality controls, photon statistics, exposure and background maps, and sky coverage. Source catalogue construction and multi-associations are briefly described. This material will be the basis for the calculation of the cluster and AGN selection functions, critical elements of the cosmological and science analyses. Conclusions. The XXL multi- data set will have a unique lasting legacy value for cosmological and extragalactic studies and will serve asa calibration resource for future dark energy studies with clusters and other X-ray selected sources. With the present article, we release the XMM XXL photon and smoothed images along with the corresponding exposure maps.

X-rays: general – large-scale structure of Unive↗
NASA NTRS1974

Numerical evaluation of electromagnetic fields due to dipole antennas in the presence of stratified media

Two numerical methods are used to evaluate the integrals that express the em fields due to dipole antennas radiating in the presence of a stratified medium. The first method is a direct integration by means of Simpson's rule. The second method is indirect and approximates the kernel of the integral by means of the fast Fourier transform. In contrast to previous analytical methods that applied only to two-layer cases the numerical methods can be used for any arbitrary number of layers with general properties.

Tsang, L.↗
DOE OSTI2023

Exact Gaussian processes for massive datasets via non-stationary sparsity-discovering kernels

Abstract A Gaussian Process (GP) is a prominent mathematical framework for stochastic function approximation in science and engineering applications. Its success is largely attributed to the GP’s analytical tractability, robustness, and natural inclusion of uncertainty quantification. Unfortunately, the use of exact GPs is prohibitively expensive for large datasets due to their unfavorable numerical complexity of $$O(N^3)$$ O ( N 3 ) in computation and $$O(N^2)$$ O ( N 2 ) in storage. All existing methods addressing this issue utilize some form of approximation—usually considering subsets of the full dataset or finding representative pseudo-points that render the covariance matrix well-structured and sparse. These approximate methods can lead to inaccuracies in function approximations and often limit the user’s flexibility in designing expressive kernels. Instead of inducing sparsity via data-point geometry and structure, we propose to take advantage of naturally-occurring sparsity by allowing the kernel to discover—instead of induce—sparse structure. The premise of this paper is that the data sets and physical processes modeled by GPs often exhibit natural or implicit sparsities, but commonly-used kernels do not allow us to exploit such sparsity. The core concept of exact, and at the same time sparse GPs relies on kernel definitions that provide enough flexibility to learn and encode not only non-zero but also zero covariances. This principle of ultra-flexible, compactly-supported, and non-stationary kernels, combined with HPC and constrained optimization, lets us scale exact GPs well beyond 5 million data points.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1985

Time-dependent nonlinear cosmic ray shocks confirming abstract

Numerical studies of time dependent cosmic ray shock structures in planar geometry are interesting because analytical time-independent solutions are available which include the non-linear reactions on the plasma flow. A feature of these time asymptotic solutions is that for higher Mach numbers (M approximately 5) and for a low cosmic ray upstream pressure the solution is not uniquely determined by the usual conservation laws of mass, momentum and energy. These numerical solutions clearly indicate that much work needs to be done before we understand shock acceleration as a time dependent process. The slowness of the process is possibly due to the fact that there is a diffusive flux into the downstream region in addition to the usual advective losses. Analytic investigations of thie phenomenon are required.

Dorfi, E. A.↗
NASA NTRS1990

Numerical simulation of solar wind density fluctuations and their effects on VLF radio interferometry

A mission consisting of an array of small satellites has been proposed to investigate radio sources at frequencies from approximately 2 to 20 MHz, a range unobservable from the ground because of ionospheric absorption and aberrations. Such a mission seems feasible and relatively inexpensive, but further investigation is still necessary. The solar wind is a major concern because it contains turbulence-induced random density variations that cause a fluctuating phase shift in any electromagnetic signal passing through, sometimes making radio source imaging by interferometry impossible. A numerical simulation has been developed which attempts to accurately characterize solar wind density fluctuations. Results have allowed more concrete conclusions to be drawn about observing constraints and have shown that interferometric imaging at radio frequencies below 10 MHz with baselines up to 100 km is possible at solar elongation angles of about 90 deg or more.

Williamson, Robert S., III↗
NASA NTRS1994

Time-stable boundary conditions for finite-difference schemes solving hyperbolic systems: Methodology and application to high-order compact schemes

We present a systematic method for constructing boundary conditions (numerical and physical) of the required accuracy, for compact (Pade-like) high-order finite-difference schemes for hyperbolic systems. First a proper summation-by-parts formula is found for the approximate derivative. A 'simultaneous approximation term' is then introduced to treat the boundary conditions. This procedure leads to time-stable schemes even in the system case. An explicit construction of the fourth-order compact case is given. Numerical studies are presented to verify the efficacy of the approach.

Carpenter, Mark H.↗
NASA NTRS1991

Investigation of the structure of the electromagnetic field and related phenomena, generated by the Active Satellite

The altitude dependencies of the moduli of the electric field E in the VLF and LF frequency bands (f sub B much less than F less than f sub B) and in the altitude range of the ionosphere Z equals (400 to 2500) km up to Z equals 6000 km (the bottom of the magnetosphere) were calculated by the linear theory. The amplitudes of the field have large maxima in four regions: the axis field (E sub o) close to the direction of the Earth's magnetic field line B sub o, beta approximately 0 degrees, the fields (E sub St), (E sub RevSt) and (E sub Res) in the Storey, Reversed Story and Resonance cones, beta approximately (0 approaches 20) degrees. Their maxima are very pronounced close to the low hybrid frequency F sub L. The nonlinear heating of the magnetoplasma under the action of an electric field Ee (sup iwt) was recently expanded by the macroscopic theory by the author. The velocities, collision frequencies and temperatures of all the constituents of a magnetoplasma - electrons, ions, and neutral particles - are taken into account. Formulae and numerical results are presented for the ionosphere in the frequency band F equals (1 to 10 exp 4) kHz and altitude range Z approximately (100 - 1000) km. Some results of calculations by the self consistent solution of the basis system of equations are also discussed.

Alpert, Yakov L.↗
NASA NTRS2022

A NASA-Air Force Precipitation Analysis for Near-Real-Time Operations

This article describes a new precipitation analysis algorithm developed by NASA for time-sensitive operations at the United States Air Force. Implemented as part of the Land Information System—a land modeling and data assimilation software framework—this NASA–Air Force Precipitation Analysis (NAFPA) combines numerical weather prediction model outputs with rain gauge measurements and satellite estimates to produce global, gridded 3-h accumulated precipitation fields at approximately 10-km resolution. Input observations are subjected to quality control checks before being used by the Bratseth analysis algorithm that converges to optimal interpolation. NAFPA assimilates up to 3.5 million observations without artificial data thinning or selection. To evaluate this new approach, a multiyear reanalysis is generated and intercompared with eight alternative precipitation products across the contiguous United States, Africa, and the monsoon region of eastern Asia. NAFPA yields superior accuracy and correlation over low-latency (up to 14 h) alternatives (numerical weather prediction and satellite retrievals), and often outperforms high-latency (up to 3.5 months) products, although the details for the latter vary by region and product. The development of NAFPA offers a high-quality, near-real-time product for use in meteorological, land surface, and hydrological research and applications.

Africa↗