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

Transonic flow calculations using a flux vector splitting method for the Euler equations

A study of the flux vector splitting method of Steger and Warming for the solution of the time dependent Euler equations in strong conservation law form for arbitrary two-dimensional geometries is presented. The procedure employed here differs from that of Buning and Steger in that it uses a different algorithm and employs implicit boundary conditions. Moreover, the method, as implemented here, does not contain any explicit smoothing or any adjustable parameters. Calculations were carried out for an NACA 0012 airfoil at various Mach numbers and angles of attack, and cylinders. Steady symmetric solutions were obtained for the full cylinder at a freestream Mach number of .5 without imposing a symmetry condition. In general, good agreement with other methods was obtained.

Seaford, C. M.↗
DOE OSTI2020

Two-loop radiative jet function for exclusive B-meson and Higgs decays

The rare radiative B -meson decay \( {B}^{-}\to \gamma {\mathrm{\ell}}^{-}\overline{v} \) and the radiative Higgs-boson decay h → γγ mediated by light-quark loops both receive large logarithmic corrections in QCD, which can be resummed using factorization theorems derived in soft-collinear effective theory. In these factorization theorems the same radiative jet function appears, which is a central object in the study of factorization beyond the leading order in scale ratios. We calculate this function at two-loop order both in momentum space and in a dual space, where its renormalization-group evolution takes on a simpler form. We also derive the two-loop anomalous dimension of the jet function and present the exact solution to its evolution equation at two-loop order. Another important outcome of our analysis is the explicit form of the two-loop anomalous dimension of the B -meson light-cone distribution amplitude in momentum space.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Embedded pairs for optimal explicit strong stability preserving Runge–Kutta methods

We construct a family of embedded pairs for optimal explicit strong stability preserving Runge–Kutta methods of order 2 ≤ p ≤ 4 to be used to obtain numerical solution of spatially discretized hyperbolic PDEs. In this construction, the goals include non-defective property, large stability region, and small error values as defined in Dekker and Verwer (1984) and Kennedy et al. (2000). The new family of embedded pairs offer the ability for strong stability preserving (SSP) methods to adapt by varying the step-size. Through several numerical experiments, we assess the overall effectiveness in terms of work versus precision while also taking into consideration accuracy and stability.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Respecting causality for training physics-informed neural networks

While the popularity of physics-informed neural networks (PINNs) is steadily rising, to this date PINNs have not been successful in simulating dynamical systems whose solution exhibits multi-scale, chaotic or turbulent behavior. In this work we attribute this shortcoming to the inability of existing PINNs formulations to respect the spatio-temporal causal structure that is inherent to the evolution of physical systems. We argue that this is a fundamental limitation and a key source of error that can ultimately steer PINN models to converge towards erroneous solutions. We address this pathology by proposing a simple re-formulation of PINNs loss functions that can explicitly account for physical causality during model training. We demonstrate that this simple modification alone is enough to introduce significant accuracy improvements, as well as a practical quantitative mechanism for assessing the convergence of a PINNs model. Here we provide state-of-the-art numerical results across a series of benchmarks for which existing PINNs formulations fail, including the chaotic Lorenz system, the Kuramoto–Sivashinsky equation in the chaotic regime, and the Navier–Stokes equations. To the best of our knowledge, this is the first time that PINNs have been successful in simulating such systems, introducing new opportunities for their applicability to problems of industrial complexity.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Addendum to SAND2023-09604 Xyce lumped-element transmission line model verification to support Empire-Cable cable SGEMP analyses

This report supplements the Verification of Empire-Cable SAND report by expanding on the use of Xyce to simulate the coupling to a transmission line cable model. While Empire-Cable solves its governing equations on a high-order, finite-element mesh with an an implicit-in-time formulation, Xyce must use a first order graph for the circuit and explicit-in-time approach to be compatible with non-linear electrical device models. Thus, given the different solution methodologies in Xyce as compared to Empire-Cable, the convergence rates are expected to be different but the overall quality of the solution should be the same. The original four canonical problems studied in the Empire-Cable verification report are replicated here running in Xyce using transmission line modeling parameters from the verification report. Overall, agreement between the codes is excellent with Xyce’s convergence rates being limited mostly to first order due to the circuit network approximation of a transmission line being a first order approximation.

42 ENGINEERING↗
NASA NTRS1974

Numerical calculations of two dimensional, unsteady transonic flows with circulation

The feasibility of obtaining two-dimensional, unsteady transonic aerodynamic data by numerically integrating the Euler equations is investigated. An explicit, third-order-accurate, noncentered, finite-difference scheme is used to compute unsteady flows about airfoils. Solutions for lifting and nonlifting airfoils are presented and compared with subsonic linear theory. The applicability and efficiency of the numerical indicial function method are outlined. Numerically computed subsonic and transonic oscillatory aerodynamic coefficients are presented and compared with those obtained from subsonic linear theory and transonic wind-tunnel data.

Beam, R. M.↗
NASA NTRS1989

The upper-branch stability of compressible boundary layer flows

The upper-branch linear and nonlinear stability of compressible boundary layer flows is studied using the approach of Smith and Bodonyi (1982) for a similar incompressible problem. Both pressure gradient boundary layers and Blasius flow are considered with and without heat transfer, and the neutral eigenrelations incorporating compressibility effects are obtained explicitly. The compressible nonlinear viscous critical layer equations are derived and solved numerically and the results indicate some solutions with positive phase shift across the critical layer. Various limiting cases are investigated including the case of much larger disturbance amplitudes and this indicates the structure for the strongly nonlinear critical layer of the Benney-Bergeon (1969) type. It is also shown how a match with the inviscid neutral inflexional modes arising from the generalized inflexion point criterion, is achieved.

Gajjar, J. S. B.↗
DOE OSTI2022

Specific Ion Effects in Lanthanide–Amphiphile Structures at the Air–Water Interface and Their Implications for Selective Separation

The use of surfactants to attract dissolved ions to water surfaces and interfaces is an essential step in both solvent-based and solvent-free separation processes. Here, we have studied the interactions of lanthanide ions in the aqueous subphase with monolayers of dihexadecyl phosphate at air–water interfaces. With heavier lanthanides (atomic number Z ≥ 65) in the subphase, the floating layer can be compressed to an area/molecule of about half the molecular cross section, indicating bilayer formation. X-ray fluorescence and reflectivity data support this conclusion. In the presence of lighter lanthanides (Z < 65), only monolayers are observed. Subphase-concentration-dependent studies using Er 3+ (heavier) and Nd 3+ (lighter) lanthanides show a stepwise progression, with ions attaching to the monolayer only when the solution concentration is >3 × 10 –7 M. Above ~10 –5 M, bilayers form but only in the presence of the heavier lanthanide. Grazing incidence X-ray diffraction shows evidence of lateral ion–ion correlations in the bilayer structure but not in monolayers. Explicit solvent all-atom molecular dynamics simulations confirm the elevated ion–ion correlation in the bilayer system. This bilayer structure isolates heavier lanthanides but not lighter lanthanides from an aqueous solution and is therefore a potential mechanism for the selective separation of heavier lanthanides.

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

Numerical simulation of flow through the Langley parametric scramjet engine

The numerical simulation of a three-dimensional turbulent, reacting flow through the entire Langley parametric scramjet engine has been obtained using a piecewise elliptic approach. The last section in the combustor has been analyzed using a parabolized Navier-Stokes code. The facility nozzle flow was analyzed as a first step. The outflow conditions from the nozzle were chosen as the inflow conditions of the scramjet inlet. The nozzle and the inlet simulation were accomplished by solving the three-dimensional Navier-Stokes equations with a perfect gas assumption. The inlet solution downstream of the scramjet throat was used to provide inflow conditions for the combustor region. The first two regions of the combustor were analyzed using the MacCormack's explicit scheme. However, the source terms in the species equations were solved implicitly. The finite rate chemistry was modeled using the two-step reaction model of Rogers and Chinitz. A complete reaction model was used in the PNS code to solve the last combustor region. The numerical solutions provide an insight of the flow details in a complete hydrogen-fueled scramjet engine module.

Srinivasan, Shivakumar↗
NASA NTRS1980

An idealized model of homogeneous solid propellant combustion

The idealization of solid propellant combustion which consists of a single condensed-phase reaction and a single gas-phase reaction is discussed. A limiting-case solution to the gas-phase conservation equations is developed and tested by comparison with transformed thermocouple data. Closed-form expressions are derived for the burning rate and temperature sensitivity. Both quantities are linked explicitly to the solid and gas activation energies. It is concluded that the combustion phenomenology associated with the above idealization is still not fully understood. It is suggested that a careful study of the various limiting-case solutions tempered by a small set of numerical benchmark calculations might further such understanding.

Miller, M. S.↗
DOE OSTI2021

Error estimates of finite difference methods for the Dirac equation in the massless and nonrelativistic regime

We present four frequently used finite difference methods and establish the error bounds for the discretization of the Dirac equation in the massless and nonrelativistic regime, involving a small dimensionless parameter 0 < ε &NestedLessLess; 1 inversely proportional to the speed of light. In the massless and nonrelativistic regime, the solution exhibits rapid motion in space and is highly oscillatory in time. Specifically, the wavelength of the propagating waves in time is at O(ε), while in space, it is at O(1) with the wave speed at O(ε -1 ). We adopt one leap-frog, two semi-implicit, and one conservative Crank-Nicolson finite difference methods to numerically discretize the Dirac equation in one dimension and establish rigorously the error estimates which depend explicitly on the time step τ, mesh size h, and the small parameter ε. The error bounds indicate that, to obtain the “correct” numerical solution in the massless and nonrelativistic regime, i.e., 0 < ε &NestedLessLess; 1, all these finite difference methods share the same ε-scalability as time step τ = O(ε 3/2 ) and mesh size h = O(ε 1/2 ). A large number of numerical results are reported to verify the error estimates.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Hydration Free Energies of Polypeptides from Popular Implicit Solvent Models versus All-Atom Simulation Results Based on Molecular Quasichemical Theory

Calculating the hydration free energy of a macromolecule in all-atom simulations has long remained a challenge, necessitating the use of models wherein the effect of the solvent is captured without explicit account of solvent degrees of freedom. This situation has changed with developments in the molecular quasi-chemical theory (QCT)-an approach that enables calculation of the hydration free energy of macromolecules within all-atom simulations at the same resolution as is possible for small molecular solutes. The theory also provides a rigorous and physically transparent framework to conceptualize and model interactions in molecular solutions and thus provides a convenient framework to investigate the assumptions in implicit solvent models. In this study, we compare the results using molecular QCT versus predictions from EEF1, ABSINTH, and GB/SA implicit solvent models for polyglycine and polyalanine solutes covering a range of chain lengths and conformations. The hydration free energies or the differences in hydration free energies between conformers obtained from the implicit solvent models do not agree with explicit solvent results, with the deviations being largest for the group additive EEF1 and ABSINTH models. GB/SA does better in capturing the qualitative trends seen in explicit solvent results. Finally, analysis founded on QCT reveals the critical importance of the cooperativity of hydration that is inherent in the hydrophilic and hydrophobic contributions to hydration-physics that is not well captured in additive models but somewhat better accounted for by means of a dielectric in the GB/SA approach.

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

Modelling of Trapped Radiation Near Jupiter

Energetic (62 to approx. 130 MeV) proton fluxes measured with the Galileo Probe inside Jupiter's main ring (radii 122,500 to 128,940 km) show a modest increase as Jupiter is approached. Solutions of a reduced (1-d) equatorial diffusion equation with-constant losses described by a lifetime tau match those data for tau greater than or equal to 10(exp 9) s for diffusion coefficient D(sub LL) = 10(exp -9) L(exp 4)/s, 10(exp -9) L(exp -3)/s, or 10(exp -10) L(exp 4)/s (this particle population may be undergoing nearly loss-free inward radial diffusion). Tau would increase as Jupiter is approached if its value were determined by pitch angle scattering due to EM wave-particle interactions. if the absolute amplitudes of the waves' magnetic fluctuations did not vary with radial distance. Exploration of numerical solutions of the diffusive transport equation for Jupiter's magnetospheric ions at die planetary ring and closer to Jupiter has been done. Explicit range-energy relationships for energy loss in ring matter (SiO2), modeled as a continuous disk, are incorporated. A spatial resolution of approx. 1000 km in the radial direction is used in the vicinity of the main ring; off-equatorial ion fluxes have been included (assuming the atmosphere is a perfect absorber). For protons the maximum effect of energy loss in the microscopic size range of ring particles is expected at energies approx. 0.1 Mev. Presumably there could be sufficiently large, unobserved bodies within the ring that would have to be modeled as discrete objects, as planetary satellites would be. Within limitations of the model and computational resources it's planned to characterize solutions for transport of magnetospheric ions past the planetary ring; the procedures that are used are fairly well known and may be applied generally to the simpler magnetospheric models.

Mihalov, John D.↗
DOE OSTI2023

Closing in on Hydrologic Predictive Accuracy: Combining the Strengths of High-Fidelity and Physics-Agnostic Models

Applications of process-based models (PBM) for predictions are confounded by multiple uncertainties and computational burdens, resulting in appreciable errors. A novel modeling framework combining a high-fidelity PBM with surrogate and machine learning (ML) models is developed to tackle these challenges and applied for streamflow prediction. A surrogate model permits high computational efficiency of a PBM solution at a minimum loss of its accuracy. A novel probabilistic ML model partitions the PBM-surrogate prediction errors into reducible and irreducible types, quantifying their distributions that arise due to both explicitly perceived uncertainties (such as parametric) or those that are entirely hidden to the modeler (not included or unexpected). Using this approach, we demonstrate a substantial improvement of streamflow predictive accuracy for a case study urbanized watershed. Such a framework provides an efficient solution combining the strengths of high-fidelity and physics-agnostic models for a wide range of prediction problems in geosciences.

58 GEOSCIENCES↗
NASA NTRS1982

Mixed time integration methods for transient thermal analysis of structures

The computational methods used to predict and optimize the thermal structural behavior of aerospace vehicle structures are reviewed. In general, two classes of algorithms, implicit and explicit, are used in transient thermal analysis of structures. Each of these two methods has its own merits. Due to the different time scales of the mechanical and thermal responses, the selection of a time integration method can be a different yet critical factor in the efficient solution of such problems. Therefore mixed time integration methods for transient thermal analysis of structures are being developed. The computer implementation aspects and numerical evaluation of these mixed time implicit-explicit algorithms in thermal analysis of structures are presented. A computationally useful method of estimating the critical time step for linear quadrilateral element is also given. Numerical tests confirm the stability criterion and accuracy characteristics of the methods. The superiority of these mixed time methods to the fully implicit method or the fully explicit method is also demonstrated.

Liu, W. K.↗
NASA NTRS1983

Mixed time integration methods for transient thermal analysis of structures

The computational methods used to predict and optimize the thermal-structural behavior of aerospace vehicle structures are reviewed. In general, two classes of algorithms, implicit and explicit, are used in transient thermal analysis of structures. Each of these two methods has its own merits. Due to the different time scales of the mechanical and thermal responses, the selection of a time integration method can be a difficult yet critical factor in the efficient solution of such problems. Therefore mixed time integration methods for transient thermal analysis of structures are being developed. The computer implementation aspects and numerical evaluation of these mixed time implicit-explicit algorithms in thermal analysis of structures are presented. A computationally-useful method of estimating the critical time step for linear quadrilateral element is also given. Numerical tests confirm the stability criterion and accuracy characteristics of the methods. The superiority of these mixed time methods to the fully implicit method or the fully explicit method is also demonstrated.

Liu, W. K.↗
NASA NTRS2008

High-Order Energy Stable WENO Schemes

A new third-order Energy Stable Weighted Essentially NonOscillatory (ESWENO) finite difference scheme for scalar and vector linear hyperbolic equations with piecewise continuous initial conditions is developed. The new scheme is proven to be stable in the energy norm for both continuous and discontinuous solutions. In contrast to the existing high-resolution shock-capturing schemes, no assumption that the reconstruction should be total variation bounded (TVB) is explicitly required to prove stability of the new scheme. A rigorous truncation error analysis is presented showing that the accuracy of the 3rd-order ESWENO scheme is drastically improved if the tuning parameters of the weight functions satisfy certain criteria. Numerical results show that the new ESWENO scheme is stable and significantly outperforms the conventional third-order WENO finite difference scheme of Jiang and Shu in terms of accuracy, while providing essentially nonoscillatory solutions near strong discontinuities.

Yamaleev, Nail K.↗
NASA NTRS1985

Numerical solution of the vertical structure equation in the normal mode method

In the present model of multilayered stability stratification, aimed at obtaining the analytic eigensolutions of the vertical structure equation, each layer is characterized by its own static stability value. By requiring continuity of pressure and vertical velocity across each interface level, and by imposing suitable upper and lower boundary conditions, matching eigensolutions are obtained in terms of the Bessel functions. Attention is given to an explicit example of a double-layered stratified atmosphere which demonstrates the mathematical manipulations involved; the resultant vertical structure functions are used to check the accuracy of the numerical solutions by the finite difference and finite element methods.

Sasaki, Y. K.↗