Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Smooth 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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
NASA NTRS1977

Landsat-D thematic mapper simulation using aircraft multispectral scanner data

A simulation of imagery from the upcoming Landsat-D Thematic Mapper was accomplished by using selected channels of aircraft 24-channel multispectral scanner data. The purpose was to simulate Thematic Mapper 30-meter resolution imagery, to compare its spectral quality with the original aircraft MSS data, and to determine changes in thematic classification accuracy for the simulated imagery. The original resolution of approximately 7.5 meters IFOV and simulated resolution of 15, 30, and 60 meters were used to indicate the trend of spectral quality and classification accuracy. The study was based in a 6.5 square kilometer area of urban Los Angeles having a diversity of land use. The original imagery was reduced in resolution by two related methods: pixel matrix averaging, and matrix smoothing with a unity box filter, followed by matrix averaging. Thematic land use classification using training sites and a Bayesian maximum-likelihood algorithm was performed at three levels of standard deviation - 1.0, 2.0, and 3.0 sigma. Plots of relative standard deviation showed that for larger training sites with a normal distribution of data, as the resolution decreased, the distribution range of density values decreased. Also, the classification accuracies for three levels of standard deviation increased as resolution decreased. However, the indication is that a point of diminishing returns had been reached, and 30 meters IFOV should be the best for multispectral classification of urban scenes.

Clark, J.↗
NASA NTRS1991

Impact cratering calculations

In the course of carrying out finite difference calculations, it was discovered that for large craters, a previously unrecognized type of crater (diameter) growth occurred which was called lip wave propagation. This type of growth is illustrated for an impact of a 1000 km (2a) silicate bolide at 12 km/sec (U) onto a silicate half-space at earth gravity (1 g). The von Misses crustal strength is 2.4 kbar. The motion at the crater lip associated with this wave type phenomena is up, outward, and then down, similar to the particle motion of a surface wave. It is shown that the crater diameter has grown d/a of approximately 25 to d/a of approximately 4 via lip propagation from Ut/a = 5.56 to 17.0 during the time when rebound occurs. A new code is being used to study partitioning of energy and momentum and cratering efficiency with self gravity for finite-sized objects rather than the previously discussed planetary half-space problems. These are important and fundamental subjects which can be addressed with smoothed particle hydrodynamic (SPH) codes. The SPH method was used to model various problems in astrophysics and planetary physics. The initial work demonstrates that the energy budget for normal and oblique impacts are distinctly different than earlier calculations for silicate projectile impact on a silicate half space. Motivated by the first striking radar images of Venus obtained by Magellan, the effect of the atmosphere on impact cratering was studied. In order the further quantify the processes of meteor break-up and trajectory scattering upon break-up, the reentry physics of meteors striking Venus' atmosphere versus that of the Earth were studied.

Ahrens, Thomas J.↗
NASA NTRS1996

Resolving Magnetic Flux Patches at the Surface of the Core

The geomagnetic field at a given epoch can be used to partition the surface of the liquid outer core into a finite number of contiguous regions in which the radial component of the magnetic flux density, B (sub r), is of one sign. These flux patches are instrumental in providing detail to surface fluid flows inferred from the changing geomagnetic field and in evaluating the validity of the frozen-flux approximation on which such inferences rely. Most of the flux patches in models of the modem field are small and enclose little flux compared to the total unsigned flux emanating from the core. To demonstrate that such patches are not required to explain the most spatially complete and accurate data presently available, those from the Magsat mission, I have constructed a smooth core field model that fits the Magsat data but does not possess small flux patches. I conclude that our present knowledge of the geomagnetic field does not allow us to resolve these features reliably at the core-mantle boundary; thus we possess less information about core flow than previously believed.

OBrien, Michael S.↗
DOE OSTI2021

A primal–dual algorithm for risk minimization

In this paper, we develop an algorithm to efficiently solve risk-averse optimization problems posed in reflexive Banach space. Such problems often arise in many practical applications as, e.g., optimization problems constrained by partial differential equations with uncertain inputs. Unfortunately, for many popular risk models including the coherent risk measures, the resulting risk-averse objective function is nonsmooth. Here, this lack of differentiability complicates the numerical approximation of the objective function as well as the numerical solution of the optimization problem. To address these challenges, we propose a primal–dual algorithm for solving large-scale nonsmooth risk-averse optimization problems. This algorithm is motivated by the classical method of multipliers and by epigraphical regularization of risk measures. As a result, the algorithm solves a sequence of smooth optimization problems using derivative-based methods. We prove convergence of the algorithm even when the subproblems are solved inexactly and conclude with numerical examples demonstrating the efficiency of our method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Advancement of hybrid fluid-kinetic modeling for HEDP and ICF science

We report on the development progress of a hybrid fluid-kinetic code for simulating fluids and plasmas in a wide range of environments, such as laser–matter interactions, inertial confinement fusion, magnetic confinement fusion, and pulsed power. The suite of numerical tools under development utilizes heterogeneous computer architectures and leverages the benefits of particle–based simulation techniques. By working to combine the kinetic particle-in-cell (PIC) model with a particle-based fluid simulation technique, such as smoothed particle hydrodynamics, we are developing a flexible framework capable of accurately modeling complex flows within and between kinetic and fluid regimes. The TriForce code is under development as a C++ framework for parallel, 3D, particle-based, hybrid fluid-kinetic plasma simulations. The fluid half of TriForce will be based upon the meshless smoothed-particle-hydrodynamics (SPH) approach, well-suited for shear, mixing, and turbulence, whereas the kinetic half resembles a traditional particle-in-cell (PIC) code; other particle-based approaches to fluid modeling that do use a mesh are also possible to use and are under investigation. Maxwell’s electromagnetic field equations are solved either via explicit or implicit algorithms or approximated via resistive magnetohydrodynamics (MHD) using an Ohm’s law and resulting induction equation (extended MHD is under development). A primary goal of enabling direct comparisons, from the same code, between results from the variants of MHD and implicit electromagnetic solutions is to improve our fundamental understanding of systems with magnetic fields. The code is under development to recover results from both radiation-MHD and fully kinetic codes in those limits, and is continuing to be developed from other follow-on grants to operate in between where both descriptions may co-exist and interact. For certain applications, it is desired for a simulation to contain fluid ions and electrons as well as kinetic ions and electrons. Typically, it is too computationally intensive to model a full-scale ICF or HEDP experiment fully kinetically since many cycles are expended with very small time steps on modeling the fluid part of a material that is well treated by the fluid approximation. In this case, many traditional PIC particles can be replaced with a single fluid particle representing the thermal part of the distribution function, and there are fewer needed kinetic particles, which describe the non-thermal part and can be sub-cycled relative to the fluid particle advance. Furthermore, a pure fluid code may, depending on the problem, simply lack many physically important details that are beyond the scope of its reduced approximations and assumptions. In this report, we summarize the objectives achieved in the development of the collisional and kinetic half of the code, and the physics problems to which the code has been applied in the areas of advanced and innovative fusion concepts, pulsed power, and magneto-inertial fusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Cholesky-based experimental design for Gaussian process and kernel-based emulation and calibration.

Gaussian processes and other kernel-based methods are used extensively to construct approximations of multivariate data sets. The accuracy of these approximations is dependent on the data used. This paper presents a computationally efficient algorithm to greedily select training samples that minimize the weighted L p error of kernel-based approximations for a given number of data. The method successively generates nested samples, with the goal of minimizing the error in high probability regions of densities specified by users. The algorithm presented is extremely simple and can be implemented using existing pivoted Cholesky factorization methods. Training samples are generated in batches which allows training data to be evaluated (labeled) in parallel. For smooth kernels, the algorithm performs comparably with the greedy integrated variance design but has significantly lower complexity. Numerical experiments demonstrate the efficacy of the approach for bounded, unbounded, multi-modal and non-tensor product densities. We also show how to use the proposed algorithm to efficiently generate surrogates for inferring unknown model parameters from data using Bayesian inference.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1933

The Aerodynamic Characteristics of Airfoils as Affected by Surface Roughness

The effect on airfoil characteristics of surface roughness of varying degrees and types at different locations on an airfoil was investigated at high values of the Reynolds number in a variable density wind tunnel. Tests were made on a number of National Advisory Committee for Aeronautics (NACA) 0012 airfoil models on which the nature of the surface was varied from a rough to a very smooth finish. The effect on the airfoil characteristics of varying the location of a rough area in the region of the leading edge was also investigated. Airfoils with surfaces simulating lap joints were also tested. Measurable adverse effects were found to be caused by small irregularities in airfoil surfaces which might ordinarily be overlooked. The flow is sensitive to small irregularities of approximately 0.0002c in depth near the leading edge. The tests made on the surfaces simulating lap joints indicated that such surfaces cause small adverse effects. Additional data from earlier tests of another symmetrical airfoil are also included to indicate the variation of the maximum lift coefficient with the Reynolds number for an airfoil with a polished surface and with a very rough one.

HOCKER RAY W↗
NASA NTRS1999

Gaussian Decomposition of Laser Altimeter Waveforms

We develop a method to decompose a laser altimeter return waveform into its Gaussian components assuming that the position of each Gaussian within the waveform can be used to calculate the mean elevation of a specific reflecting surface within the laser footprint. We estimate the number of Gaussian components from the number of inflection points of a smoothed copy of the laser waveform, and obtain initial estimates of the Gaussian half-widths and positions from the positions of its consecutive inflection points. Initial amplitude estimates are obtained using a non-negative least-squares method. To reduce the likelihood of fitting the background noise within the waveform and to minimize the number of Gaussians needed in the approximation, we rank the "importance" of each Gaussian in the decomposition using its initial half-width and amplitude estimates. The initial parameter estimates of all Gaussians ranked "important" are optimized using the Levenburg-Marquardt method. If the sum of the Gaussians does not approximate the return waveform to a prescribed accuracy, then additional Gaussians are included in the optimization procedure. The Gaussian decomposition method is demonstrated on data collected by the airborne Laser Vegetation Imaging Sensor (LVIS) in October 1997 over the Sequoia National Forest, California.

Hofton, Michelle A.↗
DOE OSTI2020

Toward a systematic improvement of the fixed-node approximation in diffusion Monte Carlo for solids—A case study in diamond

While Diffusion Monte Carlo (DMC) is in principle an exact stochastic method for ab initio electronic structure calculations, in practice, the fermionic sign problem necessitates the use of the fixed-node approximation and trial wavefunctions with approximate nodes (or zeros). This approximation introduces a variational error in the energy that potentially can be tested and systematically improved. Here, we present a computational method that produces trial wavefunctions with systematically improvable nodes for DMC calculations of periodic solids. These trial wavefunctions are efficiently generated with the configuration interaction using a perturbative selection made iteratively (CIPSI) method. A simple protocol in which both exact and approximate results for finite supercells are used to extrapolate to the thermodynamic limit is introduced. This approach is illustrated in the case of the carbon diamond using Slater–Jastrow trial wavefunctions including up to one million Slater determinants. Fixed-node DMC energies obtained with such large expansions are much improved, and the fixed-node error is found to decrease monotonically and smoothly as a function of the number of determinants in the trial wavefunction, a property opening the way to a better control of this error. The cohesive energy extrapolated to the thermodynamic limit is in close agreement with the estimated experimental value. Interestingly, this is also the case at the single-determinant level, thus, indicating a very good error cancellation in carbon diamond between the bulk and atomic total fixed-node energies when using single-determinant nodes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Global Microscopic Description of Nucleon-Nucleus Scattering with Quantified Uncertainties

We develop for the first time a microscopic global nucleon-nucleus optical potential with quantified uncertainties suitable for analyzing nuclear reaction experiments at next-generation rare-isotope beam facilities. Within the improved local density approximation and without any adjustable parameters, we begin by computing proton-nucleus and neutron-nucleus optical potentials from a set of five nuclear forces from chiral effective field theory for 1800 target nuclei in the mass range 12 ≤ A ≤ 242 for energies between 0 MeV < E ≲ 150 MeV. We then parameterize a global optical potential for each chiral force that depends smoothly on the projectile energy as well as the target nucleus mass number and isospin asymmetry. Uncertainty bands for elastic scattering observables are generated from a full covariance analysis of the parameters entering in the description of our global optical potential and benchmarked against existing experimental data for stable target nuclei. Since our approach is purely microscopic, we anticipate a similar quality of the model for nucleon scattering on unstable isotopes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

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↗
NASA NTRS1972

Delta-90 Interplanetary Monitoring Platform-H (IMP-H) flash flight report

The Delta-90 launch vehicle and the IMP-H spacecraft were successfully launched from Pad B, Complex 17, Cape Kennedy Air Force Station, Florida, at 2120:00.559 EDT on September 22, 1972. The countdown proceeded smoothly to liftoff with no major difficulties or unscheduled holds. The Delta-90/IMP-H were launched on a pad azimuth of 115 degrees, the vehicle ten rolled to 95 degrees from the north placing the spacecraft in a highly elliptical transfer orbit. Firing the spacecraft kickmotor at 1136 EDT, September 25, 1972, injected the spacecraft into its final desirable near-circular orbit approximately half way between the planet earth and its moon. Vehicle performance of all stages appeared nominal with all sequenced events occurring at the expected times. Data acquisition from all range stations was very good. Damage to the launch pad caused by liftoff was nominal.

Source record↗
NASA NTRS1973

An experimental evaluation of metallic diaphragms for positive fuel expulsion in the Atmosphere Explorer, hydrazine propulsion subsystem

Four Arde conospheroid metallic diaphragms were tested at NASA to evaluate their capability for use in the orbit adjust propulsion subsystem (OAPS) of the Atmosphere Explorer spacecraft (AE's C, D, and E). The diaphragms will be used for positive propellant expulsion and spacecraft center of mass (c.m.) control. A leak-free cycle life capability of nine reversals was demonstrated. The diaphragms rolled smoothly from ring to ring in a predictable manner on the first reversal. Varying amounts of diaphragm cocking and ring skipping were observed on subsequent reversals. The diaphragm pressure differential did not exceed 10 psid during any reversal. Cycle life capability, reversal mode, and pressure differential were not affected by sudden reversals, environmental tests, or 18,000 partial reversals. An expulsion efficiency of approximately 97 percent was demonstrated. The results of these tests show that metallic diaphragms can be used as an effective means of positive fuel expulsion; however, to achieve spacecraft c.m. control, the diaphragm must not be reversed prior to flight.

Woodruff, W. L.↗
NASA NTRS1975

Magnetospheric protons and electrons encountered by the moon in the plasma sheet

Results are summarized for a study of plasma-sheet particles (and their flow) encountered by the moon during its passage through the geomagnetic tail. The study is based on analysis of data obtained with a modulated Faraday cup on board the lunar-anchored spacecraft Explorer 35. It is shown that the electrons have a rapidly fluctuating non-Maxwellian energy distribution with a mean energy of several hundred electron volts and a density of approximately 0.2 per cu cm. Protons with energies of the order of 1 keV were usually detected above the instrument background when flowing towards earth at about 200 km/sec. It is suggested that strong terrestrial polar winds during the early history of the earth-moon system could have caused some erosion of the front side of the moon and that the relative smoothness of the same side could be explained in terms of gravitational shielding by the earth from the interplanetary rock flux.

Prakash, A.↗
NASA NTRS1986

Spherical harmonic expansions of the Earth's gravitational potential to degree 360 using 30' mean anomalies

Two potential coefficient fields that are complete to degree and order 360 have been computed. One field (OSU86E) excludes geophysically predicted anomalies while the other (OSU86F) includes such anomalies. These fields were computed using a set of 30' mean gravity anomalies derived from satellite altimetry in the ocean areas and from land measurements in North America, Europe, Australia, Japan and a few other areas. Where no 30' data existed, 1 deg x 1 deg mean anomaly estimates were used if available. No rigorous combination of satellite and terrestrial data was carried out. Instead advantage was taken of the adjusted anomalies and potential coefficients from a rigorous combination of the GEML2' potential coefficient set and 1 deg x 1 deg mean gravity anomalies. The two new fields were computed using a quadrature procedure with de-smoothing factors. The spectra of the new fields agree well with the spectra of the fields with 1 deg x 1 deg data out to degree 180. Above degree 180 the new fields have more power. The fields have been tested through comparison of Doppler station geoid undulations with undulations from various geopotential models. The agreement between the two types of undulations is approximately + or - 1.6 m. The use of a 360 field over a 180 field does not significantly improve the comparison. Instead it allows the comparison to be done at some stations where high frequency effects are important. In addition maps made in areas of high frequency information (such as trench areas) clearly reveal the signal in the new fields from degree 181 to 360.

Rapp, Richard H.↗
NASA NTRS1994

Observations of hysteresis in solar cycle variations among seven solar activity indicators

We show that smoothed time series of 7 indices of solar activity exhibit significant solar cycle dependent differences in their relative variations during the past 20 years. In some cases these observed hysteresis patterns start to repeat over more than one solar cycle, giving evidence that this is a normal feature of solar variability. Among the indices we study, we find that the hysteresis effects are approximately simple phase shifts, and we quantify these phase shifts in terms of lag times behind the leading index, the International Sunspot Number. Our measured lag times range from less than one month to greater than four months and can be much larger than lag times estimated from short-term variations of these same activity indices during the emergence and decay of major active regions. We argue that hysteresis represents a real delay in the onset and decline of solar activity and is an important clue in the search for physical processes responsible for changing solar emission at various wavelengths.

Bachmann, Kurt T.↗
NASA NTRS2004

Chandra Observations of the X-ray Environs of SN 1998bw/GRB 980425

We report X-ray studies of the environs of SN 1998bw and GRB 980425 using the Chandra X-Ray Observatory 1281 days after the gamma-ray burst (GRB). Eight X-ray point sources were localized, three and five each in the original error boxes, S1 and S2, assigned for variable X-ray counteparts to the GRB by BeppoSAX. The sum of the discrete X-ray sources plus continuous emission in S2 observed by Chandra on day 1281 is within a factor of 1.5 of the maximum and the upper limits seen by BeppoSAX. We conclude that S2 is the sum of several variable sources that have not disappeared and therefore is not associated with the GRB. Within S1, clear evidence is seen for a decline of approximately a factor of 12 between day 200 and day 1281. One of the sources in S 1, S 1 a, is coincident with the well-determined radio location of SN 1998bw and is certainly the remnant of that explosion. The nature of the other sources is also discussed. Combining our observation of the supernova with others of the GRB afterglow, a smooth X-ray light curve, spanning approx. 1400 days, is obtained by assuming that the burst and supernova were coincident at 35.6 Mpc. When this X-ray light curve is compared with those of the X-ray af "erglows" of ordinary GRBs, X-ray flashes, and ordinary supernovae, evidence emerges for at least two classes of light curves, perhaps bounding a continuum. By 3-10 yr, all these phenomena seem to converge on a common X-ray luminosity, possibly indicative of the supernova underlying them all. This convergence strengthens the conclusion that SN 1998 bw aid GRB 980425 took place in the same object.One possible explanation for the two classes is that a (nearly) standard GRB was observed at different angles, in which case X-ray afterglows with intermediate luminosities should eventually be discovered. Finally, we comment on the contribution of GRB afterglows to the ultraluminous X-ray source population.

Kouveliotou, C.↗
DOE OSTI2023

Insights into the failure mode of the Chelyabinsk meteor from high-fidelity simulation

In 2013, an asteroid with a diameter of approximate 20 m broke up in the atmosphere above the Russian city of Chelyabinsk producing what is known as an airburst event. The event was recorded on video from multiple angles and a 500-kg fragment was recovered from Lake Chebarkul shortly thereafter. These aspects, coupled with the asteroid’s large size, distinguish the Chelyabinsk event as unique. In this study, we use the smoothed particle hydrodynamics method to examine the breakup mode of a Chelyabinsk-sized monolithic asteroid and probe whether a monolithic structure is consistent with observations from the Chelyabinsk event. We model the object as a homogeneous sphere with properties consistent with Chelyabinsk meteorites. We use planar-2D simulations to test the sensitivity of results to variations in the material model and to qualitatively inspect the cause of crack nucleation and fragmentation. We then compare results from a 3D simulation to observations from the event. Our 3D simulation predicts a burst height within 3 km of the observed peak luminosity suggesting the responsible asteroid may have been monolithic.

79 ASTRONOMY AND ASTROPHYSICS↗