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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

High-resolution Y-axis readout for delay-line microchannel anodes

We have devised a method to obtain a high-resolution Y-axis event position determination from microchannel plate delay-line detectors. The method is based on the double-delay-line wedge-wedge charge partition principle of Lampton et al. (1990), where the X axis is read out by a pair of side-by-side delay lines and the Y-axis coordinate is determined by comparing the charges on the two parallel delay lines. However, our new method abandons the common field return path of the two delay lines, and splits them into two adjacent but independent transformer coupled timing circuits having no common ground connection except via the charge-measurement amplifiers. Thus we eliminate the large electrostatic capacitance of the delay lines from the input of the charge-measuring system. Over its limited working field, the spatial resolution in the Y coordinate can easily equal the resolution provided in the X coordinate by the delay line. Applications to extreme and FUV photon-counting spectroscopy are envisioned.

Raffanti, R.↗
NASA NTRS1993

The piecewise-linear predictor-corrector code - A Lagrangian-remap method for astrophysical flows

We describe a time-explicit finite-difference algorithm for solving the nonlinear fluid equations. The method is similar to existing Eulerian schemes in its use of operator-splitting and artificial viscosity, except that we solve the Lagrangian equations of motion with a predictor-corrector and then remap onto a fixed Eulerian grid. The remap is formulated to eliminate errors associated with coordinate singularities, with a general prescription for remaps of arbitrary order. We perform a comprehensive series of tests on standard problems. Self-convergence tests show that the code has a second-order rate of convergence in smooth, two-dimensional flow, with pressure forces, gravity, and curvilinear geometry included. While not as accurate on idealized problems as high-order Riemann-solving schemes, the predictor-corrector Lagrangian-remap code has great flexibility for application to a variety of astrophysical problems.

Lufkin, Eric A.↗
NASA NTRS1996

A technique for determining cloud free vs cloud contaminated pixels in satellite imagery

Since the first earth orbiting satellite sent pictures of the earth back to them, atmospheric scientists have been focused on the possibilities of using that information as both a forecasting tool and as a meteorological research tool. With the latest generation of Geostationary Operational Environmental Satellites (GOES) now entering service, that view of the earth yields views at a frequency and resolution never before available. These satellites have imagers with a five band multi-spectral capability with high spatial resolution. In addition, the sounder has eighteen thermal infrared (IR) channels plus one low-resolution visible band. With a resolution as small as one kilometer, GOES provides scientists with a powerful eye on the atmosphere. Menzel and Purdom (1994) detail both the imager and sounder capability as well as other systems on the GOES satellites. Immediately apparent in the visible channel are the patterns of clouds swirling over both oceans and continents. These clouds range in size from huge planetary systems covering thousands of kilometers to puffy fair weather cumulus clouds on the order of half a kilometer in size. With the IR sensors temperature patterns are observed. High clouds appear very cold, while low stratus field show temperatures near that of the surface. The surface, in turn, generally appears warmer than the clouds. It would seem then a simple manner to determine cloud and surface temperature from the imagery, but such is not the case. While most of the atmospheric constituents are well mixed and homogeneous, water vapor is not. The water molecule, because of its unique structure and vibration modes, affects the transmittance of the atmosphere most notably in the infrared regions. There are regions of the IR spectrum where water vapor acts as a strong absorber, and at others it is nearly transparent. The transparent wavelengths are called windows, and one such window occurs at 11.2 microns. Adjacent to this window at 12.7 microns which is strongly absorbed by water vapor. These two wavelengths form what is known as a split window, the utility of which was used. Using the linearized form of the radiative transfer equation, they were able to use the split window to determine the amount of water vapor present in the atmosphere. Jedlovec developed the physical split-window (PSW) technique which determines the integrated water content (IWC). The PSW method using Visible Infrared Spin Scan Radiometer (VISSR) Atmospheric Sounder (VAS) found on the older versions of the GOES satellites was used. Recently, Jedlovec and colleagues have been attempting to apply the PSW method using full disk IR imagery obtained by the new generation of GOES satellites. IWC is essential for improved analysis and prediction of convective storms which have been observed to develop in regions of both strong and rapidly evolving moisture gradients. It has also been used in the prediction of clouds and precipitation.

Wohlman, Richard A.↗
NASA NTRS2009

Constraining the Depth of a Martian Magma Ocean through Metal-Silicate Partitioning Experiments: The Role of Different Datasets and the Range of Pressure and Temperature Conditions

Mars accretion is known to be fast compared to Earth. Basaltic samples provide a probe into the interior and allow reconstruction of siderophile element contents of the mantle. These estimates can be used to estimate conditions of core formation, as for Earth. Although many assume that Mars went through a magma ocean stage, and possibly even complete melting, the siderophile element content of Mars mantle is consistent with relatively low pressure and temperature (PT) conditions, implying only shallow melting, near 7 GPa and 2073 K. This is a pressure range where some have proposed a change in siderophile element partitioning behavior. We will examine the databases used for parameterization and split them into a low and higher pressure regime to see if the methods used to reach this conclusion agree for the two sets of data.

Righter, K.↗
NASA NTRS1999

On Multi-Dimensional Unstructured Mesh Adaption

Anisotropic unstructured mesh adaption is developed for a truly multi-dimensional upwind fluctuation splitting scheme, as applied to scalar advection-diffusion. The adaption is performed locally using edge swapping, point insertion/deletion, and nodal displacements. Comparisons are made versus the current state of the art for aggressive anisotropic unstructured adaption, which is based on a posteriori error estimates. Demonstration of both schemes to model problems, with features representative of compressible gas dynamics, show the present method to be superior to the a posteriori adaption for linear advection. The performance of the two methods is more similar when applied to nonlinear advection, with a difference in the treatment of shocks. The a posteriori adaption can excessively cluster points to a shock, while the present multi-dimensional scheme tends to merely align with a shock, using fewer nodes. As a consequence of this alignment tendency, an implementation of eigenvalue limiting for the suppression of expansion shocks is developed for the multi-dimensional distribution scheme. The differences in the treatment of shocks by the adaption schemes, along with the inherently low levels of artificial dissipation in the fluctuation splitting solver, suggest the present method is a strong candidate for applications to compressible gas dynamics.

Wood, William A.↗
NASA NTRS1993

Extension of state-vector splitting to the Navier-Stokes equations

A multi-dimensional upwind scheme for solving the Euler equations has been extended to solve the Navier-Stokes equations. These methods share the property of not requiring body-fitted computational grids. A time advance scheme which is second-order accurate in time and space has been developed for use in conjunction with state-vector splitting to accurately resolve regions with significant gradients in flow quantities. The scheme has been tested on several simple viscous flow problems for which analytic solutions are known; the agreement of computation with theory is excellent.

Gooch, Carl F.↗
DOE OSTI2025

Baryon fraction from the BAO amplitude: a consistent approach to parameterizing perturbation growth

Galaxy clustering constrains the baryon fraction Omega_b/Omega_m through the amplitude of baryon acoustic oscillations and the suppression of perturbations entering the horizon before recombination. This produces a different pre-recombination distribution of baryons and dark matter. After recombination, the gravitational potential responds to both components in proportion to their mass, allowing robust measurement of the baryon fraction. This is independent of new-physics scenarios altering the recombination background (e.g. Early Dark Energy). The accuracy of such measurements does, however, depend on how baryons and CDM are modeled in the power spectrum. Previous template-based splitting relied on approximate transfer functions that neglected part of information. We present a new method that embeds an extra parameter controlling the balance between baryons and dark matter in the growth terms of the perturbation equations in the CAMB Boltzmann solver. This approach captures the baryonic suppression of CDM prior to recombination, avoids inconsistencies, and yields a clean parametrization of the baryon fraction in the linear power spectrum, separating out the simple physics of growth due to the combined matter potential. We implement this framework in an analysis pipeline using Effective Field Theory of Large-Scale Structure with HOD-informed priors and validate it against noiseless LCDM and EDE cosmologies with DESI-like errors. The new scheme achieves comparable precision to previous splitting while reducing systematic biases, providing a more robust way to baryon-fraction measurements. In combination with BBN constraints on the baryon density and Alcock-Paczynski estimates of the matter density, these results strengthen the use of baryon fraction measurements to derive a Hubble constant from energy densities, with future DESI and Euclid data expected to deliver competitive constraints.

Crespi, Andrea [U. Waterloo (main); Waterloo U., I↗
NASA NTRS1978

A split-recoupled-semidirect computational technique applied to transonic flow over lifting airfoils

A new version of the semidirect iterative method eliminates significant restrictions of previous versions of the method. A semidirect method solves finite-difference equations by a rapid globally implicit iterative process driven by a fast direct elliptic solver. The new approach can treat complex systems of equations in an efficient 'correction form', and allows the use of general, nonorthogonal, boundary-fitted coordinate transformations. These features are expected to lead to significant practical applications with conservation-equation systems in either two or three dimensions. The present application to the full potential equations for steady transonic flow over an airfoil at angle of attack illustrates the utility of the technique.

Martin, E. D.↗
DOE OSTI2025

Pooled PPIseq: Screening the SARS-CoV-2 and human interface with a scalable multiplexed protein-protein interaction assay platform

Protein-Protein Interactions (PPIs) are a key interface between virus and host, and these interactions are important to both viral reprogramming of the host and to host restriction of viral infection. In particular, viral-host PPI networks can be used to further our understanding of the molecular mechanisms of tissue specificity, host range, and virulence. At higher scales, viral-host PPI screening could also be used to screen for small-molecule antivirals that interfere with essential viral-host interactions, or to explore how the PPI networks between interacting viral and host genomes co-evolve. Current high-throughput PPI assays have screened entire viral-host PPI networks. However, these studies are time consuming, often require specialized equipment, and are difficult to further scale. Here, we develop methods that make larger-scale viral-host PPI screening more accessible. This approach combines the mDHFR split-tag reporter with the iSeq2 interaction-barcoding system to permit massively-multiplexed PPI quantification by simple pooled engineering of barcoded constructs, integration of these constructs into budding yeast, and fitness measurements by pooled cell competitions and barcode-sequencing. We applied this method to screen for PPIs between SARS-CoV-2 proteins and human proteins, screening in triplicate >180,000 ORF-ORF combinations represented by >1,000,000 barcoded lineages. Our results complement previous screens by identifying 74 putative PPIs, including interactions between ORF7A with the taste receptors TAS2R41 and TAS2R7, and between NSP4 with the transmembrane KDELR2 and KDELR3. We show that this PPI screening method is highly scalable, enabling larger studies aimed at generating a broad understanding of how viral effector proteins converge on cellular targets to effect replication.

60 APPLIED LIFE SCIENCES↗
NASA NTRS2003

Toward Immersed Boundary Simulation of High Reynolds Number Flows

In the immersed boundary (IB) method, the surface of an object is reconstructed with forcing terms in the underlying flow field equations. The surface may split a computational cell removing the constraint of the near wall gridlines to be aligned with the surface. This feature greatly simplifies the grid generation process which is cumbersome and expensive in particular for structured grids and complex geometries. The IB method is ideally suited for Cartesian flow solvers. The flow equations written in Cartesian coordinates appear in a very simple form and several numerical algorithms can be used for an efficient solution of the equations. In addition, the accuracy of numerical algorithms is dependent on the underlying grid and it usually deteriorates when the grid deviates from a Cartesian mesh. The challenge for the IB method lies in the representation of the wall boundaries and in providing an adequate near wall flow field resolution. The issue of enforcing no-slip boundary conditions at the immersed surface has been addressed by several authors by imposing a local reconstruction of the solution. Initial work by Verzicco et al. was based on a simple linear, one-dimensional operator and this approach proved to be accurate for boundaries largely aligned with the grid lines. Majumdar et al. used various multidimensional and high order polynomial interpolations schemes. These high order schemes, however, are keen to introduce wiggles and spurious extrema. Iaccarino & Verzicco and Kalitzin & Iaccarino proposed a tri-linear reconstruction for the velocity components and the turbulent scalars. A modified implementation that has proven to be more robust is reported in this paper. The issue of adequate near wall resolution in a Cartesian framework can initially be addressed by using a non-uniform mesh which is stretched near the surface. In this paper, we investigate an unstructured approach for local grid refinement that utilizes Cartesian mesh features. The computation of high Reynolds number wall bounded flows is particularly challenging as it requires the consideration of thin turbulent boundary layers, i.e. near wall regions with large gradients of the flow field variables. For such flows, the representation of the wall boundary has a large impact on the accuracy of the computation. It is also critical for the robustness and convergence of the flow solver.

Kalitzin, Georgi↗
DOE OSTI2024

Performance of explicit and IMEX MRI multirate methods on complex reactive flow problems within modern parallel adaptive structured grid frameworks

Large-scale multiphysics simulations are computationally challenging due to the coupling of multiple processes with widely disparate time scales. The advent of exascale computing systems exacerbates these challenges since these systems enable ever-increasing size and complexity. In recent years, there has been renewed interest in developing multirate methods as a means to handle the large range of time scales, as these methods may afford greater accuracy and efficiency than more traditional approaches of using implicit-explicit (IMEX) and low-order operator splitting schemes. However, to date there have been few performance studies that compare different classes of multirate integrators on complex application problems. In this work, we study the performance of several newly developed multirate infinitesimal (MRI) methods, implemented in the SUNDIALS solver package, on two reacting flow model problems built on structured mesh frameworks. The first model revisits prior work on a compressible reacting flow problem with complex chemistry that is implemented using BoxLib but where we now include comparisons between a new explicit MRI scheme with the multirate spectral deferred correction (SDC) methods in the original paper. The second problem uses the same complex chemistry as the first problem, combined with a simplified flow model, but runs at a large spatial scale where explicit methods become infeasible due to stability constraints. Two recently developed IMEX MRI multirate methods are tested. These methods rely on advanced features of the AMReX framework on which the model is built, such as multilevel grids and multilevel preconditioners. The results from these two problems show that MRI multirate methods can offer significant performance benefits on complex multiphysics application problems and that these methods may be combined with advanced spatial discretization to compound the advantages of both.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Automated Signal Timing Plan Reconstruction Using High-Resolution Event-Based Controller Data for Digital Twins

Transportation digital twins are essential tools for evaluating emerging technologies such as connected and automated vehicles, adaptive traffic signal control, and mobility optimization strategies. Realistic digital twins require accurate emulation of real-world signal controllers and detailed signal timing plans. However, signal timing plans are often unavailable or difficult to access, forcing researchers and modelers to rely on assumed fixed timings or halt their analysis. To overcome this challenge, we present a method that directly estimates signal timing plan parameters using high-resolution, event-based data from traffic signal controllers. The proposed method extracts key parameters, including cycle length, offset, phase sequence, coordinated phases, phase-specific minimum and maximum green durations, vehicle extensions, and splits under coordination. A rule-based deterministic signal timing reconstruction algorithm based on traffic signal operation rules, such as those outlined in the Signal Timing Manual, is developed and validated. We evaluate this method, which uses high-resolution controller event logs and verified signal timing plans, on 94 signalized intersections in Nashville, Tennessee, demonstrating their ability to generate accurate, simulation-ready signal timing plans for tools such as SUMO and Vissim.

Saroj, Abhilasha [ORNL] (ORCID:0000000191178063)↗
NASA NTRS1979

A new traffic control design method for large networks with signalized intersections

The paper presents a traffic control design technique for application to large traffic networks with signalized intersections. It is shown that the design method adopts a macroscopic viewpoint to establish a new traffic modelling procedure in which vehicle platoons are subdivided into main stream queues and turning queues. Optimization of the signal splits minimizes queue lengths in the steady state condition and improves traffic flow conditions, from the viewpoint of the traveling public. Finally, an application of the design method to a traffic network with thirty-three signalized intersections is used to demonstrate the effectiveness of the proposed technique.

Leininger, G. G.↗
DOE OSTI2021

Joint ptycho-tomography with deep generative priors

Abstract Joint ptycho-tomography is a powerful computational imaging framework to recover the refractive properties of a 3D object while relaxing the requirements for probe overlap that is common in conventional phase retrieval. We use an augmented Lagrangian scheme for formulating the constrained optimization problem and employ an alternating direction method of multipliers (ADMM) for the joint solution. ADMM allows the problem to be split into smaller and computationally more efficient subproblems: ptychographic phase retrieval, tomographic reconstruction, and regularization of the solution. We extend our ADMM framework with plug-and-play (PnP) denoisers by replacing the regularization subproblem with a general denoising operator based on machine learning. While the PnP framework enables integrating such learned priors as denoising operators, tuning of the denoiser prior remains challenging. To overcome this challenge, we propose a denoiser parameter to control the effect of the denoiser and to accelerate the solution. In our simulations, we demonstrate that our proposed framework with parameter tuning and learned priors generates high-quality reconstructions under limited and noisy measurement data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Constraints on the Hubble constant from supernova Refsdal’s reappearance

The Hubble constant measures the expansion rate of the Universe, but different methods give inconsistent values. Kelly et al. studied a supernova that had its light split into multiple images by the gravitational lensing of a foreground galaxy cluster. They combined measurements of the time delay between the multiple images with predictions made by multiple lensing models of the cluster. Furthermore, this approach allowed the authors to make a blinded measurement of the Hubble constant, finding a value that is more consistent with that derived from the cosmic microwave background than with the cosmic distance ladder.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Reflection Optical Imaging to Study Oxygen Evolution Reactions

Discovery of high-efficiency water-splitting catalysts is critical for clean hydrogen fuel generation. Quantitative and reliable chemical imaging methods can play a crucial role in providing critical catalyst efficiency information for materials development and screening. Here, we developed a reflection optical imaging method to map the performance of transition metal oxide catalysts for oxygen evolution reaction (OER). The generated high concentration of oxygen in the electrolyte will decrease the refractive index of the solutions near the electrode and change the reflectivity of optical images. We successfully imaged and measured the OER on cobalt oxide, and the localized current density and overpotential distribution maps have been extracted. The oxidation and reduction of cobalt oxide can be imaged with our reflection optical imaging method as well. These maps showed that the catalytic activities are quite different at different locations, and some regions showed much higher activities. This is due to the non-uniformity created by the catalyst synthesis process and demonstrates the importance of using the optical imaging method to map the localized catalytic efficiency. In addition, the NiFe alloy’s OER catalytic activities have been studied as well. The proposed reflection optical imaging method could find many applications in catalyst activity measurement and screening.

Electrochemistry↗
NASA NTRS1972

Development of automatic through-insulation welding for microelectric interconnections

The capability to automatically route, remove insulation from, and weld small-diameter solid conductor wire is presented. This would facilitate the economical small-quantity production of complex miniature electronic assemblies. An engineering model of equipment having this capability was developed and evaluated. Whereas early work in the use of welded magnet wire interconnections was concentrated on opposed electrode systems, and generally used heat to melt the wire insulation, the present method is based on a concentric electrode system and a wire feed system which splits the insulation by application of pressure prior to welding. The work deals with the design, fabrication, and evaluation testing of an improved version of this concentric electrode system. Two different approaches to feeding the wire to the concentric electrodes were investigated. It was concluded that the process is feasible for the interconnection of complex miniature electronic assemblies.

Arnett, J. C.↗
NASA NTRS1977

Surface photovoltage in exciton absorption range in CdS

The high resolution, intrinsic spectra of surface photovoltage are reported for semiconducting n-type CdS single crystals. At reduced temperatures (120-160 K) the spectra exhibit three sharp maxima due to A, B and C free exciton transitions. Energy positions of these lines and valence band parameters (spin-orbit and crystal field splittings) estimated from surface photovoltage are in good agreement with values obtained by other methods. The excitonic transitions are very sensitive to surface treatment, i.e. polishing, etching, background illumination and surface doping. The mechanism of direct interaction of free excitons with surface states is proposed to explain exciton lines in surface photovoltage.

Morawski, A.↗