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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

Application of a self-adaptive grid method to complex flows

A directional-split, modular, user-friendly grid point distribution code is applied to several test problems. The code is self-adaptive in the sense that grid point spacing is determined by user-specified constants denoting maximum and minimum grid spacings and constants relating the relative influence of smoothness and orthogonality. Estimates of truncation error, in terms of flow-field gradients and/or geometric features, are used to determine the point distribution. Points are redistributed along grid lines in a specified direction in an elliptic manner over a user-specified subdomain, while orthogonality and smoothness are controlled in a parabolic (marching) manner in the remaining directions. Multidirectional adaption is achieved by sequential application of the method in each coordinate direction. The flow-field solution is redistributed onto the newly distributed grid points after each unidirectional adaption by a simple one-dimensional interpolation scheme. For time-accurate schemes such interpolation is not necessary and time-dependent metrics are carried in the fluid dynamic equations to account for grid movement.

Deiwert, G. S.↗
DOE OSTI2025

A Parameter-masked Mock Data Challenge for Beyond-two-point Galaxy Clustering Statistics

The past few years have seen the emergence of a wide array of novel techniques for analyzing high-precision data from upcoming galaxy surveys, which aim to extend the statistical analysis of galaxy clustering data beyond the linear regime and the canonical two-point (2pt) statistics. We test and benchmark some of these new techniques in a community data challenge named “Beyond-2pt,” initiated during the Aspen 2022 Summer Program “Large-Scale Structure Cosmology beyond 2-Point Statistics,” whose first round of results we present here. The challenge data set consists of high-precision mock galaxy catalogs for clustering in real space, in redshift space, and on a light cone. Participants in the challenge have developed end-to-end pipelines to analyze mock catalogs and extract unknown (“masked”) cosmological parameters of the underlying ΛCDM models with their methods. The methods represented are density-split clustering, nearest neighbor statistics, BACCO power spectrum emulator, void statistics, LEFTfield field-level inference using effective field theory (EFT), and joint power spectrum and bispectrum analyses using both EFT and simulation-based inference. In this work, we review the results of the challenge, focusing on problems solved, lessons learned, and future research needed to perfect the emerging beyond-2pt approaches. The unbiased parameter recovery demonstrated in this challenge by multiple statistics and the associated modeling and inference frameworks supports the credibility of cosmology constraints from these methods. The challenge data set is publicly available, and we welcome future submissions from methods that are not yet represented.

Krause, Elisabeth [Univ. of Arizona, Tucson, AZ (U↗
DOE OSTI2024

Two splits, three ways: advances in double splitting quenches

In this work we introduce a method for calculating holographic duals of BCFTs with more than two boundaries. We apply it to calculating the dynamics of entanglement entropy in a 1+1d CFT that is instantaneously split into multiple segments and calculate the entanglement entropy as a function of time for the case of two splits, showing that our approach reproduces earlier results for the double split case. Our manuscript lays the groundwork for future calculations of the entanglement entropy for more than two splits and systems at nonzero temperature.

AdS-CFT correspondence↗
NASA NTRS1989

Numerical Simulation of Flow Through an Artificial Heart

A solution procedure was developed that solves the unsteady, incompressible Navier-Stokes equations, and was used to numerically simulate viscous incompressible flow through a model of the Pennsylvania State artificial heart. The solution algorithm is based on the artificial compressibility method, and uses flux-difference splitting to upwind the convective terms; a line-relaxation scheme is used to solve the equations. The time-accuracy of the method is obtained by iteratively solving the equations at each physical time step. The artificial heart geometry involves a piston-type action with a moving solid wall. A single H-grid is fit inside the heart chamber. The grid is continuously compressed and expanded with a constant number of grid points to accommodate the moving piston. The computational domain ends at the valve openings where nonreflective boundary conditions based on the method of characteristics are applied. Although a number of simplifing assumptions were made regarding the geometry, the computational results agreed reasonably well with an experimental picture. The computer time requirements for this flow simulation, however, are quite extensive. Computational study of this type of geometry would benefit greatly from improvements in computer hardware speed and algorithm efficiency enhancements.

Rogers, Stuart E.↗
DOE OSTI2021

Computation-guided optimization of split protein systems

Splitting bioactive proteins into conditionally reconstituting fragments is a powerful strategy for building tools to study and control biological systems. However, split proteins often exhibit a high propensity to reconstitute, even without the conditional trigger, limiting their utility. Current approaches for tuning reconstitution propensity are laborious, context-specific or often ineffective. Here, we report a computational design strategy grounded in fundamental protein biophysics to guide experimental evaluation of a sparse set of mutants to identify an optimal functional window. We hypothesized that testing a limited set of mutants would direct subsequent mutagenesis efforts by predicting desirable mutant combinations from a vast mutational landscape. This strategy varies the degree of interfacial destabilization while preserving stability and catalytic activity. We validate our method by solving two distinct split protein design challenges, generating both design and mechanistic insights. Furthermore, this new technology will streamline the generation and use of split protein systems for diverse applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1993

A multigrid method for steady Euler equations on unstructured adaptive grids

A flux-difference splitting type algorithm is formulated for the steady Euler equations on unstructured grids. The polynomial flux-difference splitting technique is used. A vertex-centered finite volume method is employed on a triangular mesh. The multigrid method is in defect-correction form. A relaxation procedure with a first order accurate inner iteration and a second-order correction performed only on the finest grid, is used. A multi-stage Jacobi relaxation method is employed as a smoother. Since the grid is unstructured a Jacobi type is chosen. The multi-staging is necessary to provide sufficient smoothing properties. The domain is discretized using a Delaunay triangular mesh generator. Three grids with more or less uniform distribution of nodes but with different resolution are generated by successive refinement of the coarsest grid. Nodes of coarser grids appear in the finer grids. The multigrid method is started on these grids. As soon as the residual drops below a threshold value, an adaptive refinement is started. The solution on the adaptively refined grid is accelerated by a multigrid procedure. The coarser multigrid grids are generated by successive coarsening through point removement. The adaption cycle is repeated a few times. Results are given for the transonic flow over a NACA-0012 airfoil.

Riemslagh, Kris↗
NASA NTRS1979

A hybridized method for computing high-Reynolds-number hypersonic flow about blunt bodies

A hybridized method for computing the flow about blunt bodies is presented. In this method the flow field is split into its viscid and inviscid parts. The forebody flow field about a parabolic body is computed. For the viscous solution, the Navier-Stokes equations are solved on orthogonal parabolic coordinates using explicit finite differencing. The inviscid flow is determined by using a Moretti type scheme in which the Euler equations are solved, using explicit finite differences, on a nonorthogonal coordinate system which uses the bow shock as an outer boundary. The two solutions are coupled along a common data line and are marched together in time until a converged solution is obtained. Computed results, when compared with experimental and analytical results, indicate the method works well over a wide range of Reynolds numbers and Mach numbers.

Weilmuenster, K. J.↗
DOE OSTI2021

A Hybrid Finite-Volume, Discontinuous Galerkin Discretization for the Radiative Transport Equation

In this report we propose a hybrid spatial discretization for the radiative transport equation that combines a second-order discontinuous Galerkin (DG) method and a second-order finite-volume (FV) method. The strategy relies on a simple operator splitting that has been used previously to combine different angular discretizations. Unlike standard FV methods with upwind fluxes, the hybrid approach is able to accurately simulate problems in scattering dominated regimes. However, it requires less memory and yields a faster computational time than a uniform DG discretization. In addition, the underlying splitting allows naturally for hybridization in both space and angle. Numerical results are given to demonstrate the efficiency of the hybrid approach in the context of discrete ordinate angular discretizations and Cartesian spatial grids.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2004

Improved Phase-Mask Fabrication of Fiber Bragg Gratings

An improved method of fabrication of Bragg gratings in optical fibers combines the best features of two prior methods: one that involves the use of a phase mask and one that involves interference between the two coherent laser beams. The improved method affords flexibility for tailoring Bragg wavelengths and bandwidths over wide ranges. A Bragg grating in an optical fiber is a periodic longitudinal variation in the index of refraction of the fiber core. The spatial period (Bragg wavelength) is chosen to obtain enhanced reflection of light of a given wavelength that would otherwise propagate relatively unimpeded along the core. Optionally, the spatial period of the index modulation can be made to vary gradually along the grating (such a grating is said to be chirped ) in order to obtain enhanced reflection across a wavelength band, the width of which is determined by the difference between the maximum and minimum Bragg wavelengths. In the present method as in both prior methods, a Bragg grating is formed by exposing an optical fiber to an ultraviolet-light interference field. The Bragg grating coincides with the pattern of exposure of the fiber core to ultraviolet light; in other words, the Bragg grating coincides with the interference fringes. Hence, the problem of tailoring the Bragg wavelength and bandwidth is largely one of tailoring the interference pattern and the placement of the fiber in the interference pattern. In the prior two-beam interferometric method, a single laser beam is split into two beams, which are subsequently recombined to produce an interference pattern at the location of an optical fiber. In the prior phase-mask method, a phase mask is used to diffract a laser beam mainly into two first orders, the interference between which creates the pattern to which an optical fiber is exposed. The prior two-beam interferometric method offers the advantage that the period of the interference pattern can be adjusted to produce gratings over a wide range of Bragg wavelengths, but offers the disadvantage that success depends on precise alignment and high mechanical stability. The prior phase-mask method affords the advantages of compactness of equipment and relative insensitivity to both misalignment and vibration, but does not afford adjustability of the Bragg wavelength. The present method affords both the flexibility of the prior two-beam interferometric method and the compactness and stability of the prior phase-mask method. In this method (see figure), a laser beam propagating along the x axis is normally incident on a phase mask that lies in the (y,z) plane. The phase of light propagating through the mask is modulated with a spatial periodicity, p, along the y axis chosen to diffract the laser light primarily to first order at the angle . (The zero-order laser light propagating along the x axis can be used for alignment and thereafter suppressed during exposure of the fiber.) The diffracted light passes through a concave cylindrical lens, which converts the flat diffracted wave fronts to cylindrical ones, as though the light emanated from a line source. Then two parallel flat mirrors recombine the diffracted beams to form an interference field equivalent to that of two coherent line sources at positions A and B (virtual sources). The interference pattern is a known function of the parameters of the apparatus and of position (x,y) in the interference field. Hence, the tilt, wavelength, and chirp of the Bragg grating can be chosen through suitable adjustments of the apparatus and/or of the position and orientation of the optical fiber. In particular, the Bragg wavelength can be adjusted by moving the fiber along the x axis, and the bandwidth can be modified over a wide range by changing the fiber tilt angle or by moving the phase mask and/or the fiber. Alignment is easy because the zero-order beam defines the x axis. The interference is relatively stable and insensitive to the mechanical vibration because of the gh symmetry and compactness of the apparatus, the fixed positions of the mirrors and lens, and the consequent fixed positions of the two virtual line sources, which are independent of the translations of the phase mask and the laser relative to the lens.

Grant, Joseph↗
NASA NTRS1989

Application of a self-adaptive grid method to complex flows

A directional-split modular user-friendly grid-point distribution code is applied to several test problems. Grid-point spacing is determined by user-specified constants denoting maximum and minimum grid spacings and relating the relative influence of smoothness and orthogonality. Estimates of truncation error, in terms of flowfield gradients and/or geometric features, are used to determine the point distribution.

Deiwert, G. S.↗
NASA NTRS1938

Calculation of the chordwise load distribution over airfoil sections with plain, split, or serially hinged trailing-edge flaps

A method is presented for the rapid calculation of the incremental chordwise normal-force distribution over an airfoil section due to the deflection of a plain flap or tab, a split flap, or a serially hinged flap. This report is intended as a supplement to NACA Report no. 631, wherein a method is presented for the calculation of the chordwise normal-force distribution over an airfoil without a flap or, as it may be considered, an airfoil with flap (or flaps) neutral. The method enables the determination of the form and magnitude of the incremental normal-force distribution to be made for an airfoil-flap combination for which the section characteristics have been determined. A method is included for the calculation of the flap normal-force and hinge-moment coefficients without necessitating a determination of the normal-force distribution.

Allen, H Julian↗
NASA NTRS1988

Choice of implicit and explicit operators for the upwind differencing method

The flux-vector and flux-difference splittings of Steger-Warming, van Leer and Roe are tested in all possible combinations on the implicit and explicit operators that can be distinguished in implicit relaxation methods for the steady Euler and Navier-Stokes equations. The tests include one-dimensional inviscid nozzle flow, and two-dimensional inviscid and viscous shock reflection. Roe's splitting, as anticipated, is found to uniformly yield the most accurate results. On the other hand, an approximate Roe splitting of the implicit operator (the complete Roe splitting is too complicated for practical use) proves to be the least robust with regard to convergence to the steady state. In this respect, the Steger-Warming splitting is the most robust; it leads to convergence when combined with any of the splittings in the explicit operator, although not necessarily in the most efficient way.

Liou, Meng-Sing↗
NASA NTRS1988

Choice of implicit and explicit operators for the upwind differencing method

The flux-vector and flux-difference splittings of Steger-Warming, Van Leer and Roe are tested in all possible combinations in the implicit and explicit operators that can be distinguished in implicit relaxation methods for the steady Euler and Navier-Stokes equations. The tests include one-dimensional inviscid nozzle flow, and two-dimensional inviscid and viscous shock reflection. Roe's splitting, as anticipated, is found to uniformly yield the most accurate results. On the other hand, an approximate Roe splitting of the implicit operator (the complete Roe splitting is too complicated for practical use) proves to be the least robust with regard to convergence to the steady state. In this respect, the Steger-Warming splitting is the most robust: it leads to convergence when combined with any of the splittings in the explicit operator, although not necessarily in the most efficient way.

Liou, Meng-Sing↗
NASA NTRS1978

The CLASSY clustering algorithm: Description, evaluation, and comparison with the iterative self-organizing clustering system (ISOCLS)

A clustering method, CLASSY, was developed, which alternates maximum likelihood iteration with a procedure for splitting, combining, and eliminating the resulting statistics. The method maximizes the fit of a mixture of normal distributions to the observed first through fourth central moments of the data and produces an estimate of the proportions, means, and covariances in this mixture. The mathematical model which is the basic for CLASSY and the actual operation of the algorithm is described. Data comparing the performances of CLASSY and ISOCLS on simulated and actual LACIE data are presented.

Lennington, R. K.↗
NASA NTRS1987

An implicit flux-split algorithm for the compressible Navier-Stokes equations

An implicit upwind scheme for the compressible Navier-Stokes equations is described and applied to the internal flow in a dual-throat nozzle. The method is second-order accurate spatially and naturally dissipative. A spatially-split approximate factorization method is used to obtain efficient steady-state solutions on the NASA Langley VPS-32 (CYBER 205) supercomputer.

Thomas, James L.↗
DOE OSTI2024

The 1 → 3 massive splitting functions from QCD factorization and SCET

Splitting functions are universal functions describing the collinear dynamics of gauge theories, and as such are crucial ingredients for a wide variety of calculations in perturbative QCD. We present analytic results for the triple collinear splitting functions in QCD with a single massive parton. We derive the splitting functions using two distinct methods; first by expanding the squared matrix elements in the collinear limit, and secondly by using soft-collinear effective theory with massive quarks. We find agreement between these two approaches, providing a strong check of our results. Additionally, we also check all iterated and soft limits of our results, finding agreement with predictions from factorization. Our results provide an important ingredient for higher order perturbative calculations involving massive partons, and for the description of the collinear dynamics of heavy flavor jets.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

QTRAJ 1.0: A Lindblad equation solver for heavy-quarkonium dynamics

Here we introduce an open-source package called QTraj that solves the Lindblad equation for heavy-quarkonium dynamics using the quantum trajectories algorithm. The package allows users to simulate the suppression of heavy-quarkonium states using externally-supplied input from 3+1D hydrodynamics simulations. The code uses a split-step pseudo-spectral method for updating the wave-function between jumps, which is implemented using the open-source multi-threaded FFTW3 package. This allows one to have manifestly unitary evolution when using real-valued potentials. In this paper, we provide detailed documentation of QTraj 1.0, installation instructions, and present various tests and benchmarks of the code.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗