Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “PHASE VELOCITY”

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 397 records · Page 22

Advanced Multi-Dimensional Capacitance Sensors Based Subsea Multiphase Mass Flow Meter to Measure and Monitor Offshore Enhanced Oil

This innovation is based on advanced multi-dimensional extensions of ECVT sensors that involve ECVT, Displacement Current Phase Tomography (DCPT), Maxwell-Wagner-Sillars polarization effect (MWS), and velocimetry which exploit the variation of electric properties between the oil, water, and gas coming out of the well. Capacitance sensors are embedded on the inside of a pipe spool and placed in-line with the subsea oil line piping on the extraction end of the well. The difference in dielectric, dielectric loss, surface polarization, and velocity of each phase are used in the multi-dimensional algorithm to measure the volume fraction, distribution, velocity, mass flow rate, and flow regime of the mixture moving through the pipe.

02 PETROLEUM↗

Phase transitions in a reaction-diffusion model on a line with boundaries

A one-dimensional model on a line of length L is investigated, which involves particle diffusion as well as single particle annihilation. There are also creation and annihilation at the boundaries. The static and dynamical behaviors of the system are studied. It is seen that the system could exhibit a dynamical phase transition. For small drift velocities, the relaxation time does not depend on the absorption rates at the boundaries. This is the fast phase. For large velocities, the smaller of the absorption rates at boundaries enter the relaxation rate and makes it longer. This is the slow phase. Finally, the effect of a random particle creation in the bulk is also investigated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Forward modeling of Doppler velocity interferometer system for improved shockwave measurements

We report that when near-instantaneous shocks are recorded by using a Doppler velocity interferometer velocity interferometer system for any reflector (VISAR), they typically exceed the detector’s ability to react and “skipped fringes” result, where its visibility briefly reduces. Traditionally, replacing skipped fringes required guesswork in analysis, which increased arrival time errors. Second, the use of long but velocity-sensitive interferometer delays with fast detectors which can resolve the delay has traditionally been avoided because of the fear of confusing the arrival time signal. However, shorter delays produce smaller fringe phase shifts per velocity and, thus, can decrease velocity precision. We realize that while some loss of fringe information occurs at shock events, this is often just a partial loss and the residual fringe information can still hold valuable information. We describe a forward model (FM) of the interferometer action and detector blurring that assists with VISAR fringe analysis at skipping events: (1) more precise shock arrival times, (2) arrival time precision not limited by long delays, and (3) improved ghost subtraction, which improves accuracy over a broad time region. Lastly, we demonstrate the utility of a FM on National Ignition Facility and Laboratory for Laser Energetics Omega shots and discuss data process methods to increase the high time resolution of VISAR systems despite their limitations and complexities.

47 OTHER INSTRUMENTATION↗

Radial and poloidal correlation reflectometry on Experimental Advanced Superconducting Tokamak

An X-mode polarized V band (50 GHz–75 GHz) radial and poloidal correlation reflectometry is designed and installed on Experimental Advanced Superconducting Tokamak (EAST). Two frequency synthesizers (12 GHz–19 GHz) are used as sources. Signals from the sources are up-converted to V band using active quadruplers and then coupled together for launching through one single pyramidal antenna. Two poloidally separated antennae are installed to receive the reflected waves from plasma. This reflectometry system can be used for radial and poloidal correlation measurement of the electron density fluctuation. In ohmically heated plasma, the radial correlation length is about 1.5 cm measured by the system. The poloidal correlation analysis provides a means to estimate the fluctuation velocity perpendicular to the main magnetic field. In the present paper, the distance between two poloidal probing points is calculated with ray-tracing code and the propagation time is deduced from cross-phase spectrum. Fluctuation velocity perpendicular to the main magnetic field in the core of ohmically heated plasma is about from −1 km/s to −3 km/s.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Impact of hohlraum cooling on ignition metrics for inertial fusion implosions

This paper extends the evaluation of ignition metrics to include the impact of hohlraum cooling before peak implosion velocity in radiation driven implosions. First, we provide an extension of the results for the key hot spot stagnation quantities from the 2018 paper. The modified analytic expressions presented here match the Hydra results for these National Ignition Facility scale implosions both with and without hohlraum cooling before peak velocity if the effective ablation pressure P abl (effective) = P abl (t pv - 0.5 ns) is used in the analytic formulas, where t pv is the time of peak implosion velocity. Second, we provide an analysis that enables a comparison of the Hydra radiation hydrodynamics code calculations utilized here with the predictions of the analytic piston model of an ICF implosion, which focused on sensitivity to time duration of the hohlraum cooling phase before peak velocity (often called the “coast time”) and the shell radius at peak velocity R pv . Third, we provide a set of ignition metrics that are valid across a wide range of capsule designs valid for implosions both with and without hohlraum cooling before peak implosion velocity is reached.

42 ENGINEERING↗

Plate motion in sheared granular fault system

Plate motion near the fault gouge layer, and the elastic interplay between the gouge layer and the plate under stick-slip conditions, is key to understanding the dynamics of sheared granular fault systems. Here, a two-dimensional implementation of the combined finite-discrete element method (FDEM), which merges the finite element method (FEM) and the discrete element method (DEM), is used to explicitly simulate a sheared granular gouge fault system. In this work we focus on investigating the influence of normal load, driving shear velocity and plate stiffness on the velocities and displacements in the direction parallel to the shear direction (x-direction) measured at locations on the upper and lower plates just adjacent to the gouge. The simulations show that during slip phases the magnitudes of the measured velocities on the upper and lower plates are proportional to the normal load and may be inversely proportional to the square root of the plate's shear modulus. Whereas, the driving shear velocity does not show distinct influence on the measured velocities. Additionally, large slip velocities are generally associated with large macroscopic friction coefficient drops. For the models subjected to smaller normal loads, larger shear velocities and with stiffer shear plates, the same magnitude of slip velocity could cause a larger drop of macroscopic friction coefficient. During stick phases, the velocities of the upper and lower plates are respectively slightly greater and slightly smaller than half of the driving shear velocity and are both in the same direction of shear. The shear strain rate of the gouge is calculated from this velocity difference between the upper and lower plate during stick phases and thus the gouge effective shear modulus can be calculated. The results show that the gouge effective shear modulus increases proportionally with normal load, while the influence of shear velocity and plate stiffness on gouge effective shear modulus is minor. The simulations address the dynamics of a laboratory-scale fault gouge system and may aid in revealing the complexities of earthquake frictional dynamics.

58 GEOSCIENCES↗

A new algorithm for three‐dimensional joint inversion of body wave and surface wave data and its application to the Southern California plate boundary region

Abstract We introduce a new algorithm for joint inversion of body wave and surface wave data to get better 3‐D P wave ( V p ) and S wave ( V s ) velocity models by taking advantage of the complementary strengths of each data set. Our joint inversion algorithm uses a one‐step inversion of surface wave traveltime measurements at different periods for 3‐D V s and V p models without constructing the intermediate phase or group velocity maps. This allows a more straightforward modeling of surface wave traveltime data with the body wave arrival times. We take into consideration the sensitivity of surface wave data with respect to V p in addition to its large sensitivity to V s , which means both models are constrained by two different data types. The method is applied to determine 3‐D crustal V p and V s models using body wave and Rayleigh wave data in the Southern California plate boundary region, which has previously been studied with both double‐difference tomography method using body wave arrival times and ambient noise tomography method with Rayleigh and Love wave group velocity dispersion measurements. Our approach creates self‐consistent and unique models with no prominent gaps, with Rayleigh wave data resolving shallow and large‐scale features and body wave data constraining relatively deeper structures where their ray coverage is good. The velocity model from the joint inversion is consistent with local geological structures and produces better fits to observed seismic waveforms than the current Southern California Earthquake Center (SCEC) model.

Fang, Hongjian↗

An enhanced AUSM + -up scheme for high-speed compressible two-phase flows on hybrid grids

Here in this work, an enhanced AUSM + -up scheme is presented for high-speed compressible two-phase flows using a six-equation two-fluid single-pressure model. Based on the observation that the AUSM + -up flux function does not take into account relative velocity between the two phases and thus is not stable and robust for computation of two-phase flows involving interaction of strong shock waves and material interfaces, the enhancement is in the form of a volume fraction coupling term and a modification of the velocity diffusion term, both proportional to the relative velocity between the two phases. These modifications in the flux function obviate the need to employ the exact Riemann solver, leading to a significantly less expensive yet robust flux scheme. Furthermore, the Tangent of Hyperbola for INterface Capturing (THINC) scheme is used in order to provide a sharp resolution for material interfaces. A number of benchmark test cases are presented to assess the performance and robustness of the enhanced AUSM + -up scheme for compressible two-phase flows on hybrid unstructured grids. Lastly, the numerical experiments demonstrate that the enhanced AUSM + -up scheme along with THINC scheme can efficiently compute high-speed two-fluid flows such as shock–bubble interactions, while accurately capturing material interfaces.

42 ENGINEERING↗

Spatially resolved properties of supernova host galaxies in SDSS-IV MaNGA

We crossmatch galaxies from Mapping Nearby Galaxies at Apache Point Observatory with the Open Supernova Catalog, obtaining a total of 132 SNe within MaNGA bundle. These 132 SNe can be classified into 67 Type Ia and 65 Type CC. We study the global and local properties of supernova host galaxies statistically. Type Ia SNe are distributed in both star-forming galaxies and quiescent galaxies, while Type CC SNe are all distributed along the star-forming main sequence. As the stellar mass increases, the Type Ia/CC number ratio increases. We find: (1) there is no obvious difference in the interaction possibilities and environments between Type Ia SN hosts and a control sample of galaxies with similar stellar mass and SFR distributions, except that Type Ia SNe tend to appear in galaxies which are more bulge-dominated than their controls. For Type CC SNe, there is no difference between their hosts and the control galaxies in galaxy morphology, interaction possibilities as well as environments; (2) compared to galaxy centers, the SN locations have smaller velocity dispersion, lower metallicity, and younger stellar population. This is a natural result of radius gradients for all these parameters. In conclusion, the SN location and its symmetrical position relative to the galaxy center, as well as regions with similar effective radii have very similar [Mg/Fe], gas-phase metallicity, gas velocity dispersion and stellar population age.

79 ASTRONOMY AND ASTROPHYSICS↗

On the Preconditioning of a Newton-Krylov Solver for a High-Order reconstructed Discontinuous Galerkin Discretization of All-Speed Compressible Flow with Phase Change for Application in Laser-Based Additive Manufacturing

This dissertation focuses on the development of a fully-implicit, high-order compressible ow solver with phase change. The work is motivated by laser-induced phase change applications, particularly by the need to develop large-scale multi-physics simulations of the selective laser melting (SLM) process in metal additive manufacturing (3D printing). Simulations of the SLM process require precise tracking of multi-material solid-liquid-gas interfaces, due to laser-induced melting/ solidi cation and evaporation/condensation of metal powder in an ambient gas. These rapid density variations and phase change processes tightly couple the governing equations, requiring a fully compressible framework to robustly capture the rapid density variations of the ambient gas and the melting/evaporation of the metal powder. For non-isothermal phase change, the velocity is gradually suppressed through the mushy region by a variable viscosity and Darcy source term model. The governing equations are discretized up to 4th-order accuracy with our reconstructed Discontinuous Galerkin spatial discretization scheme and up to 5th-order accuracy with L-stable fully implicit time discretization schemes (BDF2 and ESDIRK3-5). The resulting set of non-linear equations is solved using a robust Newton-Krylov method, with the Jacobian-free version of the GMRES solver for linear iterations. Due to the sti nes associated with the acoustic waves and thermal and viscous/material strength e ects, preconditioning the GMRES solver is essential. A robust and scalable approximate block factorization preconditioner was developed, which utilizes the velocity-pressure (vP) and velocity-temperature (vT) Schur complement systems. This multigrid block reduction preconditioning technique converges for high CFL/Fourier numbers and exhibits excellent parallel and algorithmic scalability on classic benchmark problems in uid dynamics (lid-driven cavity ow and natural convection heat transfer) as well as for laser-induced phase change problems in 2D and 3D.

42 ENGINEERING↗

Orbit-based analysis of nonlinear energetic ion dynamics in tokamaks. II. Mechanisms for rapid chirping and convective amplification

The nonlinear interactions between shear Alfvén modes and tangentially injected beam ions in the 150–400 keV range are studied numerically in realistic geometry for a JT-60U tokamak scenario. In Paper I, which was reported in the companion paper, the recently developed orbit-based resonance analysis method was used to track the resonant frequency of fast ions during their nonlinear evolution subject to large magnetic and electric drifts. Here, that method is applied to map the wave-particle power transfer from the canonical guiding center phase space into the frequency-radius plane, where it can be directly compared with the evolution of the fluctuation spectra of fast-ion-driven modes. Using this technique, we study the nonlinear dynamics of strongly driven shear Alfvén modes with low toroidal mode numbers n = 1 and n = 3. In the n = 3 case, both chirping and convective amplification can be attributed to the mode following the resonant frequency of the radially displaced particles, i.e., the usual one-dimensional phase locking process. In the n = 1 case, a new chirping mechanism is found, which involves multiple dimensions, namely, wave-particle trapping in the radial direction and phase mixing across velocity coordinates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

EXPLAINING THE OBSERVED VELOCITY DISPERSION OF DWARF GALAXIES BY BARYONIC MASS LOSS DURING THE FIRST COLLAPSE

In the widely adopted {Lambda} cold dark matter ({Lambda}CDM) scenario for galaxy formation, dwarf galaxies are the building blocks of larger galaxies. Since they formed at relatively early epochs when the background density was relatively high, they are expected to retain their integrity as satellite galaxies when they merge to form larger entities. Although many dwarf spheroidal galaxies are found in the galactic halo around the Milky Way, their phase-space density (or velocity dispersion) appears to be significantly smaller than that expected for satellite dwarf galaxies in the {Lambda}CDM scenario. In order to account for this discrepancy, we consider the possibility that they may have lost a significant fraction of their baryonic matter content during the first infall at the Hubble expansion turnaround. Such mass loss arises naturally due to the feedback by relatively massive stars that formed in their centers briefly before the maximum contraction. Through a series of N-body simulations, we show that the timely loss of a significant fraction of the dSphs initial baryonic matter content can have profound effects on their asymptotic half-mass radius, velocity dispersion, phase-space density, and the mass fraction between residual baryonic and dark matter.

79 ASTRONOMY AND ASTROPHYSICS↗

RELAP-7 Application and Enhancement for FLEX Strategies and ATF Behavior under Extended Loss of AC Power Conditions

This report summarizes the results of a three-year research project sponsored by the U.S. Department of Energy (DOE) Nuclear Energy University Program (NEUP) to enhance and apply the RELAP-7 code by adding and improving several important components (e.g., a mechanistic Reactor Core Isolation Cooling (RCIC) system model) for thermal hydraulic studies of LWRs under ELAP conditions and evaluating the time available for transition to portable FLEX equipment. The project team included University of Massachusetts–Lowell, The Ohio State University, Texas A&M University, Idaho National Laboratory and Oak Ridge National Laboratory. In the Fukushima accident, it was found that the RCIC system played a crucial role in delaying core meltdown by almost three days in Fukushima Daiichi Unit 2, because of self-regulated operation of the steam driven RCIC turbine-pump injection system. Steam flow in the convergent-divergent nozzles of the RCIC Terry turbine is two-phase non-equilibrium transonic flow with homogenous nucleation condensation. To more accurately predict the dynamic process and behavior of the transonic compressible steam flow, a one-dimensional transient two-phase analytical model is presented. A simplified four-fluid model was employed in the present work with the consideration of four separate fluid fields: vapor, liquid film, entrained droplets and condensed droplets. The mass, momentum and energy interactions between the fluids were considered and modeled. An extended seven-equation non-equilibrium critical flow model was developed to obtain the critical pressure and velocities of each phase at the nozzle throat. To predict the wetness in the divergent section, a mechanistic nucleation condensation model was integrated in the nozzle analysis model, considering the generation and consequent growth of droplets. The governing differential equations on a staggered grid were discretized using the second-order Lax-Wendroff scheme with a flux limiter, and the Semi-Implicit Method for Pressure-Linked Equation (SIMPLE) algorithm was employed to solve the discrete linear system. To demonstrate the predictability and reliability of the physical models and the numerical method proposed in the present work, three representative nozzles were modeled and simulated. The results show good agreement with the available experimental data, even for condensation shock. Then, the 1D nozzle model was employed to obtain nozzle flow tables of the Terry turbine nozzle for different working pressures which can cover the operation pressure range of the RCIC system. A mechanistic RCIC turbine-pump system model was developed and implemented in the system code TRACE to simulate dynamic responses of the RCIC system under Beyond Design Basis Accident (BDBA) conditions. The turbine-pump governing equations are based on the control volume approach of the angular momentum balance. The physics based mechanistic RCIC model was developed using the TRACE control system components (i.e., signal variables, control blocks, and tables), and incorporated into a TRACE boiling water reactor (BWR) model. The TRACE model in this report has a detailed nodalization of the reactor pressure vessel (RPV), and all of the major flow paths and system components, including the safety relief valves (SRVs) and the containment suppression pool and drywell. Based on the nozzle flow tables generated from the 1D nozzle model developed, the turbine drive torque can be calculated from table lookup. Since the detailed specifications of the RCIC pump are unavailable, the homologous curves for a Bingham pump were used in the current pump component. A station black-out (SBO) accident test problem was selected to demonstrate the TRACE RCIC model. The short-term SBO simulations were performed for two cladding materials: Zircaloy and FeCrAl, to demonstrate the effect of the accident tolerant fuel cladding on fuel heat-up under BDBA accident conditions. The wetwell plays a vital safety role in SBO and other BWR accident scenarios in that it can reduce containment pressure and supply additional core make-up water. The suppression pool temperature distribution has a very large impact on both RPV and containment pressure. Thus, another novel contribution of the project comes mainly from an improved, systems-level wetwell model which can capture buoyancy-induced thermal stratification effects due to steam injection and condensation. A two-zone stratified wetwell model has been implemented in RELAP-7 and some results from that model are presented. This wetwell model is capable of simulating thermal stratification due to a low steam mass injection rate. With a low mass flow rate, the model assumes that all the steam condenses within the pipe and the resulting plume can be approximated with a purely buoyant, heat-source driven model. The wetwell model developed with these assumptions is adequate to simulate slow transients such as extended SBO transients.

42 ENGINEERING↗

How important are electron beams in driving chromospheric evaporation in the 2014 March 29 flare?

We present high spatial resolution observations of chromospheric evaporation in the flare SOL2014-03-29T17:48. Interface Region Imaging Spectrograph observations of the Fe xxi λ1354.1 line indicate evaporating plasma at a temperature of 10 MK along the flare ribbon during the flare peak and several minutes into the decay phase with upflow velocities between 30 and 200 km s{sup −1}. Hard X-ray (HXR) footpoints were observed by the Ramaty High Energy Solar Spectroscopic Imager for two minutes during the peak of the flare. Their locations coincided with the locations of the upflows in parts of the southern flare ribbon but the HXR footpoint source preceded the observation of upflows in Fe xxi by 30–75 s. However, in other parts of the southern ribbon and in the northern ribbon, the observed upflows were not coincident with an HXR source in time or space, most prominently during the decay phase. In this case evaporation is likely caused by energy input via a conductive flux that is established between the hot (25 MK) coronal source, which is present during the whole observed time-interval, and the chromosphere. The presented observations suggest that conduction may drive evaporation not only during the decay phase but also during the flare peak. Electron beam heating may only play a role in driving evaporation during the initial phases of the flare.

79 ASTRONOMY AND ASTROPHYSICS↗

The physics of spontaneous parity-time symmetry breaking in the Kelvin–Helmholtz instability

Abstract We show that the dynamical system of an inviscid fluid with velocity shear admits parity-time (PT) symmetry, which provides a physical explanation of the well-known observation that the spectrum of the perturbation eigenmodes of the system is symmetric with respect to the real axis. It is found that the Kelvin–Helmholtz instability is triggered when and only when the PT symmetry is spontaneously broken. The analysis of PT symmetry also reveals that the relative phase between parallel velocity and pressure perturbations needs to be locked at π /2 when the instability is suppressed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Elucidating the Microscale Behavior and Phase Separation Kinetics of Thermally Responsive Ionic Liquid–Water Mixtures

Thermally responsive ionic liquids (ILs) exhibit liquid-liquid phase separation into a water-rich (WR) and ionic-liquid-rich (ILR) phase when heated above a lower critical solution temperature (LCST). This phase behavior has been leveraged for applications ranging from forward osmosis (FO) desalination, where the IL acts as a draw solute, to refrigeration and dehumidification cycles, where the IL acts as a liquid desiccant. While significant effort has been devoted to characterizing the thermodynamic and thermophysical properties of LCST ILs, their phase separation kinetics have not been investigated. In this work, we describe the macroscale phase separation kinetics (phase separation time) by gleaning insight into the microscale colloidal behavior of aqueous mixtures of four different materials, P 4444 TFA (tetrabutylphosphonium-2,4-trifluoroacetate), P 4444 DMBS (tetrabutylphosphonium-2,4-dimethyl-benzenesulfonate), N 4444 Sal (tetrabutylammonium salicylate), and P 4444 Sal (tetrabutylphosphonium salicylate) as a function of IL concentration at a separation temperature of 70 °C. We report the discontinuous microscale size distributions for each material and correlate their theoretical settling velocities to experimental phase separation times. The results indicate that a simple Stokes' law model can predict the phase separation time within reasonable accuracy. Overall, this work lays the foundation for understanding the micro- to macroscale phase separation behavior and kinetics of LCST ILs for various water-energy applications.

LCST↗

Sequentially prepared Mo-V-Based SCR catalyst for simultaneous Hg 0 oxidation and NO reduction

Molybdenum (Mo)-vanadium (V)-based selective catalytic reduction (SCR) catalyst synthesized by the sequential impregnation of Mo and W followed by V was investigated for simultaneous elemental mercury (Hg 0 ) oxidation and nitrogen oxide (NO) reduction in an existing SCR unit with respect to different TiO 2 phases, calcination temperatures, flue gas constituents, gas velocities, and reaction temperatures. Anatase phase TiO 2 and the calcination temperatures of 400 and 500 degrees C resulted in similar to 69% Hg 0 oxidation at 10 ppmv HCl and 350 degrees C. The high calcination temperature of 700 degrees C resulted in TiO 2 phase transformation from anatase to rutile and agglomeration. The modified SCR catalyst prepared with the impregnation sequence of Mo and W followed by V using anatase TiO 2 and calcination temperature 500 degrees C showed similar to 99% Hg 0 oxidation and 87% NO reduction conversions at an NH 3 /NO molar ratio of 0.9 under 350 degrees C and 5,000 hr -1 space velocity in typical sub-bituminous and lignite coal simulated flue gases. Finally, the effects of these parameters and conditions were further investigated using various characterization techniques including BET, TEM, XRD, NH 3 TPD and XAFS.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗