The application of state-variable techniques to communication and radar problems
Conference proceedings on applying differential equation or state variable techniques to radar and communication system problems
RADAR↗
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.
Multiframe X-ray diffraction on the OMEGA EP laser
A method of measuring multiple temporally resolved powder x-ray diffraction measurements during a single dynamic compression laser experiment has been demonstrated on the OMEGA EP laser. Up to four time-gated snapshots of the x-ray diffraction pattern are collected on an ultrafast x-ray framing camera with a 2θ scattering angle coverage of 48°–65° and 400-ps temporal resolution with the potential for 10°–160° and 200-ps resolution. X-ray diffraction snapshots capture the transformation to a high-pressure β phase in ramp-compressed Zr, constraining the atomic structure and pressure–density equation-of-state in a single shot. Here, this platform will be used to constrain transition pressures, understand the role of kinetics, and verify the existence of hysteresis and metastability on phase transformations at the lattice level under single-shot dynamic compression.
47 OTHER INSTRUMENTATION↗
Deviations from the Isobaric Multiplet Mass Equation due to Threshold States
Recent studies have completed the A=16 isospin quintets for states with J π =0 + and 2 + . The dependence of their masses as a function of isospin projection shows evidence for deviations from quadratic behavior indicating isospin violation beyond the expectation from two-body forces. The deviation is most pronounced for the 2 + states. Predictions from the shell model embedded in the continuum (SMEC) allow us to explain that this isospin violation is associated with a modification of the nuclear structure due to the open-quantum-system nature of the proton-rich members of the quintet. In particular, the 0 + and 2 + states in 16 Ne and the 2 + state in 16 F are threshold resonances located just above a proton-decay threshold where s-wave coupling to the continuum is expected. The measured deviations of these threshold states from the quadratic behavior of the remaining members of the multiplets makes it possible to obtain information on the magnitude and the energy dependence of the continuum-coupling energy correction. Finally, continuum coupling is also indicated for the ground state of 8 C, but this time through p-wave coupling
coulomb energies & analogue states↗
A finite element model of conduction, convection, and phase change near a solid/melt interface
Detailed understanding of heat transfer and fluid flow is required for many aerospace thermal systems. These systems often include phase change and operate over a range of accelerations or effective gravitational fields. An approach to analyzing such systems is presented which requires the simultaneous solution of the conservation laws of energy, momentum, and mass, as well as an equation of state. The variable property form of the governing equations are developed in two-dimensional Cartesian coordinates for a Newtonian fluid. A numerical procedure for solving the governing equations is presented and implemented in a computer program. The Galerkin form of the finite element method is used to solve the spatial variation of the field variables, along with the implicit Crank-Nicolson time marching algorithm. Quadratic Langrangian elements are used for the internal energy and the two components of velocity. Linear Lagrangian elements are used for the pressure. The location of the solid/liquid interface as well as the temperatures are determined form the calculated internal energy and pressure. This approach is quite general in that it can describe heat transfer without phase change, phase change with a sharp interface, and phase change without an interface. Analytical results from this model are compared to those of other researchers studying transient conduction, convection, and phase change and are found to be in good agreement. The numerical procedure presented requires significant computer resources, but this is not unusual when compared to similar studies by other researchers. Several methods are suggested to reduce the computational times.
Viterna, Larry A.↗
The computation of standard solar models
Procedures for calculating standard solar models with the usual simplifying approximations of spherical symmetry, no mixing except in the surface convection zone, no mass loss or gain during the solar lifetime, and no separation of elements by diffusion are described. The standard network of nuclear reactions among the light elements is discussed including rates, energy production and abundance changes. Several of the equation of state and opacity formulations required for the basic equations of mass, momentum and energy conservation are presented. The usual mixing-length convection theory is used for these results. Numerical procedures for calculating the solar evolution, and current evolution and oscillation frequency results for the present sun by some recent authors are given.
Ulrich, Roger K.↗
Optimization methods, flux conserving methods for steady state Navier-Stokes equation
Navier-Stokes equation as discretized by new flux conserving method proposed by Chang and Scott results in the system: vector F(vector x) = 0, where F is a vector valued function. The Optimization method we use is based on Quasi-Newton methods: given a nonlinear function vector F(vector x) = 0, we solve, Delta(vector x) = -BF(vector x), where Delta(vector x) is the correction term and B is the inverse Jacobian of F(x). Then, iteratively, vector(x(sub (i+1))) = vector(x (sub i)) + alpha.Delta(vector x(sub i)), where alpha is a line search correction term determined by a line search routine. We use the BFCG's update the Jacobian matrix B(sub k) at each iteration. It is well known that B(sub k) approaches B(*) at the solution X(*). This algorithm has several advantages over the Newton-Raphson method. For example, we do not need to calculate the Jacobian matrix at each iteration which is computationally very expensive.
Adeyeye, John↗
A Second Law Based Unstructured Finite Volume Procedure for Generalized Flow Simulation
An unstructured finite volume procedure has been developed for steady and transient thermo-fluid dynamic analysis of fluid systems and components. The procedure is applicable for a flow network consisting of pipes and various fittings where flow is assumed to be one dimensional. It can also be used to simulate flow in a component by modeling a multi-dimensional flow using the same numerical scheme. The flow domain is discretized into a number of interconnected control volumes located arbitrarily in space. The conservation equations for each control volume account for the transport of mass, momentum and entropy from the neighboring control volumes. In addition, they also include the sources of each conserved variable and time dependent terms. The source term of entropy equation contains entropy generation due to heat transfer and fluid friction. Thermodynamic properties are computed from the equation of state of a real fluid. The system of equations is solved by a hybrid numerical method which is a combination of simultaneous Newton-Raphson and successive substitution schemes. The paper also describes the application and verification of the procedure by comparing its predictions with the analytical and numerical solution of several benchmark problems.
Majumdar, Alok↗
Integration of the three dimensional heat conduction equation
Numerical integration of steady state heat conduction equations representing diurnal variations in upper atmosphere
Harris, I.↗
Hydrodynamic expansion and near-infrared absorption of x-ray heated aluminum plasmas
We use x-ray pulses from dense argon plasmas at the Z Machine (Sandia National Laboratories) to generate hypersonic aluminum plasmas akin to material ejecta during proposed planetary defense missions, fusion reactor wall excursions, and other high-energy density processes. Near-infrared absorption is used to diagnose the controlled expansion of the plasmas through cylindrical cavities following their generation from x-ray heating of solid aluminum 7075 alloy. The data are compared to multidimensional radiation hydrodynamics simulations utilizing the ALEGRA multiphysics code, accounting for the dynamics of radiation scattering, material phase change, plasma expansion, thermal re-irradiation, and interactions with the cavity and with the infrared beams. To allow for accurate simulation, density functional theory is used to apply the Hagen–Rubens relation for the far-infrared and is adjoined with a detailed configuration accounting model using the Propaceos code, producing opacities spanning 10 −1 –10 4 eV photon energy for aluminum 7075 alloy, and in comparison with pure aluminum. The model is found to agree with experimental data in the higher-fluence regime when the Hagen–Rubens relation is applied. The ejected material, which is observed to travel up to 55 km/s, is comprised of a strongly ionized, non-LTE plasma front at ∼10 eV temperature followed by a weakly ionized LTE gas at higher density. The present findings lend some confidence to the broad-range equation of state and infrared opacity models for weakly ionized aluminum plasmas while demonstrating an approach to their future refinement, with potential application to astrophysical plasmas and other extreme processes.
Adiabatic process↗
Strong gravity extruding peaks in speed of sound profiles of massive neutron stars
The speed of sound squared (SSS) s 2 in massive neutron stars (NSs) characterizes not only the stiffness of supradense neutron-rich matter within but also equivalently properties of the curved geometry due to the strong-field gravity and matter-geometry coupling. A peaked density or radius profile of s 2 has been predicted for massive NSs using various NS equation of state (EOS) models. However, the nature, cause, location and size of the peak in s 2 profiles are still very EOS model dependent. In this work, we investigate systematically s 2 profiles in massive NSs in a new approach that is independent of the nuclear EOS model and without any presumption about the NS structure and/or composition. In terms of the small quantities (reduced radius, the energy density and pressure scaled by their central values), we perform double-element perturbative expansions in solving perturbatively the scaled Tolman-Oppenheimer-Volkoff (TOV) equations and analyzing s 2 profiles from the Newtonian limit to the general relativistic (GR) case. The GR term in the TOV equations plays a twofold role: it compresses NS matter and modifies the pressure/energy density ratio from small values in Newtonian stars showing no s 2 peak to large ones for massive NSs possessing a peak in their s 2 profiles, and eventually takes away the peak in extremely compact/massive NSs approaching the causality limit. In particular, the peaked behavior in s 2 is expected to emerge near the center of massive NSs like PSR J0740+6620, while a sharp phase transition is unlikely to occur there. Furthermore, these features revealed from our analyses are universal as they are intrinsic properties of the GR stellar structure equations independent of the still very uncertain EOS of supradense neutron-rich matter in NSs.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Reversible hydrogen storage in multilayer graphane: Lattice dynamics, compressibility, and heat capacity studies
Multilayer graphane (hydride of graphite) is a crystalline hydrocarbon of composition CH, which can be synthesized from graphite and molecular hydrogen at pressures above 2GPa [V.E. Antonov et al. Carbon 100 (2016) 465]. Using X-ray diffraction, this compound was tentatively identified as the “graphane II” phase of 3D-graphane predicted by ab initio calculations [X.-D. Wen et al. PNAS 108 (2011) 6833] and consisting of layers of 2D-graphane in the “chair” conformation. When heated in a vacuum, the compound does not form any intermediate hydrocarbons and reversibly decomposes back into graphite and hydrogen at 770–920 K. In the present work, almost single-phase samples of graphite hydride and deuteride were synthesized at 7.4 GPa and 870 K. Their investigation by inelastic neutron scattering supplemented by ab initio calculations gave spectra g(E) of the phonon density of states with a gap of about 15 meV at approx. 100 meV, which is a unique identifier for the chair form of graphane. The equation of state V(P) of the hydride was studied at room temperature and hydrogen pressures up to 53 GPa by synchrotron X-ray diffraction in a diamond anvil cell. Further, the graphane II phase did not react with the surrounding hydrogen and did not undergo any phase transformations upon the compression and after heating to 1500 K at 53GPa. The high thermal and pressure stability of this exotic phase makes it an important part of the C–H system. The obtained g(E) spectra of graphite hydride and deuteride were used to calculate temperature dependences of their heat capacity. Measurements of the heat capacity at temperatures 120–673 K confirmed the good accuracy of these calculations.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Theory of elastic shells in the reference state.
Nonlinear field equations of elastic shell in terms of reference state subjected to strain displacement
SHELL THEORY↗
Absolute Hugoniot measurements in low-density plastic foams on the NIKE shock compression platform
The propagation of Mbar-range shock waves in low-density foam materials is of interest to target design in all approaches to inertial confinement fusion, high energy density physics, and laboratory astrophysics. We report absolute Hugoniot measurements for CH 1.72 N 0.086 O 0.37 plastic foams with 73 and 94 mg/cm 3 densities in the 32–107 km/s shock velocity range. The experiments were performed on the shock compression platform developed on the NIKE KrF laser facility at the U.S. Naval Research Laboratory. NIKE's 4 or 8-ns long flat-top laser pulses drive steady shock waves into foam targets at the ablative pressures of 1–7 Mbar. The propagation of the ablation and shock fronts is tracked continuously in time using streaked side-on monochromatic x-ray imaging radiography. The straight x−t trajectories of the shock and ablation fronts in the recorded streak images confirmed their steadiness. The SESAME tabulated equation of state predictions generally agrees with our Hugoniot data within a 95% confidence band. The experimental uncertainty of the evaluated shock density compression ratios remains large, indicating the need for more experiments and improved theoretical understanding of the strong shock propagation mechanisms in dry foams.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Spin-dependent interactions and heavy-quark transport in the quark-gluon plasma
We extend a previously constructed T-matrix approach to the quark-gluon plasma (QGP) to include the effects of spin-dependent interactions between partons. Following earlier work within the relativistic quark model, the spin-dependent interactions figure as relativistic corrections to the Cornell potential. When applied to the vacuum spectroscopy of quarkonia, in particular their mass splittings in S- and P-wave states, the issue of the Lorentz structure of the confining potential arises. We confirm that a significant admixture of a vector interaction (to the previously assumed scalar interaction) improves the description of the experimental mass splittings. The temperature corrections to the in-medium potential are constrained by results from thermal lattice-QCD for the equation of state (EoS) and heavy-quark (HQ) free energy in a self-consistent set-up for heavy- and light-parton spectral functions in the QGP. We then deploy the refined in-medium heavy-light T-matrix to compute the charm-quark transport coefficients in the QGP. The vector component of the confining potential, through its relativistic corrections, enhances the friction coefficient for charm quarks in the QGP over previous calculations by tens of percent at low momenta and temperatures, and more at higher momenta. Our results are promising for improving the current phenomenology of open heavy-flavor observables at Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC).
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Quasi-Newton methods for parameter estimation in functional differential equations
A state-space approach to parameter estimation in linear functional differential equations is developed using the theory of linear evolution equations. A locally convergent quasi-Newton type algorithm is applied to distributed systems with particular emphasis on parameters that induce unbounded perturbations of the state. The algorithm is computationally implemented on several functional differential equations, including coefficient and delay estimation in linear delay-differential equations.
Brewer, Dennis W.↗
Application of Scalar Filtered Density Function to Turbulent Flows Under Supercritical Condition
Abstract The scalar filtered density function (FDF) methodology is extended and employed for large eddy simulation (LES) of turbulent flows under supercritical condition. To describe real fluid behavior, the extended methodology incorporates the generalized heat and mass diffusion models along with real fluid thermodynamic relations which are derived using the cubic Peng–Robinson equation of state. These models are implemented within the stochastic differential equations comprising the scalar FDF transport. Simulations are conducted of a temporally developing mixing layer under supercritical condition and the results are assessed by comparing with data generated by direct numerical simulation (DNS) of the same layer. The consistency of the proposed FDF methodology is assessed. The LES-FDF predictions are shown to agree favorably with the DNS data and exhibit several key features pertaining to supercritical turbulent flows.
Energy & Fuels↗
Real-gas effects 1: Simulation of ideal gas flow by cryogenic nitrogen and other selected gases
The thermodynamic properties of nitrogen gas do not thermodynamically approximate an ideal, diatomic gas at cryogenic temperatures. Choice of a suitable equation of state to model its behavior is discussed and the equation of Beattie and Bridgeman is selected as best meeting the needs for cryogenic wind tunnel use. The real gas behavior of nitrogen gas is compared to an ideal, diatomic gas for the following flow processes: isentropic expansion; normal shocks; boundary layers; and shock wave boundary layer interactions. The only differences in predicted pressure ratio between nitrogen and an ideal gas that may limit the minimum operating temperatures of transonic cryogenic wind tunnels seem to occur at total pressures approaching 9atmospheres and total temperatures 10 K below the corresponding saturation temperature, where the differences approach 1 percent for both isentropic expansions and normal shocks. Several alternative cryogenic test gases - air, helium, and hydrogen - are also analyzed. Differences in air from an ideal, diatomic gas are similar in magnitude to those of nitrogen. Differences for helium and hydrogen are over an order of magnitude greater than those for nitrogen or air. Helium and hydrogen do not approximate the compressible flow of an ideal, diatomic gas.
Hall, R. M.↗
Unstructured Finite Volume Computational Thermo-Fluid Dynamic Method for Multi-Disciplinary Analysis and Design Optimization
This paper describes a finite volume computational thermo-fluid dynamics method to solve for Navier-Stokes equations in conjunction with energy equation and thermodynamic equation of state in an unstructured coordinate system. The system of equations have been solved by a simultaneous Newton-Raphson method and compared with several benchmark solutions. Excellent agreements have been obtained in each case and the method has been found to be significantly faster than conventional Computational Fluid Dynamic(CFD) methods and therefore has the potential for implementation in Multi-Disciplinary analysis and design optimization in fluid and thermal systems. The paper also describes an algorithm of design optimization based on Newton-Raphson method which has been recently tested in a turbomachinery application.
Majumdar, Alok↗