SEARCH · Engineering Papers

Results for “simulated waves”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,513 records · Page 84

All sources NASA NTRS DOE OSTI

NASA NTRS1993

Interacting scales and energy transfer in isotropic turbulence

The dependence of the energy transfer process on the disparity of the interacting scales is investigated in the inertial and far-dissipation ranges of isotropic turbulence. The strategy for generating the simulated flow fields and the choice of a disparity parameter to characterize the scaling of the interactions is discussed. The inertial range is found to be dominated by relatively local interactions, in agreement with the Kolmogorov assumption. The far-dissipation is found to be dominated by relatively non-local interactions, supporting the classical notion that the far-dissipation range is slaved to the Kolmogorov scales. The measured energy transfer is compared with the classical models of Heisenberg, Obukhov, and the more detailed analysis of Tennekes and Lumley. The energy transfer statistics measured in the numerically simulated flows are found to be nearly self-similar for wave numbers in the inertial range. Using the self-similar form measured within the limited scale range of the simulation, an 'ideal' energy transfer function and the corresponding energy flux rate for an inertial range of infinite extent are constructed. From this flux rate, the Kolmogorov constant is calculated to be 1.5, in excellent agreement with experiments.

Zhou, YE↗

NASA NTRS1994

Effect of high-latitude ionospheric convection on Sun-aligned polar caps

A coupled magnetospheric-ionospheric (M-I) magnetohydrodynamic (MHD) model has been used to simulate the formation of Sun-aligned polar cap arcs for a variety of interplanetary magnetic field (IMF) dependent polar cap convection fields. The formation process involves launching an Alfven shear wave from the magnetosphere to the ionosphere where the ionospheric conductance can react self-consistently to changes in the upward currents. We assume that the initial Alfven shear wave is the result of solar wind-magnetosphere interactions. The simulations show how the E region density is affected by the changes in the electron precipitation that are associated with the upward currents. These changes in conductance lead to both a modified Alfven wave reflection at the ionosphere and the generation of secondary Alfven waves in the ionosphere. The ensuing bouncing of the Alfven waves between the ionosphere and magnetosphere is followed until an asymptotic solution is obtained. At the magnetosphere the Alfven waves reflect at a fixed boundary. The coupled M-I Sun-aligned polar cap arc model of Zhu et al.(1993a) is used to carry out the simulations. This study focuses on the dependence of the polar cap arc formation on the background (global) convection pattern. Since the polar cap arcs occur for northward and strong B(sub y) IMF conditions, a variety of background convection patterns can exist when the arcs are present. The study shows that polar cap arcs can be formed for all these convection patterns; however, the arc features are dramatically different for the different patterns. For weak sunward convection a relatively confined single pair of current sheets is associated with the imposed Alfven shear wave structure. However, when the electric field exceeds a threshold, the arc structure intensifies, and the conductance increases as does the local Joule heating rate. These increases are faster than a linear dependence on the background electric field strength. Furthermore, above the threshold, the single current sheet pair splits into multiple current sheet pairs. For the fixed initial ionospheric and magnetospheric conditions used in this study, the separation distance between the current pairs was found to be almost independent of the background electric field strength. For either three-cell or distorted two-cell background convection patterns the arc formation favored the positive B(sub y) case in the northern hemisphere.

Sojka, J. J.↗

DOE OSTI2023

Flow Strength Measurements of Wrought and AM SS304L via Pressure Shear Plate Impact Experiments

Pressure-shear plate impact experiments were performed to quantify flow strength of wrought, as-built additively manufactured (AM), and heat-treated and recrystallized AM 304 L stainless steel (SS304L) under combined loading. Impact velocities spanned between 0.03 and 0.24 mm/μs, resulting in corresponding pressures of 0.62–5.93 GPa. Flow strength measurements are comparable for the sample variants across the studied loading conditions; however, shear wave structures significantly differ between sample type. Microstructurally aware simulations indicate local strain differences attributed to anisotropic elastic constants of large grains (~1 mm) in the as-built and heat-treated AM may impede the ability to uniformly transmit a shear wave.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Bosonic field digitization for quantum computers

Quantum simulation of quantum field theory is a flagship application of quantum computers that promises to deliver capabilities beyond classical computing. The realization of quantum advantage will require methods that can accurately predict error scaling as a function of the resolution and parameters of the model and that can be implemented efficiently on quantum hardware. In this paper, we address the representation of lattice bosonic fields in a discretized field amplitude basis, develop methods to predict error scaling, and present efficient qubit implementation strategies. A low-energy subspace of the bosonic Hilbert space, defined by a boson occupation number cutoff, can be represented with exponentially good accuracy by a low-energy subspace of a finite-size Hilbert space. The finite representation construction and the associated errors are directly related to the accuracy of the Nyquist-Shannon sampling and the finite Fourier transforms of the boson number states in the field and the conjugate-field bases. We analyze the relation between the boson mass, the discretization parameters used for wave function sampling, and the finite representation size. Numerical simulations of small size Φ 4 problems demonstrate that the boson mass optimizing the sampling of the ground state wave function is a good approximation to the optimal boson mass yielding the minimum low-energy subspace size. However, we find that accurate sampling of general wave functions does not necessarily result in accurate representation. Finally, we develop methods for validating and adjusting the discretization parameters to achieve more accurate simulations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS1985

Acoustic waves in the solar atmosphere. VII - Non-grey, non-LTE H(-) models

The propagation and shock formation of radiatively damped acoustic waves in the solar chromosphere are studied under the assumption that H(-) is the only absorber; the opacity is non-grey. Deviations from local thermodynamic equilibrium (LTE) are permitted. The results of numerical simulations show the depth dependence of the heating by the acoustic waves to be insensitive to the mean state of the atmosphere. After the waves have developed into shocks, their energy flux decays exponentially with a constant damping length of about 1.4 times the pressure scale height, independent of initial flux and wave period. Departures from LTE have a strong influence on the mean temperature structure in dynamical chromosphere models; this is even more pronounced in models with reduced particle density - simulating conditions in magnetic flux tubes - which show significantly increased temperatures in response to mechanical heating. When the energy dissipation of the waves is sufficiently large to dissociate most of the H(-) ions, a strong temperature rise is found that is reminiscent of the temperature structure in the transition zone between chromosphere and corona; the energy flux remaining in the waves then drives mass motions.

Schmitz, F.↗

DOE OSTI2021

Dispersion and the speed-limited particle-in-cell algorithm

This paper discusses temporally continuous and discrete forms of the speed-limited particle-in-cell (SLPIC) method first treated by Werner et al. [Phys. Plasmas 25, 123512 (2018)]. The dispersion relation for a 1D1V electrostatic plasma whose fast particles are speed-limited is derived and analyzed. By examining the normal modes of this dispersion relation, we show that the imposed speed-limiting substantially reduces the frequency of fast electron plasma oscillations while preserving the correct physics of lower-frequency plasma dynamics (e.g. ion acoustic wave dispersion and damping). We then demonstrate how the timestep constraints of conventional electrostatic particle-in-cell methods are relaxed by the speed-limiting approach, thus enabling larger timesteps and faster simulations. Here, these results indicate that the SLPIC method is a fast, accurate, and powerful technique for modeling plasmas wherein electron kinetic behavior is nontrivial (such that a fluid/Boltzmann representation for electrons is inadequate) but evolution is on ion timescales.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1999

Low-Frequency Oscillations and Transport Processes Induced by Multiscale Transverse Structures in the Polar Wind Outflow: A Three-Dimensional Simulation

Multiscale transverse structures in the magnetic-field-aligned flows have been frequently observed in the auroral region by FAST and Freja satellites. A number of multiscale processes, such as broadband low-frequency oscillations and various cross-field transport effects are well correlated with these structures. To study these effects, we have used our three-dimensional multifluid model with multiscale transverse inhomogeneities in the initial velocity profile. Self-consistent-frequency mode driven by local transverse gradients in the generation of the low field-aligned ion flow and associated transport processes were simulated. Effects of particle interaction with the self-consistent time-dependent three-dimensional wave potential have been modeled using a distribution of test particles. For typical polar wind conditions it has been found that even large-scale (approximately 50 - 100 km) transverse inhomogeneities in the flow can generate low-frequency oscillations that lead to significant flow modifications, cross-field particle diffusion, and other transport effects. It has also been shown that even small-amplitude (approximately 10 - 20%) short-scale (approximately 10 km) modulations of the original large-scale flow profile significantly increases low-frequency mode generation and associated cross-field transport, not only at the local spatial scales imposed by the modulations but also on global scales. Note that this wave-induced cross-field transport is not included in any of the global numerical models of the ionosphere, ionosphere-thermosphere, or ionosphere-polar wind. The simulation results indicate that the wave-induced cross-field transport not only affects the ion outflow rates but also leads to a significant broadening of particle phase-space distribution and transverse particle diffusion.

Ganguli, Supriya B.↗

DOE OSTI2022

Velocity-dependent annihilation radiation from dark matter subhalos in cosmological simulations

Here, we use the suite of Milky Way-like galaxies in the Auriga simulations to determine the contribution to annihilation radiation from dark matter subhalos in three velocity-dependent dark matter annihilation models: Sommerfeld, p-wave, and d-wave models. We compare these to the corresponding distribution in the velocity-independent s-wave annihilation model. For both the hydrodynamical and dark-matter-only simulations, only in the case of the Sommerfeld-enhanced annihilation does the total annihilation flux from subhalos exceed the total annihilation flux from the smooth halo component within the virial radius of the halo. Progressing from Sommerfeld to the s, p, and d-wave models, the contribution from the smooth component of the halo becomes more dominant, implying that for the p-wave and d-wave models the smooth component is by far the dominant contribution to the radiation. Comparing to the Galactic center excess observed by Fermi-LAT, for all simulated halos the emission is dominated by the smooth halo contribution. However, it is possible that for Sommerfeld models, extrapolation down to mass scales below the current resolution limit of the simulation would imply a non-negligible contribution to the gamma-ray emission from the Galactic Center region.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Ion density waves driving the formation of filamentary dust structures

The PlasmaKristall-4 (PK-4) experiment on the International Space Station allows for the study of the 3-dimensional interaction between plasma and dust particles. Previous simulations of the PK-4 environment have discovered fast moving ionization waves in the dc discharge. These ionization waves vary the plasma parameters by up to an order of magnitude, which may affect the mechanisms responsible for the self-organization of chains seen in the PK-4 experiment. Here, we adapt a molecular dynamics simulation to employ temporally varying plasma conditions in order to investigate the effect on the dust charging and electrostatic potential. In order to describe the differences between the average of the plasma conditions and the time-varying plasma condition, we present a model to reproduce the interaction that takes into account the negative potential from the dust grain and the positive potential from the ion wake.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2007

Receptivity of Hypersonic Boundary Layers Due to Acoustic Disturbances over Blunt Cone

The transition process induced by the interaction of acoustic disturbances in the free-stream with boundary layers over a 5-degree straight cone and a wedge with blunt tips is numerically investigated at a free-stream Mach number of 6.0. To compute the shock and the interaction of shock with the instability waves the Navier-Stokes equations are solved in axisymmetric coordinates. The governing equations are solved using the 5th -order accurate weighted essentially non-oscillatory (WENO) scheme for space discretization and using third-order total-variation-diminishing (TVD) Runge-Kutta scheme for time integration. After the mean flow field is computed, acoustic disturbances are introduced at the outer boundary of the computational domain and unsteady simulations are performed. Generation and evolution of instability waves and the receptivity of boundary layer to slow and fast acoustic waves are investigated. The mean flow data are compared with the experimental results. The results show that the instability waves are generated near the leading edge and the non-parallel effects are stronger near the nose region for the flow over the cone than that over a wedge. It is also found that the boundary layer is much more receptive to slow acoustic wave (by almost a factor of 67) as compared to the fast wave.

Kara, K.↗

NASA NTRS2005

The Planetesimal Bow Shock Model for Chondrule Formation: More Detailed Simulations in the Near Vicinity of the Planetesimal

Gas dynamic shock waves in a low temperature nebula have been considered to be a leading candidate mechanism for providing the repetitive, short-duration heating events that are believed to have been responsible for the formation of chondrules in chondrites. It has been found, for example, that shocks with Mach numbers greater than 4 or 5 would be capable of rapidly melting 0.1-1 mm sized silicate particles as required by meteoritic data. Near the nebula midplane where chondrite parent bodies are believed to have formed, possible energy sources for generating multiple shocks include mass concentrations in a gravitationally unstable nebula, tidal interactions of proto-Jupiter with the nebula, and bow waves upstream of planetesimals scattered gravitationally into eccentric orbits by proto- Jupiter. In a recent study, we have found that chondrule precursors that are melted following passage through a planetesimal bow shock would likely cool at rates that are too rapid to be consistent with meteoritic evidence. However, that study was limited to the bowshock exterior to about 1.5 planetesimal radii (measured perpendicular to the symmetry axis) to avoid complications interior to this distance where large pressure gradients and lateral flow occur as the gas flows around the planetesimal. In this paper, we reconsider the planetesimal bow shock model and report more detailed numerical simulations of chondrule precursor heating, cooling, and dynamical histories in the near vicinity of a representative planetesimal.

Hood, Lon L.↗

DOE OSTI2020

Seismic Anisotropy of Temperate Ice in Polar Ice Sheets

Here we present a series of simple shear numerical simulations of dynamic recrystallization of two-phase nonlinear viscous materials that represent temperate ice. First, we investigate the effect of the presence of water on the resulting microstructures and, second, how water influences on P wave (V p ) and fast S wave (V s ) velocities. Regardless the water percentage, all simulations evolve from a random fabric to a vertical single maximum. For a purely solid aggregate, the highest V p quickly aligns with the maximum c-axis orientation. At the same time, the maximum c-axis development reduces V s in this orientation. When water is present, the developed maximum c-axis orientation is less intense, which results in lower V p and V s . At high percentage of water, V p does not align with the maximum c-axis orientation. If the bulk modulus of ice is assumed for the water phase (i.e., implying that water is at high pressure), we find a remarkable decrease of V s while V p remains close to the value for purely solid ice. These results suggest that the decrease in V s observed at the base of the ice sheets could be explained by the presence of water at elevated pressure, which would reside in isolated pockets at grain triple junctions. Under these conditions water would not favor sliding between ice grains. However, if we consider that deformation dominates over recrystallization, water pockets get continuously stretched, allowing water films to be located at grain boundaries. This configuration would modify and even overprint the maximum c-axis-dependent orientation and the magnitude of seismic anisotropy.

58 GEOSCIENCES↗

NASA NTRS1977

Performance of twin two-dimensional wedge nozzles including thrust vectoring and reversing effects at speeds up to Mach 2.20

Transonic tunnel and supersonic pressure tunnel tests were reformed to determine the performance characteristics of twin nonaxisymmetric or two-dimensional nozzles with fixed shrouds and variable-geometry wedges. The effects of thrust vectoring, reversing, and installation of various tails were also studied. The investigation was conducted statically and at flight speeds up to a Mach number of 2.20. The total pressure ratio of the simulated jet exhaust was varied up to approximately 26 depending on Mach number. The Reynolds number per meter varied up to 13.20 x 1 million. An analytical study was made to determine the effect on calculated wave drag by varying the mathematical model used to simulate nozzle jet-exhaust plume.

Capone, F. J.↗

NASA NTRS1991

Experiment and numerical simulation of RF heating in the Tandem Mirror plasma propulsion device

Recent results of experimental and theoretical studies of ICRF heating of the plasma in the Tendem Mirror rocket are presented. The radial and axial profile of the magnetic field of the wave have been measured, and the data agree with the results from numerical simulation. One very important new finding is that the wave damped when it approached the resonance plane. This is a strong indication of beach heating effect and that RF power is absorbed by the ions in the plasma as expected. This power absorption phenomenon was also further confirmed by the refined analytical study of the wave propagation in a slab model. The electron density in the central cell has been measured by a microwave interferometer. The existence of two types of discharge similar to those of H-alpha emission have also been observed from the density measurements. A preliminary design of a flight system for a 100-day Mars transit has been undertaken. The specific weight for a 10 MW rocket is only 0.04 kg/kW.

Yang, T. F.↗

NASA NTRS2015

Tetrahedral-Mesh Simulation of Turbulent Flows with the Space-Time Conservative Schemes

Direct numerical simulations of turbulent flows are predominantly carried out using structured, hexahedral meshes despite decades of development in unstructured mesh methods. Tetrahedral meshes offer ease of mesh generation around complex geometries and the potential of an orientation free grid that would provide un-biased small-scale dissipation and more accurate intermediate scale solutions. However, due to the lack of consistent multi-dimensional numerical formulations in conventional schemes for triangular and tetrahedral meshes at the cell interfaces, numerical issues exist when flow discontinuities or stagnation regions are present. The space-time conservative conservation element solution element (CESE) method - due to its Riemann-solver-free shock capturing capabilities, non-dissipative baseline schemes, and flux conservation in time as well as space - has the potential to more accurately simulate turbulent flows using unstructured tetrahedral meshes. To pave the way towards accurate simulation of shock/turbulent boundary-layer interaction, a series of wave and shock interaction benchmark problems that increase in complexity, are computed in this paper with triangular/tetrahedral meshes. Preliminary computations for the normal shock/turbulence interactions are carried out with a relatively coarse mesh, by direct numerical simulations standards, in order to assess other effects such as boundary conditions and the necessity of a buffer domain. The results indicate that qualitative agreement with previous studies can be obtained for flows where, strong shocks co-exist along with unsteady waves that display a broad range of scales, with a relatively compact computational domain and less stringent requirements for grid clustering near the shock. With the space-time conservation properties, stable solutions without any spurious wave reflections can be obtained without a need for buffer domains near the outflow/farfield boundaries. Computational results for the isotropic turbulent flow decay, at a relatively high turbulent Mach number, show a nicely behaved spectral decay rate for medium to high wave numbers. The high-order CESE schemes offer very robust solutions even with the presence of strong shocks or widespread shocklets. The explicit formulation in conjunction with a close to unity theoretical upper Courant number bound has the potential to offer an efficient numerical framework for general compressible turbulent flow simulations with unstructured meshes.

Chang, Chau-Lyan↗

NASA NTRS1997

A Study of Fundamental Shock Noise Mechanisms

This paper investigates two mechanisms fundamental to sound generation in shocked flows: shock motion and shock deformation. Shock motion is modeled numerically by examining the interaction of a sound wave with a shock. This numerical approach is validated by comparison with results obtained by linear theory for a small-disturbance case. Analysis of the perturbation energy with Myers' energy corollary demonstrates that acoustic energy is generated by the interaction of acoustic disturbances with shocks. This analysis suggests that shock motion generates acoustic and entropy disturbance energy. Shock deformation is modeled numerically by examining the interaction of a vortex ring with a shock. These numerical simulations demonstrate the generation of both an acoustic wave and contact surfaces. The acoustic wave spreads cylindrically. The sound intensity is highly directional and the sound pressure increases with increasing shock strength. The numerically determined relationship between the sound pressure and the Mach number is found to be consistent with experimental observations of shock noise. This consistency implies that a dominant physical process in the generation of shock noise is modeled in this study.

Meadows, Kristine R.↗

NASA NTRS1989

Oxygen acoustic solitary waves in a magnetized plasma

Ion-acoustic solitary waves in a magnetized plasma containing an arbitrary mixture of H(+) and O(+) ions are studied. A nonlinear wave equation has been derived from the Poisson-Vlasov equations, including a uniform magnetic field and dissipation due to reflected electrons. When dissipation is ignored, the equation has soliton solutions associated with both oxygen and hydrogen acoustic modes, which can be either rarefactive or compressive depending on the ion concentrations and the electron/ion temperature ratio and, more weakly, on the bulk drifts of the species. If electron reflection is included, the solitary wave can be intensified. Under somewhat restrictive conditions the oxygen solitary wave is rarefactive and propagates with a velocity comparable to that observed by the Viking satellite. The three-dimensional solitons obey a relation of scales parallel to the magnetic field and in the transverse direction. Computer simulations of one-dimensional versions of the nonlinear wave equation are presented.

Qian, S.↗

NASA NTRS1990

Project summaries

Lunar base projects, including a reconfigurable lunar cargo launcher, a thermal and micrometeorite protection system, a versatile lifting machine with robotic capabilities, a cargo transport system, the design of a road construction system for a lunar base, and the design of a device for removing lunar dust from material surfaces, are discussed. The emphasis on the Gulf of Mexico project was on the development of a computer simulation model for predicting vessel station keeping requirements. An existing code, used in predicting station keeping requirements for oil drilling platforms operating in North Shore (Alaska) waters was used as a basis for the computer simulation. Modifications were made to the existing code. The input into the model consists of satellite altimeter readings and water velocity readings from buoys stationed in the Gulf of Mexico. The satellite data consists of altimeter readings (wave height) taken during the spring of 1989. The simulation model predicts water velocity and direction, and wind velocity.

Source record↗