Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “wave modeling”

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,081 records · Page 60

Standing shocks in a two-fluid solar wind

We present a numerical study of the formation of standing shocks in the solar wind using a two-fluid time-dependent model in the presence of Alfven waves. Included in this model is the adiabatic cooling and thermal conduction of both electrons and protons. In this study, standing shocks develop in the flow when additional critical points form as a result of either localized momentum addition or rapid expansion of the flow tube below the existing sonic point. While the flow speed and density exhibit the same characteristics as found in earlier studies of the formation of standing shocks, the inclusion of electron and proton heat conduction produces different signatures in the electron and proton temperature profiles across the shock layer. Owing to the strong heat conduction, the electron temperature is nearly continuous across the shock, but its gradient has a negative jump across it, thus producing a net heat flux out of the shock layer. The proton temperature exhibits the same characteristics for shocks produced by momentum addition but behaves differently when the shock is formed by the rapid divergence of the flow tube. The adiabatic cooling in a rapidly diverging flow tube reduces the proton temperature so substantially that the proton heat conduction becomes negligible in the vicinity of the shock. As a result, protons experience a positive jump in temperature across the shock. While Alfven waves do not affect the formation of standing shocks, they contribute to the change of the mmomentum and energy balance across them. We also find that for this solar wind model the inclusion of thermal conduction and adiabatic cooling for the elctrons and protons increases significantly the range of parameters characterizing the formation of standing shocks over those previously found for isothermal and polytropic models.

Habbal, Shadia R.↗

Molecular shock waves in the BN-KL region of Orion

The H2, CO, and OH IR emission-line intensities in the BN-KL region of Orion are modeled by shock waves traveling at 36 km/s with respect to an ambient gas of density 200,000 per cu cm, ionization fraction of less than 3 x 10 to the -7th, CO abundance of approximately 0.0003, and transverse magnetic field of 0.45 milligauss. A global model is considered in which the shock waves are driven by strong winds emanating from KL. The mass-loss rate has been not less than 0.003 solar mass per year for approximately 1000 years. The total wind energy released from the star is not less than 5 x 10 to the 47th ergs. The total ejected momentum is approximately 10 to the 41st g cm/S, sufficient to disrupt the dense cloud surrounding BN-KL.

Chernoff, D. F.↗

Polyimide dynamically compressed to decomposition pressures: Two-wave structures captured by velocimetry and modeling

We performed a series of six plate impact experiments on polyimide and modeled them using new reactant and products equations of state combined with an Arrhenius rate model. The first experiment was diagnosed with embedded electromagnetic velocity gauges through which we directly observed attenuation of the lead shock to an approximately constant state over a propagation distance of roughly 4 mm. Simulated gauge profiles were in excellent qualitative agreement with experiment and suggested a sluggish chemical reaction that did not proceed to completion. The remaining five experiments were conducted in a transmission geometry and diagnosed velocimetrically at the sample/window interface. All five of these yielded profiles with a sharp shock followed by a more gradual approach to maximum interface velocity that was “rounded” to varying degree. These profiles proved difficult to interpret unambiguously due to the convolution of the reactive wave upon first shock with reflection of the lead wave and reshock or release by the window. Comparison with thermochemical calculations strongly suggests that the point of maximum interface velocity corresponds to the equilibrium reshock or release locus. We discuss the implications of this point for the practice of impedance matching based on the reflected Hugoniot of reactive materials such as polymers. The reactant and thermochemical products equations of state are developmental SESAME tables 97710 and 97720, respectively.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Breakdown of vertically propagating two-dimensional gravity waves forced by orography

The propagation of orographic gravity waves into an atmosphere with exponentially decreasing density is simulated with a two-dimensional, nonlinear, time-dependent numerical model. After the stationary wave is established over the mountain, the model predicts that wave breaking causes a large reduction of the vertical momentum flux in the flow, not only at levels where wave breaking is present, but also far below the lowest occurrence of overturning. More than half of the decrease in momentum flux is explained by the presence of large amplitude, downward propagating waves, which are generated in regions of wave breaking. The downward propagating waves appear almost simultaneously with overturning, and have nonzero phase speeds, suggesting a strongly nonlinear generation mechanism that depends on local wave properties. The generation of these downward propagating waves is a robust process, insensitive to mountain height, mountain width, or density scale height. These results have important implications for observational studies of orographically generated waves as well as for schemes that seek to parameterize the effects of orogoraphy in large-scale models.

Bacmeister, Julio T.↗

Parametric raytracing modeling for NSTX-U scenario development with high harmonic fast waves and neutral beam injection

High harmonic fast waves (HHFW) are a versatile heating and current drive tool for scenario development. Extensive modeling scans were performed to find optimal parameters for different uses of HHFW in National Spherical Tokamak Experiment (NSTX-U). Scans of plasma density, temperature, magnetic field, and antenna phasing were performed both with and without neutral beam injection. For speed of calculation, the ray-tracing code GENRAY coupled to the quasilinear Fokker–Planck code CQL3D was used. CQL3D allows for a more accurate description of the fast ion population, as well as for quasilinear effects such as HHFW-induced modifications of the distribution function. Best current drive results are obtained at elevated electron temperatures and with the lowest k φ phasing. Adding neutral beams however typically strongly reduces the HHFW current drive efficiency at the low density cases due to HHFW absorption on beam ions. Results of this parametric study will feed into scenario development and predict-first whole-shot modeling of NSTX-U discharges.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A Self-Consistent Model of the Interacting Ring Current Ions with Electromagnetic ICWs

Initial results from a newly developed model of the interacting ring current ions and ion cyclotron waves are presented. The model is based on the system of two bound kinetic equations: one equation describes the ring current ion dynamics, and another equation describes wave evolution. The system gives a self-consistent description of ring current ions and ion cyclotron waves in a quasilinear approach. These two equations were solved on a global scale under non steady-state conditions during the May 2-5, 1998 storm. The structure and dynamics of the ring current proton precipitating flux regions and the wave active zones at three time cuts around initial, main, and late recovery phases of the May 4, 1998 storm phase are presented and discussed in detail. Comparisons of the model wave-ion data with the Polar/HYDRA and Polar/MFE instruments results are presented..

Khazanov, G. V.↗

Multi-Level Analysis of Pulsed Detonation Engines

The present study explores some issues concerning the operational performance of pulsed detonation engines. Zero-, one- and two-dimensional, transient models are employed in a synergistic manner to elucidate the various characteristics that can be expected from each level of analysis. The zero-dimensional model provides rapid parametric trends that help to identify the global characteristics of pulsed detonation engines. The one-dimensional model adds key wave propagation issues that are omitted in the zero-dimensional model and helps to assess its limitations. Finally, the two-dimensional model allows estimates of the first-order multi-dimensional effects and provides an initial multi-dimensional end-correction for the one-dimensional model. The zero-dimensional results indicate that the pulsed detonation engine is competitive with a rocket engine when exhausting to vacuum conditions. At finite back pressures, the PDE out-performs the rocket if the combustion pressure rise from the detonation is added to the chamber pressure in the rocket. If the two peak pressures are the same, the rocket performance is higher. Two-dimensional corrections added to the one-dimensional model result in a modest improvement in predicted specific impulse over the constant pressure boundary condition.

Ebrahimi, Houshang B.↗

Convectively Generated Gravity Waves In The Tropical Stratosphere: Case Studies And Importance For The Circulation Of The Middle Atmosphere

The advent of high altitude aircraft measurements in the stratosphere over tropical convective systems has made it possible to observe the mesoscale disturbances in the temperature field that these systems excite. Such measurements show that these disturbances have horizontal scales comparable to those of the underlying anvils (about 50-100 km) with peak to peak theta surface variations of about 300-400 meters. Moreover, correlative wind measurements from the tropical phase of the Stratosphere-Troposphere Exchange Project (STEP) clearly show that these disturbances are gravity waves. We present two case studies of anvil-scale gravity waves over convective systems. Using steady and time-dependent linear models of gravity wave propagation in the stratosphere, we show: (1) that the underlying convective systems are indeed the source of the observed phenomena; and (2) that their generating mechanism can be crudely represented as flow over a time-dependent mountain. We will then discuss the effects gravity waves of the observed amplitudes have on the circulation of the middle atmosphere, particularly the quasi-biennial, and semiannual oscillations.

Pfister, Leonhard↗

ERS-1 and Seasat scatterometer measurements of ocean winds: Model functions and the directional distribution of short waves

Calculation of accurate vector winds from scatterometers requires knowledge of the relationship between backscatter cross-section and the geophysical variable of interest. As the detailed dynamics of wind generation of centimetric waves and radar-sea surface scattering at moderate incidence angles are not well known, empirical scatterometer model functions relating backscatter to winds must be developed. Less well appreciated is the fact that, given an accurate model function and some knowledge of the dominant scattering mechanisms, significant information on the amplitudes and directional distributions of centimetric roughness elements on the sea surface can be inferred. accurate scatterometer model functions can thus be used to investigate wind generation of short waves under realistic conditions. The present investigation involves developing an empirical model function for the C-band (5.3 GHz) ERS-1 scatterometer and comparing Ku-band model functions with the C-band model to infer information on the two-dimensional spectrum of centimetric roughness elements in the ocean. The C-band model function development is based on collocations of global backscatter measurements with operational surface analyses produced by meteorological agencies. Strengths and limitations of the method are discussed, and the resulting model function is validated in part through comparison with the actual distributions of backscatter cross-section triplets. Details of the directional modulation as well as the wind speed sensitivity at C-band are investigated. Analysis of persistent outliers in the data is used to infer the magnitudes of non-wind effects (such as atmospheric stratification, swell, etc.). The ERS-1 C-band instrument and the Seasat Ku-band (14.6 GHz) scatterometer both imaged waves of approximately 3.4 cm wavelength assuming that Bragg scattering is the dominant mechanism. Comparisons of the C-band and Ku-band model functions are used both to test the validity of the postulated Bragg mechanism and to investigate the directional distribution of the imaged waves under a variety of conditions where Bragg scatter is dominant.

Freilich, Michael H.↗

Strong plasma turbulence in the earth's electron foreshock

A quantitative model is developed to account for the distribution in magnitude and location of the intense plasma waves observed in the earth's electron foreshock given the observed rms levels of waves. In this model, nonlinear strong-turbulence effects cause solitonlike coherent wave packets to form and decouple from incoherent background beam-excited weak turbulence, after which they convect downstream with the solar wind while collapsing to scales as short as 100 m and fields as high as 2 V/m. The existence of waves with energy densities above the strong-turbulence wave-collapse threshold is inferred from observations from IMP 6 and ISEE 1 and quantitative agreement is found between the predicted distribution of fields in an ensemble of such wave packets and the actual field distribution observed in situ by IMP 6. Predictions for the polarization of plasma waves and the bandwidth of ion-sound waves are also consistent with the observations. It is shown that strong-turbulence effects must be incorporated in any comprehensive theory of the propagation and evolution of electron beams in the foreshock. Previous arguments against the existence of strong turbulence in the foreshock are refuted.

Robinson, P. A.↗

Electric and magnetic observations of the structure of standing waves in the magnetosphere

Electric and magnetic field instruments on the DE 1 spacecraft provided observations of toroidal standing wave oscillations of magnetic shells in 1981 and 1982. The amplitude and nodal structure of these waves, as a function of increasing magnetic latitude, was studied and compared with model calculations of these waves. Eight pulsation events were selected for study. The pulsations in these events ranged in period from 20 to 300 s. They were observed at L values between 2.8 and 9.7 and at magnetic latitudes from the equator to 48 deg. The duration of these pulsation events was 20 to 30 min, and they usually consisted of a 5- to 10-min growth in amplitude, a 5- to 10-min interval of fluctuating amplitude, and a 5- to 20-min interval of decreasing amplitude. The relative E and B amplitudes and the phase relations between the E and B waves could be qualitatively explained by the nodal structure of the waves predicted by model calculations. Detailed quantitative comparisons of E and B amplitudes were in agreement with more recent model calculations. In a few cases, electron densities were also available as derived from plasma frequency observations; in these cases the standing wave periods calculated were reasonably close to, but somewhat less than, the observed periods.

Cahill, L. J., Jr.↗

Electromagnetic Modeling of Distributed-Source-Excitation of Coplanar Waveguides: Applications to Traveling-Wave Photomixers

In this work an electromagnetic model and subsequent design is presented for a traveling-wave, coplanar waveguide (CPW) based source that will operate in the THz frequency regime. The radio frequency (RF) driving current is a result of photoexcitation of a thin GaAs membrane using two frequency-offset lasers. The GaAs film is grown by molecular-beam-epitaxy (MBE) and displays sub-ps carrier lifetimes which enable the material conductivity to be modulated at a very high rate. The RF current flows between electrodes deposited on the GaAs membrane which are biased with a DC voltage source. The electrodes form a CPW and are terminated with a double slot antenna that couples the power to a quasi-optical system. The membrane is suspended above a metallic reflector to launch all radiation in one direction. The theoretical investigation and consequent design is performed in two steps. The first step consists of a direct evaluation of the magnetic current distribution on an infinitely extended coplanar waveguide excited by an impressed electric current distributed over a finite area. The result of the analysis is the difference between the incident angle of the laser beams and the length of the excited area that maximizes the RF power coupled to the CPW. The optimal values for both parameters are found as functions of the CPW and membrane dimensions as well as the dielectric constants of the layers. In the second step, a design is presented of a double slot antenna that matches the CPW characteristic impedance and gives good overall performance. The design is presently being implemented and measurements will soon be available.

Pasqualini, Davide↗

Electromagnetic scattering theory

Electromagnetic scattering theory is discussed with emphasis on the general stochastic variational principle (SVP) and its applications. The stochastic version of the Schwinger-type variational principle is presented, and explicit expressions for its integrals are considered. Results are summarized for scalar wave scattering from a classic rough-surface model and for vector wave scattering from a random dielectric-body model. Also considered are the selection of trial functions and the variational improvement of the Kirchhoff short-wave approximation appropriate to large size-parameters. Other applications of vector field theory discussed include a general vision theory and the analysis of hydromagnetism induced by ocean motion across the geomagnetic field. Levitational force-torque in the magnetic suspension of the disturbance compensation system (DISCOS), now deployed in NOVA satellites, is also analyzed using the developed theory.

Bird, J. F.↗

Using AI for Wave-front Estimation with the Rubin Observatory Active Optics System

Abstract The Vera C. Rubin Observatory will, over a period of 10 yr, repeatedly survey the southern sky. To ensure that images generated by Rubin meet the quality requirements for precision science, the observatory will use an active-optics system (AOS) to correct for alignment and mirror surface perturbations introduced by gravity and temperature gradients in the optical system. To accomplish this, Rubin will use out-of-focus images from sensors located at the edge of the focal plane to learn and correct for perturbations to the wave front. We have designed and integrated a deep-learning (DL) model for wave-front estimation into the AOS pipeline. In this paper, we compare the performance of this DL approach to Rubin’s baseline algorithm when applied to images from two different simulations of the Rubin optical system. We show the DL approach is faster and more accurate, achieving the atmospheric error floor both for high-quality images and low-quality images with heavy blending and vignetting. Compared to the baseline algorithm, the DL model is 40× faster, the median error 2× better under ideal conditions, 5× better in the presence of vignetting by the Rubin camera, and 14× better in the presence of blending in crowded fields. In addition, the DL model surpasses the required optical quality in simulations of the AOS closed loop. This system promises to increase the survey area useful for precision science by up to 8%. We discuss how this system might be deployed when commissioning and operating Rubin.

79 ASTRONOMY AND ASTROPHYSICS↗

Simulation of Physics-Based 0-10Hz Strong Motion Using High Performance Computing Supporting Refinements to Regional Ground Motion Models for the Central Eastern US

In collaboration with the U.S. Nuclear Regulatory Commission (NRC) the LLNL has developed a computationally efficient simulation platform designed to perform physics-based ground motion simulations for crustal earthquakes in the Stable Continental Regions of Central and Eastern US (CEUS), using high-performance computing. The main objective of the earthquake simulations was to use synthetic ground motion to provide constrains to refinements of existing ergodic Ground Motion Models (GMMs), for large magnitude earthquakes and near-fault distances, for which these models are less reliable. Physics-based broadband (0-10Hz) ground motion simulations were used to estimate the near-fault ground motion amplitudes and within event and between-event variabilities associated with fault rupture characteristics. In our simulations we used a 3D regional velocity model that was based on Saikia’s 1D velocity model (1994). In simulations performed during the first stage of this project the Saikia’s velocity model demonstrated better performance in modelling high frequency regional wave propagation for the CEUS region recorded during the Mw5.0 November 7, 2016, Cushing Oklahoma (Taylor et al., 2017), and Mw5.8 September 3, 2016, Pawnee Oklahoma earthquakes. The proposed regional 3D model includes random perturbations to the 1D background model using the stochastic scheme of Pitarka and Mellors (2021). In addition, validation analysis of the rupture generator and regional wave propagation models, using comparisons with different GMMs for Mw6.5 and Mw7.0 scenario earthquakes in the CEUS region resulted in a very good match between the simulated and empirical ground motion models. For the purposes of seismic hazard assessment at the existing and planned nuclear power plants, NRC is interested in studies aimed at improving the current ground motion models (GMM) for both Stable Continental Regions (SCR) in the Central and Eastern US and Active Crustal Regions (ACR) in the Western US. Due to lack of recorded data, these improvements require synthetic data for short fault distances and large magnitude earthquakes for which the existing recorded data is not enough to uniquely constrain the GMMs. The need for simulations and strong motion data is especially critical for the CEUS region where we do not have recorded data from potentially large damaging earthquakes with moment magnitudes 6.0 and higher. In this the project, we focused on 10Hz simulations of Mw7.0 scenario earthquakes with strike slip and thrust faulting mechanisms. We used more than 50 Mw7.0 earthquake rupture scenarios to investigate the ground motion uncertainty due to unknown earthquake rupture parameters, in particular, the slip distribution, rupture velocity, and faulting mechanism, and their implication on ground motion amplification due to forward rupture directivity effects.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

On error detection in the dynamics part of primitive equation models

Rossby-Haurwitz waves have been commonly employed as initial conditions for testing the dynamics part of global atmospheric models. Difficulties arise in connection with the primitive equation (PE) model. Chao and Geller (1982) have suggested that it would be better if normal mode solutions of the linearized primitive equations would be used as initial conditions. However, tests involving normal mode conditions are not adequate for the detection of errors in the nonlinear terms. The present investigation has the objective to review the study conducted by Chao and Geller and to provide tests for the nonlinear terms.

Chao, W. C.↗

Fluctuating charge-density-wave correlations in the three-band Hubbard model

The high-temperature superconducting cuprates host unidirectional spin- and charge-density-wave orders that can intertwine with superconductivity in nontrivial ways. While the charge components of these stripes have now been observed in nearly all cuprate families, their detailed evolution with doping varies across different materials and at high and low temperatures. We address this problem using nonperturbative determinant quantum Monte Carlo calculations for the three-band Hubbard model. Using an efficient implementation, we can resolve the model’s fluctuating spin and charge modulations and map their evolution as a function of the charge transfer energy and doping. We find that the incommensurability of the charge modulations is decoupled from the spin modulations and decreases with hole doping, consistent with experimental measurements at high temperatures. These findings support the proposal that the high-temperature charge correlations are distinct from the intertwined stripe order observed at low-temperature and in the single-band Hubbard model.

42 ENGINEERING↗