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.

SourceResearch connectionFrom yearThrough year

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

At least 1,045 records · Page 58

All sources NASA NTRS DOE OSTI

NASA NTRS1985

Solution of partial differential equations on vector and parallel computers

The present status of numerical methods for partial differential equations on vector and parallel computers was reviewed. The relevant aspects of these computers are discussed and a brief review of their development is included, with particular attention paid to those characteristics that influence algorithm selection. Both direct and iterative methods are given for elliptic equations as well as explicit and implicit methods for initial boundary value problems. The intent is to point out attractive methods as well as areas where this class of computer architecture cannot be fully utilized because of either hardware restrictions or the lack of adequate algorithms. Application areas utilizing these computers are briefly discussed.

Ortega, J. M.↗

NASA NTRS1985

Solution of partial differential equations on vector and parallel computers

Ortega, J. M.↗

NASA NTRS1998

Functional and structural mapping of human cerebral cortex: solutions are in the surfaces

The human cerebral cortex is notorious for the depth and irregularity of its convolutions and for its variability from one individual to the next. These complexities of cortical geography have been a chronic impediment to studies of functional specialization in the cortex. In this report, we discuss ways to compensate for the convolutions by using a combination of strategies whose common denominator involves explicit reconstructions of the cortical surface. Surface-based visualization involves reconstructing cortical surfaces and displaying them, along with associated experimental data, in various complementary formats (including three-dimensional native configurations, two-dimensional slices, extensively smoothed surfaces, ellipsoidal representations, and cortical flat maps). Generating these representations for the cortex of the Visible Man leads to a surface-based atlas that has important advantages over conventional stereotaxic atlases as a substrate for displaying and analyzing large amounts of experimental data. We illustrate this by showing the relationship between functionally specialized regions and topographically organized areas in human visual cortex. Surface-based warping allows data to be mapped from individual hemispheres to a surface-based atlas while respecting surface topology, improving registration of identifiable landmarks, and minimizing unwanted distortions. Surface-based warping also can aid in comparisons between species, which we illustrate by warping a macaque flat map to match the shape of a human flat map. Collectively, these approaches will allow more refined analyses of commonalities as well as individual differences in the functional organization of primate cerebral cortex.

Review↗

DOE OSTI2021

Spatially Quasi-Periodic Water Waves of Infinite Depth

Abstract We formulate the two-dimensional gravity-capillary water wave equations in a spatially quasi-periodic setting and present a numerical study of solutions of the initial value problem. We propose a Fourier pseudo-spectral discretization of the equations of motion in which one-dimensional quasi-periodic functions are represented by two-dimensional periodic functions on a torus. We adopt a conformal mapping formulation and employ a quasi-periodic version of the Hilbert transform to determine the normal velocity of the free surface. Two methods of time-stepping the initial value problem are proposed, an explicit Runge–Kutta (ERK) method and an exponential time-differencing (ETD) scheme. The ETD approach makes use of the small-scale decomposition to eliminate stiffness due to surface tension. We perform a convergence study to compare the accuracy and efficiency of the methods on a traveling wave test problem. We also present an example of a periodic wave profile containing vertical tangent lines that is set in motion with a quasi-periodic velocity potential. As time evolves, each wave peak evolves differently, and only some of them overturn. Beyond water waves, we argue that spatial quasi-periodicity is a natural setting to study the dynamics of linear and nonlinear waves, offering a third option to the usual modeling assumption that solutions either evolve on a periodic domain or decay at infinity.

Wilkening, Jon (ORCID:0000000327827596)↗

DOE OSTI2020

A compact mechanical power take-off for wave energy converters: Design, analysis, and test verification

Power Take-off (PTO) is one of the most critical components in the ocean wave conversion. This paper introduces an innovative compact mechanical PTO using a ball screw mechanism and a Mechanical Motion Rectifier (MMR) to improve the energy conversion performance. The proposed MMR PTO converts bidirectional oscillatory motion of the ocean waves into unidirectional rotation. The detailed design and working principle of the MMR PTO are explicitly explained. Dynamic analysis is presented for the three most important factors that influence the performance, including the equivalent mass, the equivalent damping, and the excitation frequency. The analytical relationship between these factors and the disengaged freewheeling motion caused by the MMR mechanism is derived, and the analytical solution is verified through dry lab testing. Experimental results in bench tests show that compared with the traditional linear PTO which use the ball screw to drive the generator directly, the MMR PTO is more efficient because of the unique freewheeling motion, achieving a maximum 81.2% energy transfer efficiency. After integrating the prototype PTO with a single-body point absorber type of wave energy converter, testing in a water tank yields total PTO efficiency up to 62.4%, showing a very promising potential for real applications.

16 TIDAL AND WAVE POWER↗

DOE OSTI2021

High-salinity brine desalination with amine-based temperature swing solvent extraction: A molecular dynamics study

Solvent extraction desalination (SED) is a promising technique for efficiently treating high salinity brines. However, only a few studies in the literature have investigated the influence of the molecular characteristics of the organic solvents on the emergent performance of the separation process. Thus, in this work we investigate the influence of the molecular structure of four different secondary amines on the desalination process of a concentrated aqueous solution of NaCl. Molecular dynamics simulations and free energy calculations are used to obtain the thermodynamic and structural properties of the pure solvent bulk phase, as well as a homogeneous solvent-water bulk mixture. Furthermore, solvent-water and solvent-brine interfacial systems are used to explicitly analyze the water partitioning behavior in the relevant temperature range for SED. The water solubilization in these secondary amines is found to be an exothermic process. Additionally, although the hydrogen-bonding interaction between the water and solvent molecules are beneficial, this interaction is not necessarily the dominant factor for ranking the water solubility performance of the different solvents. Finally, the excess volume of mixing of the bulk solvent-water mixture, as well as the electrostatic surface features of the solvents are also found to provide important clues for screening solvent behavior.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2019

Role of electrons in collision cascades in solids. I. Dissipative model

In this article, we present a detailed model for the nonadiabatic coupling between ions and electrons in energetic ion-solid interactions over a wide range of energies in concentrated solid-solution fcc alloys of the 3d transition metals Ni, Co, Fe, and Cr. The model is based on general statistical mechanical principles and results in a stochastic modification of the classical nuclei motion which is parameterized by the first-principles calculation of a dissipation function produced by explicit time-dependent electronic evolution. This model provides a full picture of an entire collision process, from the ballistic to the thermal phases of a cascade, giving a detailed description of the energy exchange between ions and electrons till their final thermalization, removing in this way some ad hoc assumptions used in the state-of-the-art atomistic two-temperature models. Here, this work is separated in two papers; in the present Part I, we report on the ab initio methodology used to translate stopping power and electron-phonon interaction into a parameterized dissipation function; Part II, to be published, addresses the nonadiabatic ion dynamics using our modified Langevin dynamics [Tamm et al. Phys. Rev. Lett. 120, 185501 (2018)] applying the dissipation functions developed here to specific collision cascade events.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

Datasets for Custom-trained Machine-learning Interatomic Potentials: Nitric Acid Aqueous Solution

This dataset was generated using an iterative active learning strategy with the ArcaNN software package (https://github.com/arcann-chem/arcann_training) to train machine-learning interatomic potentials (MLIPs) for aqueous nitric acid. Each active-learning cycle consisted of three stages: (1) training, (2) exploration, and (3) labeling. The initial training set comprised approximately 800 randomly selected configurations from a previous study by Lewis et al. (https://doi.org/10.1021/jp205510q), which investigated nitric acid solutions at 2, 3, 4, and 5 mol/L. For all configurations, single-point calculations of atomic forces and total energies were performed at the quantum density functional theory BLYP-D2 and PBE-D3 levels of theory using the CP2K Quickstep module. Valence electrons were treated explicitly, while core electrons on all atoms were represented by norm-conserving Goedecker–Teter–Hutter (GTH) pseudopotentials. Long-range dispersion interactions were accounted for using Grimme dispersion corrections. Wave functions were expanded in a mixed Gaussian-and-plane-wave scheme using TZV2P-MOLOPT basis sets for all elements and an 800 Ry auxiliary plane-wave cutoff for the electron density. Self-consistent field convergence was accelerated using orbital transformation and Direct Inversion in the Iterative Subspace, with a convergence threshold of 10^{-6}. All single-point calculations were carried out in periodic orthorhombic cells whose dimensions match those of the molecular configurations sampled from earlier trajectories. The CELL_REF keyword in CP2K was used to define a fixed reference cell, ensuring consistency in the reference data used for MLIP training, particularly when cell fluctuations are present in NpT simulations. The resulting high-fidelity energies and forces constitute the ground-truth labels used to train the MLIPs contained in this dataset.

Dinpajooh, Mohammadhasan [Pacific Northwest Nation↗

DOE OSTI2026

The Fate of Nucleated Black Holes in de Sitter Quantum Gravity

The Euclidean Nariai geometry has long been proposed as the instanton describing the nucleation of maximal-mass black holes in de Sitter space. We place this interpretation on firmer footing by showing that, once an observer is included, the gravitational path integral produces the imaginary phase required for a transition rate. As a warmup, we revisit the Hawking-Moss instanton and, as a byproduct, find that scalar fields can enhance black-hole nucleation, suggesting a quantum-gravity bound on scalar potentials with de Sitter solutions. We then study the subsequent semiclassical evolution of the nucleated black hole. We show that the previously claimed "anti-evaporation" channel is unphysical, arising from a quantum state with singular horizons. In a smooth state, the black hole instead undergoes standard thermal Hawking evaporation. We verify explicit agreement with the no-boundary state and argue that this evaporation is not subject to large quantum-gravity corrections. The nucleated black hole thus evaporates completely back to the maximally-entropic empty de Sitter vacuum, making the full process a Boltzmann fluctuation.

FOS: Physical sciences↗

NASA NTRS1980

Discrete-continuous variable structural synthesis using dual methods

Approximation concepts and dual methods are extended to solve structural synthesis problems involving a mix of discrete and continuous sizing type of design variables. Pure discrete and pure continuous variable problems can be handled as special cases. The basic mathematical programming statement of the structural synthesis problem is converted into a sequence of explicit approximate primal problems of separable form. These problems are solved by constructing continuous explicit dual functions, which are maximized subject to simple nonnegativity constraints on the dual variables. A newly devised gradient projection type of algorithm called DUAL 1, which includes special features for handling dual function gradient discontinuities that arise from the discrete primal variables, is used to find the solution of each dual problem. Computational implementation is accomplished by incorporating the DUAL 1 algorithm into the ACCESS 3 program as a new optimizer option. The power of the method set forth is demonstrated by presenting numerical results for several example problems, including a pure discrete variable treatment of a metallic swept wing and a mixed discrete-continuous variable solution for a thin delta wing with fiber composite skins.

Schmit, L. A.↗

NASA NTRS1981

Some observations on a new numerical method for solving Navier-Stokes equations

An explicit-implicit technique for solving Navier-Stokes equations is described which, is much less complex than other implicit methods. It is used to solve a complex, two-dimensional, steady-state, supersonic-flow problem. The computational efficiency of the method and the quality of the solution obtained from it at high Courant-Friedrich-Lewy (CFL) numbers are discussed. Modifications are discussed and certain observations are made about the method which may be helpful in using it successfully.

Kumar, A.↗

NASA NTRS1991

Transient modeling/analysis of hyperbolic heat conduction problems employing mixed implicit-explicit alpha method

This paper describes the evaluation of mixed implicit-explicit finite element formulations for hyperbolic heat conduction problems involving non-Fourier effects. In particular, mixed implicit-explicit formulations employing the alpha method proposed by Hughes et al. (1987, 1990) are described for the numerical simulation of hyperbolic heat conduction models, which involves time-dependent relaxation effects. Existing analytical approaches for modeling/analysis of such models involve complex mathematical formulations for obtaining closed-form solutions, while in certain numerical formulations the difficulties include severe oscillatory solution behavior (which often disguises the true response) in the vicinity of the thermal disturbances, which propagate with finite velocities. In view of these factors, the alpha method is evaluated to assess the control of the amount of numerical dissipation for predicting the transient propagating thermal disturbances. Numerical test models are presented, and pertinent conclusions are drawn for the mixed-time integration simulation of hyperbolic heat conduction models involving non-Fourier effects.

Tamma, Kumar K.↗

NASA NTRS1994

Fourth order difference methods for hyperbolic IBVP's

Fourth order difference approximations of initial-boundary value problems for hyperbolic partial differential equations are considered. We use the method of lines approach with both explicit and compact implicit difference operators in space. The explicit operator satisfies an energy estimate leading to strict stability. For the implicit operator we develop boundary conditions and give a complete proof of strong stability using the Laplace transform technique. We also present numerical experiments for the linear advection equation and Burgers' equation with discontinuities in the solution or in its derivative. The first equation is used for modeling contact discontinuities in fluid dynamics, the second one for modeling shocks and rarefaction waves. The time discretization is done with a third order Runge-Kutta TVD method. For solutions with discontinuities in the solution itself we add a filter based on second order viscosity. In case of the non-linear Burger's equation we use a flux splitting technique that results in an energy estimate for certain different approximations, in which case also an entropy condition is fulfilled. In particular we shall demonstrate that the unsplit conservative form produces a non-physical shock instead of the physically correct rarefaction wave. In the numerical experiments we compare our fourth order methods with a standard second order one and with a third order TVD-method. The results show that the fourth order methods are the only ones that give good results for all the considered test problems.

Gustafsson, Bertil↗

NASA NTRS1959

Modified Matrix Method for Calculating Steady-State Span Loading on Flexible Wings in Subsonic Flight

A method is presented for shortening the computations required to determine the steady-state span loading on flexible wings in subsonic flight. The method makes use of tables of downwash factors to find the necessary aerodynamic-influence coefficients for the application of lifting-line theory. Explicit matrix equations of equilibrium are converted into a matrix power series with a finite number of terms by utilizing certain characteristic properties of matrices. The number of terms in the series is determined by a trial-and-error process dependent upon the required accuracy of the solution. Spanwise distributions of angle of attack, airload, shear, bending moment, and pitching moment are readily obtained as functions of qm(sub R) where q denotes the dynamic pressure and mR denotes the lift-curve slope of a rigid wing. This method is intended primarily to make it practical to solve steady-state aeroelastic problems on the ordinary manually operated desk calculators, but the method is also readily adaptable to automatic computing equipment.

Gainer, Patrick A.↗

DOE OSTI2024

Double winding condensate CFT

We conjecture a two-dimensional conformal field theory built out of a linear dilaton and a compact β γ system deformed by winding condensates on each of the compact cycles. In particular, explicit expressions of the residues of the correlation functions are provided. As a worldsheet theory, it describes a stringy black hole in AdS 3 (supported by the Kalb-Ramond flux) at the Hawking-Page temperature. It is expected to be connected to the Horowitz-Polchinski-like solution near the Hagedorn temperature of thermal AdS 3 and Euclidean BTZ blackhole as we vary the temperature.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2019

Solving scattering problems in the half-line using methods developed for scattering in the full line

We reduce the solution of the scattering problem defined on the half-line [0,∞) by a real or complex potential v(x) and a general homogeneous boundary condition at x=0 to that of the extension of v(x) to the full line that vanishes for x<0. We find an explicit expression for the reflection amplitude of the former problem in terms of the reflection and transmission amplitudes of the latter. We obtain a set of conditions on these amplitudes under which the potential in the half-line develops bound states, spectral singularities, and time-reversed spectral singularities where the potential acts as a perfect absorber. We examine the application of these results in the study of the scattering properties of a δ-function potential and a finite barrier potential defined in [0,∞), discuss optical systems modeled by these potentials, and explore the configurations in which these systems act as a laser or perfect absorber. In particular, we derive an explicit formula for the laser threshold condition for a slab laser with a single mirror and establish the surprising fact that a nearly perfect mirror gives rise to a lower threshold gain than a perfect mirror. We also offer a nonlinear extension of our approach which allows for utilizing a recently developed nonlinear transfer matrix method in the full line to deal with finite-range nonlinear scattering interactions defined in the half-line.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1974

Numerical computation of viscous blunt body flows with a planar impinging shock

Two- and three-dimensional, viscous blunt body flows with planar impinging shocks are computed using an explicit, time-dependent, finite-difference method to solve the complete set of Navier-Stokes equations. The bow shock is treated as a discontinuity, while all interior shock layer detail such as shear layers, shock waves, jets and the wall boundary layer are automatically captured in the solution. Numerical results are presented for cases in which planar shock waves of different strengths and orientations are allowed to impinge on the flow field surroundings an infinite cylinder resulting in two- and three-dimensional shock interference patterns. The numerical results are compared with experiment.

Holst, T. L.↗

NASA NTRS1977

Inversion methods in temperature and aerosol remote sounding: Their commonality and differences, and some unexplored approaches

The two remote sensing problems of temperature profiling and aerosol characterization (complex refractive index, size distribution) are considered. These problems differ only in the explicit form of the source function which, for aerosols, includes contributions from both single and multiple scattering processes. When the observables are the spectral extinction or the single scattering of the source radiation, the associated problem is completely analogous to the linearized temperature inversion problem. Methods for obtaining the solution of the linear problem are classified following three main categories: (1) derivation of properties that all solutions satisfy, which must then be properties of the actual solution; (2) regularization of the ill-posed problem; and (3) data changes within their domain of uncertainty in order to avoid the basic instability. A number of unexplored methods are indicated.

Fymat, A. L.↗