Engineering Papers⌕ Search

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.

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

At least 271 records · Page 15

A Comparison of Three-Dimensional Simulations of Traveling-Wave Tube Cold-Test Characteristics Using CST MICROWAVE STUDIO and MAFIA

Previously, it was shown that MAFIA (solutions of Maxwell's equations by the Finite Integration Algorithm), a three-dimensional simulation code, can be used to produce accurate cold-test characteristics including frequency-phase dispersion, interaction impedance, and attenuation for traveling-wave tube (TWT) slow-wave structures. In an effort to improve user-friendliness and simulation time, a model was developed to compute the cold-test parameters using the electromagnetic field simulation software package CST MICROWAVE STUDIO (MWS). Cold-test parameters were calculated for several slow-wave circuits including a ferruled coupled-cavity, a folded waveguide, and a novel finned-ladder circuit using both MWS and MAFIA. Comparisons indicate that MWS provides more accurate cold-test data with significantly reduced simulation times. Both MAFIA and MWS are based on the finite integration (FI) method; however, MWS has several advantages over MAFIA. First, it has a Windows based interface for PC operation, making it very user-friendly, whereas MAFIA is UNIX based. MWS uses a new Perfect Boundary Approximation (PBA), which increases the accuracy of the simulations by avoiding stair step approximations associated with MAFIA's representation of structures. Finally, MWS includes a Visual Basic for Applications (VBA) compatible macro language that enables the simulation process to be automated and allows for the optimization of user-defined goal functions, such as interaction impedance.

Chevalier, C. T.↗

Two-dimensional simulations of internal gravity waves in the radiation zones of intermediate-mass stars

Intermediate-mass main-sequence stars have large radiative envelopes overlying convective cores. This configuration allows internal gravity waves (IGWs) generated at the convective–radiative interface to propagate towards the stellar surface. The signatures of these waves can be observed in the photometric and spectroscopic data from stars. Here, we have studied the propagation of these IGWs using two-dimensional (2D) fully non-linear hydrodynamical simulations with realistic stellar reference states from the 1D stellar evolution code, Modules for Stellar Astrophysics (MESA). When a single wave is forced, we observe wave self-interaction. When two waves are forced, we observe non-linear interaction (i.e. triadic interaction) between these waves forming waves at different wavelengths and frequencies. When a spectrum of waves similar to that found in numerical simulations is forced, we find that the surface IGW frequency slope is consistent with recent observations. This power law is similar to that predicted by linear theory for the wave propagation, with small deviations that can be an effect of non-linearities. When the same generation spectrum is applied to 3 M ⊙ models at different stellar rotation and ages, the surface IGW spectrum slope is very similar to the generation spectrum slope.

79 ASTRONOMY AND ASTROPHYSICS↗

Classical-quantum simulation of non-equilibrium Marshak waves

In the radiation hydrodynamic simulations used to design inertial confinement fusion (ICF) and pulsed power experiments, nonlinear radiation diffusion tends to dominate CPU time. This raises the interesting question of whether a quantum algorithm can be found for nonlinear radiation diffusion which provides a quantum speedup. Recently, such a quantum algorithm was introduced based on a quantum algorithm for solving systems of nonlinear partial differential equations (PDEs) which provides a quadratic quantum speedup. Here, we apply this quantum PDE (QPDE) algorithm to the problem of a non-equilibrium Marshak wave propagating through a cold, semi-infinite, optically thick target, where the radiation and matter fields are not assumed to be in local thermodynamic equilibrium. The dynamics is governed by a coupled pair of nonlinear PDEs which are solved using the QPDE algorithm, as well as two standard PDE solvers: (i) Python's py-pde solver; and (ii) the KULL ICF simulation code developed at Lawrence-Livermore National Laboratory. We compare the simulation results obtained using the QPDE algorithm and the standard PDE solvers and find excellent agreement.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Hybrid Simulation of the Shock Wave Trailing the Moon

A standing shock wave behind the Moon was predicted by Michel (1967) but never observed nor simulated. We use 1D hybrid code in order to simulate the collapse of the plasma-free cavity behind the Moon and for the first time to model the formation of this shock. Starting immediately downstream of the obstacle we consider the evolution of plasma expansion into the cavity in the frame of reference moving along with the solar wind. Well-known effects as electric charging of the cavity affecting the plasma flow and counterstreaming ion beams in the wake are reproduced. Near the apex of the inner Mach cone where the plasma flows from the opposite sides of the obstacle meet, a shock wave arises. We expect the shock to be produced at periods of high electron temperature solar wind streams (T(sub i) much less than T(sub e) approximately 100 eV). The shock is produced by the interaction of oppositely directed proton beams in the plane containing solar wind velocity and interplanetary magnetic field vectors. In the direction across the magnetic field and the solar wind velocity, the shock results from the interaction of the plasma flow with the region of the enhanced magnetic field inside the cavity that plays the role of the magnetic barrier. The appearance of the standing shock wave is expected at the distance of approximately 7R(sub M) downstream of the Moon.

Moon↗

Application of the extended boundary condition method to Monte Carlo simulations of scattering of waves by two-dimensional random rough surfaces

The extended boundary condition method is applied to Monte Carlo simulations of two-dimensional random rough surface scattering. The numerical results are compared with one-dimensional random rough surfaces obtained from the finite-element method. It is found that the mean scattered intensity from two-dimensional rough surfaces differs from that of one dimension for rough surfaces with large slopes.

Tsang, L.↗

Numerical simulations of small amplitude acoustic wave propagation in a converging-diverging nozzle

Numerical simulations of the small amplitude acoustic wave propagation in a converging-diverging nozzle proposed in Category 5 are performed and presented here. The quasi 1-D unsteady flow equations in conservative form are discretized by the DRP scheme with the artificial damping developed by Tam and Webb. Characteristic boundary conditions are used at both the inlet and outlet of the nozzle. The effect of different numerical implementations of the subsonic inflow boundary conditions on the convergence of the solution to the steady state is studied. In the case of a subsonic outflow in which a shock is formed in the nozzle, the interaction between acoustics and the shock is also investigated.

Dong, Thomas Z.↗

Subcritical dispersive shock waves upstream of planetary bow shocks and at Comet Giacobini-Zinner

The nonlinear evolution of ULF (magnetosonic) waves is studied using electromagnetic simulations. The waves were observed upstream of the planetary bow shocks at Comet Giacobini-Zinner. It is shown that as ULF waves generated by the resonant ion beam instability steepen, their polarization becomes linear. A high-frequency wave packet is generated by the steepening process. The steepened wave has a structure and a behavior similar to a subcritical dispersive shock. It is suggested that Comet Giacobini-Zinner did not have a single bow shock. It is concluded that the comet had a transition region consisting of a series of shocks which were convected by the solar wind.

Omidi, N.↗

Two-dimensional simulations of internal gravity waves in a $\mathrm{5M_⊙}$ zero-age-main-sequence model

Main-sequence intermediate-mass stars present a radiative envelope that supports internal gravity waves (IGWs). Excited at the boundary with the convective core, IGWs propagate towards the stellar surface and are suspected to impact physical processes such as rotation and chemical mixing. Using the fully compressible time-implicit code MUSIC, we study IGWs in two-dimensional simulations of a zero-age-main-sequence 5 solar mass star model up to 91 percent of the stellar radius with different luminosity and radiative diffusivity enhancements. Our results show that low-frequency waves excited by core convection are strongly impacted by radiative effects as they propagate. This impact depends on the radial profile of radiative diffusivity which increases by almost 5 orders of magnitude between the centre of the star and the top of the simulation domain. In the upper layers of the simulation domain, we observe an increase of the temperature. Our study suggests that this is due to heat added in these layers by IGWs damped by radiative diffusion. We show that non-linear effects linked to large amplitude IGWs may be relevant just above the convective core. Both these effects are intensified by the artificial enhancement of the luminosity and radiative diffusivity, with enhancement factors up to 10 4 times the realistic values. Our results also highlight that direct comparison between numerical simulations with enhanced luminosity and observations must be made with caution. Finally, our work suggests that thermal effects linked to the damping of IGWs could have a non-negligible impact on stellar structure.

79 ASTRONOMY AND ASTROPHYSICS↗

Quantum simulation of molecules without fermionic encoding of the wave function

Abstract Molecular simulations generally require fermionic encoding in which fermion statistics are encoded into the qubit representation of the wave function. Recent calculations suggest that fermionic encoding of the wave function can be bypassed, leading to more efficient quantum computations. Here we show that the two-electron reduced density matrix (2-RDM) can be expressed as a unique functional of the unencoded N -qubit-particle wave function without approximation, and hence, the energy can be expressed as a functional of the 2-RDM without fermionic encoding of the wave function. In contrast to current hardware-efficient methods, the derived functional has a unique, one-to-one (and onto) mapping between the qubit-particle wave functions and 2-RDMs, which avoids the over-parametrization that can lead to optimization difficulties such as barren plateaus. An application to computing the ground-state energy and 2-RDM of H 4 is presented.

74 ATOMIC AND MOLECULAR PHYSICS↗

Simulation of glancing shock wave and boundary layer interaction

Shock waves generated by sharp fins, glancing across a laminar boundary layer growing over a flat plate, are simulated numerically. Several basic issues concerning the resultant three-dimensional flow separation are studied. Using the same number of grid points, different grid spacings are employed to investigate the effects of grid resolution on the origin of the line of separation. Various shock strengths (generated by different fin angles) are used to study the so-called separated and unseparated boundary layer and to establish the existence or absence of the secondary separation. The usual interpretations of the flow field from previous studies and new interpretations arising from the present simulation are discussed.

Hung, Ching-Mao↗

Simulation of glancing shock wave and boundary layer interaction

Shock waves generated by sharp fins, glancing across a laminar boundary layer growing over a flat plate, are simulated numerically. Several basic issues concerning the resultant three-dimensional flow separation are studied. Using the same number of grid points, different grid spacings are employed to investigate the effects of grid resolution on the origin of the line of separation. Various shock strengths (generated by different fin angles) are used to study the so-called separated and unseparated boundary layer and to establish the existence or absence of the secondary separation. The usual interpretations of the flow field from previous studies and new interpretations arising from the present simulation are discussed.

Hung, Ching-Mao↗

The Atmospheric Effects of HSCT Emissions Simulated by a 3-Wave Interactive Model

An interactive model which couples a semi-spectral dynamical model, a radiative transfer code and a two-dimensional chemistry transport model (2-D CTM), is used to assess the atmospheric effects of the High-Speed Civil Transport (HSCT) engine emissions. The residual mean meridional circulation, the zonal-mean temperature and the eddy diffusion coefficients are calculated using zonal means and three longest zonal waves of dynamical variables integrated in the semi-spectral dynamical model. They are used in the 2-D CTM to simulate the distribution of trace gases in the atmosphere. The simulated ozone is sent to the radiative transfer code to calculate the heating rates, which drive the dynamics. This radiative coupling connects the dynamical and photochemical processes and creates feedback when the atmosphere is perturbed. It is found that in most areas the ozone depletion caused by HSCT emissions calculated using the 3-wave model has the features similar to, but with significantly larger magnitude than that calculated by the AER 2-D CTM with prescribed transport parameters and temperature. The difference is mostly due to the differences in the circulation in the two models. The radiative feedback effects are investigated by comparing the ozone depletion calculated with the baseline dynamics and with the dynamics perturbed by the HSCT emissions. The feedback through changes in the residual mean meridional circulation and the eddy diffusion coefficients has moderate effects on the simulated ozone depletion. It reduced the ozone depletion by 20-30% in northern mid and high-latitudes. However, the feedback through changes in the zonal- mean temperature is negligible.

Shia, R.-L.↗

Molecular dynamics simulation of a piston driven shock wave in a hard sphere gas

Molecular dynamics simulation is used to study the piston driven shock wave at Mach 1.5, 3, and 10. A shock tube, whose shape is a circular cylinder, is filled with hard sphere molecules having a Maxwellian thermal velocity distribution and zero mean velocity. The piston moves and a shock wave is generated. All collisions are specular, including those between the molecules and the computational boundaries, so that the shock development is entirely causal, with no imposed statistics. The structure of the generated shock is examined in detail, and the wave speed; profiles of density, velocity, and temperature; and shock thickness are determined. The results are compared with published results of other methods, especially the direct simulation Monte-Carlo method. Property profiles are similar to those generated by direct simulation Monte-Carlo method. The shock wave thicknesses are smaller than the direct simulation Monte-Carlo results, but larger than those of the other methods. Simulation of a shock wave, which is one-dimensional, is a severe test of the molecular dynamics method, which is always three-dimensional. A major challenge of the thesis is to examine the capability of the molecular dynamics methods by choosing a difficult task.

Woo, Myeung-Jouh↗

Similarity rules in gravity jitter-related spacecraft liquid propellant slosh waves excitation

The dynamical behavior of fluids, in particular the effect of surface tension on partially filled rotating fluids in a full-scale prototype Gravity Probe-B Spacecraft propellant tank and various 10 percent subscale containers with identical values of similarity parameters such as Bond number, dynamical capillary number, rotational Reynolds number, and Weber number, as well as imposed gravity jitters have been investigated. It is shown that the Bond number can be used to simulate the wave characteristics of slosh wave excitation, whereas the Weber number can be used to simulate the wave amplitude of slosh-mode excitation. It is shown that a dynamical capillary number can be used to simulate the induced perturbation of the fluid stress distribution exerted on the wall. This distribution is governed by the interaction between surface tension (slosh-wave excitation along the liquid-vapor interface) and viscous (fluid stress exerted on the wall) forces.

Hung, R. J.↗

Reconstructing the exit wave of 2D materials in high-resolution transmission electron microscopy using machine learning

Reconstruction of the exit wave function is an important route to interpreting high-resolution transmission electron microscopy (HRTEM) images. Here we demonstrate that convolutional neural networks can be used to reconstruct the exit wave from a short focal series of HRTEM images, with a fidelity comparable to conventional exit wave reconstruction. We use a fully convolutional neural network based on the U-Net architecture, and demonstrate that we can train it on simulated exit waves and simulated HRTEM images of graphene-supported molybdenum disulphide (an industrial desulfurization catalyst). We then apply the trained network to analyse experimentally obtained images from similar samples, and obtain exit waves that clearly show the atomically resolved structure of both the MoS 2 nanoparticles and the graphene support. We also show that it is possible to successfully train the neural networks to reconstruct exit waves for 3400 different two-dimensional materials taken from the Computational 2D Materials Database of known and proposed two-dimensional materials.

2D materials↗

Numerical simulation of particle-wave interaction in boundary layers

The effects of wall injection and particle motion on the spatial stability of two-dimensional plane channel flow are investigated. For this purpose, an accurate Navier-Stokes solver to simulate the space-time evolution of disturbances in three-dimensional flows has been developed. The code is operational on the NASA Langley CRAY2 and can be ported to any other supercomputer. The code has been tested extensively in tracking the spatial evolution of two-dimensional disturbances in plane channel flow and provided excellent agreement with the linear theory including at the inflow/outflow boundaries. Preliminary calculations have been performed to investigate the effects of stationary and moving sources of vortical disturbances simulating a particle traveling in the flow field. Results suggest that even at very low amplitudes, vortical disturbances act as amplifiers on the Tollmien-Schlichting waves promoting rapid instability. It is also found that slow moving particles are more dangerous than both stationary and fast moving particles for the same disturbance levels.

Biringen, S.↗

Simulating Responses of Gravitational-Wave Instrumentation

Synthetic LISA is a computer program for simulating the responses of the instrumentation of the NASA/ESA Laser Interferometer Space Antenna (LISA) mission, the purpose of which is to detect and study gravitational waves. Synthetic LISA generates synthetic time series of the LISA fundamental noises, as filtered through all the time-delay-interferometry (TDI) observables. (TDI is a method of canceling phase noise in temporally varying unequal-arm interferometers.) Synthetic LISA provides a streamlined module to compute the TDI responses to gravitational waves, according to a full model of TDI (including the motion of the LISA array and the temporal and directional dependence of the arm lengths). Synthetic LISA is written in the C++ programming language as a modular package that accommodates the addition of code for specific gravitational wave sources or for new noise models. In addition, time series for waves and noises can be easily loaded from disk storage or electronic memory. The package includes a Python-language interface for easy, interactive steering and scripting. Through Python, Synthetic LISA can read and write data files in Flexible Image Transport System (FITS), which is a commonly used astronomical data format.

Armstrong, John↗