Appraisal of an iterative method for bound states
Radial Schrodinger equation bound state eigenvalues and properties by iterative method, calculating Coulomb potentials
Anderson, N.↗
SEARCH · Engineering Papers
Results for “equations of state”
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.
On the Existence of Steady-State Solutions to the Equations Governing Fluid Flow in Networks
The steady-state solution of fluid flow in pipeline infrastructure networks driven by junction/node potentials is a crucial ingredient in various decision-support tools for system design and operation. While the nonlinear system is known to have a unique solution (when one exists), the absence of a definite result on the existence of solutions hobbles the development of computational algorithms, for it is not possible to distinguish between algorithm failure and non-existence of a solution. In this letter, we show that for any fluid whose equation of state is a scaled monomial, a unique solution exists for such nonlinear systems if the term solution is interpreted in terms of potentials and flows rather than pressures and flows. However, for gases following the CNGA equation of state, while the question of existence remains open, we construct an alternative system that always has a unique solution and show that the solution to this system is a good approximant of the true solution. Further, the existence result for flow of natural gas in networks also applies to other fluid flow networks such as water distribution networks or networks that transport carbon dioxide in carbon capture and sequestration. Most importantly, our result enables correct diagnosis of algorithmic failure, problem stiffness, and non-convergence in computational algorithms.
42 ENGINEERING↗
Stringent constraints on neutron-star radii from multimessenger observations and nuclear theory
The properties of neutron stars are determined by the nature of the matter that they contain. These properties can be constrained by measurements of the star’s size. We obtain stringent constraints on neutron-star radii by combining multimessenger observations of the binary neutron-star merger GW170817 with nuclear theory that best accounts for density-dependent uncertainties in the equation of state. We construct equations of state constrained by chiral effective field theory and marginalize over these using the gravitational-wave observations. Combining this with the electromagnetic observations of the merger remnant that imply the presence of a short-lived hypermassive neutron star, we find that the radius of a 1.4 M⊙ neutron star is $R_{1.4M_⊙} = 11.0^{+0.9}_{-0.6} km$ (90% credible interval). Finally, using this constraint, we show that neutron stars are unlikely to be disrupted in neutron star–black hole mergers; subsequently, such events will not produce observable electromagnetic emission.
79 ASTRONOMY AND ASTROPHYSICS↗
On the dynamics of agglomerating ensembles of clouds
It is shown that for dynamical time scales longer than 10 million years the ensemble of interstellar cold diffuse clouds may be described by equations of mass and momentum conservation where the pressure has the form p/v-squared/ where /v-squared/ is the mass averaged square velocity dispersion. Equilibrium cloud spectra are derived and an equation of state obtained. The equation of state undergoes a phase transition at a critical density close to the local interstellar gas density. At higher densities, more massive clouds are formed and the cloud fluid softens. It is suggested that this mechanism may naturally regulate the interstellar gas density if star formation is very efficient in the cooler high-density systems.
Cowie, L. L.↗
Numerical solution of a coupled pair of elliptic equations from solid state electronics
Iterative methods are considered for the solution of a coupled pair of second order elliptic partial differential equations which arise in the field of solid state electronics. A finite difference scheme is used which retains the conservative form of the differential equations. Numerical solutions are obtained in two ways, by multigrid and dynamic alternating direction implicit methods. Numerical results are presented which show the multigrid method to be an efficient way of solving this problem.
Phillips, T. N.↗
Numerical solution of a coupled pair of elliptic equations from solid state electronics
Iterative methods are considered for the solution of a coupled pair of second order elliptic partial differential equations which arise in the field of solid state electronics. A finite difference scheme is used which retains the conservative form of the differential equations. Numerical solutions are obtained in two ways, by multigrid and dynamic alternating direction implicit methods. Numerical results are presented which show the multigrid method to be an efficient way of solving this problem. Previously announced in STAR as N83-30109
Phillips, T. N.↗
Monoenergetic-source solutions of the steady-state cosmic-ray equation of transport
Analytic, steady state, solutions of the cosmic-ray equation of transport for the solar cavity are given for monoenergetic particles released continuously at a specified heliocentric radius. The solutions are used to illustrate energy change processes. A new feature is an apparent buildup of particles at very low energies.
Webb, G. M.↗
Simulation of Methane Oxycombustion in Supercritical Carbon Dioxide
The Allam-Fetvedt Cycle (AFC) offers the potential for electricity generation with minimal or even negative carbon emissions by using oxycombustion in supercritical carbon dioxide (sCO2). Because of the extreme high temperatures and pressures encountered in these systems, simulation of oxycombustion requires use of a complex equation of state such as the modified Soave-Redlich-Kwong (SRK) equation, coupled with equally complex correlations for transport properties. The current work simulates a model combustor using both an ideal gas equation of state and the SRK equation of state. The model combustor features a jet of CH4 mixing with and reacting with O2 in the sCO2 ambient. Two jet configurations are considered: one of pure CH4, and one of CH4 mixed with sCO2. Comparison of results shows that key output properties such as combustor temperature and flow distribution are impacted by the choice of equations of state. Inspection of simulation results shows that mixing and reaction rates, as well as the uniformity of the flow, are impacted by the choice of equation of state. These results show that even though use of a complex equation of state greatly increases the computational costs of these simulations, it is critical to accurately capturing the physics of these systems.
computational fluid dynamics↗
A Posteriori Study of a DNS Database Describing Super critical Binary-Species Mixing
Currently, the modeling of supercritical-pressure flows through Large Eddy Simulation (LES) uses models derived for atmospheric-pressure flows. Those atmospheric-pressure flows do not exhibit the particularities of high densitygradient magnitude features observed both in experiments and simulations of supercritical-pressure flows in the case of two species mixing. To assess whether the current LES modeling is appropriate and if found not appropriate to propose higher-fidelity models, a LES a posteriori study has been conducted for a mixing layer that initially contains different species in the lower and upper streams, and where the initial pressure is larger than the critical pressure of either species. An initially-imposed vorticity perturbation promotes roll-up and a double pairing of four initial span-wise vortices into an ultimate vortex that reaches a transitional state. The LES equations consist of the differential conservation equations coupled with a real-gas equation of state, and the equation set uses transport properties depending on the thermodynamic variables. Unlike all LES models to date, the differential equations contain, additional to the subgrid scale (SGS) fluxes, a new SGS term that is a pressure correction in the momentum equation. This additional term results from filtering of Direct Numerical Simulation (DNS) equations, and represents the gradient of the difference between the filtered pressure and the pressure computed from the filtered flow field. A previous a priori analysis, using a DNS database for the same configuration, found this term to be of leading order in the momentum equation, a fact traced to the existence of high-densitygradient magnitude regions that populated the entire flow; in the study, models were proposed for the SGS fluxes as well as this new term. In the present study, the previously proposed constantcoefficient SGS-flux models of the a priori investigation are tested a posteriori in LES, devoid of or including, the SGS pressure correction term. The present pressure-correction model is different from, and more accurate as well as less computationally intensive than that of the a priori study. The constant-coefficient SGS-flux models encompass the Smagorinsky (SMC), in conjunction with the Yoshizawa (YO) model for the trace, the Gradient (GRC) and the Scale Similarity (SSC) models, all exercised with the a priori study constant coefficients calibrated at the transitional state. The LES comparison is performed with the filtered- and-coarsened (FC) DNS, which represents an ideal LES solution. Expectably, an LES model devoid of SGS terms is shown to be considerably inferior to models containing SGS effects. Among models containing SGS effects, those including the pressure-correction term are substantially superior to those devoid of it. The sensitivity of the predictions to the initial conditions and grid size are also investigated. Thus, it has been discovered that, additional to the atmospheric-pressure models currently used, a new model is necessary to simulate supercritical-pressure flows. This model depends on the thermodynamic characteristics of the chemical species involved.
Bellan, Josette↗
Solving rational matrix equations in the state space with applications to computer-aided control-system design
A method of solving a class of linear matrix equations over various rings is proposed, using results from linear geometric control theory. An algorithm, successfully implemented, is presented, along with non-trivial numerical examples. Applications of the method to the algebraic control system design methodology are discussed.
Packard, A. K.↗
Axion baryogenesis puts a new spin on the Hubble tension
We show that a rotating axion field that makes a transition from a matterlike equation of state to a kinationlike equation of state around the epoch of recombination can significantly ameliorate the Hubble tension, i.e., the discrepancy between the determinations of the present-day expansion rate H 0 from observations of the cosmic microwave background on one hand and type Ia supernovae on the other. We consider a specific, UV-complete model of such a rotating axion and find that it can relax the Hubble tension without exacerbating tensions in determinations of other cosmological parameters, in particular the amplitude of matter fluctuations S 8 . We subsequently demonstrate how this rotating axion model can also generate the baryon asymmetry of our Universe, by introducing a coupling of the axion field to right-handed neutrinos. This baryogenesis model predicts heavy neutral leptons that are most naturally within reach of future lepton colliders, but in finely tuned regions of parameter space may also be accessible at the high-luminosity LHC and the beam dump experiment SHiP. Published by the American Physical Society 2024
Co, Raymond T. (ORCID:0000000283957056)↗
ARES v1.x - Performance Portable Tool to Simulate Supernovae based on Parthenon Framework
Historically, codes for simulating supernovae (such as Arepo, FLASH or LEAFS) have been at the forefront of scientific high-performance computing to the immense computational resources required for full 3D simulations. However, given the shift towards heterogenous HPC architectures, many current-generation codes are at the risk of losing their competitiveness as they are only designed to run on homogeneous CPU-only systems. There exist several efforts to enable these codes for GPU’s, however, these efforts only consider specific architectures or vendors (e.g., implement only CUDA or HIP), limiting themselves to a small range of exascale computing systems. Frameworks such as Kokkos aim to provide a framework which is agnostic of the targeted architecture, enabling the development of performant and portable code. In the Ares code, we develop a performance portable tool to simulate supernovae based on the Parthenon Framework, which in turn uses Kokkos in the background. Here, the Parthenon Framework provides an interface to the underlying mesh-refinement routines, which form the backbone of our code. In addition, we incorporate the already existing Singularity-EOS toolkit to provide us with various equations of state, primarily the Helmholtz equation of state. We also include the JINA Reaclib as a basis for our nuclear network solver. Finally, we implement a gravity solver to complete the required physics. This setup will provide us with a minimal code base to simulate supernova in a similar style to the tried-and-tested Arepo code, but in a futureproof performance portable framework.
Lim, Hyun↗
Thermodynamic properties of oxygen and nitrogen III
The final equation for nitrogen was determined. In the work on the equation of state for nitrogen, coefficients were determined by constraining the critical point to selected critical point parameters. Comparisons of this equation with all the P-density-T data were made, as well as comparisons to all other thermodynamic data reported in the literature. The extrapolation of the equation of state was studied for vapor to higher temperatures and lower temperatures, and for the liquid surface to the saturated liquid and the fusion lines. A new vapor pressure equation was also determined which was constrained to the same critical temperature, pressure, and slope (dP/dT) as the equation of state. Work on the equation of state for oxygen included studies for improving the equation at the critical point. Comparisons of velocity of sound data for oxygen were also made between values calculated with a preliminary equation of state and experimental data. Functions for the calculation of the derived thermodynamic properties using the equation of state are given, together with the derivative and integral functions for the calculation of the thermodynamic properties using the equations of state. Summary tables of the thermodynamic properties of nitrogen and oxygen are also included to serve as a check for those preparing computer programs using the equations of state.
Stewart, R. B.↗
Use of rate equations to describe laser excitation in flames
The region of validity of conventional rate equation formulations used to describe the excitation of flame gases for spectroscopic diagnostic purposes is examined by direct comparison with the density matrix formulation. The density matrix equations are formulated for an ensemble of atoms or molecules in the gas state. These equations are reduced to the steady-state limit for a two-level system under the effect of single-mode and multimode optical fields. The steady-state equations are shown to be valid only when the laser pulse rises slowly with respect to the collisional dephasing rate for single-mode excitation, and for certain special cases under multimode excitation.
Daily, J. W.↗
A non‐intrusive domain‐decomposition model reduction method for linear steady‐state partial differential equations with random coefficients
Abstract Domain decomposition methods have been proved to be an effective strategy to reduce the dimension of parametric partial differential equations (PDEs). However, existing domain decomposition methods for parametric PDEs are usually intrusive, which means domain decomposition based solvers need to be implemented from scratch for each target parametric PDE. To address this issue, we develop a new non‐intrusive domain‐decomposition model reduction method for linear steady‐state PDEs with random‐field coefficients. As a variant of our previous work by Mu and Zhang, the new method only needs access to the final linear system, that is, the global stiffness matrix and the right hand side, of a deterministic PDE solver, in order to build a domain‐decomposition‐based reduced model without intrusive implementation from scratch. The key idea is to remove the interface condition between sub‐domains and rely on the correlation between columns of the linear system to couple the sub‐domains. The non‐intrusive feature enables the applicability of the proposed method to a broader class of uncertainty quantification problems, where many legacy codes/solvers can be fully reused by our method. Two numerical examples including diffusion equations with random diffusivity and convection‐dominated transport with random velocity, are provided to demonstrate the effectiveness and efficiency of our method.
Zhang, Guannan↗
Marangoni Effects of a Drop in an Extensional Flow: The Role of Surfactant Physical Chemistry
While the changes in stresses caused by surfactant adsorption on non-deforming interfaces have been fairly well established, prior to this work, there were few studies addressing how surfactants alter stresses on strongly deforming interfaces. We chose the model problem of a drop in a uniaxial extensional flow to study these stress conditions To model surfactant effects at fluid interfaces, a proper description of the dependence of the surface tension on surface concentration, the surface equation of state, is required. We have adopted a surface equation of state that accounts for the maximum coverage limit; that is, because surfactants have a finite cross sectional area, there is an upper bound to the amount of surfactant that can adsorb in a monolayer. The surface tension reduces strongly only when this maximum coverage is approached. Since the Marangoni stresses go as the derivative of the surface equation of state times the surface concentration gradient, the non-linear equation of state determines both the effect of surfactants in the normal stress jump, (which is balanced by the product of the mean curvature of the interface times the surface tension), and the tangential stress jump, which is balanced by Marangoni stresses. First, the effects of surface coverage and intermolecular interactions among surfactants which drive aggregation of surfactants in the interface were studied. (see Pawar and Stebe, Physics of Fluids).
Stebe, Kathleen J.↗
Determination of rotor derivatives and rotor hub force and moment derivatives from flight measurements with the RSRA compound helicopter
A case study of the application of an integrated rotorcraft identification method to the linear modeling of rotor system dynamics and rotor hub loads is presented. Applying the method to flight data obtained from the RSRA compound helicopter, the rotor derivatives in the rotor state dynamic equation are identified along with blade equivalent damping and spring periodic coefficients. A rigid blade flapping equation of motion can be derived from the identified rotor state equation. It is shown that the concept of rotor hub load derivatives is useful for small maneuvering loads. The hub load derivatives can be used to relate the hub dynamic load to rotor/fuselage motion and applied inputs.
Wang, Ji C.↗
Partition function corrections
The vibration-rotation interaction coupling in diatomic molecules is analyzed and first-order expressions for the energy levels are deduced. Corrections to the partition function to account for this coupling are derived. The effects of long-range intermolecular forces due to induced dipole-induced dipole interactions, dipole-induced dipole interactions, charge-dipole interactions, and charge-charge interactions are considered. London's quantized oscillator model of the induced dipole-induced dipole interaction is developed. The partition function including effects of such dispersion forces is derived and is shown to lead to the virial equation of state. For a model combining the hard sphere and long-range attractive type potentials, this is found to lead to the van der Waals equation of state. The equation of state near the critical point is presented in reduced form (thermodynamic variables expressed in units of critical point values). Finally, approximate corrections to the electronic partition functions are derived to account for perturbations of excited electronic states induced by neighboring neutral particles, heavy ions, and fast electron collisions.
Source record↗