Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density model”

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,459 records · Page 81

BHR equations re-derived with immiscible particle effects

Compressible and variable density turbulent flows with dispersed phase effects are found in many applications ranging from combustion to cloud formation. These types of flows are among the most challenging to simulate. While the exact equations governing a system of particles and fluid are known, computational resources limit the scale and detail that can be simulated in this type of problem. Therefore, a common method is to simulate averaged versions of the flow equations, which still capture salient physics and is relatively less computationally expensive. Besnard et al. [1] developed such a model for variable density miscible turbulence, where ensemble-averaging was applied to the flow equations to yield a set of filtered equations. Besnard et al. further derived transport equations for the Reynolds stresses, the turbulent mass flux, and the density-specific volume covariance, to help close the filtered momentum and continuity equations. We re-derive the exact BHR closure equations which include integral terms owing to immiscible effects. Physical interpretations of the additional terms are proposed along with simple models. The goal of this work is to extend the BHR model to allow for the simulation of turbulent flows where an immiscible dispersed phase is non-trivially coupled with the carrier phase.

20 FOSSIL-FUELED POWER PLANTS↗

Effects of angular-momentum-changing collisions on dielectronic satellite spectra

A generalization of the low-density corona-model theory of the intensities of satellite lines of the resonance lines emitted by multiply charged helium-like ions in the X-ray spectra of high-temperature plasmas is developed in order to take account of electron collisional transitions between close autoionizing levels. A formalism for determining the population densities of the autoionizing levels and the satellite line intensities is presented, and the quantum mechanical calculation of autoionization, radiative decay and electron-induced collisional transitions is discussed. Generalized emission rate coefficients corresponding to dielectronic recombination and inner-shell-electron collisional excitation are tabulated as functions of electron density and temperature. The procedure is then used to illustrate the effects of angular-momentum-changing electron collisions on the intensities of the satellites of the 1s2p 1P to 1s2 1S resonance line of the ions Al XII and Si XIII. It is concluded that the collisional transitions have their most significant effect on the intensities of satellite transitions from metastable autoionizing states, which may be useful in determining plasma densities.

Jacobs, V. L.↗

Cosmology with decaying cosmological constant—exact solutions and model testing

We study dynamics of Λ(t) cosmological models which are a natural generalization of the standard cosmological model (the ΛCDM model). We consider a class of models: the ones with a prescribed form of Λ(t)=Λ{sub bare}+α{sup 2}/t{sup 2}. This type of a Λ(t) parametrization is motivated by different cosmological approaches. We interpret the model with running Lambda (Λ(t)) as a special model of an interacting cosmology with the interaction term −dΛ(t)/dt in which energy transfer is between dark matter and dark energy sectors. For the Λ(t) cosmology with a prescribed form of Λ(t) we have found the exact solution in the form of Bessel functions. Our model shows that fractional density of dark energy Ω{sub e} is constant and close to zero during the early evolution of the universe. We have also constrained the model parameters for this class of models using the astronomical data such as SNIa data, BAO, CMB, measurements of H(z) and the Alcock-Paczyński test. In this context we formulate a simple criterion of variability of Λ with respect to t in terms of variability of the jerk or sign of estimator (1−Ω{sub m},0−Ω{sub Λ,0}). The case study of our model enable us to find an upper limit α{sup 2} < 0.012 (2σ C.L.) describing the variation from the cosmological constant while the LCDM model seems to be consistent with various data.

79 ASTRONOMY AND ASTROPHYSICS↗

Analytical model of foreshock ion interaction with a discontinuity: A statistical study

When a solar wind discontinuity interacts with foreshock ions, foreshock transients such as hot flow anomalies and foreshock bubbles can form. These create significant dynamic pressure perturbations disturbing the bow shock, magnetopause, and magnetosphere-ionosphere system. However, presently these phenomena are not predictable. In the accompanying paper, we derived analytical equations of foreshock ion partial gyration around a discontinuity and the resultant current density. In this study, we utilize the derived current density strength to model the energy conversion from the foreshock ions, which drives the outward motion or expansion of the solar wind plasma away from the discontinuity. We show that the model expansion speeds match those from local hybrid simulations for varying foreshock ion parameters. Using MMS, we conduct a statistical study showing that the model expansion speeds are moderately correlated with the magnetic field strength variations and the dynamic pressure decreases around discontinuities with correlation coefficients larger than 0.5. We use conjunctions between ARTEMIS and MMS to show that the model expansion speeds are typically large for those already-formed foreshock transients. Our results show that our model can be reasonably successful in predicting significant dynamic pressure disturbances caused by foreshock ion-discontinuity interactions. We discuss ways to improve the model in the future.

Terry Z. Liu↗

Interior Noise Predictions in the Preliminary Design of the Large Civil Tiltrotor (LCTR2)

A prediction scheme was established to compute sound pressure levels in the interior of a simplified cabin model of the second generation Large Civil Tiltrotor (LCTR2) during cruise conditions, while being excited by turbulent boundary layer flow over the fuselage, or by tiltrotor blade loading and thickness noise. Finite element models of the cabin structure, interior acoustic space, and acoustically absorbent (poro-elastic) materials in the fuselage were generated and combined into a coupled structural-acoustic model. Fluctuating power spectral densities were computed according to the Efimtsov turbulent boundary layer excitation model. Noise associated with the tiltrotor blades was predicted in the time domain as fluctuating surface pressures and converted to power spectral densities at the fuselage skin finite element nodes. A hybrid finite element (FE) approach was used to compute the low frequency acoustic cabin response over the frequency range 6-141 Hz with a 1 Hz bandwidth, and the Statistical Energy Analysis (SEA) approach was used to predict the interior noise for the 125-8000 Hz one-third octave bands.

Grosveld, Ferdinand W.↗

Model for the electrical conductivity in dense plasma mixtures

A new density functional theory, average atom based model for the electrical conductivity of dense plasmas with a mixture of ion species, containing no adjustable parameters, is presented herein. The model takes the temperature, mass density and relative abundances of the species as input. It takes into account partial ionization, ionic structure, and core-valence orthogonality, and uses quantum mechanical calculations of cross sections. Comparison to an existing high fidelity but computationally expensive method reveals good agreement. The new model is computationally efficient and can reach high temperatures. A new mixing rule is also presented that gives reasonably accurate conductivities for high temperature plasma mixtures.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Properties of planetary fluids at high pressures and temperatures

Observational data obtained by the Voyager space probes to the giant planets Jupiter, Saturn, Uranus, and Neptune have provided valuable information, which is used to refine the picture of the nature of the interiors of these planets. Major results from the Voyager missions include observations of substantial magnetic fields and improved models of internal density distributions. The goal is to obtain equations of state and electrical conductivity data for planetary gases (H2 and He) and the ices (H2O, CH4, and NH3, and their mixtures), which are considered to be the major constituents of the giant planets. These data are needed to test theoretical data bases used to construct models of the chemical composition of planetary interiors, models which are consistent with observables such as mass, diameter, gravitational moments, rotation rate, and magnetic field. The 100 GPa (1 Mbar) pressures and several 1000 K temperatures in the giant planets can be achieved in the lab by the shock compression of liquid specimens. Results are briefly examined.

Nellis, W. J.↗

Effect of Particle Spreading Dynamics on Powder Bed Quality in Metal Additive Manufacturing

Powder spreading precedes creation of every new layer in powder bed additive manufacturing (AM). The powder spreading process can lead to powder layer defects such as porosity, poor surface roughness and particle segregation. Therefore, the creation of homogeneous layers is the first task for optimal part printing. Discrete element methods (DEM) powder spreading simulations are typically limited to a single layer and/or small number of particles. Therefore, results from such model configurations may not be generalized to multiple layer processes. In this study, a computationally efficient multi-layer powder spreading DEM simulation model is proposed. The model is calibrated experimentally using static Angle of Repose measurements. The adhesion model parameter, cohesive energy density is related to adhesive surface energy and strain energy release rate parameters. The model results show that interaction between particle and the powder spreading rake leads to noticeable variation in packing density, surface roughness, dynamic angle of repose (AOR), particle size distribution, and particle segregation. Finally, the powder model is experimentally validated using a recoater spreading rig to measure the dynamic AOR at spreading speeds consistent with recoating speeds and layer heights used in AM processes.

36 MATERIALS SCIENCE↗

Ionoshperic effects of the March 1990 magnetic storm: Comparison of theory and measurement

This paper presents a comparison of the measured and modeled ionospheric response to magnetic storms at Millstone Hill and Arecibo during March 16-23, 1990. Magnetic activity was low until midday Universal Time (UT) on day 18 when Kp reached 6, days 19 and 20 were quiet, but a large storm occurred around midnight UT on day 20 (Kp=7) ad it was a moderately disturbed (Kp=4) for the remainder of the study period. At Millstone Hill, the daytime peak density (Nm F2) showed only a modest 30% decrease in response to the first storm and recovered to prestorm values before the onset of the second storm. The model reproduces the daytime peak electron density well for this period. However, the severe storm on March 20 caused a factor of 4 deplection in electron density, while the model densities were not greatly affected. the inclusion of vibrationally excited nitrogen N2(sup *) in the model was unable to account for the observed large electron density depletions afterward March 20. The storm did not appear to affect the overall magnitude of the electron density at Arecibo very much, but did cause unusual wavelike structure in the peak density and peak height following the storm. The model reproduces the daytime Nm F2 very well for Arecibo, but after sunset the model densities decay too rapidly. This study indicates that successful modeling of severe ionospheric stroms will require better definition of the storm time inputs, especially of the neutral atmosphere.

Richards, P. G.↗

Density and pressure variability in the mesosphere and thermosphere

In an effort to isolate the essential physics of the mesosphere and the thermosphere, a steady one-dimensional density and pressure model has been developed in support of related NASA activities, i.e., projects such as the AOTV and the Space Station. The model incorporates a zeroth order basic state including both the three-dimensional wind field and its associated shear structure, etc. A first order wave field is also incorporated in period bands ranging from about one second to one day. Both basic state and perturbation quantities satsify the combined constraints of mass, linear momentum and energy conservation on the midlatitude beta plane. A numerical (iterative) technique is used to solve for the vertical wind which is coupled to the density and pressure fields. The temperature structure from 1 to 1000 km and the lower boundary conditions are specified using the U.S. Standard Atmosphere 1976. Vertical winds are initialized at the top of the Planetary Boundary Layer using Ekman pumping values over flat terrain. The model also allows for the generation of waves during the geostrophic adjustment process and incorporates wave nonlinearity effects.

Davis, T. M.↗

Photoionization Models for the Inner Gaseous Disks of Herbig Be Stars: Evidence against Magnetospheric Accretion?

We investigate the physical properties of the inner gaseous disks of three hot Herbig B2e stars, HD 76534, HD 114981, and HD 216629, by modeling CFHT-ESPaDOns spectra using non-LTE radiative transfer codes. We assume that the emission lines are produced in a circumstellar disk heated solely by photospheric radiation from the central star in order to test whether the optical and near-infrared emission lines can be reproduced without invoking magnetospheric accretion. The inner gaseous disk density was assumed to follow a simple power-law in the equatorial plane, and we searched for models that could reproduce observed lines of H i (H α and H β ), He i, Ca ii, and Fe ii. For the three stars, good matches were found for all emission line profiles individually; however, no density model based on a single power-law was able to reproduce all of the observed emission lines. Among the single power-law models, the one with the gas density varying as ∼10{sup −10}( R {sub *}/ R ){sup 3} g cm{sup −3} in the equatorial plane of a 25 R {sub *} (0.78 au) disk did the best overall job of representing the optical emission lines of the three stars. This model implies a mass for the H α -emitting portion of the inner gaseous disk of ∼10{sup −9} M {sub *}. We conclude that the optical emission line spectra of these HBe stars can be qualitatively reproduced by a ≈1 au, geometrically thin, circumstellar disk of negligible mass compared to the central star in Keplerian rotation and radiative equilibrium.

79 ASTRONOMY AND ASTROPHYSICS↗

Investigating the Role of Copper Oxide in Electrochemical CO 2 Reduction in Real Time

Copper oxides have been of considerable interest as electrocatalysts for CO 2 reduction (CO2R) in aqueous electrolytes. Yet, their role as an active catalyst in reducing the required overpotential and improving the selectivity of reaction compared with that of polycrystalline copper remains controversial. In this work, we present the use of selected-ion flow tube mass spectrometry, in concert with chronopotentiometry, in situ Raman spectroscopy, and computational modeling, to investigate CO2R on C u 2O nanoneedles, Cu 2 O nanocrystals, and Cu 2 O nanoparticles. We show experimentally that the selective formation of gaseous C2 products (i.e., ethylene) in CO 2 R is preceded by the reduction of the copper oxide (Cu 2 OR) surface to metallic copper. On the basis of density functional theory modeling, CO2R products are not formed as long as Cu 2 O is present at the surface because Cu 2 OR is kinetically and energetically more favorable than CO2R.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Analytical Modeling of a Double-Sided Flux Concentrating E-Core Transverse Flux Machine with Pole Windings

In this paper, a nonlinear analytical model based on the Magnetic Equivalent Circuit (MEC) method is developed for a double-sided E-Core Transverse Flux Machine (TFM). The proposed TFM has a cylindrical rotor, sandwiched between E-core stators on both sides. Ferrite magnets are used in the rotor with flux concentrating design to attain high airgap flux density, better magnet utilization, and higher torque density. The MEC model was developed using a series-parallel combination of flux tubes to estimate the reluctance network for different parts of the machine including air gaps, permanent magnets, and the stator and rotor ferromagnetic materials, in a two-dimensional (2-D) frame. An iterative Gauss-Siedel method is integrated with the MEC model to capture the effects of magnetic saturation. A single phase, 1 kW, 400 rpm E-Core TFM is analytically modeled and its results for flux linkage, no-load EMF, and generated torque, are verified with Finite Element Analysis (FEA). The analytical model significantly reduces the computation time while estimating results with less than 10 percent error.

30 DIRECT ENERGY CONVERSION↗

An Analytical Tool to Evaluate Defect Thermodynamics of (La,Ba)Fe1-xMxO3-δ Perovskites for Solid-Oxide Cell Applications

A modeling tool of the defect thermodynamics of (La,Ba)Fe1-xMxO3-δ perovskites which includes energetic information about oxygen vacancy formation, hydration, hydride formation, and charge disproportionation reactions has been developed1. This tool incorporates defect energies and entropies expressed as sixth order polynomial functions to allow refinements of the defect reaction equilibrium constants in the thermodynamic analysis. Calculation of (La,Ba)Fe1-xMxO3-δ Brouwer diagrams as a function of pO2/pH2O and pH2/H2O in a range of temperatures of interest is facilitated by this modeling tool. The results obtained can provide direct guidance how the electronic and ionic defect concentrations of the triple conducting perovskite materials can be used to optimize performance of solid oxide cells for energy applications. The impact of magnetic and electronic structures of the perovskites on the defect reaction energies and entropies as obtained from density function theory modeling and the role played by hydride defect species will also be discussed.

Lee, Yueh-Lin↗

Data-driven model for cross ventilation potential in high-density cities based on coupled CFD simulation and machine learning

Effective urban ventilation through decent urban planning and building design can alleviate the deterioration of the urban built environment. Furthermore, natural ventilation requirements and guidelines in current building codes and standards are either qualitative or quantitative but subject to an absolute indoor airspeed threshold without considering the outdoor wind environment. To fill this gap, this paper develops an urban-scale coupled indoor and outdoor computational fluid dynamics (CFD) model and defines a novel ventilation index to assess natural ventilation potential. The index considers wind environments of both indoor and outdoor spaces.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Variational study of fermionic and bosonic systems with non-Gaussian states: Theory and applications

We present a new variational method for investigating the ground state and out of equilibrium dynamics of quantum many-body bosonic and fermionic systems. Our approach is based on constructing variational wavefunctions which extend Gaussian states by including generalized canonical transformations between the fields. The key advantage of such states compared to simple Gaussian states is presence of non-factorizable correlations and the possibility of describing states with strong entanglement between particles. In contrast to the commonly used canonical transformations, such as the polaron or Lang–Firsov transformations, we allow parameters of the transformations to be time dependent, which extends their regions of applicability. We derive equations of motion for the parameters characterizing the states both in real and imaginary time using the differential structure of the variational manifold. The ground state can be found by following the imaginary time evolution until it converges to a steady state. Collective excitations in the system can be obtained by linearizing the real-time equations of motion in the vicinity of the imaginary time steady-state solution. Our formalism allows us not only to determine the energy spectrum of quasiparticles and their lifetime, but to obtain the complete spectral functions and to explore far out of equilibrium dynamics such as coherent evolution following a quantum quench. We illustrate and benchmark this framework with several examples: a single polaron in the Holstein and Su–Schrieffer–Heeger models, non-equilibrium dynamics in the spin-boson and Kondo models, the superconducting to charge density wave phase transitions in the Holstein model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Landmark-Warped Emulators for Models with Misaligned Functional Response

Many computer models output functional data, and in some cases, these functional data have similar, but misaligned, shape characteristics. In this paper, we introduce a general approach for building emulators for computer models that output misaligned functional data when key values in the functional response (landmarks) can be easily identified. This approach has two main parts: modeling the aligned (using the landmarks) functional data, and modeling the functions that map the misaligned data to the aligned space (warping functions). As the warping functions are required to be monotonic, we give special attention to modeling monotonic functional response data. We discuss how our approach can be easily applied for a variety of typical emulators, such as Gaussian processes, Bayesian multivariate adaptive regression splines, and Bayesian additive regression trees, and how sensitivity analysis can be performed. We demonstrate our approach by building emulators for two applications: (1) a high-energy-density physics computer model used to simulate inertial confinement fusion ignition experiments, where model outputs are highly misaligned, and (2) a multiphysics continuum hydrocode used to simulate high-velocity impact experiments, where model outputs are only slightly misaligned. In case (1) traditional methods cannot be applied, while in (2) they can be applied, but the proposed method performs significantly better.

97 MATHEMATICS AND COMPUTING↗