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 1,369 records · Page 76
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Influence of non-adiabatic effects on linear absorption spectra in the condensed phase: Methylene blue

Modeling linear absorption spectra of solvated chromophores is highly challenging as contributions are present both from coupling of the electronic states to nuclear vibrations and from solute–solvent interactions. In systems where excited states intersect in the Condon region, significant non-adiabatic contributions to absorption line shapes can also be observed. Here, we introduce a robust approach to model linear absorption spectra accounting for both environmental and non-adiabatic effects from first principles. This model parameterizes a linear vibronic coupling (LVC) Hamiltonian directly from energy gap fluctuations calculated along molecular dynamics (MD) trajectories of the chromophore in solution, accounting for both anharmonicity in the potential and direct solute–solvent interactions. The resulting system dynamics described by the LVC Hamiltonian are solved exactly using the thermalized time-evolving density operator with orthogonal polynomials algorithm (T-TEDOPA). The approach is applied to the linear absorption spectrum of methylene blue in water. We show that the strong shoulder in the experimental spectrum is caused by vibrationally driven population transfer between the bright S1 and the dark S2 states. The treatment of the solvent environment is one of many factors that strongly influence the population transfer and line shape; accurate modeling can only be achieved through the use of explicit quantum mechanical solvation. The efficiency of T-TEDOPA, combined with LVC Hamiltonian parameterizations from MD, leads to an attractive method for describing a large variety of systems in complex environments from first principles.

Dunnett, Angus J. (ORCID:0000000285797826)↗
NASA NTRS1982

Algebraic grid generation

The numerical solution of partial differential equations about irregular geometries and with varying characteristic scales has created the need for coordinate systems and associated transformations which reflect both geometric and physical requirements. The process of finding coordinate transformations in discrete representations is called 'grid generation'. The present investigation is concerned with three algebraic grid generation methods. The methods include transfinite interpolation, the multisurface method, and the two-boundary technique. Interpolation formulas, in terms of homotopic mappings and constraints in terms of point positions and/or derivatives are the essential elements of the techniques. The methods are relatively simple to understand, they are explicit and do not require extensive computational effort, and they have a high degree of generality.

Smith, R. E.↗
NASA NTRS1986

Supercomputer implementation of finite element algorithms for high speed compressible flows

Prediction of compressible flow phenomena using the finite element method is of recent origin and considerable interest. Two shock capturing finite element formulations for high speed compressible flows are described. A Taylor-Galerkin formulation uses a Taylor series expansion in time coupled with a Galerkin weighted residual statement. The Taylor-Galerkin algorithms use explicit artificial dissipation, and the performance of three dissipation models are compared. A Petrov-Galerkin algorithm has as its basis the concepts of streamline upwinding. Vectorization strategies are developed to implement the finite element formulations on the NASA Langley VPS-32. The vectorization scheme results in finite element programs that use vectors of length of the order of the number of nodes or elements. The use of the vectorization procedure speeds up processing rates by over two orders of magnitude. The Taylor-Galerkin and Petrov-Galerkin algorithms are evaluated for 2D inviscid flows on criteria such as solution accuracy, shock resolution, computational speed and storage requirements. The convergence rates for both algorithms are enhanced by local time-stepping schemes. Extension of the vectorization procedure for predicting 2D viscous and 3D inviscid flows are demonstrated. Conclusions are drawn regarding the applicability of the finite element procedures for realistic problems that require hundreds of thousands of nodes.

Thornton, E. A.↗
NASA NTRS1992

Numerical simulation of turbomachinery flows with advanced turbulence models

A three dimensional full Navier-Stokes (FNS) code is used to simulate complex turbomachinery flows. The code incorporates an explicit multistep scheme and solves a conservative form of the density averaged continuity, momentum, and energy equations. A compressible low Reynolds number form of the k-epsilon turbulence model, and a q-omega model and an algebraic Reynolds stress model have been incorporated in a fully coupled manner to approximate Reynolds stresses. The code is used to predict the viscous flow field in a backswept transonic centrifugal compressor for which laser two focus data is available. The code is also used to simulate the tip clearance flow in a cascade. The code has been extended to include unsteady Euler solutions for predicting the unsteady flow through a cascade due to incoming wakes, simulating rotor-stator interactions.

Lakshminarayana, B.↗
DOE OSTI2021

Schrödinger evolution of two-dimensional black holes

This paper systematically treats the evolving quantum state for two-dimensional black holes, with particular focus on the CGHS model, but also elucidating features generalizing to higher dimensions. This is done in Schrödinger picture(s), to exhibit the dynamic evolution of the state at intermediate times. After a review of classical solutions, also connecting to descriptions of higher-dimensional black holes, it overviews the canonical quantum treatment of the full evolution, including gravitational dynamics. Derived in an approximation to this, following conversion to “perturbation picture”, is the evolution of the quantum matter on the background geometry. Features of the evolving matter state are described, based on choice of a time slicing to put the evolution into ADM form. The choices of slicing as well as coordinates on the slices result in different quantum “pictures” for treating the evolution. If such a description is based on smooth trans-horizon slices, that avoids explicit reference to ultra-planckian modes familiar from traditional treatments, and exhibits the Hawking excitations as emerging from a “quantum atmosphere” with thickness comparable to the inverse temperature. Detailed study of the state exhibits the entanglement structure between Hawking quanta and the partner excitations inside the black hole, and the corresponding “missing information”. This explicit description also allows direct study of the evolution and features, e.g. as seen by infalling observers, of these partner excitations, helping to address various puzzles with them. Explicit treatment of the evolving state, and its extension to higher dimensions, provides further connections to information theory and a starting point for study of corrections that can unitarize evolution, arising from new quantum gravity effects — whether wormholes or something entirely different.

2D Gravity↗
NASA NTRS1981

Numerical simulation of reentry flow around the Shuttle Orbiter including real gas effects

The flight environment of a reentry vehicle is predicted from the numerical solution of fluid dynamics equations for the region between the blunt nose and the swept wings. The inviscid portion of the shock layer is modeled by the Euler equations, but the laminar viscous flow adjacent to the wall is modeled by the approximate parabolic Navier-Stokes equations. The approximations made to the axial gradients of pressure and diffusive fluxes enable the coupled inviscid and viscous equations to be solved efficiently along the body axis. The equilibrium air aftbody code contains significant improvements over its predecessor which only considers a perfect gas model and noncircular configurations. The inclusion of numerical damping either explicitly or implicitly has extended its capabilities for predicting flow field around a winged configuration at higher Machs and greater angles of attack. The results are obtained on the cylindrical coordinates and satisfactory for the Shuttle Orbiter at a free-stream Mach number of 22 and an angle of attack of 40 deg. Also discussed are the inviscid formulation and its application for wind-tunnel conditions.

Li, C. P.↗
DOE OSTI2020

Faster Johnson–Lindenstrauss transforms via Kronecker products

The Kronecker product is an important matrix operation with a wide range of applications in signal processing, graph theory, quantum computing and deep learning. In this work, we introduce a generalization of the fast Johnson–Lindenstrauss projection for embedding vectors with Kronecker product structure, the Kronecker fast Johnson–Lindenstrauss transform (KFJLT). The KFJLT reduces the embedding cost by an exponential factor of the standard fast Johnson–Lindenstrauss transform’s cost when applied to vectors with Kronecker structure, by avoiding explicitly forming the full Kronecker products. Here, we prove that this computational gain comes with only a small price in embedding power: consider a finite set of $p$ points in a tensor product of $d$ constituent Euclidean spaces $\bigotimes _{k=d}^{1}{\mathbb{R}}^{n_k}$, and let $N = \prod _{k=1}^{d}n_k$. With high probability, a random KFJLT matrix of dimension $m \times N$ embeds the set of points up to multiplicative distortion $(1\pm \varepsilon )$ provided $m \gtrsim \varepsilon ^{-2} \, \log ^{2d - 1} (p) \, \log N$. We conclude by describing a direct application of the KFJLT to the efficient solution of large-scale Kronecker-structured least squares problems for fitting the CP tensor decomposition.

Kronecker structure↗
NASA NTRS2002

High-Order Semi-Discrete Central-Upwind Schemes for Multi-Dimensional Hamilton-Jacobi Equations

We present high-order semi-discrete central-upwind numerical schemes for approximating solutions of multi-dimensional Hamilton-Jacobi (HJ) equations. This scheme is based on the use of fifth-order central interpolants like those developed in [1], in fluxes presented in [3]. These interpolants use the weighted essentially nonoscillatory (WENO) approach to avoid spurious oscillations near singularities, and become "central-upwind" in the semi-discrete limit. This scheme provides numerical approximations whose error is as much as an order of magnitude smaller than those in previous WENO-based fifth-order methods [2, 1]. Thee results are discussed via examples in one, two and three dimensions. We also pregnant explicit N-dimensional formulas for the fluxes, discuss their monotonicity and tl!e connection between this method and that in [2].

Bryson, Steve↗
DOE OSTI2022

Revisiting Multi-Material Composite Structures with Homogenized Composite Properties

Composite structures inherently develop residual stresses during their curing process. Driven predominately by mismatched thermal strains between differing materials or ply orientations, but also affected by curing process phenomena like polymer shrinkage, these residual stresses can lead to failure within composite structures. There are several methods varying in complexity that can be used to model the development of residual stresses, all of which are capable of capturing sufficient detail to understand the residual stress state at the ply level. However, explicitly modeling all plies of a layup in a composite structure can be prohibitively expensive based on the number of plies, structure size, and required element size. The computational cost can be reduced through the homogenization of the composite layup without losing much fidelity of the overall response of the structure. The homogenization process reduces the many plies of a laminate to a single lamina that reduces complexity and increases the mesh size where a single element can span multiple plies. This report focuses on verification and validation efforts for a homogenization process using a suite of finite element simulations rather than an analytic solution derived from classical laminate theory. Initial verification using representative element volumes indicated there was minimal error in the homogenization process; however, this compounded to a small, but acceptable error in strip and split ring experimental composite structures. The error does under predict the residual stress state in the strip and split ring and should be accounted for when simulating composite structures with homogenized properties.

36 MATERIALS SCIENCE↗
NASA NTRS1996

Further Development of a New, Flux-Conserving Newton Scheme for the Navier-Stokes Equations

This paper is one of a series of papers describing the development of a new numerical approach for solving the steady Navier-Stokes equations. The key features in the current development are (1) the discrete representation of the dependent variables by way of high order polynomial expansions, (2) the retention of all derivatives in the expansions as unknowns to be explicitly solved for, (3) the automatic balancing of fluxes at cell interfaces, and (4) the discrete simulation of both the integral and differential forms of the governing equations. The main purpose of this paper is, first, to provide a systematic and rigorous derivation of the conditions that are used to simulate the differential form of the Navier-Stokes equations, and second, to extend our previously-presented internal flow scheme to external flows and nonuniform grids. Numerical results are presented for high Reynolds number flow (Re = 100,000) around a finite flat plate, and detailed comparisons are made with the Blasius flat plate solution and Goldstein wake solution. It is shown that the error in the streamwise velocity decreases like r(sup alpha)(Delta)y(exp 2), where alpha approx. 0.25 and r = delta(y)/delta(x) is the grid aspect ratio.

Scott, James R.↗
DOE OSTI2021

CHARMM-GUI Drude prepper for molecular dynamics simulation using the classical Drude polarizable force field

Explicit treatment of electronic polarizability in empirical force fields (FFs) represents an extension over a traditional additive or pairwise FF and provides a more realistic model of the variations in electronic structure in condensed phase, macromolecular simulations. To facilitate utilization of the polarizable FF based on the classical Drude oscillator model, Drude Prepper has been developed in CHARMM-GUI. Drude Prepper ingests additive CHARMM protein structures file (PSF) and pre-equilibrated coordinates in CHARMM, PDB, or NAMD format, from which the molecular components of the system are identified. These include all residues and patches connecting those residues along with water, ions, and other solute molecules. This information is then used to construct the Drude FF-based PSF using molecular generation capabilities in CHARMM, followed by minimization and equilibration. In addition, inputs are generated for molecular dynamics (MD) simulations using CHARMM, GROMACS, NAMD, and OpenMM. Validation of the Drude Prepper protocol and inputs is performed through conversion and MD simulations of various heterogeneous systems that include proteins, nucleic acids, lipids, polysaccharides, and atomic ions using the aforementioned simulation packages. Stable simulations are obtained in all studied systems, including 5 μs simulation of ubiquitin, verifying the integrity of the generated Drude PSFs. Additionally, the ability of the Drude FF to model variations in electronic structure is shown through dipole moment analysis in selected systems. Finally, the capabilities and availability of Drude Prepper in CHARMM-GUI is anticipated to greatly facilitate the application of the Drude FF to a range of condensed phase, macromolecular systems.

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

Analytic evaluation of two-center molecular integrals

By using the Fourier-transform technique, the explicit expressions for the one-electron - two-center overlap integrals of Slater-type atomic orbitals up to 3d are derived. The final expressions are analytic, simple, and independent of local coordinates. Furthermore, they do not contain the nonclosed-form of exponential integrals which were presented in expressions given in earlier work. It is shown that the two-electron - two-center Coulomb integrals, as well as the hybrid integrals, can simply be expressed in terms of these integrals. The numerical instability arising from the situation in which the exponents of the two orbitals are almost equal is discussed, and a solution for this problem based on a Taylor-series expansion of the integral is suggested.

Tai, H.↗
DOE OSTI2023

The gallium anomaly

In order to test the end-to-end operations of gallium solar neutrino experiments, intense electron-capture sources were fabricated to measure the responses of the radiochemical SAGE and GALLEX/GNO detectors to known fluxes of low-energy neutrinos. Such tests were viewed at the time as a cross-check, given the many tests of 71 Ge recovery and counting that had been routinely performed, with excellent results. However, the four 51 Cr and 37 Ar source experiments yielded rates below expectations, a result commonly known as the Ga anomaly. As the intensity of the electron-capture sources can be measured to high precision, the neutrino lines they produce are fixed by known atomic and nuclear rates, and the neutrino absorption cross section on 71 Ga is tightly constrained by the lifetime of 71 Ge, no simple explanation for the anomaly has been found. To check these calibration experiments, a dedicated experiment BEST was performed, utilizing a neutrino source of unprecedented intensity and a detector optimized to increase statistics while providing some information on counting rate as a function of distance from the source. The results BEST obtained are consistent with the earlier solar neutrino calibration experiments, and when combined with those measurements, yield a Ga anomaly with a significance of approximately 4𝜎, under conservative assumptions. But BEST found no evidence of distance dependence and thus no explicit indication of new physics. In this review we describe the extensive campaigns carried out by SAGE, GALLEX/GNO, and BEST to demonstrate the reliability and precision of their experimental procedures, including 71 Ge recovery, counting, and analysis. We also describe efforts to define uncertainties in the neutrino capture cross section, which now include estimates of effects at the ≲ 0.5% level such as radiative corrections and weak magnetism. With the results from BEST, an anomaly remains even if one retains only the transition to the 71 Ge ground state, whose strength is fixed by the known lifetime of 71 Ge. We then consider the new-physics solution most commonly suggested to resolve the Ga anomaly, oscillations into a sterile fourth neutrino, 𝜈 𝑒 → 𝜈 𝑠 . We find such a solution generates substantial tension with several null experiments, owing to the large mixing angle required. While this does not exclude such solutions – the sterile sector might include multiple neutrinos as well as new interactions – it shows the need for more experimental constraints, if we are to make progress in resolving the Ga and other low-energy neutrino anomalies. We conclude by consider the role future low-energy electron-capture sources could play in this effort.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

The IsoGenie database: an interdisciplinary data management solution for ecosystems biology and environmental research

Modern microbial and ecosystem sciences require diverse interdisciplinary teams that are often challenged in “speaking” to one another due to different languages and data product types. Here we introduce the IsoGenie Database, a de novo developed data management and exploration platform, as a solution to this challenge of accurately representing and integrating heterogenous environmental and microbial data across ecosystem scales. The IsoGenieDB is a public and private data infrastructure designed to store and query data generated by the IsoGenie Project, a ~10 year DOE-funded project focused on discovering ecosystem climate feedbacks in a thawing permafrost landscape. The IsoGenieDB provides (i) a platform for IsoGenie Project members to explore the project’s interdisciplinary datasets across scales through the inherent relationships among data entities, (ii) a framework to consolidate and harmonize the datasets needed by the team’s modelers, and (iii) a public venue that leverages the same spatially explicit, disciplinarily integrated data structure to share published datasets. The IsoGenieDB is also being expanded to cover the NASA-funded Archaea to Atmosphere (A2A) project, which scales the findings of IsoGenie to a broader suite of Arctic peatlands, via the umbrella A2A Database (A2A-DB). The IsoGenieDB’s expandability and flexible architecture allow it to serve as an example ecosystems database.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1994

Tensor-GMRES method for large sparse systems of nonlinear equations

This paper introduces a tensor-Krylov method, the tensor-GMRES method, for large sparse systems of nonlinear equations. This method is a coupling of tensor model formation and solution techniques for nonlinear equations with Krylov subspace projection techniques for unsymmetric systems of linear equations. Traditional tensor methods for nonlinear equations are based on a quadratic model of the nonlinear function, a standard linear model augmented by a simple second order term. These methods are shown to be significantly more efficient than standard methods both on nonsingular problems and on problems where the Jacobian matrix at the solution is singular. A major disadvantage of the traditional tensor methods is that the solution of the tensor model requires the factorization of the Jacobian matrix, which may not be suitable for problems where the Jacobian matrix is large and has a 'bad' sparsity structure for an efficient factorization. We overcome this difficulty by forming and solving the tensor model using an extension of a Newton-GMRES scheme. Like traditional tensor methods, we show that the new tensor method has significant computational advantages over the analogous Newton counterpart. Consistent with Krylov subspace based methods, the new tensor method does not depend on the factorization of the Jacobian matrix. As a matter of fact, the Jacobian matrix is never needed explicitly.

Feng, Dan↗
NASA NTRS1998

User Manual for Beta Version of TURBO-GRD: A Software System for Interactive Two-Dimensional Boundary/ Field Grid Generation, Modification, and Refinement

TURBO-GRD is a software system for interactive two-dimensional boundary/field grid generation. modification, and refinement. Its features allow users to explicitly control grid quality locally and globally. The grid control can be achieved interactively by using control points that the user picks and moves on the workstation monitor or by direct stretching and refining. The techniques used in the code are the control point form of algebraic grid generation, a damped cubic spline for edge meshing and parametric mapping between physical and computational domains. It also performs elliptic grid smoothing and free-form boundary control for boundary geometry manipulation. Internal block boundaries are constructed and shaped by using Bezier curve. Because TURBO-GRD is a highly interactive code, users can read in an initial solution, display its solution contour in the background of the grid and control net, and exercise grid modification using the solution contour as a guide. This process can be called an interactive solution-adaptive grid generation.

Choo, Yung K.↗
NASA NTRS2002

Development of Computational Aeroacoustics Code for Jet Noise and Flow Prediction

Accurate prediction of jet fan and exhaust plume flow and noise generation and propagation is very important in developing advanced aircraft engines that will pass current and future noise regulations. In jet fan flows as well as exhaust plumes, two major sources of noise are present: large-scale, coherent instabilities and small-scale turbulent eddies. In previous work for the NASA Glenn Research Center, three strategies have been explored in an effort to computationally predict the noise radiation from supersonic jet exhaust plumes. In order from the least expensive computationally to the most expensive computationally, these are: 1) Linearized Euler equations (LEE). 2) Very Large Eddy Simulations (VLES). 3) Large Eddy Simulations (LES). The first method solves the linearized Euler equations (LEE). These equations are obtained by linearizing about a given mean flow and the neglecting viscous effects. In this way, the noise from large-scale instabilities can be found for a given mean flow. The linearized Euler equations are computationally inexpensive, and have produced good noise results for supersonic jets where the large-scale instability noise dominates, as well as for the tone noise from a jet engine blade row. However, these linear equations do not predict the absolute magnitude of the noise; instead, only the relative magnitude is predicted. Also, the predicted disturbances do not modify the mean flow, removing a physical mechanism by which the amplitude of the disturbance may be controlled. Recent research for isolated airfoils' indicates that this may not affect the solution greatly at low frequencies. The second method addresses some of the concerns raised by the LEE method. In this approach, called Very Large Eddy Simulation (VLES), the unsteady Reynolds averaged Navier-Stokes equations are solved directly using a high-accuracy computational aeroacoustics numerical scheme. With the addition of a two-equation turbulence model and the use of a relatively coarse grid, the numerical solution is effectively filtered into a directly calculated mean flow with the small-scale turbulence being modeled, and an unsteady large-scale component that is also being directly calculated. In this way, the unsteady disturbances are calculated in a nonlinear way, with a direct effect on the mean flow. This method is not as fast as the LEE approach, but does have many advantages to recommend it; however, like the LEE approach, only the effect of the largest unsteady structures will be captured. An initial calculation was performed on a supersonic jet exhaust plume, with promising results, but the calculation was hampered by the explicit time marching scheme that was employed. This explicit scheme required a very small time step to resolve the nozzle boundary layer, which caused a long run time. Current work is focused on testing a lower-order implicit time marching method to combat this problem.

Keith, Theo G., Jr.↗
NASA NTRS2004

Fiber-Optic Communication Links Suitable for On-Board Use in Modern Aircraft

The role of the Advanced Air Transportation Technologies program undertaken at the NASA Glenn Research Centers has been focused mainly on the improvement of air transportation safety, with particular emphasis on air transportation communication systems in on-board aircraft. The conventional solutions for digital optical communications systems specifically designed for local/metro area networks are, unfortunately, not capable of transporting the microwave and millimeter RF signals used in avionics systems. Optical networks capable of transporting RF signals are substantially different from the standard digital optical communications systems. The objective of this paper is to identify a number of different communication link architectures for RF/fiber optic transmission using a single backbone fiber for carrying VHF and UHF RF signals in the aircraft. To support these architectures, two approaches derived from both hybrid RF-optical and all-optical processing methodologies are discussed with single and multiple antennas for explicitly transporting VHF and UHF signals, while the relative merits and demerits of each architecture are also addressed. Furthermore, the experimental results of wavelength division multiplexing (WDM) link architecture from our test-bed platform, configured for aircraft environment to support simultaneous transmission of multiple RF signals over a single optical fiber, exhibit no appreciable signal degradation at wavelengths of both 1330 and 1550 nm, respectively. Our measurements of signal to noise ratio carried out for the transmission of FM and AM analog modulated signals at these wavelengths indicate that WDM is a fiber optic technology which is potentially suitable for avionics applications.

Nguyen, Hung↗