Search NASA⌕ Search

SEARCH · Search NASA

Results for “smoothing methods”

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,279 records · Page 71

Quickbird Satellite in-orbit Modulation Transfer Function (MTF) Measurement Using Edge, Pulse and Impulse Methods for Summer 2003

The spatial characteristics of an imaging system cannot be expressed by a single number or simple statement. However, the Modulation Transfer Function (MTF) is one approach to measure the spatial quality of an imaging system. Basically, MTF is the normalized spatial frequency response of an imaging system. The frequency response of the system can be evaluated by applying an impulse input. The resulting impulse response is termed the Point Spread function (PSF). This function is a measure of the amount of blurring present in the imaging system and is itself a useful measure of spatial quality. An underlying assumption is that the imaging system is linear and shift-independent. The Fourier transform of the PSF is called the Optical Transfer Function (OTF) and the normalized magnitude of the OTF is the MTF. In addition to using an impulse input, a knife-edge in technique has also been used in this project. The sharp edge exercises an imaging system at all spatial frequencies. The profile of an edge response from an imaging system is called an Edge Spread Function (ESF). Differentiation of the ESF results in a one-dimensional version of the Point Spread Function (PSF). Finally, MTF can be calculated through use of Fourier transform of the PSF as stated previously. Every image includes noise in some degree which makes MTF of PSF estimation more difficult. To avoid the noise effects, many MTF estimation approaches use smooth numerical models. Historically, Gaussian models and Fermi functions were applied to reduce the random noise in the output profiles. The pulse-input method was used to measure the MTF of the Landsat Thematic Mapper (TM) using 8th order even functions over the San Mateo Bridge in San Francisco, California. Because the bridge width was smaller than the 30-meter ground sample distance (GSD) of the TM, the Nyquist frequency was located before the first zero-crossing point of the sinc function from the Fourier transformation of the bridge pulse. To avoid the zero-crossing points in the frequency domain from a pulse, the pulse width should be less than the width of two pixels (or 2 GSD's), but the short extent of the pulse results in a poor signal-to-noise ratio. Similarly, for a high-resolution satellite imaging system such as Quickbird, the input pulse width was critical because of the zero crossing points and noise present in the background area. It is important, therefore, that the width of the input pulse be appropriately sized. Finally, the MTF was calculated by taking ratio between Fourier transform of output and Fourier transform of input. Regardless of whether the edge, pulse and impulse target method is used, the orientation of the targets is critical in order to obtain uniformly spaced sub-pixel data points. When the orientation is incorrect, sample data points tend to be located in clusters that result in poor reconstruction of the edge or pulse profiles. Thus, a compromise orientation must be selected so that all spectral bands can be accommodated. This report continues by outlining the objectives in Section 2, procedures followed in Section 3, descriptions of the field campaigns in Section 4, results in Section 5, and a brief summary in Section 6.

Helder, Dennis↗

Efficient generalized boundary detection

Fast, efficient, and robust compression-based methods for detecting boundaries in arbitrary datasets, including sequences (1D datasets), are desired. The methods, each employing three simple algorithms, approximate the information distance between two adjacent sliding windows within a dataset. One of the algorithms calculates an initial ordered list of subsequences; while a second algorithm updates the ordered list of subsequences by dropping a first entry and appending a last entry rather than calculating completely new ordered lists with each iteration. Large values in the distance metric are indicative of boundary locations. A smoothed z-score or a wavelet-based algorithm may then be used to locate peaks in the distance metric, thereby identifying boundary locations. An adaptive version of the method employs a collection of window sizes and corresponding weighting functions, making it more amenable to real datasets with unknown, complex, and changing structures.

Ting, Christina↗

Method development for determining the hydraulic conductivity of fractured porous media

Plausible, but unvalidated, theoretical model constructs for unsaturated hydraulic conductivity of fractured porous media are currently used in Performance Assessment (PA) modeling for cracked saltstone and concrete (Flach 2011). The Nuclear Regulatory Commission (NRC) has expressed concern about the lack of model support for these assumed Moisture Characteristic Curves (MCC) data, as noted in Requests for Additional Information (RAIs) PA-8 and SP-4 (Savannah River Remediation, LLC, 2011). The objective of this task was to advance PA model support by developing an experimental method for determining the hydraulic conductivity of fractured cementitious materials under unsaturated conditions, and to demonstrate the technique on fractured saltstone samples. The task was requested through Task Technical Request (TTR) HLW-SSF-TTR-2012-0016 and conducted in accordance with Task Technical & Quality Assurance Plan (TTQAP) SRNL-TR-2012-00090. Preliminary method development previously conducted by Kohn et al. (2012) identified transient outflow extraction as the most promising method for characterizing the unsaturated properties of fractured porous media. While the research conducted by Kohn et al. (2012) focused on fractured media analogs such as stacked glass slides, the current task focused directly on fractured saltstone. For this task, four sample types with differing fracture geometries were considered: 1) intact saltstone, 2) intact saltstone with a single saw cut, smooth surface fracture, 3) micro-fractured saltstone (induced by oven drying), and 4) micro-fractured saltstone with a single, fully-penetrating, rough-surface fracture. Each sample type was tested initially for saturated hydraulic conductivity following method ASTM D 5084 using a flexible wall permeameter. Samples were subsequently tested using the transient outflow extraction method to determine cumulative outflow as a function of time and applied pressure. Of the four sample types tested, two yielded datasets suitable for analysis (sample types 3 and 4). The intact saltstone sample (sample type 1) did not yield any measureable outflow over the pressure range of the outflow test (0-1000 cm H 2 O). This was expected because the estimated air entry pressure for intact saltstone is on the order of 100,000 cm H 2 O (Dixon et al., 2009). The intact saltstone sample with a single saw cut smooth surface fracture (sample type 2) did not produce useable data because the fracture completely drained at less than 10 cm H 2 O applied pressure. The cumulative outflow data from sample types 3 and 4 were analyzed using an inverse solution of the Richard’s equation for water flow in variably saturated porous media. This technique was implemented using the computer code Hydrus-1D (Šimunek et al., 2008) and the resulting output included the van Genuchten-Mualem water retention and relative permeability parameters and predicted saturated hydraulic conductivity (Van Genuchten, 1980; Van Genuchten et al., 1991). Estimations of relative permeability and saturated conductivity are possible because the transient response of the sample to pressure changes is recorded during the multi-step outflow extraction test. Characteristic curves were developed for sample types 3 and 4 based on the results of the transient outflow method and compared to that of intact saltstone previously reported by Dixon et al. (2009). The overall results of this study indicate that the outflow extraction method is suitable for measuring the hydraulic properties of micro-fractured porous media. The resulting cumulative outflow data can be analyzed using the computer code Hydrus-1D to generate the van Genuchten curve fitting parameters that adequately describe fracture drainage. The resulting characteristic curves are consistent with blended characteristic curves that combine the behaviors of low pressure drainage associated with fracture flow with high pressure drainage from the bulk saltstone matrix.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Epitaxial two dimensional aluminum films on silicon (111) by ultra-fast thermal deposition

Aluminum thin films are known for their extremely rough surface, which is detrimental for applications such as molecular electronics and photonics, where protrusions cause electrical shorts or strong scattering. We achieved atomically flat Al films using a highly non-equilibrium approach. Ultra-fast thermal deposition (UFTD), at rates >10 nm/s, yields RMS roughness of 0.4 to 0.8 nm for 30-50 nm thick Al films on variety of substrates. For UFTD on Si(111) substrates, the top surface follows closely the substrate topography (etch pits), indicating a 2D, layer-by-layer growth. The Al film is a mixture of (100) and (111) grains, where the latter are commensurate with the in-plane orientation of the underlying Si (epitaxy). We show the use of these ultra-smooth Al films for highly reproducible charge-transport measurements across a monolayer of alkyl phosphonic acid as well as for plasmonics applications by directly patterning them by focused ion beam to form a long-range ordered array of holes. UFTD is a one-step process, with no need for annealing, peeling, or primer layers. It is conceptually opposite to high quality deposition methods, such as MBE or ALD, which are slow and near-equilibrium processes. For Al, though, we find that limited diffusion length (and good wetting) is critical for achieving ultra-smooth thin films.

36 MATERIALS SCIENCE↗

Involvement of prostatic interstitial cells of Cajal in inflammatory cytokines-elicited catecholamines production: Implications for the pathophysiology of chronic prostatitis/chronic pelvic pain syndrome (CP/CPPS)

Highlights: • Prostatic ICCs express enzymes essential for catecholamine synthesis. • Proinflammatory cytokines stimulate catecholamine synthesis in ICCs. • NF-κB, HIF-1α and HDACs pathways are required for ICCs function. In a previous work using guinea pig prostate, we have identified a novel interstitial cells of Cajal (ICCs) which possess close contacts between sympathetic nerve bundles and smooth muscle cells. The ability of prostatic ICCs in mediating excitatory neural inputs was therefore studied using isolated murine prostate ICCs by collagenase digestion combined with FACS method. RT-PCR and Western blotting analyses revealed that prostatic ICCs under a quiescent state expressed abundantly the rate-limiting enzymes essential for catecholamine synthesis. Moreover, distinct proinflammatory cytokines (e.g. IL-1β, IL-8, ICAM-1 and TNF-α) could significantly stimulate the expression levels of the rate-limiting enzymes of catecholamine production in prostate ICCs. Mechanistically, the above-mentioned stimulatory effects of proinflammatory cytokines appeared to be mediated via activation of NF-κB, HIF-1α and HDACs signaling pathways. Considering that prostatic catecholamine overactivity serves as an essential etiology of pelvic pain by indirectly stimulating the smooth muscle cell proliferation, or by directly causing muscular spasm, our results collectively suggest that targeting the NF-κB, HIF-1α and HDACs pathways in prostate ICCs be considered as a new strategy for treatment of chronic pelvic pain syndrome (CPPS) induced by chronic prostatitis (CP). Overall, the current study should shed novel light on the biology of this unique prostate ICCs.

60 APPLIED LIFE SCIENCES↗

Image filtering as an alternative to the application of a different reconstruction kernel in CT imaging: Feasibility study in lung cancer screening

Purpose: While the acquisition of projection data in a computed tomography (CT) scanner is generally carried out once, the projection data is often removed from the system, making further reconstruction with a different reconstruction filter impossible. The reconstruction kernel is one of the most important parameters. To have access to all the reconstructions, either prior reconstructions with multiple kernels must be performed or the projection data must be stored. Each of these requirements would increase the burden on data archiving. This study aimed to design an effective method to achieve similar image quality using an image filtering technique in the image space, instead of a reconstruction filter in the projection space for CT imaging. The authors evaluated the clinical feasibility of the proposed method in lung cancer screening. Methods: The proposed technique is essentially the same as common image filtering, which performs processing in the spatial-frequency domain with a filter function. However, the filter function was determined based on the quantitative analysis of the point spread functions (PSFs) measured in the system. The modulation transfer functions (MTFs) were derived from the PSFs, and the ratio of the MTFs was used as the filter function. Therefore, using an image reconstructed with a kernel, an image reconstructed with a different kernel was obtained by filtering, which used the ratio of the MTFs obtained for the two kernels. The performance of the method was evaluated by using routine clinical images obtained from CT screening for lung cancer in five subjects. Results: Filtered images for all combinations of three types of reconstruction kernels (''smooth,''''standard,'' and ''sharp'' kernels) showed good agreement with original reconstructed images regarded as the gold standard. On the filtered images, abnormal shadows suspected as being lung cancers were identical to those on the reconstructed images. The standard deviations (SDs) for the difference between filtered images and reconstructed images ranged from 1.9 to 23.5 Hounsfield units for all kernel combinations; these SDs were much smaller than the noise SDs in the reconstructed images. Conclusions: The proposed method has good performance and is clinically feasible in lung cancer screening. This method can be applied to images reconstructed on any scanner by measuring the PSFs in each system.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Combining Multiple Gyroscope Outputs for Increased Accuracy

A proposed method of processing the outputs of multiple gyroscopes to increase the accuracy of rate (that is, angular-velocity) readings has been developed theoretically and demonstrated by computer simulation. Although the method is applicable, in principle, to any gyroscopes, it is intended especially for application to gyroscopes that are parts of microelectromechanical systems (MEMS). The method is based on the concept that the collective performance of multiple, relatively inexpensive, nominally identical devices can be better than that of one of the devices considered by itself. The method would make it possible to synthesize the readings of a single, more accurate gyroscope (a virtual gyroscope) from the outputs of a large number of microscopic gyroscopes fabricated together on a single MEMS chip. The big advantage would be that the combination of the MEMS gyroscope array and the processing circuitry needed to implement the method would be smaller, lighter in weight, and less power-hungry, relative to a conventional gyroscope of equal accuracy. The method (see figure) is one of combining and filtering the digitized outputs of multiple gyroscopes to obtain minimum-variance estimates of rate. In the combining-and-filtering operations, measurement data from the gyroscopes would be weighted and smoothed with respect to each other according to the gain matrix of a minimum- variance filter. According to Kalman-filter theory, the gain matrix of the minimum-variance filter is uniquely specified by the filter covariance, which propagates according to a matrix Riccati equation. The present method incorporates an exact analytical solution of this equation.

Bayard, David S.↗

Patched Local Lunar Gravity Solutions using GRAIL Data

We present a method to determine local gravity fields for the Moon using Gravity Recovery and Interior Laboratory (GRAIL) data. We express gravity as gridded gravity anomalies on a sphere, and we estimate adjustments to a background global start model expressed in spherical harmonics. We processed GRAIL Ka-band range-rate data with a short-arc approach, using only data over the area of interest. We determine our gravity solutions using neighbor smoothing constraints. We divided the entire Moon into twelve regions and two polar caps, with a resolution of 0.15˚ x 0.15˚ (which is equivalent to degree and order 1199 in spherical harmonics), and determined the optimal smoothing parameter for each area by comparing localized correlations between gravity and topography for each solution set. Our selected areas share nodes with surrounding areas and they are overlapping. To mitigate boundary effects, we patch the solutions together by symmetrically omitting the boundary parts of overlapping solutions. Our new solution has been iterated, and it has improved correlations with topography when compared to a fully iterated global model. Our method requires fewer resources, and can easily handle regionally varying resolution or constraints. The smooth model describes small-scale features clearly, and can be used in local studies of the structure of the lunar crust.

Sander Goossens↗

Estimation of uncertainty for contour method residual stress measurements

This paper describes a methodology for the estimation of measurement uncertainty for the contour method, where the contour method is an experimental technique for measuring a two-dimensional map of residual stress over a plane. Random error sources including the error arising from noise in displacement measurements and the smoothing of the displacement surfaces are accounted for in the uncertainty analysis. The output is a two-dimensional, spatially varying uncertainty estimate such that every point on the cross-section where residual stress is determined has a corresponding uncertainty value. Both numerical and physical experiments are reported, which are used to support the usefulness of the proposed uncertainty estimator. The uncertainty estimator shows the contour method to have larger uncertainty near the perimeter of the measurement plane. For the experiments, which were performed on a quenched aluminum bar with a cross section of 51 × 76 mm, the estimated uncertainty was approximately 5 MPa (σ/E = 7 · 10⁻⁵) over the majority of the cross-section, with localized areas of higher uncertainty, up to 10 MPa (σ/E = 14 · 10⁻⁵).

42 ENGINEERING↗

Growth of Homoepitaxial ZnO Semiconducting Films

As a high temperature wide-band-gap (3.3 eV at room temperature) semiconductor, ZnO has been used for many applications such as wave-guides, solar cells, and surface acoustic wave devices, Since ZnO has a 60 meV excitonic binding energy that makes it possible to produce excitonic lasing at room temperature, a recent surge of interest is to synthesize ZnO films for UV/blue/green laser diodes. These applications require films with a smooth surface, good crystal quality, and low defect density. Thus, homoepitaxial film growth is the best choice. Homoepitaxial films have been studied in terms of morphology, crystal structure, and electrical and optical properties. ZnO single crystal substrates grown by the hydrothermal method are mechanically polished and annealed in air for four hours before the films are deposited. The annealing temperature-dependence on ZnO substrate morphology and electrical properties is investigated. Films are synthesized by off-axis reactive sputtering deposition. This produces films that have very smooth surfaces with roughness less than or equal to 5 nm on a 5 microns x 5 microns area. The full width at half maximum of film theta rocking curves measured by the x-ray diffraction is slightly larger than that of the crystal substrate. Films are also characterized by measuring resistivity, optical transmittance, and photoluminescence. The properties of ZnO films grown on (0001) ZnO and (0001) sapphire substrates will also be compared and discussed.

Zhu, Shen↗

Exploiting Power Flow Manifold to Solve AC Optimal Power Flow

AC optimal power flow has proven difficult to solve with interior point methods on GPUs. This is largely due to challenging linear algebra problems that current state of the art massively parallel linear solvers struggle with. However, the advent of Riemannian optimization techniques and the fact that the power flow equations form a smooth manifold present an alternative approach. In this talk, we present the basics of Riemannian optimization techniques in which optimization is done directly on a manifold. Then we present computational results showing that Riemannian techniques are capable of producing solutions of comparable quality as interior point methods.

AC optimal power flow↗

GalaxyFlow: upsampling hydrodynamical simulations for realistic mock stellar catalogues

ABSTRACT Cosmological N-body simulations of galaxies operate at the level of ‘star particles’ with a mass resolution on the scale of thousands of solar masses. Turning these simulations into stellar mock catalogues requires ‘upsampling’ the star particles into individual stars following the same phase-space density. In this paper, we introduce two new upsampling methods. First, we describe GalaxyFlow, a sophisticated upsampling method that utilizes normalizing flows to both estimate the stellar phase-space density and sample from it. Secondly, we improve on existing upsamplers based on adaptive kernel density estimation (KDE), using maximum likelihood estimation to fine-tune the bandwidth for such algorithms in a way that improves both the density estimation accuracy and upsampling results. We demonstrate our upsampling techniques on a neighbourhood of the Solar location in two simulated galaxies: Auriga 6 and h277. Both yield smooth stellar distributions that closely resemble the stellar densities seen in the Gaia DR3 catalogue. Furthermore, we introduce a novel multimodel classifier test to compare the accuracy of different upsampling methods quantitatively. This test confirms that GalaxyFlow more accurately estimates the density of the underlying star particles than methods based on KDE, at the cost of being more computationally intensive.

Lim, Sung Hak (ORCID:0000000330981092)↗

Adaptive AMG with Coarsening Based on Compatible Weighted Matching

We introduce a new composite adaptive Algebraic Multigrid (composite αAMG) method to solve systems of linear equations without a-priori knowledge or assumption on characteristics of near-null components of the AMG preconditioned problem referred to as algebraic smoothness. Our version of αAMG is a composite solver built through a bootstrap strategy aimed to obtain a desired convergence rate. The coarsening process employed to build each new solver component relies on a pairwise aggregation scheme based on maximum weighted matching in a graph and on principles of compatible relaxation. The latter replaces the commonly used characterization of strength of connection in both the coarse space selection and in the interpolation scheme. The goal is to design a method leading to scalable AMG for a wide class of problems that go beyond the standard elliptic Partial Differential Equations (PDEs). In the present work, we introduce the method and demonstrate its potential when applied to symmetric positive definite linear systems arising from finite element discretization of highly anisotropic elliptic PDEs on structured and unstructured meshes.

97 MATHEMATICS AND COMPUTING↗

An explicit finite-volume time-marching procedure for turbulent flow calculations

A method was developed which calculates two-dimensional, transonic, viscous flow in ducts. The finite-volume, time-marching formulation is used to obtain steady flow solutions of the Reynolds-averaged form of the Navier-Stokes equations. The entire calculation is performed in the physical domain. Control volumes are chosen so that smoothing of flow properties, typically required for stability, is not required. Different time steps are used in the different governing equations. A new pressure interpolation scheme is introduced which improves the shock capturing ability of the method. A multi-volume method for pressure changes in the boundary layer allows calculations which use very long and thin control volumes (length/height - 1000). The method is compared with two test cases. Essentially incompressible turbulent boundary layer flow in an adverse pressure gradient is calculated and the computed distributions of mean velocity and shear are in good agreement with the measurements. Transonic viscous flow in a converging diverging nozzle is calculated; the Mach number upstream of the shock is approximately 1.25. The agreement between the calculated and measured shock strength and total pressure losses is good.

Nicholson, Stephen↗

An explicit finite-volume time-marching procedure for turbulent flow calculations

A method was developed which calculates two-dimensional, transonic, viscous flow in ducts. The finite-volume, time-marching formulation is used to obtain steady flow solutions of the Reynolds-averaged form of the Navier-Stokes equations. The entire calculation is performed in the physical domain. Control volumes are chosen so that smoothing of flow properties, typically required for stability, is not required. Different time steps are used in the different governing equations. A new pressure interpolation scheme is introduced which improves the shock capturing ability of the method. A multi-volume method for pressure changes in the boundary layer allows calculations which use very long and thin control volumes (length/height - 1000). The method is compared with two test cases. Essentially incompressible turbulent boundary layer flow in an adverse pressure gradient is calculated and the computed distributions of mean velocity and shear are in good agreement with the measurements. Transonic viscous flow in a converging diverging nozzle is calculated; the Mach number upstream of the shock is approximately 1.25. The agreement between the calculated and measured shock strength and total pressure losses is good.

Nicholson, Stephen↗

SU-F-E-09: Respiratory Signal Prediction Based On Multi-Layer Perceptron Neural Network Using Adjustable Training Samples

Purpose: Dynamic tracking of moving organs, such as lung and liver tumors, under radiation therapy requires prediction of organ motions prior to delivery. The shift of moving organ may change a lot due to huge transform of respiration at different periods. This study aims to reduce the influence of that changes using adjustable training signals and multi-layer perceptron neural network (ASMLP). Methods: Respiratory signals obtained using a Real-time Position Management(RPM) device were used for this study. The ASMLP uses two multi-layer perceptron neural networks(MLPs) to infer respiration position alternately and the training sample will be updated with time. Firstly, a Savitzky-Golay finite impulse response smoothing filter was established to smooth the respiratory signal. Secondly, two same MLPs were developed to estimate respiratory position from its previous positions separately. Weights and thresholds were updated to minimize network errors according to Leverberg-Marquart optimization algorithm through backward propagation method. Finally, MLP 1 was used to predict 120∼150s respiration position using 0∼120s training signals. At the same time, MLP 2 was trained using 30∼150s training signals. Then MLP is used to predict 150∼180s training signals according to 30∼150s training signals. The respiration position is predicted as this way until it was finished. Results: In this experiment, the two methods were used to predict 2.5 minute respiratory signals. For predicting 1s ahead of response time, correlation coefficient was improved from 0.8250(MLP method) to 0.8856(ASMLP method). Besides, a 30% improvement of mean absolute error between MLP(0.1798 on average) and ASMLP(0.1267 on average) was achieved. For predicting 2s ahead of response time, correlation coefficient was improved from 0.61415 to 0.7098.Mean absolute error of MLP method(0.3111 on average) was reduced by 35% using ASMLP method(0.2020 on average). Conclusion: The preliminary results demonstrate that the ASMLP respiratory prediction method is more accurate than MLP method and can improve the respiration forecast accuracy.

60 APPLIED LIFE SCIENCES↗

Generalized boost transformations in finite volumes and application to Hamiltonian methods

The investigation of hadron interactions within lattice QCD has been facilitated by the well-known quantisation condition, linking scattering phase shifts to finite-volume energies. Additionally, the ability to utilise systems at finite total boosts has been pivotal in smoothly charting the energy-dependent behaviour of these phase shifts. The existing implementations of the quantization condition at finite boosts rely on momentum transformations between rest and moving frames, defined directly in terms of the energy eigenvalues. This energy dependence is unsuitable in the formulation of a Hamiltonian. In this work, we introduce a novel approach to generalise the three-momentum boost prescription, enabling the incorporation of energy-independent finite-volume Hamiltonians within moving frames. We demonstrate the application of our method through numerical comparisons, employing a phenomenological ππ scattering example.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Hypersonic Flow Computations on Unstructured Meshes

A method for computing inviscid hypersonic flow over complex configurations using unstructured meshes is presented. The unstructured grid solver uses an edge{based finite{volume formulation. Fluxes are computed using a flux vector splitting scheme that is capable of representing constant enthalpy solutions. Second{order accuracy in smooth flow regions is obtained by linearly reconstructing the solution, and stability near discontinuities is maintained by locally forcing the scheme to reduce to first-order accuracy. The implementation of the algorithm to parallel computers is described. Computations using the proposed method are presented for a sphere-cone configuration at Mach numbers of 5.25 and 10.6, and a complex hypersonic re-entry vehicle at Mach numbers of 4.5 and 9.8. Results are compared to experimental data and computations made with established structured grid methods. The use of the solver as a screening tool for rapid aerodynamic assessment of proposed vehicles is described.

Bibb, K. L.↗