Numerical solution for chemically generated waves in a dilute, isothermal atmosphere
Mathematical model for impulsive heat releases by chemically generated waves in atmosphere
SEARCH · Engineering Papers
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.
Mathematical model for impulsive heat releases by chemically generated waves in atmosphere
The dynamic fragmentation in impacts into solids is examined using a physical model for the formation and growth of cracks in rocks. The physical model is then inserted into a numerical model (hydrocode) of stress wave propagation and interaction, from which the outcome of a given impact event can be computed. The hydrocode model predicts fragment sizes due to impact in terms of shock waves propagating in a homogeneous elastic medium containing a distribution of crack nucleation centers known as Weibull flaws.
The Source Physics Experiment (SPE) is a series of chemical explosions at the Nevada National Security Site (NNSS) with the goal of understanding seismic-wave generation and propagation of underground explosions. To understand explosion source physics, accurate geophysical models of the SPE site are needed. Here, we utilize a large-N seismic array deployed at the SPE phase II site to generate a shallow subsurface model of shear-wave velocity. The deployment consists of 500 geophones and covers an area of, approximately, 2.5 × 2 km. The array is located in the Yucca Flat in the northeast corner of the NNSS, Nye County, Nevada. Using ambient-noise recordings throughout the large-N seismic array, we calculate horizontal-to-vertical spectral ratios (HVSRs) across the array. We obtain 2D seismic images of shear-wave velocities across the SPE phase II site for the shallow structure of the basin. In this work the results clearly image two significant seismic impedance interfaces at ~150–500 and ~350–600 m depth. The shallower interface relates to the contrast between Quaternary alluvium and Tertiary volcanic rocks. The deeper interface relates to the contrast between Tertiary volcanic rocks and the Paleozoic bedrock. The 2D subsurface models support and extend previous understanding of the structure of the SPE phase II site. This study shows that the HVSR method in conjunction with a large-N seismic array is a quick and effective method for investigating shallow structures.
Human cognitive responses, behavioral responses, and disease dynamics co-evolve over the course of any disease outbreak, and can result in complex feedbacks. We present a dynamic agent-based model that explicitly couples the spread of disease with the spread of fear surrounding the disease, implemented within the EpiCast simulation framework. EpiCast models transmission within a realistic synthetic population, capturing individual-level interactions. In our model, fear propagates through both in-person contact and broadcast media, prompting individuals to adopt protective behaviors that reduce disease spread. In order to better understand these coupled dynamics, we create and compare a range of compartmental models to ensure that introducing additional disease states does not prevent the emergence of multiple waves in these simpler models. Additionally, we compare a range of behavioral scenarios within EpiCast, varying the level and intensity of fear and behavior change. Our results show that the addition of asymptomatic, exposed, and pre-symptomatic disease states can impact both the rate at which an outbreak progresses and its overall trajectory in compartmental models. In EpiCast, the combination of non-local fear spread via broadcasters and strong behavioral responses by fearful individuals generally leads to multiple epidemic waves, an outcome that occurs only within a narrow parameter range when fear spreads purely through local contact. Accounting for the coupled spread of fear and disease is critical for understanding disease dynamics and designing timely, targeted responses to emerging infectious threats.
Abstract The stress tensor is an important property for upper crustal studies such as those that involve pore fluids and earthquake hazards. At tectonic plate scale, plate boundary forces and mantle convection are the primary drivers of the stress field. In many local settings (10–100 s of km and <10 km depth) in tectonic plate interiors, we can simplify by assuming a constant background stress field that is perturbed by local heterogeneity in density and elasticity. Local stress orientation and sometimes magnitude can be estimated from earthquake and borehole‐based observations when available. Modeling of the local stress field often involves interpolating sparse observations. We present a new method to estimate the 3D stress field in the upper crust and demonstrate it for Oklahoma. We created a 3D material model by inverting multiple types of geophysical observations simultaneously. Integrating surface‐wave dispersion, local travel times and gravity observations produces a model of P‐wave velocity, S‐wave velocity, and density. The stress field can then be modeled using finite element simulations. The simulations are performed using our simplified view of the local stress field as the sum of a constant background stress field that is perturbed by local density and elasticity heterogeneity and gravitational body forces. An orientation of N82°E, for the maximum compressive tectonic force, best agrees with previously observed stress orientations and faulting types in Oklahoma. The gravitational contribution of the horizontal stress field has a magnitude comparable to the tectonic contribution for the upper 5 km of the subsurface.
According to collisionless shock theories, the thickness of a shock front should be of the order of the characteristic lengths of the plasmas (the Debye length, the proton and Larmor radii, etc.). Chao and Lepping (1974), found, however, that 30% of the observed interplanetary shocks at 1 AU have thicknesses much larger than these characteristic lengths. It is the objective of the present paper to investigate whether the competition between nonlinear steepening and Landau damping can result in a wave of finite width that does not steepen into a shock. A heuristic model of such a wave is developed and tested by the examples of two structures that are qualitatively shocklike, but thicker than expected from theory. It is found that both events are in the process of steepening and their limiting thicknesses due to Landau damping are greater than the corresponding proton Larmor radius for both structures as observed at Mariner 5 (nearer the sun than 1 AU) but are comparable to the proton Larmor radius for Explorer (near 1 AU) observations.
Boundary-layer winds near steeply sloped topography excite a distinct type of atmospheric wave. Basic properties of the waves were determined through analysis of 7 Mars years of temperature profiles from the Mars Climate Sounder on Mars Reconnaissance Orbiter. A more complete understanding of the waves was derived from a simulation with the NASA Ames Mars General Circulation Model. This type of wave appears in the tropics in the vicinity of conspicuous surface features, such as large volcanoes and the rim of Isidis Planitia. The amplitude is typically 5–10 K, the vertical wavelength is about 40 km, the zonal wavelength is about 1400 km, and the zonal phase speed relative to the surface is about 25 m s−1. The wave forcing comes from horizontal convergence of the near-surface winds, which generates a daily surge in the vertical winds. These convergence zones were not previously recognized as an important source of atmospheric waves on Mars. Strong forcing occurs only within a limited range of local time, producing wave packets (diurnal pulses) rather than steady oscillations. Waves are present throughout the year and their properties vary with season. Near the solstices, the wave amplitude is largest in the tropics of the summer hemisphere, where the convergence zones are stronger than in the winter hemisphere. The diurnal-mean zonal winds, which vary with latitude and season, control the direction of propagation of the predominant waves through filtering by critical levels.
A flexible accurate large-signal computer program has been developed for the design of coupled-cavity traveling wave tubes. The program is written for a TSS-360 time sharing system. The beam is described by a disk model and the slow wave structure by a sequence of cavities or cells. The computational approach is arranged so that each cavity may have different geometrical or electrical parameters than its neighbors. This allows the program user to simulate a tube of almost arbitrary complexity. Input and output couplers, severs, complicated velocity tapers, and other features peculiar to one or a few cavities may be modeled by a correct choice of input data. The beam-wave interaction is handled by a new approach in which the RF fields are expanded in solutions to the TM wave equation retaining all significant space harmonics. The program was used to perform a design study of the TWT developed for the CTS satellite. Good agreement was obtained between the predictions of the program and the measured performance of the flight tube. The internal check on power balance was satisfied within plus or minus 0.2 per cent of input beam power.
This paper was directed at active control of discrete frequency noise generated by subsonic blade rows through cancellation of the blade row interaction generated propagating acoustic waves. First discrete frequency noise generated by a rotor and stator in a duct was analyzed to determine the propagating acoustic pressure waves. Then a mathematical model was developed to analyze and predict the active control of discrete frequency noise generated by subsonic blade rows through cancellation of the propagating acoustic waves, accomplished by utilizing oscillating airfoil surfaces to generate additional control propagating pressure waves. These control waves interact with the propagating acoustic waves, thereby, in principle, canceling the acoustic waves and thus, the far field discrete frequency tones. This model was then applied to a fan exit guide vane to investigate active airfoil surface techniques for control of the propagating acoustic waves, and thus the far field discrete frequency tones, generated by blade row interactions.
In earlier experimental and analytical studies, it was found that the boundary layer on an aircraft could provide significant shielding from propeller noise at typical transport airplane cruise Mach numbers. In this paper a new three-dimensional theory is described that treats the combined effects of refraction and scattering by the fuselage and boundary layer. The complete wave field is solved by matching analytical expressions for the incident and scattered waves in the outer flow to a numerical solution in the boundary layer flow. The model for the incident waves is a near-field frequency-domain propeller source theory developed previously for free field studies. Calculations for an advanced turboprop (Prop-Fan) model flight test at 0.8 Mach number show a much smaller than expected pressure amplification at the noise directivity peak, strong boundary layer shielding in the forward quadrant, and shadowing around the fuselage. Results are presented showing the difference between fuselage surface and free-space noise predictions as a function of frequency and Mach number. Comparison of calculated and measured effects obtained in a Prop-Fan model flight test show good agreement, particularly near and aft of the plane of rotation at high cruise Mach number.
In this paper a model problem is considered that simulates an atmospheric acoustic wave propagation situation that is nonlinear. The model is derived from the basic Euler equations for the atmospheric flow and from the regular perturbations for the acoustic part. The nonlinear effects are studied by obtaining two successive linear problems in which the second one involves the solution of the first problem. Well-posedness of these problems is discussed and approximations of the radiation boundary conditions that can be used in numerical simulations are presented.
It is shown that disturbances in shock and expansion tubes can be modelled as lateral acoustic waves. The ratio of sound speed across the driver-test gas interface is shown to govern the quantity of noise in the test gas. Frequency 'focusing' which is fundamental to centered unsteady expansions is discussed and displayed in centerline pitot pressure measurements.
Results of both a linear and a nonlinear study of oxygen cyclotron waves and the associated oxygen heating are presented. Linear theory predicts that oxygen cyclotron waves will have smaller growth rates than hydrogen cyclotron waves. Results of a simulation study in which the free energy source is an initial drifting electron distribution indicate that oxygen cyclotron waves only grow to small amplitudes, while the hydrogen cyclotron waves achieve larger amplitudes. In an attempt to model more realistically the continuous ionospheric outflow, a simulation model is used, in which the electron velocity distribution is maintained by a constant flow of electrons. This latter model predicts that the oxygen waves grow to amplitudes much larger than the hydrogen waves resulting in the preferential heating of the heavier ions.
A survey of various one- and two-way interactive nested grid techniques used in hydrostatic numerical weather prediction models is presented and the advantages and disadvantages of each method are discussed. The techniques for specifying the lateral boundary conditions for each nested grid scheme are described in detail. Averaging and interpolation techniques used when applying the coarse mesh grid (CMG) and fine mesh grid (FMG) interface conditions during two-way nesting are discussed separately. The survey shows that errors are commonly generated at the boundary between the CMG and FMG due to boundary formulation or specification discrepancies. Methods used to control this noise include application of smoothers, enhanced diffusion, or damping-type time integration schemes to model variables. The results from this survey provide the information needed to decide which one-way and two-way nested grid schemes merit future testing with the Mesoscale Atmospheric Simulation System (MASS) model. An analytically specified baroclinic wave will be used to conduct systematic tests of the chosen schemes since this will allow for objective determination of the interfacial noise in the kind of meteorological setting for which MASS is designed. Sample diagnostic plots from initial tests using the analytic wave are presented to illustrate how the model-generated noise is ascertained. These plots will be used to compare the accuracy of the various nesting schemes when incorporated into the MASS model.
In the last two decades there have been extensive developments in computational unsteady transonic aerodynamics. Such developments are essential since the transonic regime plays an important role in the design of modern aircraft. Therefore, there has been a large effort to develop computational tools with which to accurately perform flutter analysis at transonic speeds. In the area of Computational Fluid Dynamics (CFD), unsteady transonic aerodynamics are characterized by the feature of modeling the motion of shock waves over aerodynamic bodies, such as wings. This modeling requires the solution of nonlinear partial differential equations. Most advanced codes such as XTRAN3S use the transonic small perturbation equation. Currently, XTRAN3S is being used for generic research in unsteady aerodynamics and aeroelasticity of almost full aircraft configurations. Use of Euler/Navier Stokes equations for simple typical sections has just begun. A brief history of the development of CFD for aeroelastic applications is summarized. The development of unsteady transonic aerodynamics and aeroelasticity are also summarized.
We show how the generation of right-handed neutrino masses in Majoron models may be associated with a first-order phase transition and accompanied by the production of a stochastic background of gravitational waves (GWs). We explore different energy scales with only renormalizable operators in the effective potential. If the phase transition occurs above the electroweak scale, the signal can be tested by future interferometers. We consider two possible energy scales for phase transitions below the electroweak scale. If the phase transition occurs at a GeV, the signal can be tested at LISA and provide a complementary cosmological probe to right-handed neutrino searches at the FASER detector. If the phase transition occurs below 100 keV, we find that the peak of the GW spectrum is two or more orders of magnitude below the putative NANOGrav GW signal at low frequencies, but well within reach of the SKA and THEIA experiments. We show how searches of very low frequency GWs are motivated by solutions to the Hubble tension in which ordinary neutrinos interact with the dark sector. We also present general calculations of the phase transition temperature and Euclidean action that apply beyond Majoron models.
Accurate aerodynamic computational predictions are essential for the safety of space vehicles, but these computations are of limited accuracy when large pressure gradients are present in the flow. The goal of the current project is to improve the state of compressible turbulence modeling for high speed flows with shock wave / turbulent boundary layer interactions (SWTBLI). Emphasis will be placed on models that can accurately predict the separated region caused by the SWTBLI. These flows are classified as nonequilibrium boundary layers because of the very large and variable adverse pressure gradients caused by the shock waves. The lag model was designed to model these nonequilibrium flows by incorporating history effects. Standard one- and two-equation models (Spalart Allmaras and SST) and the lag model will be run and compared to a new lag model. This new model, the Reynolds stress tensor lag model (lagRST), will be assessed against multiple wind tunnel tests and correlations. The basis of the lag and lagRST models are to preserve the accuracy of the standard turbulence models in equilibrium turbulence, when the Reynolds stresses are linearly related to the mean strain rates, but create a lag between mean strain rate effects and turbulence when nonequilibrium effects become important, such as in large pressure gradients. The affect this lag has on the results for SWBLI and massively separated flows will be determined. These computations will be done with a modified version of the OVERFLOW code. This code solves the RANS equations on overset grids. It was used for this study for its ability to input very complex geometries into the flow solver, such as the Space Shuttle in the full stack configuration. The model was successfully implemented within two versions of the OVERFLOW code. Results show a substantial improvement over the baseline models for transonic separated flows. The results are mixed for the SWBLI assessed. Separation predictions are not as good as the baseline models, but the over prediction of the peak heat flux downstream of the reattachment shock that plagues many models is reduced.
Ozone transport is calculated for steady, dissipative planetary waves using the Eulerian, Lagrangian mean, and residual circulation. A Lagrangian model of parcel dynamics is used to interpret planetary wave-photochemistry interaction. In chemically active regions the mean field ozone changes are found to be significant only where there are large gradients in chemical sources and sinks along parcel trajectories. The largest changes in the mean field are found in the lower stratosphere and are due to the Lagrangian mean advection. When the Lagrangian mean advection is approximated by the residual circulation, errors in the transport velocities as large as 30 pct may occur.