Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fluid wave”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Coherent fluctuations in a partially magnetized small magnetron discharge operating on molecular nitrogen

We report experimental data on rotating plasma disturbances in a quasi-steady direct-current magnetron discharge operating with molecular nitrogen at pressures of 27, 40, 67, and 133 Pa. Experiments with a segmented anode measuring the local current flow reveal up to 12 coherent structures (“spokes”) propagating in the +E×B and −E×B azimuthal directions. The current–voltage characteristics are studied and found to follow a negative resistance profile at the lower discharge currents. The existence of coherent structures and the number of spokes observed depend on the discharge current and operating pressure. Lower pressures are more favorable to the retrograde −E×B spoke rotation, which is observed in the 27, 40, and 67 Pa discharges. The results are interpreted with a gradient-driven drift wave model, and using the dispersion relation, we show that the propagation direction depends on the orientation of the local electric field within the plasma. Using an electron fluid model, we show how experimental conditions can lead to field reversals that generate −E×B spoke rotation in the lower pressure discharges.

Physics↗
DOE OSTI2023

Current fluctuations and domain depinning in quasi-two-dimensional charge-density-wave 1 T -TaS 2 thin films

We investigated the temperature dependence of the current fluctuations in thin films of the quasi-two-dimensional 1T-TaS 2 van der Waals material. The current fluctuations, determined from the derivative current–voltage characteristics of two-terminal 1T-TaS 2 devices, appear prominently at the electric fields that correspond to the transitions between various charge-density-wave macroscopic quantum condensate phases and at the onset of the depinning of the charge density wave domains. The depinning threshold field, ED, monotonically increases with decreasing temperature within the nearly commensurate charge-density-wave phase. The ED value increases with the decreasing 1T-TaS 2 film thickness, revealing the surface pinning of the charge density waves. Our analysis suggests that the domain depinning is pronounced in the nearly commensurate phase. It is induced by the electric field but facilitated by local heating. Furthermore, the measured trends for ED of the domain depinning are important for understanding the physics of charge density waves in quasi-two-dimensional crystals and for developing electronic devices based on this type of quantum materials.

2D materials↗
DOE OSTI2024

Quasilinear theory: the lost ponderomotive effects and why they matter

Quasilinear theory (QLT) has been used for modeling wave–plasma interactions for decades but remains largely heuristic. Plasma inhomogeneity, ponderomotive effects, microscopic fluctuations, and collisions are not easily accommodated from first principles in QLT, and typically are ignored entirely, due to the limitations of the standard Fourier–Laplace global-mode approach. This results in inconsistencies, for example, violation of the action conservation for nonresonant waves. However, these issues can be avoided, and the theory can be substantially generalized and corrected, if QLT is formulated using more suitable analytical tools, particularly, the Weyl symbol calculus. Here, an attempt is made to deliver an accessible review of this modern formulation, provide intuitive calculations for special cases, and elaborate on the connection with the ‘oscillation-center QLT’ originally proposed by Dewar (Phys Fluids 16:1102, 1973). A Fokker–Planck equation for a ‘dressed’ distribution is derived from the Klimontovich equation and captures quasilinear diffusion, ponderomotive forces, and interactions with background fields for a generic Hamiltonian, so many known formulations of QLT for specific plasma models become corollaries of a single unifying theory. Also, waves are allowed to be off-shell (not constrained by a dispersion relation), which allows them to accommodate microscopic fluctuations. This leads to a collision integral of the Balescu–Lenard type that has all the usual properties but is not restricted to any specific plasma model. For on-shell waves, a generalized version of the classic oscillation-center QLT is obtained. Finally, combined with the wave-kinetic equation, this formulation not only conserves particles, momentum, and energy, like the classic QLT but also reinstates the action conservation for nonresonant waves.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Design and Performance Evaluation of a Resistive Control Using a Hydraulic PTO System for the TALOS Wave Energy Converter

This study is focused on developing a numerical model to evaluate the performance of a hydraulic PTO system for the TALOS Wave Energy Converter. The WEC device is described and the architecture of the hydraulic PTO system is presented with detail. The WEC is modeled using WEC-Sim, and the PTO is modeled using the Simscape Fluids library from Simulink. The hydraulic PTO is based on a constant pressure configuration that is suitable for WEC passive control. The hydraulic system is composed by a set of rectifying valves and two hydraulic accumulators that reduce the stiffness of the system and also serve as energy storage devices. One of the advantages of this hydraulic PTO architecture is the possibility of controlling the electric generator to operate around the optimal efficiency operating point. The main components of the hydraulic PTO are off-the-shelf devices that are commercially available, which will facility a future deployment of the designed system. The design variables used for this study are the accumulator size, the maximum pressure in the accumulators, the hydraulic motor maximum displacement, and the shaft speed in the electric generator. The performance of the system is evaluated individually, using sinusoidal inputs that replicates regular wave conditions. In addition to this, the numerical model of the PTO is coupled to a WEC-Sim simulation of the TALOS Wave Energy Converter with six PTOs to generate a wave-to-wire model. The main objective of this work is to present a comprehensive design methodology that could serve as a guideline for future research efforts focused on implementing control algorithms on multi degree of freedom WECs.

hydraulic systems↗
DOE OSTI2022

Quasinormal g-modes of neutron stars with quarks

Quasinormal oscillation modes of neutron stars provide a means to probe their interior composition using gravitational wave astronomy. We compute the frequencies and damping times of composition-dependent core 𝑔-modes of neutron stars containing quark matter employing linearized perturbative equations of general relativity. We find that ignoring background metric perturbations due to the oscillating fluid, as in the Cowling approximation, underestimates the 𝑔-mode frequency by up to 10% for higher mass stars, depending on the parameters of the nuclear equation of state and how the mixed phase is constructed. The 𝑔-mode frequencies are well described by a linear scaling with the central lepton (or combined lepton and quark) fraction for nucleonic (hybrid) stars. Furthermore, our findings suggest that neutron stars with and without quarks are manifestly different with regards to their quasinormal 𝑔-mode spectrum, and may thus be distinguished from one another in future observations of gravitational waves from merging neutron stars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Center of Excellence Collaboration Projects: Second Wave

The NEAMS program aims to develop an integrated multi-physics simulation capability “pellet-to-plant” for the design and analysis of future generations of nuclear power plants. In particular, the Reactor Product Line code suite's multi-resolution hierarchy is being designed to ultimately span the full range of length and time scales present in relevant reactor design and safety analyses, as well as scale from desktop to petaflop computing platforms. In particular the NEAMS program is supporting the development of novel thermal-hydraulic codes. The Center of Excellence for Thermal Fluids Application in Nuclear Energy, launched in 2018, has as its key goals to serve as a front door to industry. The Center of Excellence for Thermal Fluids Applications in Nuclear Energy has recently launched a program to start collaborative efforts between the laboratories and industry with the objective of stimulating cooperation and increasing adoption of T/H tools developed under NEAMS by the industry-at-large. In particular, two industry partners agreed to participate in a second wave of short-term collaborations aimed at demonstrating the value of NEAMS tools to their designs. These are Framatome and TerraPower who proposed projects related to LWR (Light Water Reactor) and MSFR (Molten Salt Fast Reactor) designs respectively. In this report we present the results of these two collaborations.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Online Lead/Water Heat Exchanger Sensor/System Feasibility-PNNL (Final Report)

On-line structural health corrosion monitoring in advanced lead fast reactor heat exchangers and molten salt reactor heat exchangers is desirable for detecting tube degradation prior to leaks that may allow mixing of heat exchanger fluids or release of radiological contamination beyond the design containment boundary. This program, On-Line Lead/Water Heat Exchanger Sensor/System Feasibility – PNNL 76092, demonstrates feasibility for an ultrasonic torsional wave mode sensor attached to the outside of a long (30.5-m) heat exchanger tube in the stagnant flow area where the tube joins a heat exchanger plenum and where it is possible to protect a sensor and cable from high-force flows. The sensor must be cable connected to a monitoring instrument near the heat exchanger. The sensor and cable management approach for periodic in-service inspection will be impractical to implement on existing heat exchangers; rather permanently mounted sensors must be installed in conjunction with heat exchanger fabrication. Previous work has shown low-temperature lead zirconate titanate (PZT) piezoceramic sensors are able to detect anomalies of interest in 3.0-m long tubes. These sensors have hereby been extended to a 30.5-m long tube more representative of commercial power heat exchanger designs. The program will continue to investigate higher temperature piezoelectric ceramics and long-term performance of high temperature adhesives and sealants and mechanical pressure coupling for 350-500 °C lead reactor environments and higher temperature (700 °C) molten salt environments.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Sensitivity measurements for a 250 MHz quartz shear-horizontal surface acoustic wave biosensor under liquid viscous loading

Surface acoustic wave (SAW) devices have been used in biochemical assays due to their high sensitivity. The device sensitivity is a function of changes in the density and viscosity of the liquid. Here, we studied the effect of fluid viscosity using a 250 MHz quartz shear-horizontal (SH)-SAW biosensor by monitoring different concentrations of binary aqueous/glycerol solutions. In this study, the sensitivity of the biosensor was determined by fitting the data to models derived from perturbation theory. Measurements in water were used as the reference. For a 0% to 50% glycerol solution, an 87°–204° separation in the phase shift was observed. The slope of the plot of the phase shift vs (ηρ)0.5 was used to indicate the sensor’s sensitivity. The sensitivity for our 250 MHz quartz SH-SAW sensors was calculated to be 3.7×10−3m2sKg. The corresponding mass sensitivity was determined to be 9.25 × 105m2Kg. The limit of detection was calculated to be 36 picograms (pg), while the limit of quantification or LOQ was calculated to be 109 pg. Traditionally, liquid phase measurements have been challenging for SAW devices because liquids dampen the vibrating sensors severely. This problem has been largely solved using a transverse (shear) wave instead of the more popular longitudinal or Rayleigh waves. Liquid measurements are now possible using transverse waves, also known as shear waves, because transverse waves are only minimally attenuated by liquids. Shear-horizontal SAW sensors (SH-SAW) show great promise as label-free biosensors because of their ability to handle liquid samples. However, the viscosity of the liquid still induces loading effects and can be measured when the liquid is loaded onto the SH-SAW propagating surface (delay line). When the liquid above the delay line is perturbed by physical or chemical changes, such as binding to a receptor, it alters the propagating acoustic wave. The SH-SAW device can measure these changes in liquid properties as a change in the wave’s phase compared to the original wave. The device’s phase shift was recorded as a function of the changes in the density and viscosity of the binary glycerol solution and used to determine the sensitivity in the linear dynamic range of responses.

Materials Science↗
DOE OSTI2020

Variational nonlinear WKB in the Eulerian frame

Nonlinear WKB is a multiscale technique for studying locally plane-wave solutions of nonlinear partial differential equations (PDEs). Its application comprises two steps: (1) replacement of the original PDE with an extended system separating the large scales from the small and (2) reduction of the extended system to its slow manifold. In the context of variational fluid theories with particle relabeling symmetry, nonlinear WKB in the mean Eulerian frame is known to possess a variational structure. This much has been demonstrated using, for instance, the theoretical apparatus known as the generalized Lagrangian mean. On the other hand, the variational structure of nonlinear WKB in the conventional Eulerian frame remains mysterious. By exhibiting a variational principle for the extended equations from step (1) above, here we demonstrate that nonlinear WKB in the Eulerian frame is in fact variational. Remarkably, the variational principle for the extended system admits loops of relabeling transformations as a symmetry group. Noether’s theorem therefore implies that the extended Eulerian equations possess a family of circulation invariants parameterized by S 1 . As an illustrative example, we use our results to systematically deduce a variational model of high-frequency acoustic waves interacting with a larger-scale compressible isothermal flow.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

The influence of the shock-to-reshock time on the Richtmyer–Meshkov instability in reshock

Experiments on the Richtmyer–Meshkov instability (RMI) in a dual driver vertical shock tube (DDVST) are described. An initially planar, stably stratified membraneless interface is formed by flowing air from above and sulfur hexafluoride from below the interface location using the method of Jones & Jacobs ( Phys. Fluids , vol. 9, issue 1997, 1997, pp. 3078–3085). A random three-dimensional, multi-modal initial perturbation is imposed by vertically oscillating the gas column to produce Faraday waves. The DDVST design generates two shock waves, one originating above and one below the interface, with these shocks having independently controllable strengths and interface arrival times. The shock waves have nominal strengths of $M_L=1.17$ and $M_H=1.18$ for the shock wave originating in the light and heavy gas, respectively, with these strengths chosen to result in arrested bulk interface motion following reshock. The influence of the length of the shock-to-reshock time, as well as the order of shock arrival, on the post-reshock RMI is examined. The mixing layer width grows according to $h\propto t^\theta$ , where $\theta _H=0.36\pm 0.018$ (95 %) and $\theta _L=0.38\pm 0.02$ (95 %) for heavy and light shock first experiments, respectively, indicating no strong dependence on the order of shock wave arrival. Volume integrated specific turbulent kinetic energy (TKE) in the mixing layer versus time is found to decay according to $E_{tot}/\bar {\rho }\propto t^p$ with $p_H=-0.823\pm 0.06$ (95 %) and $p_L=-1.061\pm 0.032$ (95 %) for heavy and light shock first experiments, respectively. Notably, the 95 % confidence intervals do not overlap. Analysis on the influence of the shock-to-reshock time on turbulent length scales, transition criteria, spectra and mixing layer anisotropy are also presented.

Ferguson, Kevin (ORCID:0000000166537482)↗
DOE OSTI2026

Modeling rf sheath formation in turbulent tokamak boundary plasma

During ICRF antenna operation, complex interactions between turbulent density profiles, nonlinear RF sheaths, and RF-induced convective transport are observed to alter plasma density in the tokamak edge [D’Ippolito et al., Nucl. Fusion 38, 1543 (1998)]. In this work, we explore the physics of such interactions via numerical modeling, using a nonlinear EM/plasma/sheath code (VSim) and profiles obtained from a fluid plasma turbulence code (Hermes) in a 3D slab domain containing biased side-wall limiters. RF-rectified sheath formation on antenna and limiter surfaces is observed as electromagnetic waves launched by the antenna are refracted through the turbulent density profile. On transport timescales, such sheath potentials have been shown to influence both the mean species density and its RMS fluctuation spectrum [Smithe et al., these proceedings]. On the faster RF timescales, we demonstrate that the converse is also true – regions of high plasma density near material surfaces give rise to the highest sheath potential amplitudes. When density is turbulent and spatially nonuniform, localized regions of high sheath potential (hotspots) may develop where high-density filaments intersect material surfaces. Such hotspots are of particular concern as sources of impurity sputtering, and we explore their behavior in response to changes both to the local plasma density and to antenna operating parameters and structure. Related results exploring the role of Faraday shields and/or enclosing structures in suppressing high sheath potentials for other devices (e.g. SPARC) will also be shown.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Shock equation of state experiments in MgO up to 1.5 TPa and the effects of optical depth on temperature determination

Laser-driven shock compression enables an experimental study of phase transitions at unprecedented pressures and temperatures. One example is the shock Hugoniot of magnesium oxide (MgO), which crosses the B1–B2-liquid triple point at 400–600 GPa, 10 000–13 000 K (0.86–1.12 eV). MgO is a major component within the mantles of terrestrial planets and has long been a focus of high-pressure research. Here, we combine time-resolved velocimetry and pyrometry measurements with a decaying shock platform to obtain pressure–temperature data on MgO from 300 to 1500 GPa and 9000 to 50 000 K. Pressure–temperature–density Hugoniot data are reported at 1500 GPa. These data represent the near-instantaneous response of an MgO [100] single crystal to shock compression. We report on a prominent temperature anomaly between 400 and 460 GPa, in general agreement with previous shock studies, and draw comparison with equation-of-state models. We provide a detailed analysis of the decaying shock compression platform, including a treatment of a pressure-dependent optical depth near the shock front. We show that if the optical depth of the shocked material is larger than 1 μm, treating the shock front as an optically thick gray body will lead to a noticeable overestimation of the shock temperature.

36 MATERIALS SCIENCE↗
DOE OSTI2025

An infrasound source analysis of the OSIRIS-REx sample return capsule hypersonic re-entry

The OSIRIS-REx sample return capsule's hypersonic re-entry into the atmosphere is a rare opportunity to test a variety of sonic boom source models since the projectile dimensions are well characterized. While the as-flown flight path is unknown, the predicted flight path enables a rough approximation of the source Mach number and location. Six infrasound microphones deployed in the boom carpet along the predicted flight path recorded impulsive signals from the OSIRIS-REx re-entry. Using a suite of atmosphere profiles and the geometric acoustics approximation, we estimate locations with uncertainty estimates along the flight path from which the signals were emitted. Acoustic overpressure and signal duration predictions from Whitham's far field theory, Carlson's simplified sonic boom prediction method, and a drag-dominated hypersonic model are analyzed with uncertainty estimates from the location estimate. While the Carlson simplified sonic boom prediction method could be accurate, our preference is for the drag-dominated source model. Using this source model with an inviscid Burgers's equation solver for propagation, we obtained an excellent match to the recorded data. In conclusion, these results will help better inform future sample return capsule re-entry observation campaigns as well as contribute to a better understanding of high altitude infrasonic sources.

58 GEOSCIENCES↗
DOE OSTI2023

TEAMER: Experimental performance characterization of a shrouded axial-flow turbine

Sitkana has developed a shrouded hydrokinetic turbine with a modular, low-cost design that can be scaled to meet the needs of remote communities. With technical support from the University of Washington, Sitkana sought to experimentally characterize the mechanical power and structural loads of various 1:3.3 scale rotor geometries. In all, 11 different rotor geometries were characterized with variations in height-to-diameter ratio, blade number, and blade type (foiled versus flat). All tests were conducted in Reynolds-independent flow conditions in the Alice C. Tyler Flume at the University of Washington. Results allow Sitkana to (1) refine the optimal rotor geometry, (2) validate numerical models, and (3) predict power output for a full-scale system. This project is part of the TEAMER RFTS 8 (request for technical support) program.

16 TIDAL AND WAVE POWER↗
DOE OSTI2021

Alpha particle driven Alfvénic instabilities in ITER post-disruption plasmas

Fusion-born alpha particles in ITER disruption simulations are investigated as a possible drive of Alfvénic instabilities. The ability of these waves to expel runaway electron (RE) seed particles is explored in the pursuit of a passive, inherent RE mitigation scenario. The spatiotemporal evolution of the alpha particle distribution during the disruption is calculated using the linearized Fokker–Planck solver CODION coupled to a fluid disruption simulation. These simulations are done in the limit of no alpha particle transport during the thermal quench, which can be seen as a most pessimistic situation where there is also no RE seed transport. Under these assumptions, the radial anisotropy of the resulting alpha population provides free energy to drive Alfvénic modes during the quench phase of the disruption. We use the linear gyrokinetic magnetohydrodynamic code LIGKA to calculate the Alfvén spectrum and find that the equilibrium is capable of sustaining a wide range of modes. The self-consistent evolution of the mode amplitudes and the alpha distribution is calculated utilizing the wave-particle interaction tool HAGIS. Intermediate mode number (n = 7–15, 22–26) toroidal Alfvén eigenmodes are shown to saturate at an amplitude of up to δB/B ≈ 0.1% in the spatial regimes crucial for RE seed formation. We find that the mode amplitudes are predicted to be sufficiently large to permit the possibility of significant radial transport of REs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Population balance modeling of polyurethane foam formation with pressure‐dependent growth kernel

Abstract Polyurethane foams are widely used materials often chosen for their useful characteristics such as low thermal conductivity, ease of application, and high strength‐to‐weight ratios. Computational models are needed to predict the dynamics of the flow and expansion, and the resulting material properties, to improve manufacturing processes. In this paper, a model for PMDI, a water‐blown polyurethane foam, is presented. By extending a kinetics‐based approach by adding bubble‐scale information via a population balance equation (PBE) using the quadrature method of moments, we can track bubble size distributions during foaming. We present results from a three‐dimensional computational fluid dynamics model using arbitrary Lagrangian–Eulerian interface tracking implemented in finite element software. The model compares favorably with experimental data, including dynamics, bubble distributions measured by both camera and diffusion wave spectroscopy, and post‐test bubble size from scanning electron microscopy and density measurements from x‐ray computed tomography.

Ortiz, Weston↗
DOE OSTI2021

A process to verify numerical models for seismic fluid-structure interaction in advanced reactor vessels

Seismic design and qualification of a liquid-filled advanced nuclear reactor will have to account for fluid-structure interaction (FSI). Interaction between the tank, internal components, and contained liquid will rely on analysis of numerical models that must be verified and validated. Here this study demonstrates a verification process for models of a base-supported cylindrical tank by comparing numerical predictions and analytical solutions. The numerical models are consistent with the assumptions made to derive analytical solutions, namely, either a rigid or a linear elastic tank, ideal fluid, and small-amplitude, unidirectional, horizontal inputs. One software platform is used to illustrate the process. Seismic FSI analysis is performed using the Arbitrary Lagrangian-Eulerian (ALE) and Incompressible Computational Fluid Dynamics (ICFD) solvers in LS-DYNA. Reported responses are those used for design, including hydrodynamic pressures on the tank wall, shear forces and moments at the tank base, and wave heights of the contained liquid. The accuracy of the numerical results is discussed. The numerical models are verified for calculating the pressures on the tank wall and reactions at its base. Accurate simulation of wave action is challenging for both solvers. Recommendations for modeling, code development, and steps for verification are provided. Although focused on reactor vessels and one software platform, the verification process described herein is broadly applicable to liquid-filled vessels and other finite element codes.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Theoretical Control-Centric Modeling for Precision Model-Based Sliding Mode Control of a Hydraulic Artificial Muscle Actuator

Artificial muscles (AMs) typically rely on pneumatic sources of fluid power. The use of hydraulics can increase the power and force to weight and volume ratios of AM actuators. This paper develops a control-centric third-order single-input single-output (SISO) lumped-parameter dynamic model and sliding mode position controller based on Filippov's principle of equivalent dynamics for a braided hydraulic artificial muscle (HAM) actuator. The model predicts the nonlinear behavior of the HAM free contraction and captures the fluid and actuator nonlinear dynamic interactions in addition to the braid deformation. Model simulations are compared to experimental results for quasi-static pressurization, isometric pressurization, and open-loop square wave commands at 0.25, 0.5, and 1 Hz. Experiments of sine wave tracking at 0.25, 0.5, and 1 Hz and continuous square wave tracking at 0.067 Hz are conducted using a sliding mode controller (SMC) derived from the model. The SMC achieves a steady-state error of 6 μm at multiple setpoints within the actuator's 17 mm stroke. Compared to a proportional-integral-derivative (PID) controller, the SMC root-mean-square (RMS) error, mean error, and absolute maximum error are reduced on average by 53%, 61%, and 44%, respectively, demonstrating the benefit of model-based approaches for controlling HAMs.

42 ENGINEERING↗