Engineering Papers⌕ Search

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.

At least 1,117 records · Page 62

Two-Dimensional Network Flow Modeling of No-Vent Tank Filling of A Cryogenic Tank With Thermo-Dynamic Vent System Assisted Injector

A finite volume-based network flow modeling tool, Generalized Fluid System Simulation Program (GFSSP), was used to simulate the no-vent filling of a cryogenic tank with a thermodynamic vent system assisted injector. The tank was discretized into nodes in the radial and axial directions to calculate two-dimensional axisymmetric flow with heat transfer between the solid wall and fluid. The nodes are connected by branches where momentum equations are solved to calculate the flowrates. The mass and energy conservation equations and the equation of state are solved to calculate pressure, temperature, and resident mass at the nodes. The energy conservation equations at the solid nodes were solved to calculate the temperature of the solid wall. The system of equations was solved by a combination of successive substitution and simultaneous Newton-Raphson method. The numerical model accounts for a) different regimes of pool boiling heat transfer, b) condensation of vapor around the spray droplets, c) condensation of vapor at the interface of the cooled injector and ullage, and d) condensation at the liquid-vapor interphase. The numerical predictions of tank pressure, filling rate, and wall temperatures were compared with the test data. The paper also presents the effect of nodal discretization by comparing the results of the 1D and 2D models and examines the effect of the droplet diameter on tank filling.

Nodal Model↗

An approximate factorization solution of the Navier-Stokes equations for transonic flow using body-fitted coordinates with application to NACA 64A010 airfoils

The implementation of the approximate factorization algorithm and its ability to efficiently and accurately describe transonic flow about an NACA 64A010 airfoil section is examined. The approximate factorization algorithm is developed from the nondimensional, conservative, vectorized Navier-Stokes equations expressed in curvilinear coordinates. Equations of state and transport coefficient relations appropriate to atmospheric air are appended to close the system of partial differential equations. An algebraic turbulence model is also incorporated into the equation set. This algorithm was verified by investigating the flow about an NACA 64A010 airfoil at 0, 2, and 3.5 deg angle of attack for free-stream conditions of 2,000,000 Reynolds number and 0.8 Mach number. Overall results were in good qualitative agreement with wind tunnel data sets. However, while nondimensional times of six were attained, numerical difficulties prevented any case from reaching a true steady state.

Copper, G. K.↗

Comparing matching prescriptions between pre-equilibrium and hydrodynamic models in high-energy nuclear collisions

Abstract State-of-the-art simulations of high-energy nuclear collisions rely on hybrid setups, involving in particular a pre-equilibrium stage to let the system evolve from a far-from-equilibrium initial condition towards a near-equilibrated state after which fluid dynamics can be applied meaningfully. A known issue is the mismatch between the equation of state in the fluid-dynamical evolution and the effective one in the previous stage, which leads to discontinuities at the interface between the two models. Here we introduce a new matching prescription at this interface, based on the entropy, and we compare it with the standard one relying on local energy conservation. We study the behavior of various quantities at the switching time between the models and investigate a number of final-state hadronic observables. For the latter, we show that they are not modified significantly by the choice of matching prescription, provided an appropriate normalization is chosen for the initial state. In turn, our approach reduces sizeably the ratio of bulk over thermodynamic pressure at the beginning of the fluid-dynamical stage.

Borghini, Nicolas (ORCID:0000000239063258)↗

Role of back stress in the creep behavior of particle strengthened alloys

Recent developments in the interpolation of high-temperature steady-state creep results have introduced the concept that the stress dependence of the creep rate should be in terms of the effective stress referred to as the applied stress minus a back stress. This paper reports on back stresses taken from data on a gamma-prime-strengthened wrought nickel-base superalloy, an oxide dispersion-strengthened ODS nickel-base solid solution alloy, and an ODS nickel-base superalloy. The effect of air versus vacuum environments and the effect of dynamic changes in the strengthening microstructures on the magnitude of the back stress are assessed. The role of modulus normalization and the back stress correction in determining the true creep activation energy are examined. It is shown that the high values of the apparent stress exponent 'n' of the steady-state creep equation can be easily explained through a relationship between n, the true stress exponent of steady-state creep, and the stress which when subtracted from the applied stress results in the effective driving stress acting on the mobile dislocations during creep.

Purushothaman, S.↗

Impact conditions required for formation of melt by jetting in silicates

It is demonstrated that the process of jetting which occurs when particles collide at oblique angles may produce melt at much lower velocities than are required for melt production in head-on collisions. The minimum velocities of impact required for jetting in aluminum, bronzitite, dunite, and quartz are calculated by the method of shock polars. The analysis, which depends on stated assumptions about attained pressures and the occurrence of jetting, uses shock-velocity particle-velocity equations of state with three shock regimes. The treatment indicates that jetting should arise in bronzitite, dunite, and quartz at relative velocities as low as 1-2 km/sec. At such velocities material which passes near the stagnation point in the jet-forming region is subjected to sufficiently high pressures so that it is probably melted.

Kieffer, S. W.↗

Application of constrained optimization to active control of aeroelastic response

Active control of aeroelastic response is a complex in which the designer usually tries to satisfy many criteria which are often conflicting. To further complicate the design problem, the state space equations describing this type of control problem are usually of high order, involving a large number of states to represent the flexible structure and unsteady aerodynamics. Control laws based on the standard Linear-Quadratic-Gaussian (LQG) method are of the same high order as the aeroelastic plant. To overcome this disadvantage of the LQG mode, an approach developed for designing low order optimal control laws which uses a nonlinear programming algorithm to search for the values of the control law variables that minimize a composite performance index, was extended to the constrained optimization problem. The method involves searching for the values of the control law variables that minimize a basic performance index while satisfying several inequality constraints that describe the design criteria. The method is applied to gust load alleviation of a drone aircraft.

Newsom, J. R.↗

Constrained optimization techniques for active control of aeroelastic response

Active control of aeroelastic response is a complex problem in which the designer usually tries to satisfy many design criteria which are often conflicting in nature. To further complicate the design problem, the state space equations describing this type of control problem are usually of high order, involving a large number of states to represent the flexible structure and unsteady aerodynamics. Control laws based on the standard Linear - Quadratic - Gaussian method are of the same high order as the aeroelastic plant and may be difficult to implement in the flight computer. To overcome this disadvantage a new approach was developed for designing low-order optimized robust control laws. In this approach, a nonlinear programming algorithm is used to search for the values of control law design variables that minimize a performance index while satisfying several inequality constraints that describe the design criteria on the stability robustness and responses. The method is applied to a gust load alleviation problem and a stability robustness improvement problem of a drone aircraft.

Mukhopadhyay, Vivekananda↗

Bulk-Flow Analysis of Hybrid Thrust Bearings for Advanced Cryogenic Turbopumps

A bulk-flow analysis and computer program for prediction of the static load performance and dynamic force coefficients of angled injection, orifice-compensated hydrostatic/hydrodynamic thrust bearings have been completed. The product of the research is an efficient computational tool for the design of high-speed thrust bearings for cryogenic fluid turbopumps. The study addresses the needs of a growing technology that requires of reliable fluid film bearings to provide the maximum operating life with optimum controllable rotordynamic characteristics at the lowest cost. The motion of a cryogenic fluid on the thin film lands of a thrust bearing is governed by a set of bulk-flow mass and momentum conservation and energy transport equations. Mass flow conservation and a simple model for momentum transport within the hydrostatic bearing recesses are also accounted for. The bulk-flow model includes flow turbulence with fluid inertia advection, Coriolis and centrifugal acceleration effects on the bearing recesses and film lands. The cryogenic fluid properties are obtained from realistic thermophysical equations of state. Turbulent bulk-flow shear parameters are based on Hirs' model with Moody's friction factor equations allowing a simple simulation for machined bearing surface roughness. A perturbation analysis leads to zeroth-order nonlinear equations governing the fluid flow for the thrust bearing operating at a static equilibrium position, and first-order linear equations describing the perturbed fluid flow for small amplitude shaft motions in the axial direction. Numerical solution to the zeroth-order flow field equations renders the bearing flow rate, thrust load, drag torque and power dissipation. Solution to the first-order equations determines the axial stiffness, damping and inertia force coefficients. The computational method uses well established algorithms and generic subprograms available from prior developments. The Fortran9O computer program hydrothrust runs on a Windows 95/NT personal computer. The program, help files and examples are licensed by Texas A&M University Technology License Office. The study of the static and dynamic performance of two hydrostatic/hydrodynamic bearings demonstrates the importance of centrifugal and advection fluid inertia effects for operation at high rotational speeds. The first example considers a conceptual hydrostatic thrust bearing for an advanced liquid hydrogen turbopump operating at 170,000 rpm. The large axial stiffness and damping coefficients of the bearing should provide accurate control and axial positioning of the turbopump and also allow for unshrouded impellers, therefore increasing the overall pump efficiency. The second bearing uses a refrigerant R134a, and its application in oil-free air conditioning compressors is of great technological importance and commercial value. The computed predictions reveal that the LH2 bearing load capacity and flow rate increase with the recess pressure (i.e. increasing orifice diameters). The bearing axial stiffness has a maximum for a recess pressure rati of approx. 0.55. while the axial damping coefficient decreases as the recess pressure ratio increases. The computer results from three flow models are compared. These models are a) inertialess, b) fluid inertia at recess edges only, and c) full fluid inertia at both recess edges and film lands. The full inertia model shows the lowest flow rates, axial load capacity and stiffness coefficient but on the other hand renders the largest damping coefficients and inertia coefficients. The most important findings are related to the reduction of the outflow through the inner radius and the appearance of subambient pressures. The performance of the refrigerant hybrid thrust bearing is evaluated at two operating speeds and pressure drops. The computed results are presented in dimensionless form to evidence consistent trends in the bearing performance characteristics. As the applied axial load increases, the bearing film thickness and flow rate decrease while the recess pressure increases. The axial stiffness coefficient shows a maximum for a certain intermediate load while the damping coefficient steadily increases. The computed results evidence the paramount of centrifugal fluid inertia at low recess pressures (i.e. low loads), and where there is actually an inflow through the bearing inner diameter, accompanied by subambient pressures just downstream of the bearing recess edge. These results are solely due to centrifugal fluid inertia and advection transport effects. Recommendations include the extension of the computer program to handle flexure pivot tilting pad hybrid bearings and the ability to calculate moment coefficients for shaft angular misalignments.

SanAndres, Luis↗

A Potential Function Derivation of a Constitutive Equation for Inelastic Material Response

Physical and thermodynamic concepts are used to develop a potential function for application to high temperature polycrystalline material response. Inherent in the formulation is a differential relationship between the potential function and constitutive equation in terms of the state variables. Integration of the differential relationship produces a state variable evolution equation that requires specification of the initial value of the state variable and its time derivative. It is shown that the initial loading rate, which is directly related to the initial hardening rate, can significantly influence subsequent material response. This effect is consistent with observed material behavior on the macroscopic and microscopic levels, and may explain the wide scatter in response often found in creep testing.

Stouffer, D. C.↗

A mathematical simulation model of the CH-47B helicopter, volume 2

A nonlinear simulation model of the CH-47B helicopter, was adapted for use in a simulation facility. The model represents the specific configuration of the variable stability CH-47B helicopter. Modeling of the helicopter uses a total force approach in six rigid body degrees of freedom. Rotor dynamics are simulated using the Wheatley-Bailey equations, steady state flapping dynamics and included in the model of the option for simulation of external suspension, slung load equations of motion. Validation of the model was accomplished by static and dynamic data from the original Boeing Vertol mathematical model and flight test data. The model is appropriate for use in real time piloted simulation and is implemented on the ARC Sigma IX computer where it may be operated with a digital cycle time of 0.03 sec.

Weber, J. M.↗

The effect of suspended particles on Jean's criterion for gravitational instability

The effect that the proper inclusion of suspended particles has on Jeans' criterion for the self-gravitational instability of an unbounded nonrotating adiabatic gas cloud is examined by formulating the appropriate model system, introducing particular physically plausible equations of state and constitutive relations, performing a linear stability analysis of a uniformly expanding exact solution to these governing equations, and exploiting the fact that there exists a natural small material parameter for this problem given by N sub 1/n sub 1, the ratio of the initial number density for the particles to that for the gas. The main result of this investigation is the derivation of an altered criterion which can substantially reduce Jeans' original critical wavelength for instability. It is then shown that the existing discrepancy between Jeans' theoretical prediction using and actual observational data relevant to the Andromeda nebula M31 can be accounted for by this new criterion of assuming suspended particles of a reasonable grain size and distribution to be present.

Wollkind, David J.↗

eXTP: Enhanced X-Ray Timing and Polarimetry Mission

eXTP is a science mission designed to study the state of matter under extreme conditions of density, gravity and magnetism. Primary goals are the determination of the equation of state of matter at supra-nuclear density, the measurement of QED effects in highly magnetized star, and the study of accretion in the strong-field regime of gravity. Primary targets include isolated and binary neutron stars, strong magnetic field systems like magnetars, and stellar-mass and supermassive black holes. The mission carries a unique and unprecedented suite of state-of-the-art scientific instruments enabling for the first time ever the simultaneous spectral-timing-polarimetry studies of cosmic sources in the energy range from 0.5-30 keV (and beyond). Key elements of the payload are: the Spectroscopic Focusing Array (SFA) - a set of 11 X-ray optics for a total effective area of approx. 0.9 m(exp. 2) and 0.6 m(exp. 2) at 2 keV and 6 keV respectively, equipped with Silicon Drift Detectors offering less than 180 eV spectral resolution; the Large Area Detector (LAD) - a deployable set of 640 Silicon Drift Detectors, for a total effective area of approx. 3.4 m(exp. 2), between 6 and 10 keV, and spectral resolution better than 250 eV; the Polarimetry Focusing Array (PFA) - a set of 2 X-ray telescope, for a total effective area of 250 cm(exp. 2) at 2 keV, equipped with imaging gas pixel photoelectric polarimeters; the Wide Field Monitor (WFM) - a set of 3 coded mask wide field units, equipped with position-sensitive Silicon Drift Detectors, each covering a 90 degrees x 90 degrees field of view. The eXTP international consortium includes major institutions of the Chinese Academy of Sciences and Universities in China, as well as major institutions in several European countries and the United States. The predecessor of eXTP, the XTP mission concept, has been selected and funded as one of the so-called background missions in the Strategic Priority Space Science Program of the Chinese Academy of Sciences since 2011. The strong European participation has significantly enhanced the scientific capabilities of eXTP. The planned launch date of the mission is earlier than 2025.

X-ray astronomy↗

Relativistic core–valence-separated equation-of-motion coupled-cluster singles and doubles method: Efficient implementation and benchmark calculations

An efficient implementation for the relativistic exact two-component core–valence-separated equation-of-motion coupled-cluster singles and doubles (X2C-CVS-EOM-CCSD) method is reported. The explicit exclusion of pure valence excitations in the EOM-CCSD excited-state eigenvalue equations significantly improves the efficiency for calculations of core-excited states. Benchmark relativistic CVS-EOM-CC calculations with systematic inclusion of relativistic, correlation, and basis-set effects are shown to provide highly accurate results for core ionized and excited states involving heavy atoms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The importance of few-nucleon forces in chiral effective field theory

We study the importance of few-nucleon forces in chiral effective field theory for describing many-nucleon systems. A combinatorial argument suggests that three-nucleon forces-which are conventionally regarded as next-to-next-to-leading order-should accompany the two-nucleon force already at leading order (LO) starting with mass number $A$ ≃10–20. We find that this promotion enables the first realistic description of the 16 O ground state based on a renormalization-group-invariant LO interaction. We also performed coupled-cluster calculations of the equation of state for symmetric nuclear matter and our results indicate that LO four-nucleon forces could play a crucial role for describing heavy-mass nuclei. In conclusion, the enhancement mechanism we found is very general and could be important also in other many-body problems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

R -mode Stability of GW190814's Secondary Component as a Supermassive and Superfast Pulsar

The nature of GW190814's secondary component m 2 of mass (2.50–2.67)M ⊙ in the mass gap between the currently known maximum mass of neutron stars and the minimum mass of black holes is currently under hot debate. Among the many possibilities proposed in the literature, the m 2 was suggested as a superfast pulsar while its r-mode stability against the run-away gravitational radiation through the Chandrasekhar-Friedman-Schutz mechanism is still unknown. Previously, Fortin et al. constructed a sample of 33 unified equations of state using the same nuclear interactions from the crust to the core consistently; from that sample we use those equations that fulfill all currently known astrophysical and nuclear physics constraints to compare the minimum frequency required for m 2 to rotationally sustain a mass greater than 2.50 M ⊙ with the critical frequency above which the r-mode instability occurs. Here, we use two extreme damping models assuming that the crust is either perfectly rigid or elastic. Using the stability of 19 observed low-mass X-ray binaries as an indication that the rigid crust damping of the r-mode dominates within the models studied, we find that m 2 is r-mode-stable while rotating with a frequency higher than 870.2 Hz (0.744 times its Kepler frequency of 1169.6 Hz) as long as its temperature is lower than about 3.9 × 10 7 K, further supporting the proposal that GW190814's secondary component is a supermassive and superfast pulsar.

79 ASTRONOMY AND ASTROPHYSICS↗

Magnetohydrodynamic equilibrium. III - Helically symmetric fields. IV - Nonequilibrium of nonsymmetric hydrodynamic topologies

It is pointed out that plasma confinement in stable equilibrium states constitutes a fundamental and still unresolved question in plasma astrophysics and thermonuclear fusion research. The problem has two parts related to the equilibrium states themselves and their mechanical stability. The question of the existence of general solutions of the field and fluid equations for the steady dynamical interaction of inviscid compressible fluids of high electrical conductivity with magnetic and gravity fields is considered. In the absence of fluid motions, the presented equations become the familiar equations of magnetostatics supplemented by an equation of state. Starting from this simplest case of magnetostatic equilibrium, the investigation proceeds to the more complex case of magnetohydrodynamic equilibrium. Examples of helically symmetric fields are presented to illustrate the use of the formulation for treating the dynamics of helically symmetric hydromagnetic flows.

Tsinganos, K. C.↗

A Generalized Fluid Formulation for Turbomachinery Computations

A generalized formulation of the equations of motion of an arbitrary fluid are developed for the purpose of defining a common iterative algorithm for computational procedures. The method makes use of the equations of motion in conservation form with separate pseudo-time derivatives used for defining the numerical flux for a Riemann solver and the convergence algorithm. The partial differential equations are complemented by an thermodynamic and caloric equations of state of a complexity necessary for describing the fluid. Representative solutions with a new code based on this general equation formulation are provided for three turbomachinery problems. The first uses air as a working fluid while the second uses gaseous oxygen in a regime in which real gas effects are of little importance. These nearly perfect gas computations provide a basis for comparing with existing perfect gas code computations. The third case is for the flow of liquid oxygen through a turbine where real gas effects are significant. Vortex shedding predictions with the LOX formulations reduce the discrepancy between perfect gas computations and experiment by approximately an order of magnitude, thereby verifying the real gas formulation as well as providing an effective case where its capabilities are necessary.

Merkle, Charles L.↗

Classical and quantum simulations of 1+1-dimensional ${\mathbb{Z}}_{2}$ gauge theory at finite temperature and density

Simulating strongly coupled gauge theories at finite temperature and density is a longstanding challenge in nuclear and high-energy physics with fundamental implications for condensed matter physics. Here, we simulate such systems using minimally entangled typical thermal state (METTS) approaches, which combine classical random sampling with imaginary-time evolution, implementable on either classical or quantum computers, to estimate thermal averages of observables. We study 1+1-dimensional ${\mathbb{Z}}_{2}$ gauge theory coupled to spinless fermionic matter, which maps onto a local quantum spin chain. We benchmark both a classical matrix-product-state implementation of METTS and a recently proposed adaptive variational approach for near-term quantum devices, focusing on the equation of state and measures of fermion confinement. Of particular importance is the choice of basis for METTS sampling, which impacts both the sampling overhead and quantum circuit complexity. Our work sets the stage for future studies of strongly coupled gauge theories using classical and quantum hardware.

Chen, I-Chi [Iowa State Univ., Ames, IA (United St↗