SEARCH · Search NASA

Results for “implicit method”

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,027 records · Page 57

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Uniform accuracy of implicit-explicit Runge-Kutta (IMEX-RK) schemes for hyperbolic systems with relaxation

Implicit-explicit Runge-Kutta (IMEX-RK) schemes are popular methods to treat multiscale equations that contain a stiff part and a non-stiff part, where the stiff part is characterized by a small parameter. Here, in this work, we prove rigorously the uniform stability and uniform accuracy of a class of IMEX-RK schemes for a linear hyperbolic system with stiff relaxation. The result we obtain is optimal in the sense that it holds regardless of the value of and the order of accuracy is the same as the design order of the original scheme, i.e., there is no order reduction.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1986

Self-adaptive difference method for the effective solution of computationally complex problems of boundary layer theory

An implicit difference procedure for the solution of equations for a chemically reacting hypersonic boundary layer is described. Difference forms of arbitrary error order in the x and y coordinate plane were used to derive estimates for discretization error. Computational complexity and time were minimized by the use of this difference method and the iteration of the nonlinear boundary layer equations was regulated by discretization error. Velocity and temperature profiles are presented for Mach 20.14 and Mach 18.5; variables are velocity profiles, temperature profiles, mass flow factor, Stanton number, and friction drag coefficient; three figures include numeric data.

Schoenauer, W.↗

NASA NTRS1991

A Fourier collocation time domain method for numerically solving Maxwell's equations

A new method for solving Maxwell's equations in the time domain for arbitrary values of permittivity, conductivity, and permeability is presented. Spatial derivatives are found by a Fourier transform method and time integration is performed using a second order, semi-implicit procedure. Electric and magnetic fields are collocated on the same grid points, rather than on interleaved points, as in the Finite Difference Time Domain (FDTD) method. Numerical results are presented for the propagation of a 2-D Transverse Electromagnetic (TEM) mode out of a parallel plate waveguide and into a dielectric and conducting medium.

Shebalin, John V.↗

DOE OSTI2016

Transition State Theory Methods To Measure Diffusion in Flexible Nanoporous Materials: Application to a Porous Organic Cage Crystal

Transition state theory (TST) methods are useful for predicting adsorbate diffusivities in nanoporous materials at timescales inaccessible to molecular dynamics (MD). Most TST applications treat the nanoporous framework as rigid, which is inaccurate in highly flexible materials or where adsorbate dimensions are comparable to the size of pore aperatures. In this study, we demonstrate two computationally efficient TST methods for simulating adsorbate diffusion in nanoporous materials where framework flexibility has a significant influence on diffusion. These methods are applied to light gas diffusion in porous organic cage crystal 3 (CC3), a highly flexible molecular crystal that has shown promise in gas separation applications. Diffusion in CC3 is modeled as a series of uncorrelated adsorbate hops between cage molecules and the voids between adjacent cages. The first method we applied to compute adsorbate hopping rates in CC3 is implicit ligand sampling (ILS). In ILS TST, hopping rates are calculated in an ensemble of rigid framework snapshots captured from a fully flexible MD trajectory of the empty CC3 structure. The second TST method we applied is umbrella sampling (US), where hopping rates are computed from a series of biased MD simulations. Our ILS and US TST calculations are shown to agree well with direct MD simulation of adsorbate diffusion in CC3. We anticipate that the efficient TST methods detailed here will be broadly applicable other classes of flexible nanoporous materials such as metal–organic frameworks.

36 MATERIALS SCIENCE↗

NASA NTRS1983

Improving the convergence rate of parabolic ADI methods

The rate of convergence to steady state of parabolic Alternating Direction Implicit (ADI) solvers is analyzed in terms of the L(2)-norms of the residuals. The analysis allows one to predict the number of iterations necessary for convergence as function of the Courant number, Lambda. A simple modification of existing ADI codes is devised. It improves the convergence rate substantially and is insensitive to the Courant number in a large range of Lambda.

Abarbanel, S. S.↗

NASA NTRS1989

Analysis of a second-order-accurate finite-volume method for temporally-growing compressible shear layers

A finite-volume method for solving the compressible laminar Navier-Stokes equation in two dimensions is assessed for use in the simulation of transitional flows. The method is second-order-accurate, uses alternating-direction-implicit time integration, and a total-variation-diminishing smoothing operator in the inviscid flux terms. The test problem was a free shear layer with a forced periodicity in x at a specified wavelength and with the parallel mean flow specified as a hyperbolic-tangent profile. A grid refinement investigation verified this method to be second-order-accurate with less than 4 percent error in the growth rate of isolated modes on a 32 by 64 grid. The interaction between modes of similar amplitude was small, less than 1 percent. However, in cases where there was a dominant mode, numerical phase error pumped energy into all other modes. This error results in nonphysical growth rate for modes whose energy content is several orders of magnitude below the dominant mode. In the case when the dominant mode saturates, the pumping action stops and the growth rate of the next largest growing mode regains physical significance.

Atkins, H. L.↗

NASA NTRS1991

A Fourier-Chebyshev pseudospectral method for solving steady 3-D Navier-Stokes and heat equations in cylindrical cavities

A Fourier-Chebyshev pseudospectral method for solving steady 3D Navier-Stokes equations in cylindrical cavities is presented and discussed. The general method is pseudo-unsteady and uses a semi-implicit finite difference scheme for the time discretization. The generalized ADI (Alternating Direction Implicit) procedure is then applied to reduce the problem to a successive solution of one-dimensional problems. The spatial approximation uses a Fourier-Galerkin approximation in the periodic direction and a Chebyshev-collocation approximation in the other directions. Difficulties related to the pressure are surmounted by using the artifical compressibility method. A suitable variable change was chosen to avoid the problem of singularity at the axis generated by cylindrical coordinates. The method is first tested on an advection-diffusion equation and then on the Navier-Stokes equations. Finally, the method is illustrated by a convection problem of a differentially heated fluid.

Pulicani, J. P.↗

NASA NTRS1989

On finite element implementation and computational techniques for constitutive modeling of high temperature composites

The research work performed during the past year on finite element implementation and computational techniques pertaining to high temperature composites is outlined. In the present research, two main issues are addressed: efficient geometric modeling of composite structures and expedient numerical integration techniques dealing with constitutive rate equations. In the first issue, mixed finite elements for modeling laminated plates and shells were examined in terms of numerical accuracy, locking property and computational efficiency. Element applications include (currently available) linearly elastic analysis and future extension to material nonlinearity for damage predictions and large deformations. On the material level, various integration methods to integrate nonlinear constitutive rate equations for finite element implementation were studied. These include explicit, implicit and automatic subincrementing schemes. In all cases, examples are included to illustrate the numerical characteristics of various methods that were considered.

Saleeb, A. F.↗

DOE OSTI2019

Unified Approach to Implicit and Explicit Solvent Simulations of Electrochemical Reaction Energetics

One of the major open challenges in ab initio simulations of the electrochemical interface is the determination of electrochemical barriers under a constant driving force. Existing methods to do so include extrapolation techniques based on fully explicit treatments of the electrolyte, as well as implicit solvent models which allow for a continuous variation in electrolyte charge. Emerging hybrid continuum models have the potential to revolutionize the field, since they account for the electrolyte with little computational cost while retaining some explicit electrolyte, representing a “best of both worlds” method. In this work, we present a unified approach to determine reaction energetics from fully explicit, implicit, and hybrid treatments of the electrolyte based on a new multicapacitor model of the electrochemical interface. A given electrode potential can be achieved by a variety of interfacial structures; a crucial insight from this work is that the effective surface charge gives a good proxy of the local potential, the true driving force of electrochemical processes. In contrast, we show that the traditionally considered work function gives rise to multivalued functions depending on the simulation cell size. Furthermore, we show that the reaction energetics are largely insensitive to the countercharge distribution chosen in hybrid implicit/explicit models, which means that any of the myriad implicit electrolyte models can be equivalently applied. Furthermore, this work thus paves the way for the accurate treatment of ab initio reaction energetics of general surface electrochemical processes using both implicit and explicit electrolytes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1993

The effect of Reynolds number and turbulence on airfoil aerodynamics at -90 degrees incidence

A method has been developed for calculating the viscous flow about airfoils in with and without deflected flaps at -90 deg incidence. This method provides for the solution of the unsteady incompressible Navier-Stokes equations by means of an implicit technique. The solution is calculated on a body-fitted computational mesh using a staggered grid method. The vorticity is defined at the node points, and the velocity components are defined at the mesh-cell sides. The staggered-grid orientation provides for accurate representation of vorticity at the node points and the continuity equation at the mesh-cell centers. The method provides for the direct solution of the flow field and satisfies the continuity equation to machine zero at each time-step. The method is evaluated in terms of its ability to predict two-dimensional flow about an airfoil at -90 degrees incidence for varying Reynolds number and different boundary layer models. A laminar and a turbulent boundary layer models. A laminar and a turbulent boundary layer model are considered in the evaluation of the method. The variation of the average loading and surface pressure distribution due to flap deflection, Reynolds number, and laminar or turbulent flow are presented and compared with experimental results. The comparisons indicate that the calculated drag and drag reduction caused by flap deflection and the calculated average surface pressure are in excellent agreement with the measured results at a similar Reynolds number.

Stremel, Paul M.↗

NASA NTRS2013

Design of a Variational Multiscale Method for Turbulent Compressible Flows

A spectral-element framework is presented for the simulation of subsonic compressible high-Reynolds-number flows. The focus of the work is maximizing the efficiency of the computational schemes to enable unsteady simulations with a large number of spatial and temporal degrees of freedom. A collocation scheme is combined with optimized computational kernels to provide a residual evaluation with computational cost independent of order of accuracy up to 16th order. The optimized residual routines are used to develop a low-memory implicit scheme based on a matrix-free Newton-Krylov method. A preconditioner based on the finite-difference diagonalized ADI scheme is developed which maintains the low memory of the matrix-free implicit solver, while providing improved convergence properties. Emphasis on low memory usage throughout the solver development is leveraged to implement a coupled space-time DG solver which may offer further efficiency gains through adaptivity in both space and time.

Design↗

DOE OSTI2023

Physics-based adaptivity of a spectral method for the Vlasov–Poisson equations based on the asymmetrically-weighted Hermite expansion in velocity space

We propose a spectral method for the 1D-1V Vlasov–Poisson system where the discretization in velocity space is based on asymmetrically-weighted Hermite functions, dynamically adapted via a scaling α and shifting u of the velocity variable. Specifically, at each time instant an adaptivity criterion selects new values of α and u based on the numerical solution of the discrete Vlasov–Poisson system obtained at that time step. Once the new values of the Hermite parameters α and u are fixed, the Hermite expansion is updated and the discrete system is further evolved for the next time step. The procedure is applied iteratively over the desired temporal interval. The key aspects of the adaptive algorithm are: the map between approximation spaces associated with different values of the Hermite parameters that preserves total mass, momentum and energy; and the adaptivity criterion to update α and u based on physics considerations relating the Hermite parameters to the average velocity and temperature of each plasma species. For the discretization of the spatial coordinate, we rely on Fourier functions and use the implicit midpoint rule for time stepping. The resulting numerical method possesses intrinsically the property of fluid-kinetic coupling, where the low-order terms of the expansion are akin to the fluid moments of a macroscopic description of the plasma, while kinetic physics is retained by adding more spectral terms. Moreover, the scheme features conservation of total mass, momentum and energy associated in the discrete, for periodic boundary conditions. A set of numerical experiments confirms that the adaptive method outperforms the non-adaptive one in terms of accuracy and stability of the numerical solution.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1989

A new uniformly valid asymptotic integration algorithm for elasto-plastic-creep and unified viscoplastic theories including continuum damage

A new scheme to integrate a system of stiff differential equations for both the elasto-plastic creep and the unified viscoplastic theories is presented. The method has high stability, allows large time increments, and is implicit and iterative. It is suitable for use with continuum damage theories. The scheme was incorporated into MARC, a commercial finite element code through a user subroutine called HYPELA. Results from numerical problems under complex loading histories are presented for both small and large scale analysis. To demonstrate the scheme's accuracy and efficiency, comparisons to a self-adaptive forward Euler method are made.

Chulya, A.↗

NASA NTRS1991

A new uniformly valid asymptotic integration algorithm for elasto-plastic creep and unified viscoplastic theories including continuum damage

Chulya, Abhisak↗

DOE OSTI2018

An automated framework for NMR chemical shift calculations of small organic molecules

When using nuclear magnetic resonance (NMR) to assist in chemical identification in complex samples, researchers commonly rely on databases for chemical shift spectra. However, authentic standards are typically depended upon to build libraries experimentally. Considering complex biological samples, such as blood and soil, the entirety of NMR spectra required for all possible compounds would be infeasible to ascertain due to limitations of available standards and experimental processing time. As an alternative, we introduce the in silico Chemical Library Engine (ISiCLE) NMR chemical shift module to accurately and automatically calculate NMR chemical shifts of small organic molecules through use of quantum chemical calculations. ISiCLE performs density functional theory (DFT)-based calculations for predicting chemical properties - specifically NMR chemical shifts in this manuscript - via the open source, high-performance omputational chemistry software, NWChem. ISiCLE calculates the NMR chemical shifts of sets of molecules using any available combination of DFT method, solvent, and NMR-active nuclei, using both user-selected reference compounds and/or linear regression methods. Calculated NMR chemical shifts are provided to the user for each molecule, along with comparisons with respect to a number of metrics commonly used in the literature. Here, we demonstrate ISiCLE using a set of 312 molecules, ranging in size up to 90 carbon atoms. For each, calculation of NMR chemical shifts have been performed with 8 different levels of DFT theory, and with solvation effects using the implicit solvent Conductor-like Screening Model (COSMO). The DFT method dependence of the calculated chemical shifts have been systematically investigated through benchmarking and subsequently compared to experimental data available in the literature. Furthermore, ISiCLE has been applied to a set of 80 methylcyclohexane conformers, combined via Boltzmann weighting and compared to experimental values. We demonstrate that our protocol shows promise in the automation of chemical shift calculations and, ultimately, the expansion of chemical shift libraries.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1994

Effect of Reynolds number and turbulence on airfoil aerodynamics at -90-degree incidence

A method has been developed for calculating the viscous flow about airfoils with and without deflected flaps at -90 deg incidence. This method provides for the solution of the unsteady incompressible Navier-Stokes equations by means of an implicit technique. The solution is calculated on a body-fitted computational mesh using a staggered-grid method. The vorticity is defined at the node points, and the velocity components are defined at the mesh-cell sides. The staggered-grid orientation provides for accurate representation of vorticity at the node points and the continuity equation at the mesh-cell centers. The method provides for the noniterative solution of the flowfield and satisfies the continuity equation to machine zero at each time step. The method is evaluated in terms of its stability to predict two-dimensional flow about an airfoil at -90-deg incidence for varying Reynolds number and laminar/turbulent models. The variations of the average loading and surface pressure distribution due to flap deflection, Reynolds number, and laminar or turbulent flow are presented and compared with experimental results. The comparisom indicate that the calculated drag and drag reduction caused by flap deflection and the calculated average surface pressure are in excellent agreement with the measured results at a similar Reynolds number.

Stremel, Paul M.↗

NASA NTRS1990

Improved boundary layer heat transfer calculations near a stagnation point

A thermal design of a solar receiver has been developed for the solutions of problems involving phase-change thermal energy storage and natural convection loss. Two dimensional axisymmetrical solidification and melting of materials contained between two concentric cylinders of finite length has been studied for thermal energy storage analysis. For calculation of free convection loss inside receiver cavity, two dimensional axisymmetrical, laminar, transient free convection including radiation effects has been studied using integral/finite difference method. Finite difference equations are derived for the above analysis subject to constant or variable material properties, initial conditions, and boundary conditions. The validity of the analyses has been substantiated by comparing results of the present general method with available analytic solutions or numerical results reported in the literature. Both explicit and implicit schemes are tested in phase change analysis with different number of nodes ranging from 4 to 18. The above numerical methods have been applied to the existing solar receiver analyzing computer code as additional subroutines. The results were computed for one of the proposed Brayton cycle solar receiver models running under the actual environmental conditions. Effect of thermal energy storage on the thermal behavior of the receiver has been estimated. Due to the thermal energy storage, about 65% reduction on working gas outlet temperature fluctuation has been obtained; however, maximum temperature of thermal energy storage containment has been increased about 18%. Also, effect of natural convection inside a receiver cavity on the receiver heat transfer has been analyzed. The finding indicated that thermal stratification occurs during the sun time resulting in higher receiver temperatures at the outlet section of the gas tube, and lower temperatures at the inlet section of the gas tube when compared with the results with no natural convection. Due to heat supply from the air during the shade time, minimum temperature has been increased, while maximum temperature has been reduced due to convection loss to air. Consequently, cyclic temperature fluctuation has been reduced 29% for working gas and 16% for thermal energy storage containment. On the other hand, despite the presence of the natural convection the time-averaged temperatures for receiver components were found to be similar for two cases with/without natural convection (maximum difference was 1.8%).

Ahn, Kyung Hwan↗

DOE OSTI2016

IMPLEMENTATION AND VALIDATION OF A FULLY IMPLICIT ACCUMULATOR MODEL IN RELAP-7

This paper presents the implementation and validation of an accumulator model in RELAP-7 under the framework of preconditioned Jacobian free Newton Krylov (JFNK) method, based on the similar model used in RELAP5. RELAP-7 is a new nuclear reactor system safety analysis code being developed at the Idaho National Laboratory (INL). RELAP-7 is a fully implicit system code. The JFNK and preconditioning methods used in RELAP-7 is briefly discussed. The slightly modified accumulator model is summarized for completeness. The implemented model was validated with LOFT L3-1 test and benchmarked with RELAP5 results. RELAP-7 and RELAP5 had almost identical results for the accumulator gas pressure and water level, although there were some minor difference in other parameters such as accumulator gas temperature and tank wall temperature. One advantage of the JFNK method is its easiness to maintain and modify models due to fully separation of numerical methods from physical models. It would be straightforward to extend the current RELAP-7 accumulator model to simulate the advanced accumulator design.

Accumulator↗