Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “generalized linear model”

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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
NASA NTRS2013

Wavelength Dependence of Solar Irradiance Enhancement During X-class Flares and Its Influence on the Upper Atmosphere

The wavelength dependence of solar irradiance enhancement during flare events is one of the important factors in determining how the Thermosphere-Ionosphere (TI) system responds to flares. To investigate the wavelength dependence of flare enhancement, the Flare Irradiance Spectral Model (FISM) was run for 61X-class flares. The absolute and the percentage increases of solar irradiance at flare peaks, compared to pre-flare conditions, have clear wavelength dependences. The 0-4 nm irradiance increases much more ((is) approximately 680 on average) than that in the 14-25 nm waveband ((is) approximately 65 on average), except at 24 nm ( (is) approximately 220). The average percentage increases for the 25-105 nm and 122-190 nm wave bands are approximately 120 and approximately 35, respectively. The influence of 6 different wavebands (0-14 nm, 14-25 nm, 25-105 nm, 105-120 nm, 121.56 nm,and122-175 nm) on the thermosphere was examined for the October 28th, 2003 flare (X17-class) event by coupling FISM with the National Center for Atmospheric Research (NCAR) Thermosphere-Ionosphere-Electrodynamics General Circulation Model(TIE-GCM) under geomagnetically quiet conditions (Kp=1). While the enhancement in the0-14nm waveband caused the largest enhancement of the globally integrated solar heating, the impact of solar irradiance enhancement on the thermosphere at 400 km is largest for the 25-105 nm waveband (EUV), which accounts for about 33 K of the total 45 K temperature enhancement, and approximately 7.4% of the total approximately 11.5% neutral density enhancement. The effect of 122-175 nm flare radiation on the thermosphere is rather small. The study also illustrates that the high-altitude thermospheric response to the flare radiation at 0-175 nm is almost a linear combination of the responses to the individual wavebands. The upper thermospheric temperature and density enhancements peaked 3-5 h after the maximum flare radiation.

FISM↗
NASA NTRS2013

Wavelength Dependence of Solar Irradiance Enhancement During X-Class Flares and Its Influence on the Upper Atmosphere

The wavelength dependence of solar irradiance enhancement during flare events is one of the important factors in determining how the Thermosphere-Ionosphere (T-I) system responds to flares. To investigate the wavelength dependence of flare enhancement, the Flare Irradiance Spectral Model (FISM) was run for 61 X-class flares. The absolute and the percentage increases of solar irradiance at flare peaks, compared to pre-flare conditions, have clear wavelength dependences. The 0-14 nm irradiance increases much more (approx. 680% on average) than that in the 14-25 nm waveband (approx. 65% on average), except at 24 nm (approx. 220%). The average percentage increases for the 25-105 nm and 122-190 nm wavebands are approx. 120% and approx. 35%, respectively. The influence of 6 different wavebands (0-14 nm, 14-25 nm, 25-105 nm, 105- 120 nm, 121.56 nm, and 122-175 nm) on the thermosphere was examined for the October 28th, 2003 flare (X17-class) event by coupling FISM with the National Center for Atmospheric Research (NCAR) Thermosphere-Ionosphere-Electrodynamics General Circulation Model (TIE-GCM) under geomagnetically quiet conditions (Kp=1). While the enhancement in the 0-14 nm waveband caused the largest enhancement of the globally integrated solar heating, the impact of solar irradiance enhancement on the thermosphere at 400 km is largest for the 25-105 nm waveband (EUV), which accounts for about 33 K of the total 45 K temperature enhancement, and approx. 7.4% of the total approx. 11.5% neutral density enhancement. The effect of 122-175 nm flare radiation on the thermosphere is rather small. The study also illustrates that the high-altitude thermospheric response to the flare radiation at 0-175 nm is almost a linear combination of the responses to the individual wavebands. The upper thermospheric temperature and density enhancements peaked 3-5 h after the maximum flare radiation.

TIE-GCM↗
NASA NTRS1979

A higher order panel method for linearized supersonic flow

The basic integral equations of linearized supersonic theory for an advanced supersonic panel method are derived. Methods using only linear varying source strength over each panel or only quadratic doublet strength over each panel gave good agreement with analytic solutions over cones and zero thickness cambered wings. For three dimensional bodies and wings of general shape, combined source and doublet panels with interior boundary conditions to eliminate the internal perturbations lead to a stable method providing good agreement experiment. A panel system with all edges contiguous resulted from dividing the basic four point non-planar panel into eight triangular subpanels, and the doublet strength was made continuous at all edges by a quadratic distribution over each subpanel. Superinclined panels were developed and tested on s simple nacelle and on an airplane model having engine inlets, with excellent results.

Ehlers, F. E.↗
NASA NTRS1992

Design and analysis of low boom concepts at Langley Research Center

The objective of the sonic boom research in the current High Speed Research Program is to ultimately make possible overland supersonic flight by a high speed civil transport. To accomplish this objective, it is felt that results in four areas must demonstrate that such a vehicle would be acceptable by the general public, by the airframers, and by the airlines. It should be demonstrated: (1) that some waveform shape has the possibility of being acceptable to the general public; (2) that the atmosphere would not totally destroy such a waveform during propagation; (3) that a viable airplane could be built which produces such a waveform; and (4) that any performance penalty suffered by a low boom aircraft would be counteracted by the economic benefit of overland supersonic flight. The work being done at LaRC is in support of the third element listed above--the area of configuration design. The initial part of the paper will give a review of the theory being used for configuration designs and discuss two theory validation models which were built and tested within the past two years. Discussion of the wind tunnel and theoretical results (linear theory and higher order methods) and their implications for future designs will be included.

Darden, Christine M.↗
NASA NTRS2012

Variability Selected Low-Luminosity Active Galactic Nuclei in the 4 Ms Chandra Deep Field-South

The 4 Ms Chandra Deep Field-South (CDF-S) and other deep X-ray surveys have been highly effective at selecting active galactic nuclei (AGN). However, cosmologically distant low-luminosity AGN (LLAGN) have remained a challenge to identify due to significant contribution from the host galaxy. We identify long-term X ray variability (approx. month years, observed frame) in 20 of 92 CDF-S galaxies spanning redshifts approx equals 00.8 - 1.02 that do not meet other AGN selection criteria. We show that the observed variability cannot be explained by X-ray binary populations or ultraluminous X-ray sources, so the variability is most likely caused by accretion onto a supermassive black hole. The variable galaxies are not heavily obscured in general, with a stacked effective power-law photon index of Gamma(sub Stack) approx equals 1.93 +/- 0.13, and arc therefore likely LLAGN. The LLAGN tend to lie it factor of approx equal 6-89 below the extrapolated linear variability-luminosity relation measured for luminous AGN. This may he explained by their lower accretion rates. Variability-independent black-hole mass and accretion-rate estimates for variable galaxies show that they sample a significantly different black hole mass-accretion-rate space, with masses a factor of 2.4 lower and accretion rates a factor of 22.5 lower than variable luminous AGNs at the same redshift. We find that an empirical model based on a universal broken power-law power spectral density function, where the break frequency depends on SMBH mass and accretion rate, roughly reproduces the shape, but not the normalization, of the variability-luminosity trends measured for variable galaxies and more luminous AGNs.

Young, M.↗
NASA NTRS1994

Discrete-Layer Piezoelectric Plate and Shell Models for Active Tip-Clearance Control

The objectives of this work were to develop computational tools for the analysis of active-sensory composite structures with added or embedded piezoelectric layers. The targeted application for this class of smart composite laminates and the analytical development is the accomplishment of active tip-clearance control in turbomachinery components. Two distinct theories and analytical models were developed and explored under this contract: (1) a discrete-layer plate theory and corresponding computational models, and (2) a three dimensional general discrete-layer element generated in curvilinear coordinates for modeling laminated composite piezoelectric shells. Both models were developed from the complete electromechanical constitutive relations of piezoelectric materials, and incorporate both displacements and potentials as state variables. This report describes the development and results of these models. The discrete-layer theories imply that the displacement field and electrostatic potential through-the-thickness of the laminate are described over an individual layer rather than as a smeared function over the thickness of the entire plate or shell thickness. This is especially crucial for composites with embedded piezoelectric layers, as the actuating and sensing elements within these layers are poorly represented by effective or smeared properties. Linear Lagrange interpolation polynomials were used to describe the through-thickness laminate behavior. Both analytic and finite element approximations were used in the plane or surface of the structure. In this context, theoretical developments are presented for the discrete-layer plate theory, the discrete-layer shell theory, and the formulation of an exact solution for simply-supported piezoelectric plates. Finally, evaluations and results from a number of separate examples are presented for the static and dynamic analysis of the plate geometry. Comparisons between the different approaches are provided when possible, and initial conclusions regarding the accuracy and limitations of these models are given.

Heyliger, P. R.↗
NASA NTRS2003

Kalman Filtering with Inequality Constraints for Turbofan Engine Health Estimation

Kalman filters are often used to estimate the state variables of a dynamic system. However, in the application of Kalman filters some known signal information is often either ignored or dealt with heuristically. For instance, state variable constraints (which may be based on physical considerations) are often neglected because they do not fit easily into the structure of the Kalman filter. This paper develops two analytic methods of incorporating state variable inequality constraints in the Kalman filter. The first method is a general technique of using hard constraints to enforce inequalities on the state variable estimates. The resultant filter is a combination of a standard Kalman filter and a quadratic programming problem. The second method uses soft constraints to estimate state variables that are known to vary slowly with time. (Soft constraints are constraints that are required to be approximately satisfied rather than exactly satisfied.) The incorporation of state variable constraints increases the computational effort of the filter but significantly improves its estimation accuracy. The improvement is proven theoretically and shown via simulation results. The use of the algorithm is demonstrated on a linearized simulation of a turbofan engine to estimate health parameters. The turbofan engine model contains 16 state variables, 12 measurements, and 8 component health parameters. It is shown that the new algorithms provide improved performance in this example over unconstrained Kalman filtering.

Simon, Dan↗
NASA NTRS1984

A comparison of the properties of highly polarized QSOs versus low-polarization QSOs

An optical linear-polarization survey testing the relationships between polarization and other properties of QSOs is presented. Polarimetric observations of 239 QSOs have been made. Data from the literature are used to analyze the relationships between polarization and radio, optical, and X-ray properties. The highly polarized QSOs are generally compact radio sources, are associated with radio properties such as low-frequency variability and superluminal motion, exhibit large-amplitude rapid photometric variability, have steep nonthermal optical continua, and may show excess X-ray emission. No relationship is seen between polarization and redshift, optical luminosity, or the equivalent width of emission lines. The results are discussed in the context of both isotropic and anisotropic (beaming) models. While the associations between the optical and radio properties of highly polarized QSOs provide strong motivation for the anisotropic model, the lack of associations with redshift, optical luminosity, and emission-line strength is inconsistent with this model. It is concluded that the highly polarized QSOs are probably fundamentally different in their physical properties from low-polarization QSOs.

Moore, R. L.↗
NASA NTRS1993

Spectrum formation in supernovae - Numerical techniques

The study combines several novel techniques for spectrum simulation in the Eddington computer program which solves the comoving frame equation of transfer coupled with the statistical and radiative equilibrium equations. One of these is a generalization of the accelerated lambda iteration (ALI) scheme to include an approximate frequency-derivative operator. This greatly enhances the convergence rate of ALI in optically thick, high-velocity shear flows. Another is a partial linearization technique which is capable of efficiently solving a very large number of rate equations on a moderately sized computer. An expansion opacity and emissivity approximation is derived which makes it possible to determine the effect on the transfer and statistical equilibrium of a very large number of lines not explicitly represented in the frequency grid and additionally to treat line-blanketing from species not explicitly included in the rate equations. The utility of these techniques is illustrated with models of two supernovae.

Eastman, Ronald G.↗
NASA NTRS1998

Inversion of Crater Morphometric Data to Gain Insight on the Cratering Process

In recent years, morphometric data for Venus and several outer planet satellites have been collected, so we now have observational data of complex Craters formed in a large range of target properties. We present general inversion techniques that can utilize the morphometric data to quantitatively test various models of complex crater formation. The morphometric data we use in this paper are depth of a complex crater, the diameter at which the depth-diameter ratio changes, and onset diameters for central peaks, terraces, and peak rings. We tested the roles of impactor velocities and hydrostatic pressure vs. crustal strength, and we tested the specific models of acoustic fluidization (Melosh, 1982) and nonproportional growth (Schultz, 1988). Neither the acoustic fluidization model nor the nonproportional growth in their published formulations are able to successfully reproduce the data. No dependence on impactor velocity is evident from our inversions. Most of the morphometric data is consistent with a linear dependence on the ratio of crustal strength to hydrostatic pressure on a planet, or the factor c/pg.

Herrick, Robert R.↗
NASA NTRS1998

Efficient and Optimal Attitude Determination Using Recursive Global Positioning System Signal Operations

In this paper, a new and efficient algorithm is developed for attitude determination from Global Positioning System signals. The new algorithm is derived from a generalized nonlinear predictive filter for nonlinear systems. This uses a one time-step ahead approach to propagate a simple kinematics model for attitude determination. The advantages of the new algorithm over previously developed methods include: it provides optimal attitudes even for coplanar baseline configurations; it guarantees convergence even for poor initial conditions; it is a non-iterative algorithm; and it is computationally efficient. These advantages clearly make the new algorithm well suited to on-board applications. The performance of the new algorithm is tested on a dynamic hardware simulator. Results indicate that the new algorithm accurately estimates the attitude of a moving vehicle, and provides robust attitude estimates even when other methods, such as a linearized least-squares approach, fail due to poor initial starting conditions.

Crassidis, John L.↗
NASA NTRS2009

Massively Parallel Dantzig-Wolfe Decomposition Applied to Traffic Flow Scheduling

Optimal scheduling of air traffic over the entire National Airspace System is a computationally difficult task. To speed computation, Dantzig-Wolfe decomposition is applied to a known linear integer programming approach for assigning delays to flights. The optimization model is proven to have the block-angular structure necessary for Dantzig-Wolfe decomposition. The subproblems for this decomposition are solved in parallel via independent computation threads. Experimental evidence suggests that as the number of subproblems/threads increases (and their respective sizes decrease), the solution quality, convergence, and runtime improve. A demonstration of this is provided by using one flight per subproblem, which is the finest possible decomposition. This results in thousands of subproblems and associated computation threads. This massively parallel approach is compared to one with few threads and to standard (non-decomposed) approaches in terms of solution quality and runtime. Since this method generally provides a non-integral (relaxed) solution to the original optimization problem, two heuristics are developed to generate an integral solution. Dantzig-Wolfe followed by these heuristics can provide a near-optimal (sometimes optimal) solution to the original problem hundreds of times faster than standard (non-decomposed) approaches. In addition, when massive decomposition is employed, the solution is shown to be more likely integral, which obviates the need for an integerization step. These results indicate that nationwide, real-time, high fidelity, optimal traffic flow scheduling is achievable for (at least) 3 hour planning horizons.

Rios, Joseph Lucio↗
NASA NTRS2022

X-59 Sonic Boom Test Results from the NASA Glenn 8- by 6-Foot Supersonic Wind Tunnel

A wind tunnel test was conducted to investigate near-field sonic boom pressure signatures from a model of the X-59 Low-Boom Flight Demonstrator aircraft. A 1.62%-scale model of the aircraft in the C612A configuration was fabricated for the wind tunnel test, which took place in the NASA Glenn 8- by 6-Foot Supersonic Wind Tunnel in September and October 2021. The model had provisions for two different mounting options: a swept blade strut that attached at the top of model ahead of the inlet, and rear-entry sting that was made as one piece with a dummy nacelle, and which had a 2”-long cylindrical segment aft of the nozzle exit before tapering up in size. The blade strut allowed for a clean aft end of the model for evaluation of the shocks from that region, while the sting avoided the significant distortions of the flow and shocks from the blade strut along the top of the model. Both the sting and the strut had adapters that attached to a force balance. The model had alternate parts for ±0.5° deflections of the flaps, ailerons, and stabilator, and ±1° deflections of the T-tail horizontal surface. Off-body static pressure measurements of the flow field below the model were made by use of a pressure rail which had 420 orifices along its tip. The model was positioned at various heights from the rail by vertical movement of the wind tunnel strut, and at various longitudinal stations relative to the rail by means of a linear actuator mounted between the tunnel strut and the balance. The longitudinal positioning allowed multiple pressure signatures to be obtained along different portions of the rail. These signatures were aligned by accounting for the model longitudinal movement and then averaged to take out the effects of tunnel flow distortions and the interference of the rail flow field and shocks on the model pressure signatures. The test was run at approximate Mach numbers of 1.36, 1.4, and 1.47, and the model was set at various angles of attack and roll relative to the rail. Plots of the model signatures for all the variations of Mach number, model angles, control deflections, and height relative to the rail are provided throughout the report. Repeatability was generally very good and gave confidence in the quality of the measurements. The signatures measured at various heights from the rail provided insight into the aging of the model shocks as they propagated from 1.2 to 3 body lengths from the model. Off-track signatures up to 45° from centerline obtained by rolling the model gave indications of the shock flow fields across the width of the sonic boom carpet. The deflections of the various control surfaces allowed assessment of the boom sensitivity to the control surface movements.

Sonic boom↗
NASA NTRS2012

NASA Tech Briefs, September 2012

Topics covered include: Beat-to-Beat Blood Pressure Monitor; Measurement Techniques for Clock Jitter; Lightweight, Miniature Inertial Measurement System; Optical Density Analysis of X-Rays Utilizing Calibration Tooling to Estimate Thickness of Parts; Fuel Cell/Electrochemical Cell Voltage Monitor; Anomaly Detection Techniques with Real Test Data from a Spinning Turbine Engine-Like Rotor; Measuring Air Leaks into the Vacuum Space of Large Liquid Hydrogen Tanks; Antenna Calibration and Measurement Equipment; Glass Solder Approach for Robust, Low-Loss, Fiber-to-Waveguide Coupling; Lightweight Metal Matrix Composite Segmented for Manufacturing High-Precision Mirrors; Plasma Treatment to Remove Carbon from Indium UV Filters; Telerobotics Workstation (TRWS) for Deep Space Habitats; Single-Pole Double-Throw MMIC Switches for a Microwave Radiometer; On Shaft Data Acquisition System (OSDAS); ASIC Readout Circuit Architecture for Large Geiger Photodiode Arrays; Flexible Architecture for FPGAs in Embedded Systems; Polyurea-Based Aerogel Monoliths and Composites; Resin-Impregnated Carbon Ablator: A New Ablative Material for Hyperbolic Entry Speeds; Self-Cleaning Particulate Prefilter Media; Modular, Rapid Propellant Loading System/Cryogenic Testbed; Compact, Low-Force, Low-Noise Linear Actuator; Loop Heat Pipe with Thermal Control Valve as a Variable Thermal Link; Process for Measuring Over-Center Distances; Hands-Free Transcranial Color Doppler Probe; Improving Balance Function Using Low Levels of Electrical Stimulation of the Balance Organs; Developing Physiologic Models for Emergency Medical Procedures Under Microgravity; PMA-Linked Fluorescence for Rapid Detection of Viable Bacterial Endospores; Portable Intravenous Fluid Production Device for Ground Use; Adaptation of a Filter Assembly to Assess Microbial Bioburden of Pressurant Within a Propulsion System; Multiplexed Force and Deflection Sensing Shell Membranes for Robotic Manipulators; Whispering Gallery Mode Optomechanical Resonator; Vision-Aided Autonomous Landing and Ingress of Micro Aerial Vehicles; Self-Sealing Wet Chemistry Cell for Field Analysis; General MACOS Interface for Modeling and Analysis for Controlled Optical Systems; Mars Technology Rover with Arm-Mounted Percussive Coring Tool, Microimager, and Sample-Handling Encapsulation Containerization Subsystem; Fault-Tolerant, Real-Time, Multi-Core Computer System; Water Detection Based on Object Reflections; SATPLOT for Analysis of SECCHI Heliospheric Imager Data; Plug-in Plan Tool v3.0.3.1; Frequency Correction for MIRO Chirp Transformation Spectroscopy Spectrum; Nonlinear Estimation Approach to Real-Time Georegistration from Aerial Images; Optimal Force Control of Vibro-Impact Systems for Autonomous Drilling Applications; Low-Cost Telemetry System for Small/Micro Satellites; Operator Interface and Control Software for the Reconfigurable Surface System Tri-ATHLETE; and Algorithms for Determining Physical Responses of Structures Under Load.

Source record↗
NASA NTRSDate not supplied

Performance and Portability of a Linear Solver Across Emerging Architectures

A linear solver algorithm used by a large-scale unstructured-grid computational fluid dynamics application is examined for a broad range of familiar and emerging architectures. Efficient implementation of a linear solver is challenging on recent CPUs offering vector architectures. Vector loads and stores are essential to effectively utilize available memory bandwidth on CPUs, and maintaining performance across different CPUs can be difficult in the face of varying vector lengths offered by each. A similar challenge occurs on GPU architectures, where it is essential to have coalesced memory accesses to utilize memory bandwidth effectively. In this work, we demonstrate that restructuring a computation, and possibly data layout, with regard to architecture is essential to achieve optimal performance by establishing a performance benchmark for each target architecture in a low level language such as vector intrinsics or CUDA. In doing so, we demonstrate how a linear solver kernel can be mapped to Intel® Xeon™ and Xeon Phi™, Marvell® ThunderX2®, NEC® SX-Aurora™ TSUBASA Vector Engine, and NVIDIA® and AMD® GPUs. We further demonstrate that the required code restructuring can be achieved in higher level programming environments such as OpenACC, OCCA, and Intel® OneAPI™/SYCL, and that each generally results in optimal performance on the target architecture. Relative performance metrics for all implementations are shown, and subjective ratings for ease of implementation and optimization are suggested.

Programming models↗
NASA NTRS1977

A two-scale scattering model with application to the JONSWAP '75 aircraft microwave scatterometer experiment

The general problem of bistatic scattering from a two scale surface was evaluated. The treatment was entirely two-dimensional and in a vector formulation independent of any particular coordinate system. The two scale scattering model was then applied to backscattering from the sea surface. In particular, the model was used in conjunction with the JONSWAP 1975 aircraft scatterometer measurements to determine the sea surface's two scale roughness distributions, namely the probability density of the large scale surface slope and the capillary wavenumber spectrum. Best fits yield, on the average, a 0.7 dB rms difference between the model computations and the vertical polarization measurements of the normalized radar cross section. Correlations between the distribution parameters and the wind speed were established from linear, least squares regressions.

Wentz, F. J.↗
NASA NTRS1993

Computational mechanics needs study

In order to assess the needs in computational mechanics over the next decade, we formulated a questionnaire and contacted computational mechanics researchers and users in industry, government, and academia. As expected, we found a wide variety of computational mechanics usage and research. This report outlines the activity discussed with those contacts, as well as that in our own organizations. It should be noted that most of the contacts were made before the recent decline of the defense industry. Therefore, areas which are strongly defense-oriented may decrease in relative importance. In order to facilitate updating of this study, names of a few key researchers in each area are included as starting points for future literature surveys. These lists of names are not intended to represent those persons doing the best research in that area, nor are they intended to be comprehensive. They are, as previously stated, offered as starting points for future literature searches. Overall, there is currently a broad activity in computational mechanics in this country, with the breadth and depth increasing as more sophisticated software and faster computers become more available. The needs and desires of the workers in this field are as diverse as their background and organizational products. There seems to be some degree of software development in any organization (although the level of activity is highly variable from one organization to another) which has any research component in its mission. It seems, however, that there is considerable use of commercial software in almost all organizations. In most industrial research organizations, it appears that very little actual software development is contracted out, but that most is done in-house, using a mixture of funding sources. Government agencies vary widely in the ratio of in-house to out-house ratio. There is a considerable amount of experimental verification in most, but not all, organizations. Generally, the amount of experimental verification is more than we expected. Of all the survey contacts, one or two believe that the resources they are allocated are sufficient, but most do not. Some believe they have only half the resources they need. Some see their resource deficits as short-term, while others see it as a trend which will continue or perhaps worsen. The pessimism is stronger in the defense and aerospace industry. When considering only the nonlinear development efforts, there appears to be an even mix of geometric and material nonlinearity. There is not much particular emphasis in linear analysis unless it is in extension of current analysis capabilities to larger problems. The primary exception is concern about modeling of composites, where proven methodologies have trailed element and computer hardware development. Most of the people we spoke to use finite element techniques, but there is some finite difference and boundary element work ongoing. There is also some interest in multiple methods. Coupling of finite elements and boundary elements appears to be of high interest, since the two analysis types are complementary.

Griffin, O. Hayden, Jr.↗
NASA NTRS1996

Theoretical Basis for Finite Difference Extrapolation of Sonic Boom Signatures

Calculation of sonic boom signatures for aircraft has traditionally followed the methods of Whitham' and Walkden. The wave disturbance generated by the vehicle is obtained by area rule linearized supersonic flow methods, which yield a locally axisymmetric asymptotic solution. This solution is acoustic in nature, i.e., first order in disturbance quantities, and corresponds to ray acoustics. Cumulative nonlinear distortion of the signature is incorporated by using this solution to adjust propagation speed to first order, thus yielding a solution second order in disturbance quantities. The effects of atmospheric gradients are treated by Blokhintzov's method of geometrical acoustics. Both nonlinear signature evolution and ray tracing are applied as if the pressure field very close to the vehicle were actually that given by the source term (the 'F-function') of the asymptotic linearized flow solution. The viewpoint is thus that the flow solution exists at a small radius near the vehicle, and may be treated as an input to an extrapolation procedure consisting of ray tracing and nonlinear aging. The F-function is often regarded as a representation of a near-field pressure signature, and it is common for computational implementations to treat it interchangeably with the pressure signature. There is a 'matching radius' between the source function and the subsequent propagation extrapolation. This viewpoint has been supported by wind tunnel tests of simple models, and very typically yields correct results for actual flight vehicles. The assumption that the F-function and near-field signature are interchangeable is generally not correct. The flowfield of a vehicle which is not axisymmetric contains crossflow components which are very significant at small radii and less so at larger distances. From an acoustical viewpoint, the crossflow is equivalent to source diffraction portions of the wave field. Use of the F-function as a near field signature effectively assumes that the diminution of the crossflow/diffraction component may be applied all at once at the matching radius noted above. This approximation, though not rigorously validated, is responsible for the usual correct far-field results. On the other hand, if an actual near-field signature (either from wind tunnel or CFD data) is used at a starting point rather than one based on th effective source distribution, the predicted far-field signature is generally wrong.

Plotkin, Kenneth J.↗