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 1,009 records · Page 56

Simulating nonlinear optical processes on a superconducting quantum device

Simulating plasma physics on quantum computers is difficult because most problems of interest are nonlinear, but quantum computers are not naturally suitable for nonlinear operations. In weakly nonlinear regimes, plasma problems can be modelled as wave–wave interactions. In this paper, we develop a quantization approach to convert nonlinear wave–wave interaction problems to Hamiltonian simulation problems. We demonstrate our approach using two qubits on a superconducting device. Unlike a photonic device, a superconducting device does not naturally have the desired interactions in its native Hamiltonian. Nevertheless, Hamiltonian simulations can still be performed by decomposing required unitary operations into native gates. To improve experimental results, we employ a range of error-mitigation techniques. Apart from readout error mitigation, we use randomized compilation to transform undiagnosed coherent errors into well-behaved stochastic Pauli channels. Moreover, to compensate for stochastic noise, we rescale exponentially decaying probability amplitudes using rates measured from cycle benchmarking. We carefully consider how different choices of product-formula algorithms affect the overall error and show how a trade-off can be made to best utilize limited quantum resources. This study provides an example of how plasma problems may be solved on near-term quantum computing platforms.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

3D orthorhombic earth model effects on seismic source characterization

Most earth materials are anisotropic with regard to seismic wave-speeds, especially materials such as shales, or where oriented fractures are present. However, the base assumption for many numerical simulations is to treat earth materials as isotropic media. This is done for simplicity, the apparent weakness of anisotropy in the far field, and the lack of well-characterized anisotropic material properties for input into numerical simulations. One approach for addressing the higher complexity of actual geologic regions is to model the material as an orthorhombic medium. We have developed an explicit time-domain, finite-difference (FD) algorithm for simulating three-dimensional (3D) elastic wave propagation in a heterogeneous orthorhombic medium. The objective of this research is to investigate the errors and biases that result from modeling a non-isotropic medium as an isotropic medium. This is done by computing “observed data” by using synthetic, anisotropic simulations with the assumption of an orthorhombic, anisotropic earth model. Green’s functions for an assumed isotropic earth model are computed and then used an inversion designed to estimate moment tensors with the “observed” data. One specific area of interest is how shear waves, which are introduced in an anisotropic model even for an isotropic explosion, affect the characterization of seismic sources when isotropic earth assumptions are made. This work is done in support of the modeling component of the Source Physics Experiment (SPE), a series of underground chemical explosions at the Nevada National Security Site (NNSS).

58 GEOSCIENCES↗

Hybrid simulations of the saturated Farley-Buneman instability in the ionosphere

Numerical simulations of the Farley-Buneman instability in 2-1/2 dimensions using particle ions and fluid electrons show the growth, saturation and nonlinear behavior of two-stream waves. This hybrid technique models the saturated state of the instability for a much longer period of time than the pure particle codes that preceded it. While focusing principally on modeling the topside E region equatorial electrojet, many of these results apply to the auroral two-stream instability as well. The following features are seen in all our hybrid simulations: (1) wave growth at an angle offset from the electron drift direction where the angle depends on the strength of the driving electric field, (2) nonlinear coupling to waves traveling perpendicular to the propagation direction of the principal two-stream waves, (3) a saturated wave phase velocity at or above the sound speed but well below the velocity predicted by linear theory and (4) phase velocities which remain almost constant as a simulated radar sweeps from a horizontal direction to nearly vertical. The nonlinear electron motion dominates the behavior of these waves. Further, these simulations indicate that ion kinetic effects are not essential for the saturation of the instability and that electron temperature effects have a minor impact on the final saturated state.

Oppenheim, Meers↗

Rough surface scattering from moving ocean surfaces as an indicator of scattering mechanisms

Ever since Crombie used the Doppler spectrum of RF signals scattered from the ocean surface to provide conclusive evidence for Bragg resonant scattering, scattering from the ocean surface has served as a testing ground for the identification of governing scattering mechanisms in rough surface scattering. The moving ocean surface has the advantage for rough surface scattering that, being a dispersive medium, features with different spatial scales the tendency to travel at different speeds. Thus, by looking at the temporal signature of the scattered field, one can try to isolate the scattering mechanism. This feature has been exploited by various authors. To give a recent example, we cite the work of Plant, who examined the Doppler spectrum of the return signal to infer the presence of two-scale scattering in scattering at medium incidence angles. The problem with restricting oneself to the examination of the Doppler spectrum, as has been done traditionally, is that one restricts the phenomena under investigation to temporally homogeneous phenomena: transient phenomena, which are localized in time, are not localized in frequency space. There is increasing evidence that this type of phenomena, due to specular scattering or breaking waves, may also play an important role in determining the scattering mechanism in ocean-like surfaces. To overcome this problem, we introduce the use of the Wavelet Transform to study the frequency-temporal signatures of the scattered field from moving ocean waves. We calculate this signature using various analytic scattering theories and show that the Wavelet Transform provides a useful tool for separating the different scattering mechanisms operating in scattering in ocean-like surfaces. Next we simulate realistic nonlinear moving ocean waves and calculate the temporal scattered field signature by using the method of moments and the stop-start approximation: the surface is assumed stationary during a scattering event, but moves between incident EM pulses. We compute and compare the Doppler frequency spectrum and the Wavelet Transform of the resulting signal for various types of surfaces and show how different scattering mechanisms may be isolated by this technique. Finally, we compare these numerical results against the analytical predictions in order to obtain the domain of applicability of various scattering theories for moving ocean surfaces which may include wave breaking and specular glints.

Rodriguez, Ernesto↗

Modification of Mach 6 Freestream Flow by Pitot Probe Bow Shock as Measured by Line FLDI and SAFS

Measurements of the flow in front of flat-face Pitot probes in NASA Langley ResearchCenter’s 20-Inch Mach 6 Air Tunnel were made using a 35-point line focused laser differentialinterferometer (FLDI), a two-point FLDI, and a self-aligned focusing schlieren system. Thesesystems were used to make measurements of flow unsteadiness between the bow shock andPitot probe face, while also probing the freestream region just upstream of the bow shock.Shock standoff distances from the non-intrusive measurements match closely with analyticalsolutions using the probe face diameter. Acoustic wave angles and velocities in the freestreamcompare favorably with DNS simulations of a Mach 6 freestream flow, and wave velocities in thepost-shock region slow as they approach the probe face. Peak frequencies in the post-shockregion as measured by FLDI compare favorably to the estimates computed from the shockstandoff distance for most probe diameters. With increasing unit Reynolds number, fluctuationintensities in the post-shock region were attenuated below the intensities of the freestream, witha spectral, spatial, and probe size dependence. Pitot probe power spectral density estimatesare presented from the test campaign, which can be compared with FLDI and SAFS data toevaluate the applicability of the transfer functions derived from computational simulations.

Joshua M Weisberger↗

Stationary Planetary Waves in Northern Late Winter: MRO/MARCI Observations and Mars Climate Model Simulations

Mars reveals similar, yet also rather different, atmospheric circulation patterns compared to those on Earth. In both atmospheres, solar differential heating drives global Hadley circulation cells. However during solstice on Mars, its Hadley cells are hemispherically asymmetric: an intense, deep, cross-hemisphere single cell dominates with rising motion in the summer hemisphere and sinking motion in the winter hemisphere. Both planets also exhibit thermally indirect (i.e., eddy-driven) Ferrel circulation cells in middle and high latitudes. In addition, Earth and Mars exhibit distinctive large-scale orography and, in a broadly defined context, continentality. For Mars’ northern midlatitudes, Tharsis in the western hemisphere, and Arabia Terra and Elysium in the eastern hemisphere, are the primary large-scale topographic features. In the southern-midlatitudes, Tharsis and Argyre in the western hemisphere, and Hellas in the eastern hemisphere are the key topographic features which can influence large-scale circulation patterns. Such underlying orographic complexes not only cause significant latitudinal excursions of the seasonal mean westerly circumnavigating polar vortex but also significantly modulate the intensity and preferred geographic regions of traveling baroclinic weather systems.

Hollingsworth, J. L.↗

Nonlinear 3-D evolution of bounded kinetic Alfven waves due to shear flow and collisionless tearing instability

The three-dimensional nonlinear simulation of the evolution of a kinetic Alfven wave is presented for simple boundary conditions that resemble those of a discrete auroral breakup arc. The model used is a variant of the reduced megnetohydrodynamic description that includes, in addition, the important dispersive effect of electron inertia. It is found that with a 'generator' boundary condition corresponding to a current generator and a fully conducting 'load' boundary, the wave evolves three-dimensionally through a combination of shear flow and collisionless tearing mode instability.

Seyler, Charles E., Jr.↗

Incorporating the Effect of Large‐Scale Vertical Motion on Convection Through Convective Mass Flux Adjustment in E3SMv2

Abstract Recent observational studies suggest that the large‐scale dynamical forcing (vertical motion) plays important roles in deep convection development. In this study we propose a convective mass flux adjustment (MAdj) approach to represent the dynamical effects of large‐scale vertical motion on convection in the Department of Energy's Energy Exascale Earth System Model version 2 (E3SMv2). With MAdj, convection is enhanced (suppressed) when there is large‐scale ascending (descending) motion at the planetary boundary layer top. The coupling of convection with large‐scale circulation significantly improves the simulation of climate variability in E3SMv2 across multiple scales from the diurnal cycle, convectively coupled equatorial waves, to the Madden‐Julian Oscillation (MJO). The standard E3SMv2 tends to simulate overly weak diurnal amplitude of precipitation and overly weak variance of convectively coupled equatorial Kelvin and westward inertio‐gravity (WIG) waves. It also fails to simulate the essential characteristics of the MJO: continuous eastward propagation. With MAdj, the amplitude of diurnal cycle of precipitation is systematically increased and its probability density distribution is much closer to observations. The MAdj can also simulate more realistic eastward propagation of the MJO and much stronger convectively coupled Kelvin and WIG waves. Moreover, the MAdj approach slightly improves the climatology simulations in precipitation, cloud, radiation, circulation, temperature, and moisture fields, with overall root‐mean‐square error (RMSE) of major climatological fields reduced by about 2%. The MAdj approach suppresses excessive grid‐scale precipitation, reducing precipitation wet biases over South China Sea, Philippine Sea, Himalayas, and South Pacific Convergence Zone in western Pacific in summer.

54 ENVIRONMENTAL SCIENCES↗

Demonstration of a new unstructured mesh IMC x-ray transport capability in LAP codes

The Advanced Simulation and Computing (ASC) Transport project’s Jayenne Implicit Monte Carlo (IMC) transport library now includes an unstructured mesh capability and is available in a Lagrangian Applications Project (LAP) code. In this presentation, we discuss recent work by the LAP and Transport projects that provides an IMC transport capability for radiation hydrodynamics in the Lagrangian frame. Verification problems testing the new capabilities have been simulated and analyzed, i.e. Marshak wave, Mach 45, Su-Olsen, picket fence, and crooked pipe, both in one and two dimensions. We also present results on two stretch goal problems: a simplified COAX high energy density physics experiment and a supernova shock simulation. Finally, we identify current limitations and future work needed to bring a full capability to the user community.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Theory and Design Tools For Studies of Reactions to Abrupt Changes in Noise Exposure

Study plans, a pre-tested questionnaire, a sample design evaluation tool, a community publicity monitoring plan, and a theoretical framework have been developed to support combined social/acoustical surveys of residents' reactions to an abrupt change in environmental noise, Secondary analyses of more than 20 previous surveys provide estimates of three parameters of a study simulation model; within individual variability, between study wave variability, and between neighborhood variability in response to community noise. The simulation model predicts the precision of the results from social surveys of reactions to noise, including changes in noise. When the study simulation model analyzed the population distribution, noise exposure environments and feasible noise measurement program at a proposed noise change survey site, it was concluded that the site could not yield sufficient precise estimates of human reaction model to justify conducting a survey. Additional secondary analyses determined that noise reactions are affected by the season of the social survey.

Fields, James M.↗

Are Ring Current Ions Lost in Electromagnetic Ion Cyclotron Wave Dispersion Relation?

Electromagnetic ion cyclotron (EMIC) waves are widely observed in the inner and outer magnetosphere, at geostationary orbit, at high latitudes along the plasmapause, and at the ionospheric altitudes. Interaction of the Ring Current (RC) ions and EMIC waves causes ion scattering into the loss cone and leads to decay of the RC, especially during the main phase of storms when the RC decay times of about one hour or less are observed. The oblique EMIC waves damp due to Landau resonance with the thermal plasmaspheric electrons, and subsequent transport of the dissipating wave energy into the ionosphere below causes an ionosphere temperature enhancement. Induced scattering of these waves by the plasmaspheric thermal ions leads to ion temperature enhancement, and forms a so-called hot zone near the plasmapause where the temperature of core plasma ions can reach tens of thousands of degrees. Relativistic electrons in the outer radiation belt also interact well with the EMIC waves, and during the main and/or recovery phases of the storms these electrons can easily be scattered into the loss cone over a time scale from several hours to a day. The plasma density distribution in the magnetosphere and the ion content play a critical role in EMIC wave generation and propagation, but the wave dispersion relation in the known RC-EMIC wave interaction models is assumed to be determined by the thermal plasma distribution only. In these models, the modification of the EMIC wave dispersion relation caused by the RC ions is not taken into account, and the RC ions are only treated as a source of free energy in order to generate EMIC waves. At the same time, the RC ions can dominate the thermal magnetospheric content in the night MLT sector at great L shells during the main and/or recovery storm phase. In this study, using our self-consistent RC-EMIC wave model [Khazanov et al., 2006], we simulate the May 1998 storm in order to quantify the global EMIC wave redistribution caused by taking into account the RC ions in the EMIC wave dispersion relation. The dramatic wave pattern redistribution is observed in the postdusk-predawn MLT sector (night sector) for L greater than 5. We found the intense EMIC waves (about a few nT) there during the main and early recovery phases of the storm. The observed wave generation in this sector is caused by taking into account the EMIC wave dispersion change due to the RC ions. There are no waves at these locations in our model if the RC ions are taken into account in the wave growth rate only, and the wave dispersion relation is only governed by the thermal plasmaspheric model.

Khazanov, G. V.↗

Effect of Ring Current Ions on Electromagnetic Ion Cyclotron Wave Dispersion Relation

Electromagnetic ion cyclotron (EMIC) waves are widely observed in the inner and outer magnetosphere, at geostationary orbit, at high latitudes along the plasmapause, and at the ionospheric altitudes. Interaction of the Ring Current (RC) ions and EMIC waves causes ion scattering into the loss cone and leads to decay of the RC, especially during the main phase of storms when the RC decay times of about one hour or less are observed. The oblique EMIC waves damp due to Landau resonance with the thermal plasmaspheric electrons, and subsequent transport of the dissipating wave energy into the ionosphere below causes an ionosphere temperature enhancement. Induced scattering of these waves by the plasmaspheric thermal ions leads to ion temperature enhancement, and forms a so-called hot zone near the plasmapause where the temperature of core plasma ions can reach tens of thousands of degrees. Relativistic electrons in the outer radiation belt also interact well with the EMIC waves, and during the main and/or recovery phases of the storms these electrons can easily be scattered into the loss cone over a time scale from several hours to a day. The plasma density distribution in the magnetosphere and the ion content play a critical role in EMIC wave generation and propagation, but the wave dispersion relation in the known RC-EMIC wave interaction models is assumed to be determined by the thermal plasma distribution only. In these models, the modification of the EMIC wave dispersion relation caused by the RC ions is not taken into account, and the RC ions are only treated as a source of free energy in order to generate EMIC waves. At the same time, the RC ions can dominate the thermal magnetospheric content in the night MLT sector at great L shells during the main and/or recovery storm phase. In this study, using our self-consistent RC-EMIC wave model [Khazanov et al., 2006], we simulate the May 1998 storm in order to quantify the global EMIC wave redistribution caused by taking into account the RC ions in the EMIC wave dispersion relation. The dramatic wave pattern redistribution is observed in the postdusk-predawn MLT sector (night sector) for L greater than 5. We found the intense EMIC waves (about a few nT) there during the main and early recovery phases of the storm. The observed wave generation in this sector is caused by taking into account the EMIC wave dispersion change due to the RC ions. There are no waves at these locations in our model if the RC ions are taken into account in the wave growth rate only, and the wave dispersion relation is only governed by the thermal plasmaspheric model.

Gamayunov, K. V.↗

Performance characterization and modeling of an oscillating surge wave energy converter

Abstract Testing wave energy converters in the ocean could be expensive and complex, which necessitates the use of numerical modeling. However, accurately modeling the response of wave energy converters with high-fidelity simulations can be computationally intensive in the design stage where different configurations must be considered. Reduced-order models based on simplified equations of motion can be very useful in the design, optimization, or control of wave energy converters. Given the complex dynamics of wave energy converters, accurate representation, and evaluation of relative contributions by different forces are required. This effort is concerned with a performance characterization of the hydrodynamic response of an oscillating surge wave energy converter that is based on a reduced-order model. A state-space model is used to represent the radiation damping term. Morison’s representation of unsteady forces is used to account for the nonlinear damping. Wave tank tests are performed to validate simulations. A free response simulation is used to determine the coefficients of the state-space model. Torque-forced simulations are used to identify the coefficients of the nonlinear damping term for different amplitudes and wave frequencies. The impact of varying these coefficients on the response is investigated. An assessment of the capability of the model in predicting the hydrodynamic response under irregular forcing is performed. The results show that the maximum error is 3% when compared with high-fidelity simulations. It is determined that the nonlinear damping is proportional to the torque amplitude and its effects are more pronounced as the amplitude of the flap oscillations increases.

16 TIDAL AND WAVE POWER↗

Simulations of Arctic ozone depletion with current and doubled levels of CO2

Results from idealized 3-D simulations of a dynamical-radiative-photochemical model of the stratosphere are presented for the Northern Hemisphere winter and spring. For a simulation of a quiescent winter, it is found that with current levels of CO2 only modest polar ozone depletion occurs, consistent with observations. For a second simulation with the same planetary wave amplitudes in the upper troposphere but with doubled CO2, the model predicts a northern hemisphere ozone hole comparable to that observed in Antarctica with almost complete ozone destruction at 20 km. Reasons for the marked difference between the simulations are identified.

Butchart, Neal↗

MESSENGER Orbital Observations of Large-Amplitude Kelvin-Helmholtz Waves at Mercury's Magnetopause

We present a survey of Kelvi\ n-Helmholtz (KH) waves at Mercury's magnetopause during MESSENGER's first Mercury year in orb it. The waves were identified on the basis of the well-established sawtooth wave signatures that are associated with non-linear KH vortices at the magnetopause. MESSENGER frequently observed such KH waves in the dayside region of the magnetosphere where the magnetosheath flow velocity is still sub -sonic, which implies that instability growth rates at Mercury's magnetopau~ are much larger than at Earth. We attribute these greater rates to the limited wave energy dissipation in Mercury's highly resistive regolith. The wave amplitude was often on the order of ' 00 nT or more, and the wave periods were - 10- 20 s. A clear dawn-dusk asymmetry is present in the data, in that all of the observed wave events occurred in the post-noon and dusk-side sectors of the magnetopause. This asymmetry is like ly related to finite Larmor-radius effects and is in agreement with results from particle-in-cell simulations of the instability. The waves were observed almost exclusively during periods when the north-south component of the magnetosheath magnetic field was northward, a pattern similar to that for most terrestrial KH wave events. Accompanying plasma measurements show that the waves were associated with the transport of magnetosheath plasma into the magnetosphere.

Sundberg, Torbjorn↗

Emission of electromagnetic waves as a stopping mechanism for nonlinear collisionless ionization waves in a high-β regime

A high energy density plasma embedded in a neutral gas is able to launch an outward-propagating nonlinear electrostatic ionization wave that traps energetic electrons. The trapping maintains a strong sheath electric field, enabling rapid and long-lasting wave propagation aided by field ionization. Here, using 1D3V kinetic simulations, we examine the propagation of the ionization wave in the presence of a transverse MG-level magnetic field with the objective to identify qualitative changes in a regime where the initial thermal pressure of the plasma exceeds the pressure of the magnetic field (Beta > 1). Our key finding is that the magnetic field stops the propagation by causing the energetic electrons sustaining the wave to lose their energy by emitting an electromagnetic wave. The emission is accompanied by the magnetic field expulsion from the plasma and an increased electron loss from the trapping wave structure. The described effect provides a mechanism mitigating rapid plasma expansion for those applications that involve an embedded plasma, such as high-flux neutron production from laser-irradiated deuterium gas jets.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Design of a Two-Body Wave Energy Converter Featuring Controllable Geometry

While the field of wave energy has been the subject of numerical simulation, scale model testing, and precommercial project testing for decades, wave energy technologies remain in the early stages of development and must continue to prove themselves as a promising modern renewable energy field. A wave energy converter (WEC) concept currently being explored is the variable-geometry WEC (VGWEC), which aims to add an extra control option to WEC design. VGWECs attempt to incorporate controllable geometric features to adjust the floating body hydrodynamics to favor either power absorption, load shedding, or other operational goals. These variable geometry components have been proposed to be controlled on a sea-state-to-sea-state or wave-to-wave time scale depending on the force (or torque) and bandwidth limitations of the actuators required to manipulate just the controllable geometric hull features. Having control over both the WEC geometry components and the power takeoff (PTO) offers the potential to improve overall system performance and reliability if a cost-effective solution can be found for a given WEC architecture. This paper will present the recent developments and results of a VGWEC concept that incorporates variable-geometry modules into a two-body WEC. In the proposed VGWEC concept, the variable-geometry modules consist of air-inflatable bags in the surface float and a water inflatable ring in the subsurface body. The surface float is tethered directly to the subsurface body through tether lines, each connected to a separate PTO. Adjusting the geometry of both the surface and subsurface bodies along with the PTO coefficients can maximize power in design sea states while reducing motion response and PTO forces when transitioning to sea states where rated power is reached and load shedding is prioritized. The ability to transition between operating condition is expected to increase the sea state operational map and power capacity.

geometry control↗