Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite density”

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 721 records · Page 40

A microfluidic study of transient flow states in permeable media using fluorescent particle image velocimetry

Velocity fields in flow in permeable media are of great importance to many subsurface processes such as geologic storage of CO 2 , oil and gas extraction, and geothermal systems. Steady-state flow is characterized by velocity fields that do not change significantly over time. The flow field transitions to a new steady state once it experiences a disturbance such as a change in flow rate or in pressure gradient. This transition is often assumed to be instantaneous, which justifies the expression of constitutive relations as functions of instantaneous phase saturations. This work examines the evolution of velocity fields in a surrogate quasi-2D permeable medium using a microfluidic device, a microscopy system, and a high-speed camera. Tracer particles are injected into the medium along with DI water. The evolution of the velocity field is examined by tracing these particles in the captured images using the standard high-density particle image velocimetry algorithm founded on cross-correlation. The results suggest that the transition between steady states for an incompressible fluid takes a finite and non-negligible amount of time that is independent of the magnitude of the change in pressure gradient. The existence of transient states and the nature of the response during these states are readily interpreted by the principle of least action where flow gradually establishes an optimal configuration such that energy dissipation is 2 minimized. The findings provide evidence against the applicability of the assumption that flowing phases relax instantaneously to their steady states and, hence, against the accuracy of the classical multiphase extension of Darcy's law.

42 ENGINEERING↗

Models of Saturn's rings which satisfy the optical observations

A theoretical model of Saturn's rings is investigated which includes the shadowing effect and realistic anisotropic phase functions for the ring particles. The effects of multiple scattering and the finite size of the sun, including the penumbra, are rigorously included. The permissible range of the relevant parameters, including optical thickness, single scattering albedo, volume density, and phase function are investigated by comparing the theoretical results to observations of the ring brightness vs phase angle, wavelength, and elevation of the sun and earth. Anisotropic scattering by the ring particles is necessary in order to match the observations. The color dependence of the opposition effect is interpreted in terms of the albedo spectrum of the ring particles.

Kawata, Y.↗

Unsteady design-point flow phenomena in transonic compressors

High-frequency response probes which had previously been used exclusively in the MIT Blowndown Facility were successfully employed in two conventional steady state axial flow compressor facilities to investigate the unsteady flowfields of highly loaded transonic compressors at design point operation. Laser anemometry measurements taken simultaneously with the high response data were also analyzed. The time averaged high response data of static and total pressure agreed quite well with the conventional steady state instrumentation except for flow angle which showed a large spread in values at all radii regardless of the type of instrumentation used. In addition, the time resolved measurements confirmed earlier test results obtained in the MIT Blowdown Facility for the same compressor. The results of these tests have further revealed that the flowfields of highly loaded transonic compressors are heavily influenced by unsteady flow phenomena. The high response measurements exhibited large variations in the blade to blade flow and in the blade passage flow. The observed unsteadiness in the blade wakes is explained in terms of the rotor blades' shed vorticity in periodic vortex streets. The wakes were modeled as two-dimensional vortex streets with finite size cores. The model fit the data quite well as it was able to reproduce the average wake shape and bi-modal probability density distributions seen in the laser anemometry data. The presence of vortex streets in the blade wakes also explains the large blade to blade fluctuations seen by the high response probes which is simply due to the intermittent sampling of the vortex street as it is swept past a stationary probe.

Gertz, J. B.↗

Deformational models of rifting and folding on Venus

Features of presumed tectonic origin on Venus are reviewed, and lithospheric strength envelopes are derived based on laboratory measurements of the deformational properties of crustal and subcrustal rocks, extrapolated to conditions appropriate to Venus. Models for rifting and folding are developed that use this lithospheric structure and take into account both brittle and ductile yielding as well as finite elastic strength. For both rifting and folding, structures with characteristic widths and spacings are predicted whose size depends on the thickness of the lithosphere, density contrast, and elastic properties of the layer. Finally, the model predictions are compared with the widths and spacings of observed tectonic features, and it is concluded that they are consistent with a relatively strong mantle layer separated from a thin brittle surface layer by a ductile lower crust. These results allow constraints to be placed on the crustal thickness and thermal gradient on Venus.

Banerdt, W. Bruce↗

Comparative dynamics of the ionospheres of Venus and Mars at large solar zenith angles

Simple three-ion spectral models of the dynamics and ion chemistry in the upper ionospheres of Venus and Mars have been constructed by including two minor ions, H(+) and O2(+) for Venus and O(+) and CO2(+) for Mars, along with the major ion, O(+) for Venus and O2(+) for Mars. Horizontal flow velocities and ion densities have been calculated in each case and compared with available measured data. For Venus, the present calculations yield results in agreement with earlier finite-difference model studies and with experimental observations. Calculations for Mars are in agreement with the limited observational data obtained by the Viking landers.

Singhal, R. P.↗

Comparison of ALE and SPH Simulations of Vertical Drop Tests of a Composite Fuselage Section into Water

Simulation of multi-terrain impact has been identified as an important research area for improved prediction of rotorcraft crashworthiness within the NASA Subsonic Rotary Wing Aeronautics Program on Rotorcraft Crashworthiness. As part of this effort, two vertical drop tests were conducted of a 5-ft-diameter composite fuselage section into water. For the first test, the fuselage section was impacted in a baseline configuration without energy absorbers. For the second test, the fuselage section was retrofitted with a composite honeycomb energy absorber. Both tests were conducted at a nominal velocity of 25-ft/s. A detailed finite element model was developed to represent each test article and water impact was simulated using both Arbitrary Lagrangian Eulerian (ALE) and Smooth Particle Hydrodynamics (SPH) approaches in LS-DYNA, a nonlinear, explicit transient dynamic finite element code. Analytical predictions were correlated with experimental data for both test configurations. In addition, studies were performed to evaluate the influence of mesh density on test-analysis correlation.

Jackson, Karen E.↗

Confirmation of the PPLB Derivative Discontinuity: Exact Chemical Potential at Finite Temperatures of a Model System

The landmark 1982 paper of Perdew, Parr, Levy, and Balduz (often called PPLB) laid the foundation for our modern understanding of the role of the derivative discontinuity in density functional theory, which drives much development to account for its effects. A simple model for the chemical potential at vanishing temperature played a crucial role in their argument. We investigate the validity of this model in the simplest non-trivial system to which it can be applied and which can be easily solved exactly, the Hubbard dimer. We find exact agreement in the crucial zero-temperature limit, and show the model remains accurate for a significant range of temperatures. We identify how this range depends on the strength of correlations. We extend the model to approximate free energies accounting for the derivative discontinuity, a feature missing in standard semilocal approximations. Here, we provide a correction to this approximation to yield even more accurate free energies. We discuss the relevance of these results for warm dense matter.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Resolving black-hole microstructure with new momentum carriers

All known horizonless black-hole microstate geometries correspond to brane sources that acquire a finite size, and hence break the spherical symmetry of the black hole. We construct, for the first time, solutions with zero horizon area that have the same charges as a three-charge F1-NS5-P Type-IIA black hole and preserve this spherical symmetry. The momentum of these solutions is carried by longitudinal D0-D4 density fluctuations inside the NS5-branes. We argue that these solutions should be interpreted as the long-throat limit of a family of smooth, horizonless microstate geometries, called superstrata, where such geometries degenerate. The existence of these geometries indicates that a finite-size horizon does not appear even in the singular corners of the moduli space of three-charge microstate geometries.

79 ASTRONOMY AND ASTROPHYSICS↗

FDDO and DSMC analyses of rarefied gas flow through 2D nozzles

Two different approaches, the finite-difference method coupled with the discrete-ordinate method (FDDO), and the direct-simulation Monte Carlo (DSMC) method, are used in the analysis of the flow of a rarefied gas expanding through a two-dimensional nozzle and into a surrounding low-density environment. In the FDDO analysis, by employing the discrete-ordinate method, the Boltzmann equation simplified by a model collision integral is transformed to a set of partial differential equations which are continuous in physical space but are point functions in molecular velocity space. The set of partial differential equations are solved by means of a finite-difference approximation. In the DSMC analysis, the variable hard sphere model is used as a molecular model and the no time counter method is employed as a collision sampling technique. The results of both the FDDO and the DSMC methods show good agreement. The FDDO method requires less computational effort than the DSMC method by factors of 10 to 40 in CPU time, depending on the degree of rarefaction.

Chung, Chan-Hong↗

Symmetry-energy expansion with strange dense matter

The quantum chromodynamics (QCD) phase diagram at large densities and low temperatures can be probed using both neutron stars and low-energy heavy-ion collisions. Heavy-ion collisions are nearly isospin-symmetric systems, whereas neutron stars are highly isospin asymmetric since they are neutron rich. The symmetry-energy expansion is used to connect these regimes across isospin asymmetry. However, the current symmetry-energy expansion does not account for strange particles. In this work, we include finite strangeness by redefining the isospin-asymmetry parameter and the symmetry-energy expansion in a way that is consistent with QCD SU(3) flavor symmetry. Furthermore, our new symmetry energy works well beyond typical neutron star central densities and admits a skewness term in the presence of strangeness for the case of weak equilibrium.

Equations of state of nuclear matter↗

Large Kohn Anomaly and Phonon Collapse Induced by Charge Density Wave in UPt 2⁢ Si 2

Using high-energy-resolution inelastic x-ray scattering, we observe anomalous softening and damping of the transverse acoustic phonon in UPt 2 ⁢Si 2 as the system is cooled toward the charge density wave (CDW) transition temperature 𝑇 CDW . The phonon exhibits a marked Kohn-type anomaly around the CDW wave vector 𝑄 CDW and becomes overdamped within a finite momentum range already well above 𝑇 CDW . The dispersion anomaly is consistent with potential Fermi surface nesting, which together with the extended phonon collapse indicates strong electron-phonon coupling. The transition temperature estimated from the phonon softening is markedly lower than 𝑇 CDW , consistent with a primarily electronic instability rather than a phonon-driven transition. In conclusion, our results establish UPt 2⁢ Si 2 as a prime example of a strongly correlated electron CDW system with exceptionally strong electron-phonon coupling driving phonon softening and collapse.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Massively parallel axisymmetric fluid model for streamer discharges

A highly parallelizable fluid plasma simulation tool based upon the first-order drift-diffusion equations is discussed. Atmospheric pressure plasmas have densities and gradients that require small element sizes in order to accurately simulate the plasm resulting in computational meshes on the order of millions to tens of millions of elements for realistic size plasma reactors. To enable simulations of this nature, parallel computing is required and must be optimized for the particular problem. Here, a finite-volume, electrostatic drift-diffusion implementation for low-temperature plasma is discussed. The implementation is built upon the Message Passing Interface (MPI) library in C++ using Object Oriented Programming. The underlying numerical method is outlined in detail and benchmarked against simple streamer formation from other streamer codes. Electron densities, electric field, and propagation speeds are compared with the reference case and show good agreement. Convergence studies are also performed showing a minimal space step of approximately 4 μm required to reduce relative error to below 1% during early streamer simulation times and even finer space steps are required for longer times. Additionally, strong and weak scaling of the implementation are studied and demonstrate the excellent performance behavior of the implementation up to 100 million elements on 1024 processors. Lastly, different advection schemes are compared for the simple streamer problem to analyze the influence of numerical diffusion on the resulting quantities of interest.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Large Area Detection of Microstructural Defects with Multi-Mode Ultrasonic Signals

Cyclic loading or other stresses can lead to development of cracks and crack growth in mechanical structures, leading to eventual failure. While ultrasound imaging can be used for non-destructive testing of such structures, conventional ultrasound techniques are often limited by crack size, density, and areal coverage. An effective characterization of real-world, large-area structures is required at an early damage stage to prevent catastrophic failure and predict remaining life. In this study, a new nonlinear ultrasonic testing (NUT) method is proposed for large-area monitoring of practical structures with arbitrary complexity by using multiple-mode guided-wave ultrasonic signals. The proposed guided-wave NUT technique requires single-element transducers, simple electronics, and a mixed time-frequency domain signal processing. As a proof-of-concept demonstration, numerical simulations and experiments are performed on an A36 carbon steel beam assembly with previously formed microstructural defects that cause nonlinearities in ultrasonic response. The quadratic dependence of the nonlinear wave excitation on the input ultrasonic signal amplitude is shown by numerical simulations, and such a nonlinear ultrasonic response is experimentally observed in the zone with a high density of microstructural defects.

36 MATERIALS SCIENCE↗

Comparison of Several Numerical Methods for Simulation of Compressible Shear Layers

An investigation is conducted on several numerical schemes for use in the computation of two-dimensional, spatially evolving, laminar variable-density compressible shear layers. Schemes with various temporal accuracies and arbitrary spatial accuracy for both inviscid and viscous terms are presented and analyzed. All integration schemes use explicit or compact finite-difference derivative operators. Three classes of schemes are considered: an extension of MacCormack's original second-order temporally accurate method, a new third-order variant of the schemes proposed by Rusanov and by Kutier, Lomax, and Warming (RKLW), and third- and fourth-order Runge-Kutta schemes. In each scheme, stability and formal accuracy are considered for the interior operators on the convection-diffusion equation U(sub t) + aU(sub x) = alpha U(sub xx). Accuracy is also verified on the nonlinear problem, U(sub t) + F(sub x) = 0. Numerical treatments of various orders of accuracy are chosen and evaluated for asymptotic stability. Formally accurate boundary conditions are derived for several sixth- and eighth-order central-difference schemes. Damping of high wave-number data is accomplished with explicit filters of arbitrary order. Several schemes are used to compute variable-density compressible shear layers, where regions of large gradients exist.

Kennedy, Christopher A.↗

Circuit Quantum Electrodynamics in Hyperbolic Space: From Photon Bound States to Frustrated Spin Models

Circuit quantum electrodynamics is one of the most promising platforms for efficient quantum simulation and computation. In recent groundbreaking experiments, the immense flexibility of superconducting microwave resonators was utilized to realize hyperbolic lattices that emulate quantum physics in negatively curved space. Here we investigate experimentally feasible settings in which a few superconducting qubits are coupled to a bath of photons evolving on the hyperbolic lattice. We compare our numerical results for finite lattices with analytical results for continuous hyperbolic space on the Poincaré disk. We find good agreement between the two descriptions in the long-wavelength regime. We show that photon-qubit bound states have a curvature-limited size. We propose to use a qubit as a local probe of the hyperbolic bath, for example, by measuring the relaxation dynamics of the qubit. We find that, although the boundary effects strongly impact the photonic density of states, the spectral density is well described by the continuum theory. We show that interactions between qubits are mediated by photons propagating along geodesics. We demonstrate that the photonic bath can give rise to geometrically frustrated hyperbolic quantum spin models with finite-range or exponentially decaying interaction.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A fully-integrated lattice Boltzmann method for fluid–structure interaction

Here we present a fully-integrated lattice Boltzmann (LB) method for fluid–structure interaction (FSI) simulations that efficiently models deformable solids in complex suspensions and active systems. Our Eulerian method (LBRMT) couples finite-strain solids to the LB fluid on the same fixed computational grid with the reference map technique (RMT). An integral part of the LBRMT is a new LB boundary condition for moving deformable interfaces across different densities. With this fully Eulerian solid–fluid coupling, the LBRMT is well-suited for parallelization and simulating multi-body contact without remeshing or extra meshes. We validate its accuracy via a benchmark of a deformable solid in a lid-driven cavity, then showcase its versatility through examples of soft solids rotating and settling. The LBRMT achieves a spatial convergence rate between first-order and second-order for FSI simulations and is designed for low to intermediate Reynolds number flows with finite inertia at small Mach numbers. With simulations of complex suspensions mixing, we highlight the potential of the LBRMT for studying collective behavior in soft matter and biofluid dynamics.

97 MATHEMATICS AND COMPUTING↗

Parametric instabilities of circularly polarized Alfven waves in high-beta plasmas

CGL relations including the effect of finite ion Larmor radius are used to consider a class of parametric instabilities of finite-amplitude, circularly polarized Alfven waves in high-beta plasmas. The disperison relation governing the instabilities is a sixth-order polynomial which is solved numerically. There are two types of instabilities: a modulational instability at k is less than k(0) and a relatively weak and narrow bandwidth instability at k is less than approximately k(0), where k and k(0) are the wavenumbers of the unstable density fluctuation and the 'pump' wave, respectively. It is shown that these instabilities can occur for left-handed pump waves and that the modulational instability is unstable over a very broad band in k with a maximum growth rate at finite k is not equal to 0.

Hamabata, Hiromitsu↗

Power spectral density estimation by spline smoothing in the frequency domain

An approach, based on a global averaging procedure, is presented for estimating the power spectrum of a second order stationary zero-mean ergodic stochastic process from a finite length record. This estimate is derived by smoothing, with a cubic smoothing spline, the naive estimate of the spectrum obtained by applying FFT techniques to the raw data. By means of digital computer simulated results, a comparison is made between the features of the present approach and those of more classical techniques of spectral estimation.

Defigueiredo, R. J. P.↗