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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
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 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↗
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↗
DOE OSTI2024

High-Order Wall-Modeled Large-Eddy Simulation of High-Lift Configuration

This paper presents the assessment of several recent enhancements for a high-order wall-modeled large-eddy simulation (WMLES) approach and demonstrates order independence with a fixed data exchange location in the wall model. The two enhancements include the use of isotropic tetrahedral elements to improve accuracy and an explicit subgrid-scale model, the Vreman model, to improve accuracy and robustness. The [Formula: see text] study focused on the high-lift Common Research Model (HL-CRM) at the angle of attack of 19.57 deg, a benchmark problem from the 4th AIAA High-Lift Prediction Workshop. Solution polynomial orders of [Formula: see text], and 5 were used in the study. The study demonstrated [Formula: see text] independence in integrated forces, pitch moment, velocity profile in the wall-normal direction, and surface flow topology. It also showed that a [Formula: see text] order of at least 3 ([Formula: see text]) was needed to correctly predict the external inviscid flow and the surface flow topology. Thereafter, [Formula: see text] simulations over several other angles of attack demonstrated that the high-order WMLES approach can correctly predict the maximum lift and flow separation regions for HL-CRM with about 40 million degrees of freedom (DOF) compared to at least 250 million DOF required by second-order methods.

Engineering↗
DOE OSTI2024

Spectral Analysis of Regular Material Point Method and its Application to Study High Pressure Reverse Osmosis Membrane Compaction and Embossing

Material Point Method (MPM) is gaining widespread interest in applied continuum mechanics. The fact that all the continuum properties are stored on the particles (or material points) and the governing equations are solved on these material points makes MPM extremely suited to problems involving severe material deformations, such as crack propagation, soil movement, and fluid flows. Despite its popularity, only a few studies have focused on the numerical properties of MPM. This presentation introduces a global spectral analysis of the regular material point method. Contrary to previous studies, the analysis focuses on the numerical properties of the method in the spectral space. The amplification factor is derived as a function of the non- dimensional wave numbers. It provides insights into the stability and dissipative properties of the method for various CFL and Fourier numbers. The effect of the grid shape functions, number of particles per cell and their locations inside the grid cell are also analyzed. The EXAGOOP MPM solver (https://github.com/NREL/Exagoop.git) is developed at the National Renewable Energy Laboratory as a part of the NAWI UHPRO project and is based on the AMReX framework. A single-level, uniform cartesian grid is used as the background mesh, while the particle class in AMReX is used to manage the material point operations. Linear hat and B-splines are used as grid shape functions, while the time integration is performed using explicit Euler time integration. EXAGOOP is both CPU and GPU compatible and has been demonstrated to work well on multiple compute architectures. The performance of EXAGOOP on various computing architectures is presented along with its application to study compaction and embossing of high-pressure reverse osmosis membranes. The MPM solution accurately reproduces the membrane deformation. The deformed pore size and structure simulated using MPM also agree well with experimental SEM images.

material point method↗
DOE OSTI2022

Hidden Sectors from Multiple Line Bundles for the B−L$B-L$ MSSM

Abstract We give a formalism for constructing hidden sector bundles as extensions of sums of line bundles in heterotic M‐theory. Although this construction is generic, we present it within the context of the specific Schoen threefold that leads to the physically realistic MSSM model. We discuss the embedding of the line bundles, the existence of the extension bundle, and a number of necessary conditions for the resulting bundle to be slope‐stable and thus supersymmetric. An explicit example is presented, where two line bundles are embedded into the factor of the maximal subgroup of the hidden sector E 8 gauge group, and then enhanced to a non‐Abelian bundle by extension. For this example, there are in fact six inequivalent extension branches, significantly generalizing that space of solutions compared with hidden sectors constructed from a single line bundle.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Simulating pitch angle scattering using an explicitly solvable energy-conserving algorithm

Particle distribution functions evolving under the Lorentz operator can be simulated with the Langevin equation for pitch angle scattering. This approach is frequently used in particle based Monte-Carlo simulations of plasma collisions, among others. However, most numerical treatments do not guarantee energy conservation, which may lead to unphysical artifacts such as numerical heating and spectra distortions. In this paper, we present a novel structure-preserving numerical algorithm for the Langevin equation for pitch angle scattering. Similar to the well-known Boris algorithm, the proposed numerical scheme takes advantage of the structure-preserving properties of the Cayley transform when calculating the velocity-space rotations. The resulting algorithm is explicitly solvable, while preserving the norm of velocities down to machine precision. We demonstrate that the method has the same order of numerical convergence as the traditional stochastic Euler-Maruyama method. The numerical scheme is benchmarked by simulating the pitch angle scattering of a particle beam, and comparing with the analytical solution. Benchmark results show excellent agreement with theoretical predictions, showcasing the remarkable long time accuracy of the proposed algorithm.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

An efficient numerical model for predicting residual stress and strain in parts manufactured by laser powder bed fusion

Abstract Computational modeling of additively manufactured structures plays an increasingly important role in product design and optimization. For laser powder bed fusion processes, the accurate modeling of stress and distortion requires large amount of computational cost due to very localized heat input and evolving complex geometries. The current study takes advantage of a graphics processing unit accelerated explicit finite element analysis code and approximated heat conduction analysis to predict the macroscopic thermo-mechanical behavior in laser selective melting. Adjacent layers and tracks were lumped to reduce the number of time steps and elements in the finite element model. The effects of track and layer grouping on prediction accuracy and solution efficiency are investigated to provide a guidance for a cost-effective simulation. Thin-wall builds from Inconel alloy 625 (IN625) powders were simulated by applying the developed modeling approach to get the detailed residual stress and distortion at a computational speed 50 times higher than conventional approach. Under repeated heating and cooling cycles, a high tensile stress was produced near surfaces of a build due to a larger shrinkage on surface than that in central area. It is also shown that horizontal stresses concentrate near the root and top layers of the IN625 build. The predicted residual elastic strain distribution was validated by the experimental measurement using x-ray synchrotron diffraction.

36 MATERIALS SCIENCE↗
DOE OSTI2022

High level requirements for an IFE development program

This White Paper lays out a framework to help guide the formulation of an IFE technology program that explicitly takes account of the needs of the end-user power plant design, construction, licensing, commissioning and operational phases. This viewpoint seeks to ensure that near-term decisions are consistent with the long-term goals, and seeks to minimize effort spent on developments that are incompatible with an integrated system solution.

Dunne, Mike↗
DOE OSTI2022

The Schwarz alternating method for transient solid dynamics

Abstract In our earlier work, we formulated the Schwarz alternating method as a means for concurrent multiscale coupling in finite deformation solid mechanics for quasi‐static problems. Herein, we advance this method for the study of transient dynamic multiscale solid mechanics problems where information is exchanged back and forth between small and large scales. The extension to dynamics relies on the notion of a global time stepper. Within each global time step, the subdomains are coupled by the standard Schwarz iterative process. Remarkably, each subdomain can use its own time step or even its own time integrator to advance its solution in time, provided that they synchronize at each global time step. We study the performance of the Schwarz method on several examples designed for this purpose. Our numerical experiments demonstrate that the method is capable of coupling regions with different mesh resolutions, different element types, and different time integration schemes (e.g., implicit and explicit), all without introducing any artifacts that afflict other coupling methods for transient dynamics. Finally, we apply the dynamic Schwarz alternating method to the simulation of a bolted joint subjected to dynamic loading, as a demonstration of the performance of the method in a realistic scenario.

Mota, Alejandro↗
DOE OSTI2019

Economic analysis of a nuclear hybrid energy system in a stochastic environment including wind turbines in an electricity grid

The electricity market in the U.S. is moving toward an energy mix including more variable renewable energy (VRE) sources like wind and solar. This increasing VRE penetration is altering the profile of the net electricity demand. Normally, flexible generators (e.g. gas turbines) are used to absorb fluctuations in the electricity demand. Therefore, an increase in VRE also leads to an increase in installed capacity for such flexible generators which, in turn, may increase the price of electricity. Nuclear-Renewable Hybrid Energy Systems (N-R HES) are seen as a solution to this problem. The N-R HES contains a nuclear power plant electrically and/or thermally linked to an industrial process (IP). This paper demonstrates a new methodology and presents the corresponding RAVEN/Modelica-based software framework to evaluate the financial performance of N-R HES. The novelty of the method is that it explicitly incorporates the stochastic nature of N-R HES inputs such as electricity demand and wind speeds. The method is applied to a generic (rather than region-specific) N-R HES. The analysis shows that under the right market conditions, economic benefits can be achieved by adding a suitable IP to the base-load generator. In particular, there exists an optimal relationship between the IP capacity and the demand. Furthermore, global trends indicate VRE penetration raises the cost of electricity due to the added volatility, an issue only mitigated when operating with high profit margins for the IP.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

A pseudospectral implicit particle-in-cell method with exact energy and charge conservation

The standard particle-in-cell (PIC) method employs explicit finite-difference (FD) methods (e.g. the leap-frog scheme) for both spatial and temporal integrations. Here, we employ a pseudospectral method for solving the Poisson equation and a fully implicit time integration to achieve exact energy conservation. The advantage of a pseudospectral field solver is its spectral accuracy in solving field solutions. Earlier studies of implicit time integration of PIC FD equations can enforce exact energy exchange between field and particles, resulting in exact energy-conserving schemes. Here, we prove that the exact energy conservation property can be carried over to the pseudospectral scheme. Simultaneously, we provide a solution to ensure a pseudospectral charge continuity equation. We demonstrate the new scheme in a 2D electrostatic PIC code. In conclusion, theoretical results are confirmed via numerical examples.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Integrated Quantum-Classical Protocol for the Realistic Description of Solvated Multinuclear Mixed-Valence Transition-Metal Complexes and Their Solvatochromic Properties

Linear cyanide-bridged polymetallic complexes, which undergo photoinduced metal-to-metal charge transfer, represent prototypical systems for studying long-range electron-transfer reactions and understanding the role played by specific solute–solvent interactions in modulating the excited-state dynamics. Here, to tackle this problem, while achieving a statistically meaningful description of the solvent and of its relaxation, one needs a computational approach capable of handling large polynuclear transition-metal complexes, both in their ground and excited states, as well as the ability to follow their dynamics in several environments up to nanosecond time scales. Here, we present a mixed quantum classical approach, which combines large-scale molecular dynamics (MD) simulations based on an accurate quantum mechanically derived force field (QMD-FF) and self-consistent QMD polarized point charges, with IR and UV–vis spectral calculations to model the solvation dynamics and optical properties of a cyano-bridged trinuclear mixed-valence compound (trans-[(NC) 5 Fe III (μ-CN)Ru II (pyridine) 4 (μ-NC)Fe III (CN) 5 ] 4– ). We demonstrate the reliability of the QMD-FF/MD approach in sampling the solute conformational space and capturing the local solute–solvent interactions by comparing the results with higher-level quantum mechanics/molecular mechanics (QM/MM) MD reference data. The IR spectra calculated along the classical MD trajectories in different solvents correctly predict the red shift of the CN stretching band in the aprotic medium (acetonitrile) and the subtle differences measured in water and methanol, respectively. By explicitly including the solvent molecules around the cyanide ligands and calculating the thermal averaged absorption spectra using time-dependent density functional theory calculations within the Tamm–Dancoff approximation, the experimental solvatochromic shift is quantitatively reproduced going from water to methanol, while it is overestimated for acetonitrile. This discrepancy can likely be traced back to the lack of important dispersion interactions between the solvent cyano groups and the pyridine substituents in our micro solvation model. The proposed protocol is applied to the ground state in water, methanol, and acetonitrile and can be flexibly generalized to study excited-state nonequilibrium solvation dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗