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,747 records · Page 97
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Life Cycle Cost Analysis of Prestressed Concrete Poles Subjected to Wind, Surges, and Waves

Prestressed concrete (PC) poles are becoming popular choices to support coastal power transmission systems. However, the existing literature does not offer a detailed analysis of the effectiveness of PC poles in terms of long-term vulnerabilities and the direct and indirect costs. This is due to (1) lack of fragility models for PC poles and (2) lack of probabilistic wind, storm surge, and wave models in coastal settings. In this study, we address these gaps through a series of Monte Carlo simulations to estimate fragility of PC poles as a function of age and hazard (wind, surges, and waves) intensity, and the development of a probabilistic hazard model based on 10,000 years of synthetic tropical cyclone data. The probabilistic hazard model is used in conjunction with high-resolution hydrodynamic models to generate realizations of coastal wind, storm surges, and waves for the Louisiana and Mississippi coasts. A comprehensive life cycle cost analysis for a service life of 70 years considering direct and indirect losses is conducted to compare the performance of a transmission line located in Pascagoula, Mississippi, when wood poles are replaced by PC poles. Results showed that aging has a minor effect on the reliability of PC poles, highlighting the advantages of replacing wood poles with PC poles, especially in coastal areas. In addition, PC poles are significantly more cost-effective compared with wood poles over their life cycle, leading to an estimated saving of $11.55 million (68.17% reduction). The results of this study provide key insight to inform decision-making processes to keep the coastal power grids resilient and cost-effective against future storm hazards.

natural disasters↗
NASA NTRS1970

The Darwin model as a tool for electromagnetic plasma simulation

The Darwin model of electromagnetic interaction is presented as a self-consistent theory, and is shown to be an excellent approximation to the Maxwell theory for slow electromagnetic waves. Since the fast waves of the Maxwell theory are absent, it is convenient for use in the computer simulation of the electromagnetic dynamics of nonrelativistic plasma.

Kaufman, A. N.↗
NASA NTRS2001

Terahertz-Regime, Micro-VEDs: Evaluation of Micromachined TWT Conceptual Designs

Summary form only given. The Terahertz (THz) region of the electromagnetic spectrum (approx.300-3000 GHz) has enormous potential for high-data-rate communications, spectroscopy, astronomy, space research, medicine, biology, surveillance, remote sensing, industrial process control, etc. The most critical roadblock to full exploitation of the THz band is lack of coherent radiation sources that are powerful (0.01-10.0 W continuous wave), efficient (>1 %), frequency agile (instantaneously tunable over 1% bandwidths or more), reliable, and relatively inexpensive. Micro-machined Vacuum Electron Devices (micro-VEDs) represent a promising solution. We describe prospects for miniature, THz-regime TWTs fabricated using micromachining techniques. Several approx.600 GHz conceptual designs are compared. Their expected performance has been analyzed using SD, 2.51), and 3D TWT codes. A folded waveguide (FWG) TWT forward-wave amplifier design is presented based on a Northrop Grumman (NGC) optimized design procedure. This conceptual device is compared to the simulated performance of a novel, micro-VED helix TWT. Conceptual FWG TWT backward-wave amplifiers and oscillators are also discussed. A scaled (100 GHz) FWG TWT operating at a relatively low voltage (-12 kV) is under development at NGC. Also, actual-size micromachining experiments are planned to evaluate the feasibility of arrays of micro-VED TWTs. Progress and results of these efforts are described. This work was supported, in part by AFOSR, ONR, and NSF.

Booske, John H.↗
NASA NTRS1992

Wave equations and computational models for sonic boom propagation through a turbulent atmosphere

The improved simulation of sonic boom propagation through the real atmosphere requires greater understanding of how the transient acoustic pulses popularly termed sonic booms are affected by atmospheric turbulence. A nonlinear partial differential equation that can be used to simulate the effects of smaller-scale atmospheric turbulence on sonic boom waveforms is described. The equation is first order in the time derivative and involves an extension of geometrical acoustics to include diffraction phenomena. Various terms in the equation are explained in physical terms. Such terms include those representing convection at the wave speed, diffraction, molecular relaxation, classical dissipation, and nonlinear steepening. The atmospheric turbulence enters through an effective sound speed, which varies with all three spatial coordinates, and which is the sum of the local sound speed and the component of the turbulent flow velocity projected along a central ray that connects the aircraft trajectory with the listener.

Pierce, Allan D.↗
DOE OSTI2025

Absorbing boundary conditions in material point method adopting perfectly matched layer theory

This study focuses on solving the numerical challenges of imposing absorbing boundary conditions for dynamic simulations in the material point method (MPM). To attenuate elastic waves leaving the computational domain, the current work integrates the Perfectly Matched Layer (PML) theory into the implicit MPM framework. The proposed approach introduces absorbing particles surrounding the computational domain that efficiently absorb outgoing waves and reduce reflections, allowing for accurate modeling of wave propagation and its further impact on geotechnical slope stability analysis. The study also includes several benchmark tests to validate the effectiveness of the proposed method, such as several types of impulse loading and symmetric and asymmetric base shaking. The conducted numerical tests also demonstrate the ability to handle large deformation problems, including the failure of elasto-plastic soils under gravity and dynamic excitations. The findings extend the capability of MPM in simulating continuous analysis of earthquake-induced landslides, from shaking to failure.

58 GEOSCIENCES↗
DOE OSTI2022

Monitoring of mechanical changes in a pipe assembly with complex geometry using multi-mode acoustic signals

Guided wave ultrasonic testing has been widely used for nondestructive and long-range evaluation of diverse mechanical structures, such as pipes and plates. In general, a large number of guided wave modes exist even in the simplest geometries, and they exhibit highly dispersive characteristics. Thus, sophisticated excitation methods and elaborate modal analyses need to be implemented to excite, acquire, and process well-defined pure modes as part of such a measurement method. If the interrogated structure has complex geometry, for example, with cross-sectional variations or multiple joints, the utilization of pure modes becomes significantly more difficult. To overcome this difficulty, this paper proposes the use of multi-mode, multi-frequency guided waves for long-range inspection of structures with arbitrarily complex geometries. Furthermore, using broadband chirp signals and multi-mode time–frequency analysis in a prototypical complex pipe assembly that includes multiple pipe diameters, elbows, and flange connections, we experimentally show sensitive detection of mechanical changes (both material addition and loss). In addition, we demonstrate that the complexity of the structure can be modeled in the digital domain with sufficient detail to allow accurate numerical simulation of the acoustic response using this method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Floating platform effects on power generation in spar and semisubmersible wind turbines

The design and financing of commercial-scale floating offshore wind projects require a better understanding of how power generation differs between newer floating turbines and well-established fixed-bottom turbines. In floating turbines, platform mobility causes additional rotor motion that can change the time-averaged power generation. In this work, OpenFAST simulations examine the power generated by the National Renewable Energy Laboratory's 5-MW reference turbine mounted on the OC3-UMaine spar and OC4-DeepCWind semisubmersible floating platforms, subjected to extreme irregular waves and below-rated turbulent inflow wind from large-eddy simulations of a neutral atmospheric boundary layer. For these below-rated conditions, average power generation in floating turbines is most affected by two types of turbine displacements: an average rotor pitch angle that reduces power, caused by platform pitch; and rotor motion upwind-downwind that increases power, caused by platform surge and pitch. The relative balance between these two effects determines whether a floating platform causes power gains or losses compared to a fixed-bottom turbine; for example, the spar creates modest (3.1%–4.5%) power gains, whereas the semisubmersible creates insignificant (0.1%–0.2%) power gains for the simulated conditions. Furthermore, platform surge and pitch motions must be analyzed concurrently to fully capture power generation in floating turbines, which is not yet universal practice. Finally, a simple analytical model for predicting average power in floating turbines under below-rated wind speeds is proposed, incorporating effects from both the time-averaged pitch displacement and the dynamic upwind-downwind displacements.

17 WIND ENERGY↗
DOE OSTI2025

The Role of Kinetic Instabilities and Waves in Collisionless Magnetic Reconnection

Magnetic reconnection converts magnetic field energy into particle energy by breaking and reconnecting magnetic field lines. Magnetic reconnection is a kinetic process that generates a wide variety of kinetic waves via wave-particle interactions. Kinetic waves have been proposed to play an important role in magnetic reconnection in collisionless plasmas by, for example, contributing to anomalous resistivity and diffusion, particle heating, and transfer of energy between different particle populations. These waves range from below the ion cyclotron frequency to above the electron plasma frequency and from ion kinetic scales down to electron Debye length scales. This review aims to describe the progress made in understanding the relationship between magnetic reconnection and kinetic waves. We focus on the waves in different parts of the reconnection region, namely, the diffusion region, separatrices, outflow regions, and jet fronts. Particular emphasis is placed on the recent observations from the Magnetospheric Multiscale (MMS) spacecraft and numerical simulations, which have substantially increased the understanding of the interplay between kinetic waves and reconnection. Some of the ongoing questions related to waves and reconnection are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Roadmap on methods and software for electronic structure based simulations in chemistry and materials

This Roadmap article provides a succinct, comprehensive overview of the state of electronic structure methods and software for molecular and materials simulations. Seventeen distinct sections collect insights by 51 leading scientists in the field. Each contribution addresses the status of a particular area, as well as current challenges and anticipated future advances, with a particular eye towards software related aspects and providing key references for further reading. Foundational sections cover density functional theory and its implementation in real-world simulation frameworks, Green's function based many-body perturbation theory, wave-function based and stochastic electronic structure approaches, relativistic effects and semiempirical electronic structure theory approaches. Subsequent sections cover nuclear quantum effects, real-time propagation of the electronic structure, challenges for computational spectroscopy simulations, and exploration of complex potential energy surfaces. The final sections summarize practical aspects, including computational workflows for complex simulation tasks, the impact of current and future high-performance computing architectures, software engineering practices, education and training to maintain and broaden the community, as well as the status of and needs for electronic structure based modeling from the vantage point of industry environments. Overall, the field of electronic structure software and method development continues to unlock immense opportunities for future scientific discovery, based on the growing ability of computations to reveal complex phenomena, processes and properties that are determined by the make-up of matter at the atomic scale, with high precision.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Air-coupled tsunamis generated from impacts and airbursts: Our understanding before Hunga-Tonga Hunga-Ha'apai

The effort to prevent or mitigate the effects of an impact on Earth is known as planetary defense. A significant component of planetary defense research involves risk assessment. Much of our understanding of the risk from near-Earth objects comes from the geologic record in the form of impact craters, but not all asteroid impacts are crater-forming events. Small asteroids explode before reaching the surface, generating an airburst, and most impacts into the ocean do not penetrate the water to form a crater in the sea floor. The risk from these non-crater-forming ocean impacts and airbursts is difficult to quantify and represents a significant uncertainty in our assessment of the overall threat. We are currently working to better understand impact scenarios that can generate dangerous tsunamis. One of the suggested mechanisms for the production of asteroid–generated tsunamis is by direct coupling of the pressure wave to the water, analogous to the means by which a moving weather front can generate a meteotsunami. To test this hypothesis, we ran a series of airburst simulations and provided time-resolved pressure and wind profiles to use as source functions for tsunami models. We used the CTH hydrocode to model the various airburst scenarios to compare to the results of other simulations and provide time dependent boundary conditions as input to shallow-water wave propagation codes. The strongest and most destructive meteotsunamis are generated by atmospheric pressure oscillations with amplitudes of only a few hPa1 (mbar), corresponding to changes in sea level of a few cm. The resulting wave is strongest when there is a resonance between the ocean and the atmospheric forcing. A Proudman resonance takes place when the atmospheric disturbance’s translational speed (U) equals the longwave phase speed $\sqrt{gh}$ of shallow water wave. Coupling is strongest when the Froude number (Fr=U/c) is unity. A weather front propagates much slower than the speed of sound, so meteotsunamis are most common and dangerous in shallow bodies of water such as the Mediterranean Sea or Lake Michigan. By contrast, the blast wave from an airburst or crater-forming impact propagates at a speed faster than a tsunami in the deepest ocean, and a Proudman resonance cannot be achieved even though the overpressures are orders of magnitude greater. However, blast wave profiles are N-waves in which a sharp shock wave leading to overpressure is followed by a more gradual rarefaction to a much longer-duration underpressure phase. Even though the blast outruns the water wave it is forcing, the tsunami should continue to be driven by the out-of-resonance gradient associated with the suction phase, which may depend strongly on the details of the airburst or impact scenario. The open question is whether there are any conditions under which such an airburst-driven tsunami can be dangerous enough to contribute to the overall impact risk. We have also identified other potential mechanisms for airburst-generated tsunamis: 1) reaction force at the surface from the plume ejected into space, which carries significant momentum, 2) expanding toroidal vortices at the surface, which travel more slowly than the shock wave and can generate a Proudman resonance in relatively shallow ocean (such as continental shelf), and 3) steam explosion from seawater ablation by a “Type II” (Libyan Desert Glass-type) airburst in which the hot vapor jet descends to the surface. On January 15, 2022, the Hunga-Tonga Hunga-Ha’apai volcano, located approximately 60 km north of Tongatapu, the main island of Tonga, violently erupted with a powerful explosion, culminating the period of volcanic activity that started in December of 2021. This event and resulting tsunamis provided an existence proof for the air pressure wave coupling mechanism we proposed. It also suggests that it can be stronger and more significant over much greater distances than we contemplated, leading to global tsunamis associated with impact events on land as well as in the water. Large atmospheric explosions generate global Lamb waves with larger amplitudes, longer periods, and slower speeds than the local and regional blast waves we modeled prior to that event. This paper reviews our analysis and modeling of airburst-driven tsunamis prior to the 2022 Hunga-Tonga Hunga-Ha’apai tsunami, which was the subject of two presentations at the 2023 Planetary Defense Conference and is the subject of another paper currently in preparation.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1997

ENSO Mechanisms

This current leading theory of the El Nino-Southern Oscillation (ENSO) phenomenon involves reflection of equatorial waves on the western ocean boundary to shift the ENSO phase from El Nino to La Nina and vice versa. However, recent satellite altimetry and In situ observations indicate that this theory has several flaws. These include imperfect equatorial wave reflection on the western ocean boundary, a maxima of the simulated coupled wind-sea surface temperature (SST) interaction located too far into the eastern equatorial Pacific, and an excessive role of thermocline displacement on SST in this region. We show that the central equatorial Pacific SST is fundamental in driving the EIVSO ocean-atmosphere coupled system, through the east-west displacement of the eastern edge of the warm pool, perfectly related with ENSO. Within the equatorial wave-guide, the dominance of surface zonal advection in these displacements is demonstrated with four different satellite and in situ data sets and three ocean models. This demonstration is supported by the evidence of a convergence of water masses into the eastern edge of the warm pool, resulting in a well-defined salinity front. All these results lead us to propose a notable modification of the leading delayed-action oscillator theory for the oscillatory nature of ENSO. Simulations with a linearized coupled ocean-atmosphere model result in 3-6 year ENSO-like oscillations, with many Of the variable model parameters found very close to their observed values. This simple model suggests that ocean processes ignored or underestimated in the delayed action oscillator theory, such as zonal current convergence, zonal advection of sea surface temperature and equatorial wave reflection from the eastern-ocean boundary, are fundamental to the development of ENSO, in particular to its manifestations in the central equatorial Pacific.

Picaut, Joel↗
DOE OSTI2024

A protocol and analysis of year-long simulations of global storm-resolving models and beyond

We propose a protocol to evaluate and analyze year-long simulations of global storm-resolving models (GSRMs). The proposed protocol complements an earlier 40-day simulation protocol under the DYAMOND (DYnamics of the Atmospheric general circulation Modeled On Non-hydrostatic Domains) project to allow the analysis of the seasonal cycle and associated climatic relevant phenomena. This intercomparison aims to reveal how GSRMs, which can simulate mesoscale convective systems (MCSs) in the global domain, reproduce atmospheric large-scale structures related to convection beyond month-long simulations. The intercomparison for one-year simulations is conducted by either atmosphere-only models or atmosphere–ocean coupled models with atmospheric horizontal mesh sizes less than 5 km. We recommend the continuous four seasons from March 2020 to February 2021 as a target period for the intercomparison but with options for many groups to join more flexibly. The output variables are collected at 0.25° resolution, and archives of a small set of native grid variables are encouraged to analyze tropical cyclones and MCSs. Through the proposed global storm-resolving simulation, we will evaluate the climatological distributions of the atmospheric large-scale circulations, such as the Intertropical Convergence Zone (ITCZ), monsoon, midlatitude jets, their time evolution, and the upscale impacts on them. We present sample analyses from a one-year simulation using the 3.5 km mesh Nonhydrostatic Icosahedral Atmospheric Model (NICAM), revealing the realistic zonal contrast of tropical precipitation, no double ITCZ structure, the reasonable midlatitude jet position and intensity but a weak bias of storm track activities, and a warm bias over the Eurasia during boreal winter. We also clarify the cross-scale interaction, such as the effects of cold pools on mean precipitation over the Maritime Continent through the precipitation diurnal cycle and the effects of resolved gravity waves on midlatitude mean flows. The proposed one-year simulation protocol is referred to as the “Sendai Protocol.” This protocol is not unique or definite for evaluating GSRMs; we prospect a hierarchical set of experiments from short-term to multi-year simulations as GSRM intercomparisons.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1996

Three-Dimensional Electromagnetic Monte Carlo Particle-in-Cell Simulations of Critical Ionization Velocity Experiments in Space

Although the existence of the critical ionization velocity (CIV) is known from laboratory experiments, no agreement has been reached as to whether CIV exists in the natural space environment. In this paper we move towards more realistic models of CIV and present the first fully three-dimensional, electromagnetic particle-in-cell Monte-Carlo collision (PIC-MCC) simulations of typical space-based CIV experiments. In our model, the released neutral gas is taken to be a spherical cloud traveling across a magnetized ambient plasma. Simulations are performed for neutral clouds with various sizes and densities. The effects of the cloud parameters on ionization yield, wave energy growth, electron heating, momentum coupling, and the three-dimensional structure of the newly ionized plasma are discussed. The simulations suggest that the quantitative characteristics of momentum transfers among the ion beam, neutral cloud, and plasma waves is the key indicator of whether CIV can occur in space. The missing factors in space-based CIV experiments may be the conditions necessary for a continuous enhancement of the beam ion momentum. For a typical shaped charge release experiment, favorable CIV conditions may exist only in a very narrow, intermediate spatial region some distance from the release point due to the effects of the cloud density and size. When CIV does occur, the newly ionized plasma from the cloud forms a very complex structure due to the combined forces from the geomagnetic field, the motion induced emf, and the polarization. Hence the detection of CIV also critically depends on the sensor location.

Wang, J.↗
DOE OSTI2025

Interleaved Cuk Converter Wave Energy System With Advanced Control and Grid Support Functions for Seamless Integration

This paper presents an innovative wave energy conversion system that integrates an interleaved Cuk converter with advanced nonlinear control for seamless grid integration. The system efficiently extracts power on the DC side using the interleaved Cuk converter, while an inverter manages power transfer to the grid or load on the AC side. A nonlinear control architecture, based on the Lyapunov energy function, ensures stable and optimal operation under varying wave conditions, effectively addressing the challenges of variability and unpredictability inherent in wave energy. The proposed system also incorporates features to enhance power quality and minimize losses, making it a robust solution for renewable energy integration. The system's effectiveness in harnessing wave energy and achieving smooth grid integration is validated through comprehensive computer simulations in MATLAB/Simulink, with case study results demonstrating its capabilities.

HYDRO ENERGY,TIDAL AND WAVE POWER↗
NASA NTRS2007

Inter-annual and Long-term Temperature Variations in the Mesopause Region at High Latitudes Generated by the Stratospheric QBO

The Numerical Spectral Model (NSM) simulates the Quasi-biennial Oscillation (QBO) that dominates the zonal circulation of the lower stratosphere at low latitudes. In the model, the QBO is generated with parameterized small-scale gravity waves (GW), which are partially augmented in 3D with planetary waves owing to baroclinic instability. Due to GW filtering, the QBO extends into the upper mesosphere, evident in UARS zonal wind and TIMED temperature measurements. While the QBO zonal winds are confined to equatorial latitudes, even in simulations with latitude-independent wave source, the associated temperature variations extend to high latitudes. The meridional circulation redistributes some of the QBO energy to focus it partially onto the Polar Regions. The resulting QBO temperature variations away from the equator tend to increase at higher altitudes to produce inter-annual variations that can exceed 5 K in the polar mesopause region -- and our 3D model simulations show that the effect is variable from year to year and can produce large differences between the two hemispheres, presumably due to interactions involving the seasonal variations. Modeling studies with the NSM have shown that long-term variations can also be generated by the QBO interacting with the seasonal cycles through OW node-filtering. A 30-month QBO, optimally synchronized by the 6-month Semi-Annual Oscillation (SAO), thus produces a 5-year or semi-decadal (SD) oscillation -- and observational evidence for that has been provided by a recent analysis of stratospheric NCEP data. In a simulation with the 2D version of the NSM, this SD oscillation extends into the upper mesosphere, and we present results to show that the related temperature variations could contribute significantly to the long-term variations of the polar mesopause region. Quasi-decadal variations could furthermore arise from the modeled solar cycle modulations of the QBO and 12-month annual oscillation. Our numerical results are discussed in the context of the observed low summer temperatures reproduced by the model, to demonstrate that the above interannual and long-term variations could contribute significantly to the climatology of Polar Mesospheric Clouds (PMC) investigated by the Aeronomy of Ice in the Mesosphere (AIM) mission.

Mayr, Hans G.↗
DOE OSTI2026

Shakeup and shakeoff spectra in the electron capture decay of atomic 7 Be

The most stringent laboratory-based experimental limits on the existence of submegaelectronvolt sterile neutrinos are currently set by decay spectroscopy of radioactive 7 Be embedded into superconducting sensors. The systematic uncertainties are dominated by the modeling of the electron shakeup and shakeoff spectra that are not based on state-of-the-art atomic theory and do not include electron correlations or relativistic effects. We have used the multiconfiguration Dirac-Fock formalism to obtain correlated wave functions ab initio and compute all single and double shake processes in the electron capture decay of atomic 7 Be . The simulations can explain some but not all of the observed spectral features, likely because the wave functions are modified by the Ta sensor material that the 7 Be is embedded into. The new models also show that the 𝐿/𝐾 electron capture ratio of 7 Be in Ta has previously been slightly underestimated revising the previous value of 0.070(7) to a new value of 0.0756(20).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1990

TS - Dean interactions in curved channel flow

A weakly nonlinear theory is developed to study the interaction of TS waves and Dean vortices in curved channel flow. The prediction obtained from the theory agree well with results obtained from direct numerical simulations of curved channel flow, especially for low amplitude disturbances. At low Reynolds numbers the wave interaction is generally stabilizing to both disturbances, though as the Reynolds number increases, many linearly unstable TS waves are further destabilized by the presence of Dean vortices.

Singer, Bart A.↗
NASA NTRS1998

Can Kelvin-Helmholtz Instabilities of Jet-Like Structures and Plumes Cause Solar Wind Fluctuations at 1AU

During its south polar passage in 1994, Ulysses sampled the solar wind emerging from the south polar coronal hole. In analysing these data two types of density fluctuations lasting a few hours have been reported, one characterized by fluctuations in velocity gradients ("microstreams") and the other by magnetic plus thermal-pressure balance structures ("PBS"). Microstreams were both temporal and spatial in nature. At higher frequencies, MHD turbulence was observed and found to be less evolved than in the ecliptic, but essentially independent of heliographic latitude. It is argued here that microstreams, PBS, and MHD turbulence could all be the remnant of shears associated with plumes and other filamentary structures ("jets") which have been reported to exist in coronal holes. The shear between a jet and its ambient can become unstable to the MHD Kelvin-Helmholtz ("KH') instability at 5-10 solar radius and the propagating instability can cause fluctuations like those seen by Ulysses. This motivates us to simulate coronal jets using a 3D MHD ZEUS code. The first 2D results have just started to come and are promising. To study the KH instability the jet is perturbed at the boundary with a linear amplitude and fixed frequency. The jet seems to pass through various distinct phases, one of which Is apparently dominated by KH instabilities. These instabilities drive oblique shocks into the jet as the turbulent eddies contact the jet surface. It is known that KH instabilities and internal shock waves are partially suppressed by magnetic field tension. Hence, in simulating far along the jet the Instability is expected to produce Alfvenic fluctuations like those seen near 1 AU.

Parhi, Shyamsundar↗