Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “wave equation”

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 415 records · Page 23

On the equations governing the axisymmetric perturbations of the Kerr black hole

It is shown how Teukolsky's equation, governing the perturbations of the Kerr black hole, can be reduced, in the axisymmetric case, to a one-dimensional wave equation with four possible potentials. The potentials are implicitly dependent on the frequency; and besides, depending on circumstances, they can be complex. In all cases (i.e. whether or not the potentials are real or complex), the problem of the reflexion and the transmission of gravitational waves by the potential barriers can be formulated, consistently, with the known conservation laws. It is further shown that all four potentials lead to the same reflexion and transmission coefficients.

Chandrasekhar, S.↗

Nondestructive Evaluation of Concrete: Elastic Property Imaging Through Full Waveform Inversion

Concrete is a major construction material worldwide and plays a crucial role in the nuclear industry. The elastic properties of concrete are prone to change and degrade while in service, as it is often subjected to extreme operational and environmental conditions. An accurate evaluation of concrete's elastic properties is thus essential to ensure structural integrity and safety. This is especially true for concrete in nuclear power plants, where irradiation effects significantly impact concrete mechanical properties. There are various methods to assess these properties, with ultrasound-based techniques showing high potential due to their nondestructive nature, cost-effectiveness, and safety. While several nondestructive evaluation methods exist, most rely on idealizations such as assuming homogeneous material and plane wavefronts. In this work, we address these issues by introducing an ultrasound-based nondestructive method aimed at reconstructing spatially varying images of concrete mechanical properties. By accurately modeling wave physics, including scattering and reflection, we overcome several of the aforementioned idealizations and aim to utilize the full waveform for imaging material properties through depth, resulting in more reliable images. Full waveform inversion (FWI) was first introduced by geophysicists to reconstruct subsurface elastic property images. The goal is to minimize the difference between simulated and recorded wavefield signals, often through gradient-based optimization algorithms. While FWI is primarily conducted using the acoustic approximation of the wave equation, few works focus on elastic FWI, where the goal is to reconstruct images of not only the pressure wave speed but also the shear wave speed and density (or their equivalents). This work explores the potential of using elastic FWI to predict concrete mechanical properties as an initial effort for a more accurate monitoring of concrete conditions in service. Reconstructing images of different elastic parameters enables more specificity and accurate condition diagnosis. This paper will detail this approach and provide examples demonstrating the effectiveness of elastic FWI in reconstructing comprehensive maps of concrete mechanical properties.

42 - ENGINEERING↗

Computer program for analysis of coupled-cavity traveling wave tubes

A flexible, accurate, large signal computer program was developed for the design of coupled cavity traveling wave tubes. The program is written in FORTRAN IV for an IBM 360/67 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 geometrical or electrical parameters different from those of 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 an approach in which the radio frequency fields are expanded in solutions to the transverse magnetic wave equation. All significant space harmonics are retained. The program was used to perform a design study of the traveling-wave tube developed for the Communications Technology Satellite. Good agreement was obtained between the predictions of the program and the measured performance of the flight tube.

Connolly, D. J.↗

Discretized energy minimization in a wave guide with point sources

An anti-noise problem on a finite time interval is solved by minimization of a quadratic functional on the Hilbert space of square integrable controls. To this end, the one-dimensional wave equation with point sources and pointwise reflecting boundary conditions is decomposed into a system for the two propagating components of waves. Wellposedness of this system is proved for a class of data that includes piecewise linear initial conditions and piecewise constant forcing functions. It is shown that for such data the optimal piecewise constant control is the solution of a sparse linear system. Methods for its computational treatment are presented as well as examples of their applicability. The convergence of discrete approximations to the general optimization problem is demonstrated by finite element methods.

Propst, G.↗

An Analytical Comparison of the Acoustic Analogy and Kirchhoff Formulation for Moving Surfaces

The Lighthill acoustic analogy, as embodied in the Ffowcs Williams-Hawkings (FW-H) equation, is compared with the Kirchhoff formulation for moving surfaces. A comparison of the two governing equations reveals that the main Kirchhoff advantage (namely nonlinear flow effects are included in the surface integration) is also available to the FW-H method if the integration surface used in the FW-H equation is not assumed impenetrable. The FW-H equation is analytically superior for aeroacoustics because it is based upon the conservation laws of fluid mechanics rather than the wave equation. This means that the FW-H equation is valid even if the integration surface is in the nonlinear region. This is demonstrated numerically in the paper. The Kirchhoff approach can lead to substantial errors if the integration surface is not positioned in the linear region. These errors may be hard to identify. Finally, new metrics based on the Sobolev norm are introduced which may be used to compare input data for both quadrupole noise calculations and Kirchhoff noise predictions.

Brentner, Kenneth S.↗

A Study of Supersonic Surface Sources: The Ffowcs Williams-Hawkings Equation and the Kirchhoff Formula

In this paper we address the mathematical problem of noise generation from high speed moving surfaces. The problem we are solving is the linear wave equation with sources on a moving surface. The Ffowcs Williams-Hawkings (FW-H) equation as well as the govern- ing equation for deriving the Kirchhoff formula for moving surfaces are both this type of partial differential equation. We give a new exact solution of this problem here in closed form which is valid for subsonic and supersonic motion of the surface but it is particularly suitable for supersonically moving surfaces. This new solution is the simplest of all high speed formulations of Langley and is denoted formulation 4 following the tradition of numbering of our major results for the prediction of the noise of rotating blades. We show that for a smooth surface moving at supersonic speed, our solution has only removable singularities. Thus it can be used for numerical work.

Farassat, F.↗

A strongly implicit procedure for the compressible Navier-Stokes equations

The strongly implicit procedure of Stone is used to obtain numerical solutions of the compressible Navier-Stokes equations in conservative form. In contrast to the spatially split Douglas-Gunn type methods, the method is shown to be numerically stable for the three-dimensional wave equation. The method is applied to a variety of external and internal two-dimensional flow problems involving shock wave boundary-layer interaction for both laminar and turbulent flows. The results are in good agreement with other methods and/or experiments. The storage penalty associated with the method is discussed and a simple, yet effective, means of minimizing the problem is presented.

Walters, R. W.↗

Subsonic and Supersonic Effects in Bose-Einstein Condensate

A paper presents a theoretical investigation of subsonic and supersonic effects in a Bose-Einstein condensate (BEC). The BEC is represented by a time-dependent, nonlinear Schroedinger equation that includes terms for an external confining potential term and a weak interatomic repulsive potential proportional to the number density of atoms. From this model are derived Madelung equations, which relate the quantum phase with the number density, and which are used to represent excitations propagating through the BEC. These equations are shown to be analogous to the classical equations of flow of an inviscid, compressible fluid characterized by a speed of sound (g/Po)1/2, where g is the coefficient of the repulsive potential and Po is the unperturbed mass density of the BEC. The equations are used to study the effects of a region of perturbation moving through the BEC. The excitations created by a perturbation moving at subsonic speed are found to be described by a Laplace equation and to propagate at infinite speed. For a supersonically moving perturbation, the excitations are found to be described by a wave equation and to propagate at finite speed inside a Mach cone.

Zak, Michail↗

A study of shock-associated magnetohydrodynamic waves in the solar wind

Three major topics were addressed, one theoretical and two observational. The topics were: (1) an attempt to understand the evolution of the large-amplitude magnetohydrodynamic (MHD) waves in the foreshock, using a nonlinear wave equation called the Derivative Nonlinear Schrodinger equation (henceforth DNLS) as a model, (2) using the extensive set of ISE data to test for the presence of various nonlinear wave processes which might be present, and (3) a study of plasma turbulence in the interstellar medium which might be physically similar to that in the solar wind. For these investigations we used radioastronomical techniques. Good progress was made in each of these areas and a separate discussion of each is given.

Spangler, Steven R.↗

Two-step perfectly matched layer for arbitrary-order pseudo-spectral analytical time-domain methods

Numerical simulation of an electrodynamic system in empty space requires the implementation of open boundary conditions (BC) to terminate the solution of Maxwell’s equations on the boundaries of the computational domain. The Perfectly Matched Layer (PML) has become the method of choice for open BC with wave equations, as it is straightforward and relatively easy to implement, and offers very efficient and user-adjustable absorption rates. PMLs are most often employed with the Finite-Difference Time-Domain (FDTD) algorithm, which in its most common implementation offers second-order accuracy in space and time on Cartesian grids. Yet, simulations (including some class of electromagnetic Particle-In-Cell simulations) that require higher precision may resort to higher-order Maxwell solvers employing extended finite-difference stencils, or even to pseudo-spectral Maxwell solvers, for which a general, versatile and efficient formulation of the PML has been missing so far. Here we propose a novel “two-step” formulation of the PML that is simple, very versatile and can be used as is with any Maxwell solver. In particular, it is applicable to a large class of Maxwell solvers including the arbitrary-order Pseudo-Spectral Analytical Time-Domain (PSATD) solver, which offers arbitrarily low numerical dispersion when increasing solver order and becomes dispersion-free at infinite order. Analysis and numerical simulations demonstrate that the new formulation is as efficient as the standard PML formulation, both for the FDTD and the PSATD implementations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Active stability augmentation of large space structures: A stochastic control problem

A problem in SCOLE is that of slewing an offset antenna on a long flexible beam-like truss attached to the space shuttle, with rather stringent pointing accuracy requirements. The relevant methodology aspects in robust feedback-control design for stability augmentation of the beam using on-board sensors is examined. It is framed as a stochastic control problem, boundary control of a distributed parameter system described by partial differential equations. While the framework is mathematical, the emphasis is still on an engineering solution. An abstract mathematical formulation is developed as a nonlinear wave equation in a Hilbert space. That the system is controllable is shown and a feedback control law that is robust in the sense that it does not require quantitative knowledge of system parameters is developed. The stochastic control problem that arises in instrumenting this law using appropriate sensors is treated. Using an engineering first approximation which is valid for small damping, formulas for optimal choice of the control gain are developed.

Balakrishnan, A. V.↗

A Mathematical Formulation of the SCOLE Control Problem. Part 2: Optimal Compensator Design

The study initiated in Part 1 of this report is concluded and optimal feedback control (compensator) design for stability augmentation is considered, following the mathematical formulation developed in Part 1. Co-located (rate) sensors and (force and moment) actuators are assumed, and allowing for both sensor and actuator noise, stabilization is formulated as a stochastic regulator problem. Specializing the general theory developed by the author, a complete, closed form solution (believed to be new with this report) is obtained, taking advantage of the fact that the inherent structural damping is light. In particular, it is possible to solve in closed form the associated infinite-dimensional steady-state Riccati equations. The SCOLE model involves associated partial differential equations in a single space variable, but the compensator design theory developed is far more general since it is given in the abstract wave equation formulation. The results thus hold for any multibody system so long as the basic model is linear.

Balakrishnan, A. V.↗

Long time behavior of unsteady flow computations

This paper addresses a specific issue of time accuracy in the calculation of external aerodynamic problems. The class of problems discussed consists of inviscid compressible subsonic flows. These problems are governed by a convective equation. A key issue that is not understood is the long time behavior of the solution. This is important if one desires transient calculations of problems governed by the Euler equations or its derivatives such as the small disturbance equations or the potential formulations for the gust problem. Difficulties arise for two dimensional problems where the time rate decay solutions of the wave equation is slow. In concert with the above mentioned problem, exterior flows require proper modeling of the boundary conditions. In particular, this requires the truncation of infinite regions into finite regions with the aid of artificial boundaries. These boundary conditions must be consistent with the physics of the unbounded problem as well as consistent in time and space. Our treatment of the problem is discussed in detail and examples are given to verify the results.

Hariharan, S. I.↗

Radial energy transport by magnetospheric ULF waves: Effects of magnetic curvature and plasma pressure

The 'radial' transport of energy by internal ULF waves, stimulated by dayside magnetospheric boundary oscillations, is analyzed in the framework of one-fluid magnetohydrodynamics. (the term radial is used here to denote the direction orthogonal to geomagnetic flux surfaces.) The model for the inhomogeneous magnetospheric plasma and background magnetic field is axisymmetric and includes radial and parallel variations in the magnetic field, magnetic curvature, plasma density, and low but finite plasma pressure. The radial mode structure of the coupled fast and intermediate MHD waves is determined by numerical solution of the inhomogeneous wave equation; the parallel mode structure is characterized by a Wentzel-Kramer-Brillouin (WKB) approximation. Ionospheric dissipation is modeled by allowing the parallel wave number to be complex. For boudnary oscillations with frequencies in the range from 10 to 48 mHz, and using a dipole model for the background magnetic field, the combined effects of magnetic curvature and finite plasma pressure are shown to (1) enhance the amplitude of field line resonances by as much as a factor of 2 relative to values obtained in a cold plasma or box-model approximation for the dayside magnetosphere; (2) increase the energy flux delivered to a given resonance by a factor of 2-4; and (3) broaden the spectral width of the resonance by a factor of 2-3. The effects are attributed to the existence of an 'Alfven buoyancy oscillation,' which approaches the usual shear mode Alfven wave at resonance, but unlike the shear Alfven mode, it is dispersive at short perpendicular wavelengths. The form of dispersion is analogous to that of an internal atmospheric gravity wave, with the magnetic tension of the curved background field providing the restoring force and allowing radial propagation of the mode. For nominal dayside parameters, the propagation band of the Alfven buoyancy wave occurs between the location of its (field line) resonance and that of the fast mode cutoff that exists at larger radial distances.

Kouznetsov, Igor↗

On the equations governing the gravitational perturbations of the Kerr black hole

Teukolsky's radial equation governing the general, non-axisymmetric, gravitational perturbations of the Kerr black hole is reduced to the form of a one-dimensional wave equation by making use of the transformation which enables the treatment of the non-axisymmetric modes in exactly the same way as the axisymmetric modes.

Chandrasekhar, S.↗

Rigorous solutions for sound radiation from circular ducts with hyperbolic horns or infinite plane baffle

A rigorous treatment is presented of sound radiation from circular ducts with either a hyperbolic horn or an infinite plane baffle. In the analysis hyperboloidal wave functions are used, which are defined here, for the first time, as a class of eigensolutions of the wave equation for oblate spheroidal co-ordinates. The numerical results include the complex conversion (or reflection) coefficients and the radiation directivity for various incident wave modes, spinning modes as well as axisymmetric modes. The solutions are valid for the whole frequency range including frequencies above and below the cut-off frequencies of the duct modes involved.

Cho, Y. C.↗

Reduction of the equation for lower hybrid waves in a plasma to a nonlinear Schroedinger equation

Equations describing the nonlinear propagation of waves in an anisotropic plasma are rarely exactly soluble. However it is often possible to make approximations that reduce the exact equations into a simpler equation. The use of MACSYMA to make such approximations, and so reduce the equation describing lower hybrid waves into the nonlinear Schrodinger equation which is soluble by the inverse scattering method is demonstrated. MACSYMA is used at several stages in the calculation only because there is a natural division between calculations that are easiest done by hand, and those that are easiest done by machine.

Karney, C. F. F.↗

Three-Dimensional Radiative-Transfer Equation

Progress made toward interpretation of radiometric observations. Paper discusses equation of radiative transfer in three-dimensional, inhomogeneous, scattering medium illuminated from above and bounded below by laterallyinhomogeneous, reflective plane. Representation of radiation field with full three-dimensional variability derived by use of spatial Fourier transform and matrix-operator techniques developed previously for one-dimensional version of problem. Equations useful for radiometric measurements from aircraft and spacecraft. Although derivations and resulting equations complicated, use of Fourier-transform, matrix-operator approach to solve practical problems simpler than direct solution of complete three-dimensional, linear wave equations.

Martonchik, J. V.↗