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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS1984

Universal single level implicit algorithm for gasdynamics

A single level effectively explicit implicit algorithm for gasdynamics is presented. The method meets all the requirements for unconditionally stable global iteration over flows with mixed supersonic and supersonic zones including blunt body flow and boundary layer flows with strong interaction and streamwise separation. For hyperbolic (supersonic flow) regions the method is automatically equivalent to contemporary space marching methods. For elliptic (subsonic flow) regions, rapid convergence is facilitated by alternating direction solution sweeps which bring both sets of eigenvectors and the influence of both boundaries of a coordinate line equally into play. Point by point updating of the data with local iteration on the solution procedure at each spatial step as the sweeps progress not only renders the method single level in storage but, also, improves nonlinear accuracy to accelerate convergence by an order of magnitude over related two level linearized implicit methods. The method derives robust stability from the combination of an eigenvector split upwind difference method (CSCM) with diagonally dominant ADI(DDADI) approximate factorization and computed characteristic boundary approximations.

Lombard, C. K.↗
NASA NTRS2019

Entropy Stable Method for the Euler Equations Revisited: Central Differencing via Entropy Splitting and SBP

The two decades old high order central differencing via entropy splitting and summation-by-parts (SBP) difference boundary closure of Olsson & Oliger, Gerritsen & Olsson, and Yee et al. (15, 7, 37) is revisited. The objective of this paper is to prove for the first time that the entropy split scheme is an entropy stable method for central differencing with SBP operators for both periodic and non-periodic boundary conditions for nonlinear Euler equations. Standard high order spatial central differencing as well as high order central spatial DRP (dispersion relation preserving) spatial differencing is part of the entropy stable methodology framework. The proof is to replace the spatial derivatives by summation-by-parts (SBP) difference operators in the entropy split form of the equations using the physical entropy of the Euler equations. The numerical boundary closure follows directly from the SBP operator. No additional numerical boundary procedure is required. In contrast, Tadmor-type entropy conserving schemes (31) using mathematical entropies and more recently in (35], do not naturally come with a numerical boundary closure and a generalized SBP operator has to be developed (18). Long time integration of 2D and 3D test cases is included to show the comparison of this efficient entropy stable method with the Tadmor-type of entropy conservative methods. Studies also include the comparison among the three skew-symmetric splittings on their nonlinear stability and accuracy performance without added numerical dissipations for smooth flows. These are, namely, entropy splitting, Ducros et al. splitting and the Kennedy & Grubber splitting.

Sjogreen, Bjorn↗
NASA NTRS1989

Flux-vector splitting for unsteady calculations on dynamic meshes

The method of flux vector splitting used is that of Van Leer. The fluxes split in this manner have the advantage of being continuously differentiable at eigenvalue sign changes and this allows normal shocks to be captured with at most two interior zones, although in practice only one zone is usually observed. The fluxes as originally derived, however did not include the necessary terms appropriate for calculations on a dynamic mesh. The extension of the splitting to include these terms while retaining the advantages of the original splitting is the main purpose of this investigation. In addition, the use of multiple grids to reduce the computer time is investigated. A subiterative procedure to eliminate factorization and linearization error so that larger time steps can be used is also investigated.

Anderson, W. Kyle↗
NASA NTRS2020

Shape Estimation for Elongated Deformable Object using B-spline Chained Multiple Random Matrices Model

In this paper, a B-spline chained multiple random matrix models (RMMs) representation is proposed to model geometric characteristics of an elongated deformable object. The hyper degrees of freedom structure of the elongated deformable object make its shape estimation challenging. Based on the likelihood function of the proposed B-spline chained multiple RMMs, an expectation-maximization (EM) method is derived to estimate the shape of the elongated deformable object. A split and merge method based on the Euclidean minimum spanning tree (EMST) is proposed to provide initialization for the EM algorithm. The proposed algorithm is evaluated for the shape estimation of the elongated deformable objects in scenarios, such as the static rope with various configurations (including configurations with intersection), the continuous manipulation of a rope and a plastic tube, and the assembly of two plastic tubes. The execution time is computed and the accuracy of the shape estimation results is evaluated based on the comparisons between the estimated width values and its ground-truth, and the intersection over union (IoU) metric.

Gang Yao↗
DOE OSTI2024

FY24 Q4 quarterly report: LLNL's node activity within HydroGEN STCH

In this quarterly report, we Computed the driving force for the water splitting reaction as a function of reaction progress and continued Development of the numerical method for water splitting reactions of H 2 production on the particle surface within 3D porous microstructures.

08 HYDROGEN↗
NASA NTRS1981

A preliminary study of relaxation methods for the inviscid conservative gasdynamics equations using flux splitting

Plus minus flux vector split schemes are combined with the multigrid relaxation method to obtain fast rates of iterative (or steady state) convergence. On model nonlifting biconvex airfoil flow, steady state convergence is obtained using multigrid in 40 iteractions for subsonic cases and in 120 iterations for transonic cases. A 65 x 33 stretched grid is used. However, the multigrid scheme used here (and there are numerous variations) sometimes could not be made to work for transonic cases with a large number of supersonic points. Moreover, the multigrid method is only several times faster than the conventional algorithm which itself could perhaps be accelerated by other means. Nethertheless, with proper adjustment the multigrid method will ultimately lead to very fast rates of iterative convergence.

Steger, J. L.↗
NASA NTRS1989

Free-vibration characteristics and correlation of a space station split-blanket solar array

Two methods for studying the free-vibration characteristics of a large split-blanket solar array in a zero-g cantilevered configuration are presented. The zero-g configuration corrresponds to an on-orbit configuration of the Space Station solar array. The first method applies the equations of continuum mechanics to determine the natural frequencies of the array; the second uses the finite element method program, MSC/NASTRAN. The stiffness matrix from the NASTRAN solution was found to be erroneously grounded. The results from the two methods are compared. It is concluded that the grounding does not seriously compromise the solution to the elastic modes of the solar array. However, the correct rigid body modes need to be included to obtain the correct dynamic model.

Carney, Kelly S.↗
NASA NTRS1989

Free-vibration characteristics and correlation of a Space Station split-blanket solar array

Two methods for studying the free-vibration characteristics of a large split-blanket solar array in a zero-g cantilevered configuration are presented. The zero-g configuration corresponds to an on-orbit configuration of the Space Station solar array. The first method applies the equations of continuum mechanics to determine the natural frequencies of the array; the second uses the finite element method program, MSC/NASTRAN. The stiffness matrix from the NASTRAN solution was found to be erroneously grounded. The results from the two methods are compared. It is concluded that the grounding does not seriously compromise the solution to the elastic modes of the solar array. However, the correct rigid body modes need to be icluded to obtain the correct dynamic model.

Carney, Kelly S.↗
NASA NTRS1994

On Bi-Grid Local Mode Analysis of Solution Techniques for 3-D Euler and Navier-Stokes Equations

A procedure is presented for utilizing a bi-grid stability analysis as a practical tool for predicting multigrid performance in a range of numerical methods for solving Euler and Navier-Stokes equations. Model problems based on the convection, diffusion and Burger's equation are used to illustrate the superiority of the bi-grid analysis as a predictive tool for multigrid performance in comparison to the smoothing factor derived from conventional von Neumann analysis. For the Euler equations, bi-grid analysis is presented for three upwind difference based factorizations, namely Spatial, Eigenvalue and Combination splits, and two central difference based factorizations, namely LU and ADI methods. In the former, both the Steger-Warming and van Leer flux-vector splitting methods are considered. For the Navier-Stokes equations, only the Beam-Warming (ADI) central difference scheme is considered. In each case, estimates of multigrid convergence rates from the bi-grid analysis are compared to smoothing factors obtained from single-grid stability analysis. Effects of grid aspect ratio and flow skewness are examined. Both predictions are compared with practical multigrid convergence rates for 2-D Euler and Navier-Stokes solutions based on the Beam-Warming central scheme.

Ibraheem, S. O.↗
NASA NTRS1996

A Method to Analyze and Optimize the Load Sharing of Split Path Transmissions

Split-path transmissions are promising alternatives to the common planetary transmissions for rotorcraft. Heretofore, split-path designs proposed for or used in rotorcraft have featured load-sharing devices that add undesirable weight and complexity to the designs. A method was developed to analyze and optimize the load sharing in split-path transmissions without load-sharing devices. The method uses the clocking angle as a design parameter to optimize for equal load sharing. In addition, the clocking angle tolerance necessary to maintain acceptable load sharing can be calculated. The method evaluates the effects of gear-shaft twisting and bending, tooth bending, Hertzian deformations within bearings, and movement of bearing supports on load sharing. It was used to study the NASA split-path test gearbox and the U.S. Army's Comanche helicopter main rotor gearbox. Acceptable load sharing was found to be achievable and maintainable by using proven manufacturing processes. The analytical results compare favorably to available experimental data.

Krantz, Timothy L.↗
NASA NTRS2019

Entropy Stable Method for the Euler Equations Revisited: Central Differencing via Entropy Splitting and SBP

The two decades old high order central differencing via entropy splitting and summation-by-parts (SBP) difference boundary closure of Ols- son & Oliger (1994), Gerritsen & Olsson (1996), and Yee et al. (2000) is revisited. The entropy splitting is a form of skew-symmetric splitting of the nonlinear Euler flux derivatives. Central differencing applied to the entropy splitting form of the Euler flux derivatives together with SBP difference operators will, hereafter, be referred to as entropy split schemes. This study is prompted by the recent growing interest in numerical methods for which a discrete entropy conservation law holds, a discrete global entropy conservation can be proved and/or the numerical method possesses a stable entropy in the framework of SBP difference operators and L2-energy norm estimate. The objective of this paper is to recast the entropy split scheme as the re- cent definition of an entropy stable method for central differencing with SBP operators for both periodic and non-periodic boundary conditions for non- linear Euler equations. Standard high order spatial central differencing as well as high order central spatial DRP (dispersion relation preserving) spatial differencing is part of the entropy stable methodology framework. Long time integration of 2D and 3D test cases is included to show the comparison of this efficient entropy stable method with the Tadmor-type of entropy conservative methods. Studies also include the comparison among the three skew-symmetric splittings on their nonlinear stability and accuracy performance without added numerical dissipations for smooth flows. These are, namely, entropy splitting, Ducros et al. splitting and the Kennedy & Grub- ber splitting.

Sjogreen, Bjorn↗
DOE OSTI2020

HTGR Simulation Methods & International Collaborations

ART-GCR “Methods” activity is split between Experimental Validation data from the ANL NSTF and OSU HTTF (next three presentations). HTGR core simulation (this presentation). International collaboration within OECD Generation-IV (Gen-IV) and USA/Japan bi-lateral agreements (this presentation) HTGR Simulation Methods No new NE-52 funding for HTGR Methods support in FY20; ~$200K FY19 carry-over funds only. Consists of international code-to-code benchmarks (IAEA CRP on HTGR UAM and OECD/NEA MHTGR-350) and refinement of a few-group Pebble Bed Reactor (PBR) cross section (XS) generation methodology. Funding will be requested in FY21 to produce the final reports for the two benchmarks and continue the development of the PBR XS generation methodology. Additional (non-ART) HTGR-related support work at INL NEAMS: HTR-Application work package at INL Create a benchmark for the pebble shuffling and depletion algorithms being developed for NEAMS Griffin code. iFOA award with X-Energy: Develop independent Monte Carlo model of Xe-100 design. Independent design confirmatory analysis of Xe-100 design using NEAMS tools Griffin and Pronghorn. Support X-Energy design team to use their own legacy design tools (VSOP99 and MGT). Support NEAMS Griffin and Pronghorn development team for the iFOA needs (received $50K additional funding for required development).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1981

Eliminating gaps in split rings

Simple installation method allows thinner, lighter tether rings than conventional procedures, saving expensive materials. Installer inverts ring with pliers before it is slid over cable, then returns it to its original position after installation. Ring is in correct orientation, and coils are tightly compressed for high reliability fastening.

Gould, R. W.↗
DOE OSTI2024

Numerical coupling of aerosol emissions, dry removal, and turbulent mixing in the E3SM Atmosphere Model version 1 (EAMv1) – Part 2: A semi-discrete error analysis framework for assessing coupling schemes

Abstract. Part 1 (Wan et al., 2024) of this study discusses the motivation and empirical evaluation of a revision to the aerosol-related numerical process coupling in the atmosphere component of the Energy Exascale Earth System Model version 1 (EAMv1) to address the previously reported issue of strong sensitivity of the simulated dust aerosol lifetime and dry removal rate to the model's vertical resolution. This paper complements that empirical justification of the revised scheme with a mathematical justification leveraging a semi-discrete analysis framework for assessing the splitting error of process coupling methods. The framework distinguishes the error due to numerical splitting from the error due to the time integration method(s) used for each individual process. Such a distinction results in a framework that provides an intuitive understanding of the causes of the splitting error. The application of this framework to the dust life cycle in EAMv1 confirms (i) that the original EAMv1 scheme artificially strengthens the effect of dry removal processes and (ii) that the revised splitting reduces that artificial strengthening. While the error analysis framework is presented in the context of the dust life cycle in EAMv1, the framework can be broadly leveraged to evaluate process coupling schemes, both in other physical problems and for any number of processes. This framework will be particularly powerful when the various process implementations support a variety of time integration approaches. Whereas traditional local truncation error approaches require separate consideration of each combination of time integration methods, this framework enables evaluation of coupling schemes independent of particular time integration approaches for each process while still allowing for the incorporation of these specific time integration errors if so desired. The framework also explains how the splitting error terms result from (i) the integration of individual processes in isolation from other processes and (ii) the choices of input state and time step size for the isolated integration of processes. Such a perspective has the potential for the rapid development of alternative coupling approaches that utilize knowledge both about the desired accuracy and about the computational costs of individual processes.

58 GEOSCIENCES↗
NASA NTRS1988

Accurate solutions, parameter studies and comparisons for the Euler and potential flow equations

Parameter studies are conducted using the Euler and potential flow equation models for steady and unsteady flows in both two and three dimensions. The Euler code is an implicit, upwind, finite volume code which uses the Van Leer method of flux vector splitting which has been recently extended for use on dynamic meshes and maintain all the properties of the original splitting. The potential flow code is an implicit, finite difference method for solving the transonic small disturbance equations and incorporates both entropy and vorticity corrections into the solution procedures thereby extending its applicability into regimes where shock strength normally precludes its use. Parameter studies resulting in benchmark type calculations include the effects of spatial and temporal refinement, spatial order of accuracy, far field boundary conditions for steady flow, frequency of oscillation, and the use of subiterations at each time step to reduce linearization and factorization errors. Comparisons between Euler and potential flow results are made, as well as with experimental data where available.

Anderson, W. Kyle↗
NASA NTRS1981

Temporal and spatial inconsistencies of time-split finite-difference schemes

The properties of an implicit time-split algorithm, which utilizes locally one dimensional spatial steps, are examined using the two-dimensional heat conduction equation as the test problem. Both temporal and spatial inconsistencies inherent in the scheme are identified. A consistent, implicit splitting approach is developed. The relationship between this method and other time-split implicit schemes is explained, and stability problems encountered with the method in three dimensions are discussed.

Dwoyer, D. L.↗
NASA NTRS1994

A New Global Potential Energy Surface for the Hydroperoxyl Radical, HO2: Reaction Coefficients for H + O2 and Vibrational Splittings for H Atom Transfer

A new analytic global potential energy surface describing the hydroperoxyl radical system H((sup 2)S) + O2(X (sup 3)Sigma((sup -)(sub g))) (reversible reaction) HO2 ((X-tilde) (sup 2)A'') (reversible reaction) O((sup 3)P) + O H (X (sup 2)Pi) has been fitted using the ab initio complete active space SCF (self-consistent-field)/externally contracted configuration interaction (CASSCF/CCI) energy calculations of Walch and Duchovic. Results of quasiclassical trajectory studies to determine the rate coefficients of the forward and reverse reactions at combustion temperatures will be presented. In addition, vibrational energy levels were calculated using the quantum DVR-DGB (discrete variable representation-distributed Gaussian basis) method and the splitting due to H atom migration is investigated. The material of the proposed presentation was reviewed and the technical content will not reveal any information not already in the public domain and will not give any foreign industry or government a competitive advantage.

Dateo, Christopher E.↗
DOE OSTI2024

Robust wind farm layout optimization

Wake interactions in wind farms cause losses in annual energy production (AEP) on the order of 10%. Wind farm designers optimize the layout of the farm to mitigate wake losses, especially in the dominant site-specific wind directions. As wind turbines and wind farms grow in scale, optimization becomes more complex. Offshore wind farms regularly comprise more than 100 wind turbines and are characterized by complex boundaries due to shipping lanes, neighboring wind farms, and other constraints. Layout optimization methods are broadly split between gradient-based and gradient-free approaches. Gradient-based approaches can converge quickly and perform well for smaller, academic problems but are often sensitive to initial conditions and tuning parameters and require expert knowledge to use. On the other hand, gradient-free approaches can be more robust to problem complexities. We present a robust layout optimization approach based on a random search algorithm. The algorithm is intended for those who are not optimization experts and has few tuning parameters that need specification to achieve satisfactory results. Unlike off-the-shelf methods, which use generally available, non-domain-specific optimization routines that accept as inputs an optimization function and constraint definitions, this approach takes advantage of the relative computational costs of the different evaluations by evaluating cheaper computations first (boundary and minimum distance constraints) and running expensive AEP evaluations only if all other checks pass. Moreover, an outer genetic algorithm allows multiple solutions to evolve in parallel, enabling rapid solution development on high-performance computers. We discuss the relative ease of selecting necessary tuning parameters and demonstrate the efficacy of the genetic random search on a complex layout problem consisting of placing 70 turbines in a nonconvex and unconnected boundary region.

17 WIND ENERGY↗