Search NASA⌕ Search

SEARCH · Search NASA

Results for “explicit solution”

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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Entropy stable modal discontinuous Galerkin schemes and wall boundary conditions for the compressible Navier-Stokes equations

Entropy stable schemes ensure that physically meaningful numerical solutions also satisfy a semi-discrete entropy inequality under appropriate boundary conditions. Here in this work, we describe a discretization of viscous terms in the compressible Navier-Stokes equations which enables a simple and explicit imposition of entropy stable no-slip and reflective (symmetry) wall boundary conditions for discontinuous Galerkin (DG) discretizations. Specifically, we derive methods for imposing adiabatic no-slip and reflective (symmetry) boundary conditions for modal entropy stable DG formulations which preserve a semi-discrete entropy inequality. Numerical results confirm the robustness and accuracy of the proposed approaches.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Structure and Thermodynamic Stability of Zeolitic Imidazolate Framework Surfaces

Extensive efforts over the last several decades have focused on the possibility of “crystal engineering” metal–organic framework (MOF) materials with bulk properties tailored toward specific applications. However, beyond their bulk structure, MOF interfaces and surfaces can also play an important role in governing the materials properties relevant to many applications, including both interfacial mass and/or charge transport. We presently examine the diversity of stable surface structures/terminations of zeolitic imidazolate frameworks (ZIFs), a subclass of MOFs, under a variety of conditions that are characteristic of either gas-phase, “post-synthetic” conditions (relevant to many MOF applications) or solution-phase conditions (typical of those used during MOF synthesis). We construct surface phase diagrams to predict the most stable ZIF surface terminations as a function of external parameters (including temperature, adsorbate pressures, and pH), making explicit comparison against prior experimental observations, when possible. Overall, we find that the resulting phase diagrams can be used to explain the results of prior experimental studies of ZIF terminations across a variety of conditions and provide important insights into the factors that govern the structure of ZIF interfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

Evaluation of Implicit‐Explicit Additive Runge‐Kutta Integrators for the HOMME‐NH Dynamical Core

Abstract The nonhydrostatic High‐Order Method Modeling Environment (HOMME‐NH) atmospheric dynamical core supports acoustic waves that propagate significantly faster than the advective wind speed, thus greatly limiting the time step size that can be used with standard explicit time integration methods. Resolving acoustic waves is unnecessary for accurate climate and weather prediction. This numerical stiffness is addressed herein by considering implicit‐explicit additive Runge‐Kutta (ARK IMEX) methods that can treat the acoustic waves in a stable manner without requiring implicit treatment of nonstiff modes. Various ARK IMEX methods are evaluated for their efficiency in producing accurate solutions, ability to take large time step sizes, and sensitivity to grid cell length ratio. Both the gravity wave test and baroclinic instability test from the 2012 Dynamical Core Model Intercomparison Project are used to recommend 5 of the 27 ARK IMEX methods tested for use in HOMME‐NH.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Verification of an energy-conserving semi-implicit electrostatic particle-in-cell scheme for modeling high-density plasma at scale

A verification study of a semi-implicit energy-conserving electrostatic particle-in-cell algorithm is presented. The algorithm relaxes the time-step and mesh-size constraints that require resolution of the plasma period and Debye length associated with traditional explicit momentum-conserving particle-in-cell algorithms. Physical implications and applicability of using the semi-implicit scheme for modeling high-density plasmas are discussed. Where possible, numerical results are compared against analytical solutions. The simulation results indicate that the algorithm is stable at time steps larger than twice the inverse plasma frequency and cell sizes larger than the Debye length. It is found that the algorithm gives adequate results, provided that the distribution function and the spatiotemporal scales dictating the physics of the problem are resolved. As such, the algorithm may provide a robust method for kinetic modeling of high-density plasmas at scale.

Cyclotron resonance↗
DOE OSTI2023

Anti-polyelectrolyte and polyelectrolyte effects on conformations of polyzwitterionic chains in dilute aqueous solutions

Abstract Polyzwitterions (PZs) are considered as model synthetic analogs of intrinsically disordered proteins. Based on this analogy, PZs in dilute aqueous solutions are expected to attain either globular (i.e. molten, compact) or random coil conformations. Addition of salt is expected to open these conformations. To the best of our knowledge, these hypotheses about conformations of PZs have never been verified. In this study, we test these hypotheses by studying effects of added salt [potassium bromide (KBr)] on gyration and hydrodynamic radii of poly(sulfobetaine methacrylate) in dilute aqueous solutions using dynamic light scattering and small-angle X-ray scattering, respectively. Effects of zwitteration are revealed by direct comparisons of the PZs with the polymers of the same backbone but containing (1) no explicit charges on side groups such as poly(2-dimethylaminoethyl methacrylate)s and (2) explicit cationic side groups with tertiary amino bromide pendants. Zeta-potential measurements, transmission electron microscopy, and ab initio molecular dynamics simulations reveal that the PZs acquire net positive charge in near salt-free conditions due to protonation but retain coiled conformations. Added KBr leads to nonmonotonic changes exhibiting an increase followed by a decrease in radius of gyration (and hydrodynamic radius), which are called antipolyelectrolyte and polyelectrolyte effects, respectively. Charge regulation and screening of charge–charge interactions are discussed in relation to the antipolyelectrolyte and polyelectrolyte effects, respectively, which highlight the importance of salt in affecting net charge and conformations of PZs.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Implementation of Perturbation Theory and Sensitivity Capabilities in Griffin

Griffin is a Multiphysics Object-Oriented Simulation Environment (MOOSE) based reactor Multiphysics analysis application, jointly developed by Argonne and Idaho National Laboratories under the DOE-NE NEAMS program. This fiscal year, capabilities for reactivity and sensitivity evaluation using perturbation methods were implemented and verified. The First Order Perturbation Method (FOPT) was employed to compute reactivity worth resulting from small perturbations in input parameters, while the Generalized Perturbation Theory (GPT) was used to evaluate sensitivities of a range of response types, including reaction rate ratio, k-eigenvalue, neutron generation time, and effective delayed neutron fraction. These perturbation methods enable users to quantify how response quantities change due to a perturbation in a input parameter without explicitly performing an additional transport simulation for each perturbed state. In particular, the GPT formulation accounts for indirect effects arising from flux changes by solving generalized inhomogeneous equations, for which a Neumann series-based iterative solution method was developed and implemented in Griffin. The implemented reactivity and sensitivity evaluation capabilities were verified using two test problems: an infinite homogeneous system and a two-dimensional hexagonal core. The results showed excellent agreement with reference solutions obtained by a direct method based on finite difference approximation as well as GPT-based results from the PERSENT code, confirming the accuracy of both reactivity and sensitivity evaluations. Additionally, preliminary uncertainty quantification (UQ) results were obtained by combining the sensitivity values computed using GPT and external covariance data, demonstrating that the implemented sensitivity results can be reliably used for uncertainty calculations. To further demonstrate the generality and practical strength of the implementation, the sensitivity evaluation capability was successfully applied to the Empire microreactor with a geometrically complex design that poses significant modeling challenges. The results confirm that Griffin enables sensitivity evaluations even for irregular and highly heterogeneous reactor configurations, thereby establishing a foundation for UQ applications in advanced reactor designs and analyses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1973

Stability of matrix equations arising in electromagnetics

Many problems in electromagnetics can be formulated as a Fredholm equation of the first kind. In these problems it is desired to determine some function f(x) from its integral transform g(y). Except for special cases for which an explicit inverse operator is known, the equation is usually discretized by the method of moments into a matrix equation Af = g, which is then solved on a computer. A method of calculating the sensitivity of the solution of Af = g due to errors in either the vector g or the matrix A is developed. The usefulness of this method is demonstrated with problems from antenna theory, scattering theory, and remote sensing.

Klein, C.↗
NASA NTRS1983

An explicit example of Hopf bifurcation in fluid mechanics

It is observed that a complete and explicit example of Hopf bifurcation appears not to be known in fluid mechanics. Such an example is presented for the rotating Benard problem with free boundary conditions on the upper and lower faces, and horizontally periodic solutions. Normal modes are found for the linearization, and the Veronis computation of the wave numbers is modified to take into account the imposed horizontal periodicity. An invariant subspace of the phase space is found in which the hypotheses of the Joseph-Sattinger theorem are verified, thus demonstrating the Hopf bifurcation. The criticality calculations are carried through to demonstrate rigorously, that the bifurcation is subcritical for certain cases, and to demonstrate numerically that it is subcritical for all the cases in the paper.

Kloeden, P.↗
NASA NTRS1997

Fully Implicit Ablation and Thermal Response Program for Spacecraft Heatshield Analysis

A fully implicit ablation and thermal response program has been developed for the simulation of one-dimensional transient transport of thermal energy in a multilayer stack of isotropic materials and structure which can ablate from a front surface and decompose in-depth. Equations and numerical procedures for solution are described. Solutions are compared with those of Aerotherm Charring Material Thermal Response and Ablation Program, and with the arcjet data. The code is numerically more stable, and solves much wider range of problems compared with the existing explicit code. Applications of the code for the analysis of aeroshell heatshields of Stardust, Mars 2001, and Mars Microprobe using the advanced Light Weight Ceramic Ablators developed at the NASA Ames Research Center are presented and discussed in detail.

Chen, Y. K.↗
NASA NTRS2011

A Boundary Condition Relaxation Algorithm for Strongly Coupled, Ablating Flows Including Shape Change

Implementations of a model for equilibrium, steady-state ablation boundary conditions are tested for the purpose of providing strong coupling with a hypersonic flow solver. The objective is to remove correction factors or film cooling approximations that are usually applied in coupled implementations of the flow solver and the ablation response. Three test cases are considered - the IRV-2, the Galileo probe, and a notional slender, blunted cone launched at 10 km/s from the Earth's surface. A successive substitution is employed and the order of succession is varied as a function of surface temperature to obtain converged solutions. The implementation is tested on a specified trajectory for the IRV-2 to compute shape change under the approximation of steady-state ablation. Issues associated with stability of the shape change algorithm caused by explicit time step limits are also discussed.

Gnoffo, Peter A.↗
NASA NTRS1995

The space-time solution element method: A new numerical approach for the Navier-Stokes equations

This paper is one of a series of papers describing the development of a new numerical method for the Navier-Stokes equations. Unlike conventional numerical methods, the current method concentrates on the discrete simulation of both the integral and differential forms of the Navier-Stokes equations. Conservation of mass, momentum, and energy in space-time is explicitly provided for through a rigorous enforcement of both the integral and differential forms of the governing conservation laws. Using local polynomial expansions to represent the discrete primitive variables on each cell, fluxes at cell interfaces are evaluated and balanced using exact functional expressions. No interpolation or flux limiters are required. Because of the generality of the current method, it applies equally to the steady and unsteady Navier-Stokes equations. In this paper, we generalize and extend the authors' 2-D, steady state implicit scheme. A general closure methodology is presented so that all terms up through a given order in the local expansions may be retained. The scheme is also extended to nonorthogonal Cartesian grids. Numerous flow fields are computed and results are compared with known solutions. The high accuracy of the scheme is demonstrated through its ability to accurately resolve developing boundary layers on coarse grids. Finally, we discuss applications of the current method to the unsteady Navier-Stokes equations.

Scott, James R.↗
DOE OSTI2022

A new multi-model absolute difference-based sensitivity (MMADS) analysis method to screen non-influential processes under process model and parametric uncertainty

Process-based models have been widely used for hydrologic modeling, and it is a common practice to use sensitivity analysis methods for excluding non-influential hydrologic processes from further investigation and/or model improvement. This study develops a new method called multi-model absolute difference-based sensitivity (MMADS) analysis method to screen non-influential system processes and parameters. MMADS is conceptually similar to the Morris method for addressing parametric uncertainty, but has a unique feature to address both process model uncertainty (i.e., a process may be represented by multiple process models) and process model parameter uncertainty (i.e., parameters associated with a process model are random). MMADS first evaluates absolute differences of a quantity of interest (i.e., a system model output) by varying process models and/or process model parameter values, and then calculates the mean and variance of the differences for investigating process influence. The mean measures overall influence of the process on the quantity of interest, and the variance estimates influence of nonlinear effects of the process and/or its interactions with other processes. MMADS is an extension of the Morris method from a parameter space to a joint parameter-model space for explicitly addressing both process model uncertainty and model parameter uncertainty. The performance of MMADS is evaluated by using two numerical experiments. One experiment is based on Sobol’s G*-function with ten product elements, and has analytical solutions of the MMADS mean and variance of absolute differences. The other experiment is for groundwater flow modeling which considers three processes (i.e., recharge, geology, and snowmelt) that interact with each other. Finally, results indicate that MMADS is computationally efficient and can identify non-influential processes of complex hydrological systems.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1979

Comptonization of X-rays by low-temperature electrons

A method is described for calculating the spectrum that results from the Compton scattering of a monochromatic source of X-rays by low-temperature electrons, both for initial-value relaxation problems and for steady-state spatial diffusion problems. The method gives an exact solution of the inital-value problem for evolution of the spectrum in an infinite homogeneous medium if Klein-Nishina corrections to the Thomson cross section are neglected. This, together with approximate solutions for problems in which Klein-Nishina corrections are significant and/or spatial diffusion occurs, shows spectral structure near the original photon wavelength that may be used to infer physical conditions in cosmic X-ray sources. Explicit results, shown for examples of time relaxation in an infinite medium and spatial diffusion through a uniform sphere, are compared with results obtained by Monte Carlo calculations and by solving the appropriate Fokker-Planck equation.

Illarionov, A.↗
NASA NTRS1983

Analysis of a hybrid, uni-directional laminate with damage

The fracture behavior of hybrid (buffer strip) laminated composites is predicted by a novel method in terms of material properties, geometry, and internal damage. Attention is given to the cases of broken fibers in a unidirectional half-plane, adjoined half-planes having different fiber and matrix properties, and the solution of two half-planes bounding a third, distinct region of finite width. The analysis is based on a materials modeling approach using the classical shear-lag assumption to represent the stress transfer between fibers. Explicit fiber and matrix properties of the three regions are retained, and changes in the laminate behavior as a function of the material properties, buffer strip width and initial crack length are discussed.

Goree, J. G.↗
NASA NTRS1987

Analysis of structures with rotating, flexible substructures

A new methodology has been developed for the dynamic analysis of flexible structures, parts of which may be experiencing discrete motion relative to other parts. This methodology provides the capability of representing the continuum deformations typically treated using finite element methods. In addition, it provides the capability of representing the discrete motion at joints traditionally available with multibody methods. After decomposing the structure into substructures and associating a frame of reference with each substructure, the equations of motion for each substructure can be written explicitly including contributions due to the frame of reference generalized coordinates. By expanding the set of constraints to include constraints that eliminate the redundancy introduced by the frame generalized coordinates, the equations of motion become amenable to solution. The first digital computer program using this methodology, the General Rotorcraft Aeromechanical Stability Program (GRASP), was introduced in 1986. Although GRASP is limited to applications involving steady-state rotation, extension to arbitrary motions (including spin-up) can be accomplished by the selective retention of nonlinear terms in this formulation.

Hopkins, A. Stewart↗
NASA NTRS1989

Stability and dynamic analysis of a slender column with curved longitudinal stiffeners

The results of a stability design study are presented for a slender column with curved longitudinal stiffeners for large space structure applications. Linear stability analyses are performed using a link-plate representation of the stiffeners to determine stiffener local buckling stresses. Results from a set of parametric analyses are used to determine an approximate explicit expression for stiffener local buckling in terms of its geometric parameters. This expression along with other equations governing column stability and mass are assembled into a determinate system describing minimum mass stiffened column design. An iterative solution is determined to solve this system and a computer program incorporating this routine is presented. Example design problems are presented which verify the solution accuracy and illustrate the implementation of the solution routine. Also, observations are made which lead to a greatly simplified first iteration design equation relating the percent increase in column mass to the percent increase in column buckling load. From this, generalizations are drawn as to the mass savings offered by the stiffened column concept. Finally, the percent increase in fundamental column vibration frequency due to the addition of deployable stiffeners is studied.

Lake, Mark S.↗
NASA NTRS1990

Multistep matrix integrators for real-time simulation

An explicit linear multistep matrix-integration technique is presented for vector systems of ODEs which employs the stability region placement approach to permit the time-step to be chosen independently of system eigenvalues. Closed-form solutions for the general p-step method and the case where the system matrix has zero eigenvalues are given. It is shown that system mode shapes are preserved over the integration process, and that the technique remains applicable to systems with eigenvalues at their origin without need for computing a matrix inversion.

De Abreu-Garcia, J. A.↗
NASA NTRS1994

Improving the efficiency of aerodynamic shape optimization

The computational efficiency of an aerodynamic shape optimization procedure that is based on discrete sensitivity analysis is increased through the implementation of two improvements. The first improvement involves replacing a grid-point-based approach for surface representation with a Bezier-Bernstein polynomial parameterization of the surface. Explicit analytical expressions for the grid sensitivity terms are developed for both approaches. The second improvement proposes the use of Newton's method in lieu of an alternating direction implicit methodology to calculate the highly converged flow solutions that are required to compute the sensitivity coefficients. The modified design procedure is demonstrated by optimizing the shape of an internal-external nozzle configuration. Practically identical optimization results are obtained that are independent of the method used to represent the surface. A substantial factor of 8 decrease in computational time for the optimization process is achieved by implementing both of the design procedure improvements.

Burgreen, Greg W.↗