Search NASA⌕ Search

SEARCH · Search NASA

Results for “splitting methods”

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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
NASA NTRS2001

Did the Eruption of the Mt. Pinatubo Volcano Affect Cirrus Properties?

Some observations suggest that the volcanic aerosols produced by the Mt. Pinatubo eruption may have altered cirrus properties. We look for evidence that such modification of cirrus is extensive enough to be climatically significant by comparing three satellite-based cirrus datasets produced by the ISCCP (International Satellite Cloud Climatology Project) analysis, the 'split-window' analysis, and 3I analysis. Since the former two have not been compared in detail before, we conduct such a comparison here. When applied to AVHRR (Advanced Very High Resolution Radiometer) data, both the ISCCP and split-window analyses identify about 0.2 - 0.3 cirrus cloud amount in tropical latitudes; however, there are detailed differences of classification for about half of these clouds. The discrepancies are attributed to the simplified assumptions made by both methods. The latter two datasets are derived from infrared radiances, so they are much less sensitive to volcanic aerosols than the ISCCP analysis. After the Mt. Pinatubo eruption, the ISCCP results indicate a dramatic decrease of thin cirrus (cloud top pressure less than 440 mb and visible optical thickness less than 1.3) over ocean, accompanied by a comparable increase of altocumulus and cumulus clouds; over land, there are no significant changes. In contrast, results from the split-window and 3I analyses show little change in thin cirrus amount over either ocean or land that is associated with the volcanic eruption. The ISCCP results can, therefore, be understood as a misclassification of thin cirrus because the additional reflected sunlight by the volcanic aerosol makes the cirrus clouds appear to be optically thicker. Examination of the split-window signature and the infrared emissivities from 3I show no significant change in infrared emissivity (or optical thickness). These results indicate that the Mt. Pinatubo volcanic aerosol did not have a significant systematic effect on tropical cirrus properties, but rather produced only temporary, local effects. Hence, these results indicate that there was no significant climate feedback produced by aerosol-cirrus-radiative interactions.

Luo, Zhengzhao↗
DOE OSTI2024

Double diffraction imaging of x-ray induced structural dynamics in single free nanoparticles

Abstract Because of their high photon flux, x-ray free-electron lasers (FEL) allow to resolve the structure of individual nanoparticles via coherent diffractive imaging (CDI) within a single x-ray pulse. Since the inevitable rapid destruction of the sample limits the achievable resolution, a thorough understanding of the spatiotemporal evolution of matter on the nanoscale following the irradiation is crucial. We present a technique to track x-ray induced structural changes in time and space by recording two consecutive diffraction patterns of the same single, free-flying nanoparticle, acquired separately on two large-area detectors opposite to each other, thus examining both the initial and evolved particle structure. We demonstrate the method at the extreme ultraviolet (XUV) and soft x-ray Free-electron LASer in Hamburg (FLASH), investigating xenon clusters as model systems. By splitting a single XUV pulse, two diffraction patterns from the same particle can be obtained. For focus intensities of about 2 × 10 12 W cm −2 we observe still largely intact clusters even at the longest delays of up to 650 picoseconds of the second pulse, indicating that in the highly absorbing systems the damage remains confined to one side of the cluster. Instead, in case of five times higher flux, the diffraction patterns show clear signatures of disintegration, namely increased diameters and density fluctuations in the fragmenting clusters. Future improvements to the accessible range of dynamics and time resolution of the approach are discussed.

Sauppe, M. (ORCID:0000000163123102)↗
DOE OSTI2020

Separating the intrinsic alignment signal and the lensing signal using self-calibration in photo- z surveys with KiDS450 and KV450 Data

ABSTRACT To reach the full potential for the next generation of weak lensing surveys, it is necessary to mitigate the contamination of intrinsic alignments (IAs) of galaxies in the observed cosmic shear signal. The self-calibration (SC) of IAs provides an independent method to measure the IA signal from the survey data and the photometric redshift information. It operates differently from the marginalization method based on the IA modelling. In this work, we present the first application of SC to the KiDS450 data and the KV450 data, to split directly the intrinsic shape–galaxy density (Ig) correlation signal and the gravitational shear–galaxy density (Gg) correlation signal, using the information from photometric redshift (photo-z). We achieved a clear separation of the two signals and performed several validation tests. Our measured signals are found to be in general agreement with the KiDS450 cosmic shear best-fitting cosmology, for both lensing and IA measurements. For KV450, we use partial (high-z) data, and our lensing measurements are also in good agreement with KV450 cosmic shear best fit, while our IA signal suggests a larger IA amplitude for the high-z sample. We discussed the impact of photo-z quality on IA detection and several other potential systematic biases. Finally, we discuss the potential application of the information extracted for both the lensing signal and the IA signal in future surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Angular clustering properties of the DESI QSO target selection using DR9 Legacy Imaging Surveys

ABSTRACT The quasar target selection for the upcoming survey of the Dark Energy Spectroscopic Instrument (DESI) will be fixed for the next 5 yr. The aim of this work is to validate the quasar selection by studying the impact of imaging systematics as well as stellar and galactic contaminants, and to develop a procedure to mitigate them. Density fluctuations of quasar targets are found to be related to photometric properties such as seeing and depth of the Data Release 9 of the DESI Legacy Imaging Surveys. To model this complex relation, we explore machine learning algorithms (random forest and multilayer perceptron) as an alternative to the standard linear regression. Splitting the footprint of the Legacy Imaging Surveys into three regions according to photometric properties, we perform an independent analysis in each region, validating our method using extended Baryon Oscillation Spectroscopic Survey (eBOSS) EZ-mocks. The mitigation procedure is tested by comparing the angular correlation of the corrected target selection on each photometric region to the angular correlation function obtained using quasars from the Sloan Digital Sky Survey (SDSS) Data Release 16. With our procedure, we recover a similar level of correlation between DESI quasar targets and SDSS quasars in two-thirds of the total footprint and we show that the excess of correlation in the remaining area is due to a stellar contamination that should be removed with DESI spectroscopic data. We derive the Limber parameters in our three imaging regions and compare them to previous measurements from SDSS and the 2dF QSO Redshift Survey.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Defensive Islanding to Enhance the Resilience of Distribution Systems Against Cyber-Induced Failures

The extensive integration of communication, computation, and control technologies into cyber-physical power systems (CPPSs) has increased the vulnerabilities of CPPSs to cyberattacks. This calls for developing solutions that assess and reduce the impacts of cyber-induced failures on CPPSs. This paper proposes a defensive islanding strategy to isolate impacted parts of the CPPS and form self-sufficient islanded grids with an objective of minimum load curtailment. The defensive islanding aims to split a power system into smaller grids to improve its resilience against a potential extreme event. A clustering approach that leverages the hierarchical spectral clustering method is utilized for the optimal defensive islanding. The proposed approach captures the fragility behavior and loading conditions of power system components due to cyber-induced failures. A graphical-based coupling framework is used to map the impacts of cyber failures into operation of power system components. The proposed method is demonstrated on a modified 33-node distribution feeder system integrated with distributed energy resources. The amount of load curtailment and radiality constraints have been used to evaluate the performance of the proposed clustering strategies. The results show the capability of the proposed algorithm to create islands considering the cyber-induced failures for enhanced resilience.

cyber-induced failures↗
DOE OSTI2019

Schottky Anomalies in the Low-Temperature Specific Heat of YBa{sub 2}Cu{sub 3}O{sub y} HTSC

The evolution of the Schottky anomalies in the temperature dependence of the specific heat (in the temperature range T = 2–10 K and the magnetic field range H = 0–9 T) of polycrystalline YBa{sub 2}Cu{sub 3}O{sub y} HTSCs is studied when the oxygen content in coarse-crystalline samples changes and when the degree of structural disordering in optimally doped fine-crystalline samples of the same compound changes. The number of paramagnetic centers with uncompensated spin moments, which generate the Schottky anomalies, and the splitting of the lower energy levels are estimated. A correlation between these parameters and both the number of oxygen vacancies in chain planes and the method of oxygen ordering in chains is found. The experimental data are compared with the results of studying the magnetic susceptibility at T > T{sub c} in the same samples. Based on this comparison, we can conclude that the nature of the Schottky centers is related to the formation of paramagnetic centers at breaks in the Cu1–O4 copper–oxygen chains in the basal planes of the crystal structure.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Arrays of glass wedges for multi-dimensional optical diagnostics

There is a common need in the advancement of optical diagnostic techniques to increase the dimensionality of measurements. For example, point measurements could be improved to multi-point, line, planar, volumetric, or time-resolved volumetric measurements. In this work, a unique optical element is presented to enable multi-dimensional measurements, namely, an array of glass wedges. A light source is passed through the wedges, and different portions of the illumination are refracted by different amounts depending on the glass wedge angle. Subsequent optics can be used to focus the light to multiple points, lines, or planes. Basic characterization of a glass wedge array is presented. Additional wedge-array configurations are discussed, including the use of a periodic intensity mask for multi-planar measurements via structured illumination. The utility of this optical element is briefly demonstrated in (a) multi-planar flame particulate measurements, (b) multi-point femtosecond-laser electronic excitation tagging for flow velocimetry, and (c) multi-line nitric oxide molecular tagging velocimetry in a hypersonic shock-tunnel. One significant advantage of this optical component is its compatibility with high-energy laser sources, which may be a limiting factor with other beam-splitting or beam-forming elements such as some diffractive optics. Additionally, an array of glass wedges is simple and easily customizable compared to other methods for forming multiple closely spaced illumination patterns. In conclusion, suggestions for further development and applications are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1986

Efficient solution methods for the Navier-Stokes equations

Implicit finite difference schemes for solving two-dimensional and three-dimensional Euler and thin layer Navier-Stokes equations are addressed. The methods are demonstrated in fully vectorized codes for a Cray type architecture. The Beam and Warming implicit approximate factorization algorithm in generalized coordinates is used. The methods are either time accurate or accelerated non-time accurate steady state schemes. Acceleration and efficiency modifications such as matrix reduction, diagonalization, and flux split schemes are presented. Two dimensional inviscid and viscous calculations (e.g., airfoils with a deflected spoiler, circulation control airfoils, and unsteady buffeting) and of three dimensional viscous elliptical bodies, exhausting boattails, and generic oblique wing computations are discussed.

Pulliam, T. H.↗
NASA NTRS1991

Compressible flows with periodic vortical disturbances around lifting airfoils

A numerical method is developed for solving periodic, three-dimensional, vortical flows around lifting airfoils in subsonic flow. The first-order method that is presented fully accounts for the distortion effects of the nonuniform mean flow on the convected upstream vortical disturbances. The unsteady velocity is split into a vortical component which is a known function of the upstream flow conditions and the Lagrangian coordinates of the mean flow, and an irrotational field whose potential satisfies a nonconstant-coefficient, inhomogeneous, convective wave equation. Using an elliptic coordinate transformation, the unsteady boundary value problem is solved in the frequency domain on grids which are determined as a function of the Mach number and reduced frequency. The numerical scheme is validated through extensive comparisons with known solutions to unsteady vortical flow problems. In general, it is seen that the agreement between the numerical and analytical results is very good for reduced frequencies ranging from 0 to 4, and for Mach numbers ranging from .1 to .8. Numerical results are also presented for a wide variety of flow configurations for the purpose of determining the effects of airfoil thickness, angle of attack, camber, and Mach number on the unsteady lift and moment of airfoils subjected to periodic vortical gusts. It is seen that each of these parameters can have a significant effect on the unsteady airfoil response to the incident disturbances, and that the effect depends strongly upon the reduced frequency and the dimensionality of the gust. For a one-dimensional (transverse) or two-dimensional (transverse and longitudinal) gust, the results indicate that airfoil thickness increases the unsteady lift and moment at the low reduced frequencies but decreases it at the high reduced frequencies. The results show that an increase in airfoil Mach number leads to a significant increase in the unsteady lift and moment for the low reduced frequencies, but a significant decrease for the high reduced frequencies.

Scott, James R.↗
NASA NTRS1991

A finite-difference, frequency-domain numerical scheme for the solution of the linearized unsteady Euler equations

A numerical method is developed for solving periodic, three-dimensional, vortical flows around lifting airfoils in subsonic flow. The first-order method, that is presented, fully accounts for the distortion effects of the nonuniform mean flow on the convected upstream vortical disturbances. The unsteady velocity is split into a vortical component which is a known function of the upstream flow conditions and the Lagrangian coordinates of the mean flow, and an irrotational field whose potential satisfies a nonconstant-coefficient, inhomogeneous, convective wave equation. Using an elliptic coordinate transformation, the unsteady boundary value problem is solved in the frequency domain on grids which are determined as a function of the Mach number and reduced frequency. Extensive comparisons are made with known solutions to unsteady vortical flow problems, and it is seen that the agreement is generally very good for reduced frequencies ranging from 0 up to 4.

Scott, James R.↗
NASA NTRS2006

Spacecraft hazard avoidance utilizing structured light

At JPL, a <5 kg free-flying micro-inspector spacecraft is being designed for host-vehicle inspection. The spacecraft includes a hazard avoidance sensor to navigate relative to the vehicle being inspected. Structured light was selected for hazard avoidance because of its low mass and cost. Structured light is a method of remote sensing 3-dimensional structure of the proximity utilizing a laser, a grating, and a single regular APS camera. The laser beam is split into 400 different beams by a grating to form a regular spaced grid of laser beams that are projected into the field of view of an APS camera. The laser source and the APS camera are separated forming the base of a triangle. The distance to all beam intersections of the host are calculated based on triangulation.

remote sensing↗
NASA NTRS2016

Dielectric Waveguides Splitter and Hybrid/Isolator for Bidirectional Link

A system, method, device, and apparatus provide a dielectric waveguide splitter/bi-directional link. A dielectric substrate fabricated into a first Y-junction waveguide with a first port splitting into a first branch leading to a second port and a second branch leading to a third port. An angle between the first branch and the second branch is below ninety degrees (90.degree.). The dielectric waveguide splitter enables millimeter-wave (mmWave) transmission between the first port and the second port while reducing feedback of the mmWave between the second and third port. Two Y-junction waveguides may be fabricated back-to-back to provide simultaneous bidirectional mmWave transmission at a single frequency.

Tang, Adrian Joseph↗
DOE OSTI2022

Fast and scalable quantum Monte Carlo simulations of electron-phonon models

We introduce methodologies for highly scalable quantum Monte Carlo simulations of electron-phonon models, and report benchmark results for the Holstein model on the square lattice. The determinant quantum Monte Carlo (DQMC) method is a widely used tool for simulating simple electron-phonon models at finite temperatures, but incurs a computational cost that scales cubically with system size. Alternatively, near-linear scaling with system size can be achieved with the hybrid Monte Carlo (HMC) method and an integral representation of the Fermion determinant. Here, we introduce a collection of methodologies that make such simulations even faster. To combat "stiffness" arising from the bosonic action, we review how Fourier acceleration can be combined with time-step splitting. To overcome phonon sampling barriers associated with strongly-bound bipolaron formation, we design global Monte Carlo updates that approximately respect particle-hole symmetry. To accelerate the iterative linear solver, we introduce a preconditioner that becomes exact in the adiabatic limit of infinite atomic mass. Finally, we demonstrate how stochastic measurements can be accelerated using fast Fourier transforms. Here, these methods are all complementary and, combined, may produce multiple orders of magnitude speedup, depending on model details.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Quantum Monte Carlo Method in the Steady State

We present a numerically exact steady-state inchworm Monte Carlo method for nonequilibrium quantum impurity models. Rather than propagating an initial state to long times, the method is directly formulated in the steady state. This eliminates any need to traverse the transient dynamics and grants access to a much larger range of parameter regimes at vastly reduced computational costs. We benchmark the method on equilibrium Green’s functions of quantum dots in the noninteracting limit and in the unitary limit of the Kondo regime. We then consider correlated materials described with dynamical mean field theory and driven away from equilibrium by a bias voltage. We show that the response of a correlated material to a bias voltage differs qualitatively from the splitting of the Kondo resonance observed in bias-driven quantum dots.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1983

On a finite-difference method for solving transient viscous flow problems

A method has been developed to solve the unsteady, compressible Navier-Stokes equation with the property of consistency and the ability of minimizing the equation stiffness. It relies on innovative extensions of the state-of-the-art finite-difference techniques and is composed of: (1) the upwind scheme for split-flux and the central scheme for conventional flux terms in the inviscid and viscous regions, respectively; (2) the characteristic treatment of both inviscid and viscous boundaries; (3) an ADI procedure compatible with interior and boundary points; and (4) a scalar matrix coefficient including viscous terms. The performance of this method is assessed with four sample problems; namely, a standing shock in the Laval duct, a shock reflected from the wall, the shock-induced boundary-layer separation, and a transient internal nozzle flow. The results from the present method, an existing hybrid block method, and a well-known two-step explicit method are compared and discussed. It is concluded that this method has an optimal trade-off between the solution accuracy and computational economy, and other desirable properties for analyzing transient viscous flow problems.

Li, C. P.↗
NASA NTRS2018

Adjustments and Uncertainty Quantification for SLS Aerodynamic Sectional Loads

This paper presents a method for adjusting sectional loads to match target values for integrated force and moment coefficients. In a typical application, the sectional load profile for one flight condition is calculated from Computational Fluid Dynamics (CFD) while the integrated forces and moments are measured in a wind tunnel experiment. These two methods do not generally result in identical predictions, and this leads to an inherent inconsistency between different data products. This paper aims to provide a procedure to remove that inconsistency. A sectional load profile for a launch vehicle splits the rocket into slices along its length and calculates the aerodynamic loading on each slice, which leads to a one-dimensional aerodynamic load profile that is used for structural analysis. Adjusting sectional loads, also known as line loads, is a nontrivial matter due to several consistency constraints. For example, the adjusted sectional normal force profile must be consistent with both the integrated normal force and pitching moment. To avoid such inconsistency issues, this paper presents a method using a Proper Orthogonal Decomposition (POD) to generate basis functions to adjust the sectional load profiles. As a corollary, this correction method enables the creation of an uncertainty quantification for sectional loads that is consistent with the dispersed integrated force and moment database and its uncertainty quantification. Several extensions to this technique, such as applying the method to the surface pressures, are considered.

Quantification↗
DOE OSTI2020

FFT-based solver for higher-order and multi-phase-field fracture models applied to strongly anisotropic brittle materials

This paper presents the application of a fast Fourier transform (FFT) based method to solve two phase field models designed to simulate crack growth of strongly anisotropic materials in the brittle regime. By leveraging the ability of the FFT-based solver to generate solutions with higher-order and global continuities, we design two simple algorithms to capture the complex fracture patterns (e.g. sawtooth, and curved crack growth) common in materials with strongly anisotropic surface energy via the multi-phase-field and high-order phase-field frameworks. A staggered operator-split solver is used where both the balance of linear momentum and the phase field governing equations are formulated in the periodic domain. The unit phase field of the initial failure region is prescribed by the penalty method to alleviate the sharp material contrast between the initial failure region and the base material. The discrete frequency vectors are generalized to estimate the second and fourth order gradients such that the Gibbs effect near shape interfaces or jump conditions can be suppressed. Furthermore, a preconditioner is adopted to improve the convergence rate of the iterative linear solver. Three numerical experiments are used to systematically compare the performance of the FFT-based method in the multi-phase-field and high-order phase-field frameworks.

42 ENGINEERING↗
NASA NTRS1995

An Extended Ab Initio and Theoretical Thermodynamics Studies of the Bergman Reaction and the Energy Splitting of the Singlet Ortho-, Meta-, and Para-Benzynes

The autoaromatization of (Z)-hex-3-ene-1,5-diyne to the singlet biradical para-benzyne has been reinvestigated by state of the art ab initio methods. Previous CCSD(T)/6-31G(d,p) and CASPT2[0]/ANO[C(5s4p2d1f)/H(3s2p)] calculations estimated the the reaction heat at 298 K to be 8-10 and 4.9 plus or minus 3.2 kcal/mol, respectively. Recent NO- and oxygen-dependent trapping experiments and collision-induced dissociation threshold energy experiments estimate the heat of reaction to be 8.5 plus or minus 1.0 at 470 K (recomputed to 9.5 plus or minus 1.0 at 298 K) and 8.4 plus or minus 3.0 kcal/mol at 298 K, respectively. New theoretical estimates at 298 K predict the values at the basis set limit for the CCSD(T) and CASPT2(g1) methods to be 12.7 plus or minus 2.0 and 5.4 plus or minus 2.0 kcal/mol, respectively. The experimentally predicted electronic contribution to the heat of activation is 28.6 kcal/mol. This can be compared with 25.5 and 29.8 kcal/mol from the CASPT2[g1] and the CCSD(T) methods, respectively. The new study has in particular improved on the one-particle basis set for the CCSD(T) method as compared to earlier studies. For the CASPT2 investigation the better suited CASPT2[g1] approximation is utilized. The original CASPT2 method, CASPT2[0], systematically favors open shell systems relative to closed shell systems. This was previously corrected empirically. The study shows that the energy difference between CCSD(T) and CASPT2[g1] at the basis set limit is estimated to be 7 plus or minus 2 kcal/mol. The study also demonstrates that the estimated heat of reaction is very sensitive to the quality of the basis set.

Lindh, Roland↗