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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,099 records · Page 61

All sources NASA NTRS DOE OSTI

NASA NTRS1994

Preconditioned conjugate gradient methods for the Navier-Stokes equations

A preconditioned Krylov subspace method (GMRES) is used to solve the linear systems of equations formed at each time-integration step of the unsteady, two-dimensional, compressible Navier-Stokes equations of fluid flow. The Navier-Stokes equations are cast in an implicit, upwind finite-volume, flux-split formulation. Several preconditioning techniques are investigated to enhance the efficiency and convergence rate of the implicit solver based on the GMRES algorithm. The superiority of the new solver is established by comparisons with a conventional implicit solver, namely line Gauss-Seidel relaxation (LGSR). Computational test results for low-speed (incompressible flow over a backward-facing step at Mach 0.1), transonic flow (trailing edge flow in a transonic turbine cascade), and hypersonic flow (shock-on-shock interactions on a cylindrical leading edge at Mach 6.0) are presented. For the Mach 0.1 case, overall speedup factors of up to 17 (in terms of time-steps) and 15 (in terms of CPU time on a CRAY-YMP/8) are found in favor of the preconditioned GMRES solver, when compared with the LGSR solver. The corresponding speedup factors for the transonic flow case are 17 and 23, respectively. The hypersonic flow case shows slightly lower speedup factors of 9 and 13, respectively. The study of preconditioners conducted in this research reveals that a new LUSGS-type preconditioner is much more efficient than a conventional incomplete LU-type preconditioner.

Ajmani, Kumud↗

NASA NTRSDate not supplied

Discretization Error Estimation and Control for Farfield Acoustic Signatures

We investigate the utility of adjoint-based error estimates for sonic boom farfield simulations governed by solutions of the augmented Burgers’ equation. Solution of this nonlinear system uses operator splitting with a second-order finite volume discretization in space and second-order Runge-Kutta time marching, while the absorption and molecular relaxation are solved using second-order central differencing. The discretization error in selected ground sonic boom cost functionals is estimated using the method of adjoint-weighted residuals. Key elements of the implementation process are emphasized with details provided on the practical aspects as appliedto the sonic boom farfield propagation. We establish the accuracy of the adjoint solutions usingcomplex step and finite difference approaches, and examine the accuracy of the error estimates using analytical N-wave solutions. We then apply it to a pressure waveform corresponding to the X-59 research aircraft. The investigations demonstrate that the method of adjoint-weighted residuals accurately predicts the level of discretization error present in sonic boom farfield simulations while offering insight into which features of the near field signal are the primary drivers of ground noise metrics. The numerical results indicate that at sampling frequencies as low as50kHz, discretization error in the propagation is under 0.01 dB[A] for realistically complex examples.

CST↗

DOE OSTI2024

Pressure stability in explicitly coupled simulations of poromechanics with application to CO 2 sequestration

We study in detail the pressure stabilizing effects of the non-iterated fixed-stress splitting in poromechanical problems which are nearly undrained and incompressible. When applied in conjunction with a spatial discretization which does not satisfy the discrete inf–sup condition, namely a mixed piecewise linear–piecewise constant spatial discretization, the explicit fixed-stress scheme can have a pressure stabilizing effect in transient problems. This effect disappears, however, upon time step refinement or the attainment of steady state. The interpretation of the scheme as an Augmented Lagrangian method similar to Uzawa iteration for incompressible flow helps explain these results. Moreover, due to the slowly evolving solution within undrained seal regions, we show that the explicit fixed-stress scheme requires very large time steps to reveal its pressure stabilizing effect in examples of geologic CO 2 sequestration. We note that large time steps can result in large errors in drained regions, such as the aquifer or reservoir regions of these examples, and can prevent convergence of nonlinear solvers in the case of multiphase flows, which can make the explicit scheme an unreliable source of pressure stabilization. We conclude by demonstrating that pressure jump stabilization is as effective in the explicit fixed-stress setting as in the fully implicit setting for undrained problems, while maintaining the stability and convergence of the fixed-stress split for drained problems.

58 GEOSCIENCES↗

DOE OSTI2025

Asymptotic-preserving dynamical low-rank method for the stiff nonlinear Boltzmann equation

In kinetic theory, numerically solving the full Boltzmann equation is extremely expensive. This is because the Boltzmann collision operator involves a high-dimensional, nonlinear integral that must be evaluated at each spatial grid point and every time step. The challenge becomes even more pronounced in the fluid (strong collisionality) regime, where the collision operator exhibits strong stiffness, causing explicit time integrators to impose severe stability restrictions. In this paper, we propose addressing this problem through a dynamical low-rank (DLR) approximation. The resulting algorithm requires evaluating the Boltzmann collision operator only r 2 times, where r, the rank of the approximation, is much smaller than the number of spatial grid points. We propose a novel DLR integrator, called the XL integrator, which reduces the number of steps compared to the available alternatives (such as the projector splitting or basis update & Galerkin (BUG) integrator). For a class of problems including the Boltzmann collision operator which enjoys a separation property between physical and velocity space, we further propose a specialized version of the XL integrator, called the sXL integrator. This version requires solving only one differential equation to update the low-rank factors. Furthermore, the proposed low-rank schemes are asymptotic-preserving, meaning they can capture the asymptotic fluid limit in the case of strong collisionality. Our numerical experiments demonstrate the efficiency and accuracy of the proposed methods across a wide range of regimes, from non-stiff (kinetic) to stiff (fluid).

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2018

Fourier amplitude sensitivity test applied to dynamic combined finite-discrete element methods-based simulations

Fracture propagation plays a key role for a number of applications of interest to the scientific community, from dynamic fracture processes like spallation and fragmentation in metals to failure of ceramics, airplane wings, etc. Simulations of material deformation and fracture propagation rely on accurate knowledge of material characteristics such as material strength and the amount of energy being dissipated during the fracture process. Within the combined finite-discrete element method (FDEM) framework material fracture behavior is typically described through a parametrized softening curve, which defines a stress-strain relationship unique to each material. We apply the Fourier amplitude sensitivity test to explore how each of these parameters influences the simulated damage processes and to determine the key input parameters that have the most impact on the model response. We present several sensitivity numerical experiments for the simulation of a split Hopkinson pressure bar (SHPB) test for weathered granite samples using different combinations of model parameters. We validate the obtained results against SHPB experimental data. The experiments show that the model is mostly sensitive to parameters related to tensile and shear strengths, even in the presence of other parameter perturbations. Lastly, the results suggest that the specification of tensile and shear strengths at the interfaces dominate the stress-time history of the FDEM simulation of SHPB test.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Hydrogen Production and Li-Ion Battery Performance with MoS2-SiNWs-SWNTs@ZnONPs Nanocomposites

This study explores the hydrogen generation potential via water-splitting reactions under UV-vis radiation by using a synergistic assembly of ZnO nanoparticles integrated with MoS2, single-walled carbon nanotubes (SWNTs), and crystalline silicon nanowires (SiNWs) to create the MoS2-SiNWs-SWNTs@ZnONPs nanocomposites. A comparative analysis of MoS2 synthesized through chemical and physical exfoliation methods revealed that the chemically exfoliated MoS2 exhibited superior performance, thereby being selected for all subsequent measurements. The nanostructured materials demonstrated exceptional surface characteristics, with specific surface areas exceeding 300 m2 g−1. Notably, the hydrogen production rate achieved by a composite comprising 5% MoS2, 1.7% SiNWs, and 13.3% SWNTs at an 80% ZnONPs base was approximately 3909 µmol h−1g−1 under 500 nm wavelength radiation, marking a significant improvement of over 40-fold relative to pristine ZnONPs. This enhancement underscores the remarkable photocatalytic efficiency of the composites, maintaining high hydrogen production rates above 1500 µmol h−1g−1 even under radiation wavelengths exceeding 600 nm. Furthermore, the potential of these composites for energy storage and conversion applications, specifically within rechargeable lithium-ion batteries, was investigated. Composites, similar to those utilized for hydrogen production but excluding ZnONPs to address its limited theoretical capacity and electrical conductivity, were developed. The focus was on utilizing MoS2, SiNWs, and SWNTs as anode materials for Li-ion batteries. This strategic combination significantly improved the electronic conductivity and mechanical stability of the composite. Specifically, the composite with 56% MoS2, 24% SiNWs, and 20% SWNTs offered remarkable cyclic performance with high specific capacity values, achieving a complete stability of 1000 mA h g−1 after 100 cycles at 1 A g−1. These results illuminate the dual utility of the composites, not only as innovative catalysts for hydrogen production but also as advanced materials for energy storage technologies, showcasing their potential in contributing to sustainable energy solutions.

Chemistry↗

NASA NTRS1980

Nonrelativistic hyperfine splitting in muonic helium by adiabatic perturbation theory

Huang and Hughes have discussed the hyperfine splitting delta mu of muonic helium using a variational approach. In this paper, the Born-Oppenheimer approximation is used to simplify the evaluation of delta mu in the nonrelativistic limit. The first-order perturbed wave function of the electron is obtained in closed form by modifying the method used by Dalgarno and Lynn. The result is close to those of Huang and Hughes which required a large Hylleraas expansion and considerable extrapolation.

Drachman, R. J.↗

NASA NTRS1988

Reanalysis procedure for large structural systems

Global-basis-vector approximate reanalysis techniques for use in automated structural optimization schemes are developed and demonstrated. The vibrational response of a modified structure (MS) is estimated by lumping the design variables into a single tracing parameter, applying an operator splitting procedure to express the FEM equations of the MS in terms of the original-structure equations plus correction terms, and reducing the MS equations via a classical Bubnov-Galerkin scheme. The sensitivity of the vibrational response to structural modifications is evaluated, and the relationship between this method and the preconditioned conjugate-gradient method (Noor and Peters, 1988) is explored. Numerical results for linear static and free vibration problems involving beamlike lattices, double-layered hexahedral grids, and structural gridworks are presented in tables and graphs and briefly characterized.

Noor, Ahmed K.↗

NASA NTRS1992

Distributed traffic signal control using fuzzy logic

We present a distributed approach to traffic signal control, where the signal timing parameters at a given intersection are adjusted as functions of the local traffic condition and of the signal timing parameters at adjacent intersections. Thus, the signal timing parameters evolve dynamically using only local information to improve traffic flow. This distributed approach provides for a fault-tolerant, highly responsive traffic management system. The signal timing at an intersection is defined by three parameters: cycle time, phase split, and offset. We use fuzzy decision rules to adjust these three parameters based only on local information. The amount of change in the timing parameters during each cycle is limited to a small fraction of the current parameters to ensure smooth transition. We show the effectiveness of this method through simulation of the traffic flow in a network of controlled intersections.

Chiu, Stephen↗

NASA NTRS2012

Autonomous Information Unit for Fine-Grain Data Access Control and Information Protection in a Net-Centric System

As communication and networking technologies advance, networks will become highly complex and heterogeneous, interconnecting different network domains. There is a need to provide user authentication and data protection in order to further facilitate critical mission operations, especially in the tactical and mission-critical net-centric networking environment. The Autonomous Information Unit (AIU) technology was designed to provide the fine-grain data access and user control in a net-centric system-testing environment to meet these objectives. The AIU is a fundamental capability designed to enable fine-grain data access and user control in the cross-domain networking environments, where an AIU is composed of the mission data, metadata, and policy. An AIU provides a mechanism to establish trust among deployed AIUs based on recombining shared secrets, authentication and verify users with a username, X.509 certificate, enclave information, and classification level. AIU achieves data protection through (1) splitting data into multiple information pieces using the Shamir's secret sharing algorithm, (2) encrypting each individual information piece using military-grade AES-256 encryption, and (3) randomizing the position of the encrypted data based on the unbiased and memory efficient in-place Fisher-Yates shuffle method. Therefore, it becomes virtually impossible for attackers to compromise data since attackers need to obtain all distributed information as well as the encryption key and the random seeds to properly arrange the data. In addition, since policy can be associated with data in the AIU, different user access and data control strategies can be included. The AIU technology can greatly enhance information assurance and security management in the bandwidth-limited and ad hoc net-centric environments. In addition, AIU technology can be applicable to general complex network domains and applications where distributed user authentication and data protection are necessary. AIU achieves fine-grain data access and user control, reducing the security risk significantly, simplifying the complexity of various security operations, and providing the high information assurance across different network domains.

Chow, Edward T.↗

DOE OSTI2025

Experimental determination of the magnetic anisotropy in five-coordinated Co( II ) field-induced single molecule magnets

Magnetic anisotropy of the central metal atom is a crucial property of single molecule magnets (SMMs). Small structural changes can alter the magnetic properties, and accurate experimental methods to investigate magnetic anisotropy are therefore critical. Here, we investigate two five-coordinated Co( II ) SMMs, [CoCl 2 Cltpy] (1) and [CoBr 2 Cltpy] (2) (Cltpy = 4′-chloro-2,2′:6′,2′′-terpyridine), through multiple techniques. Ab initio theoretical calculations performed on the two compounds show that both possess axial magnetic anisotropy with the magnetic easy axis pointing towards one of the terminal halogen atoms. Theoretical calculations on SMMs are typically done on isolated molecular species, and to validate this approximation the magnetic anisotropy was further studied through experimental techniques. EPR measurements confirm an axial anisotropy of 1, and magnetic measurements provide experimental Zero-Field Splitting (ZFS) parameters, showing that the values from theoretical calculations are slightly overestimated. The X-ray electron density determined from 20 K single-crystal synchrotron X-ray diffraction data provides estimated d-orbital populations also suggesting axial magnetic anisotropy in both systems, and furthermore suggesting a more pronounced axiality in 1 compared to 2. This is in good agreement with the results obtained from both magnetic measurements and theoretical calculations. The magnetic anisotropy of 1 is quantified experimentally through polarized powder neutron diffraction via the site susceptibility method, confirming an axial magnetic anisotropy of the compound. A slight deviation in the easy axis direction is observed between experimental and theoretical results. This, together with the overestimation of the ZFS parameters from theoretical calculations, shows that experimental investigation of the magnetic anisotropy of SMMs is of high relevance. Magnetic anisotropy of the central metal atom is a crucial property of single molecule magnets (SMMs).

Slavensky, Hannah H. [Aarhus Univ. (Denmark)] (ORC↗

NASA NTRS1993

An upwind multigrid method for solving viscous flows on unstructured triangular meshes

A multigrid algorithm is combined with an upwind scheme for solving the two dimensional Reynolds averaged Navier-Stokes equations on triangular meshes resulting in an efficient, accurate code for solving complex flows around multiple bodies. The relaxation scheme uses a backward-Euler time difference and relaxes the resulting linear system using a red-black procedure. Roe's flux-splitting scheme is used to discretize convective and pressure terms, while a central difference is used for the diffusive terms. The multigrid scheme is demonstrated for several flows around single and multi-element airfoils, including inviscid, laminar, and turbulent flows. The results show an appreciable speed up of the scheme for inviscid and laminar flows, and dramatic increases in efficiency for turbulent cases, especially those on increasingly refined grids.

Bonhaus, Daryl Lawrence↗

DOE OSTI2023

Refining the GENEVA method for Higgs boson production via gluon fusion

We describe a number of improvements to the Geneva method for matching NNLO calculations to parton shower programs. In particular, we detail changes to the resummed calculation used in the matching procedure, including disentangling the cross section dependence on factorisation and beam scales, and an improved treatment of timelike logarithms. We also discuss modifications in the implementation of the splitting functions which serve to make the resummed calculation differential in the higher multiplicity phase space. These changes improve the stability of the numerical cancellation of the nonsingular term at small values of the resolution parameter. As a case study, we consider the gluon-initiated Higgs boson production process gg → H. We validate the NNLO accuracy of our predictions against independent calculations, and compare our showered and hadronised results with recent data taken at the ATLAS and CMS experiments in the diphoton decay channel, finding good agreement.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS1986

An isoparametric spectral element method for solution of the Navier-Stokes equations in complex geometry

High-order (p-type) finite element methods combine the geometric flexibility of standard low-order finite element (or finite volume) techniques with the rapid convergence properties of spectral methods. Various p-type schemes have been proposed for elliptic problems. However, an application of p-type methods to more complex equations, in particular the passive scalar or Navier-Stokes equations, is complicated by the nature of the equations and by efficiency considerations as regards the solution procedure. Patera (1984) has proposed a p-type method for the Navier-Stokes equations, taking into account the arising additional complications. The present study is concerned with an isoparametric spectral element discretization for the solution of flow problems in fully general, curvy geometries. Attention is given to the isoparametric formulation for elliptic equations and a time-splitting scheme for the Navier-Stokes equations.

Korczak, K. Z.↗

NASA NTRS1987

A new algorithm for the Navier-Stokes equations applied to transonic flows over wings

A new algorithm has been developed for the Navier-Stokes equations. For the convective terms, this algorithm employs flux vector splitting together with a locally rotated coordinate system that is aligned with the streamwise direction. The algorithm has been implemented into a Navier-Stokes code, the NASA Ames Research Center's TNS code, and several cases of steady flow have been calculated. The results show improvements over the original method in the code. First, in a case of separated flow, there is better agreement with the experimental results and also there is an absence of numerical oscillations that occur with the original method. Second, in a case with a strong supersonic to subsonic shock wave, comparisons show an improvement in the resolution of the shock wave and the reexpansion singularity.

Goorjian, Peter M.↗

NASA NTRS1989

A dual potential formulation of the Navier-Stokes equations

A dual potential formulation for numerically solving the Navier-Stokes equations is developed and presented. The velocity field is decomposed using a scalar and vector potential. Vorticity and dilatation are used as the dependent variables in the momentum equations. Test cases in two dimensions verify the capability to solve flows using approximations from potential flow to full Navier-Stokes simulations. A three-dimensional incompressible flow formulation is also described. An interesting feature of this approach to solving the Navier-Stokes equations is the decomposition of the velocity field into a rotational part (vector potential) and an irrotational part (scalar potential). The Helmholtz decomposition theorem allows this splitting of the velocity field. This approach has had only limited use since it increases the number of dependent variables in the solution. However, it has often been used for incompressible flows where the solution scheme is known to be fast and accurate. This research extends the usage of this method to fully compressible Navier-Stokes simulations by using the dilatation variable along with vorticity. A time-accurate, iterative algorithm is used for the uncoupled solution of the governing equations. Several levels of flow approximation are available within the framework of this method. Potential flow, Euler and full Navier-Stokes solutions are possible using the dual potential formulation. Solution efficiency can be enhanced in a straightforward way. For some flows, the vorticity and/or dilatation may be negligible in certain regions (e.g., far from a viscous boundary in an external flow). It is possible to drop the calculation of these variables then and optimize the solution speed. Also, efficient Poisson solvers are available for the potentials. The relative merits of non-primitive variables versus primitive variables for solution of the Navier-Stokes equations are also discussed.

Gegg, S. G.↗

NASA NTRSDate not supplied

Pilot Workload Rating Predictions Using Image Data and Recurrent Neural Networks

In this work, we augmented existing methods for estimating pilot workload ratings with deep neural networks trained using data from simulated flight tests in the Vertical Motion Simulator (VMS). We used an existing method, Spare Capacity Operations Estimator (SCOPE), along with a recurrent neural network and conducted comparison studies between the two methods individually, and when used together. We found that using both methods together can improve the result over using either approach alone. In our first test case, we achieved an improved linear correlation coefficient of 0.409 over that of SCOPE alone at 0.352 on the training dataset. Through cross validation, we also found that the results may be dependent on the split of training vs. validation data, and that further investigation should be conducted to understand what additional inputs to the neural network model should be made.

Image Data↗

DOE OSTI2020

Performance of position-sensitive resistive silicon detectors in the Forward Array Using Silicon Technology (FAUST)

The silicon telescopes of the Forward Array Using Silicon Technology (FAUST) have been recently upgraded with silicon detectors which use resistive charge-splitting to measure the position of charged particles. This is in addition to the measurement of the total energy that these charged particles deposit in the detector. The upgrade results in increased angular resolution with a much smaller number of signals than silicon strip detectors with similar resolution. Here, a novel method of calibration has been used to demonstrate the efficacy of these position-sensitive detectors, and the detectors have been used in a commissioning experiment.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗