SEARCH · Engineering Papers

Results for “reasonable accuracy”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 919 records · Page 51

All sources NASA NTRS DOE OSTI

NASA NTRS1989

Accuracy of telemetry signal power loss in a filter as an estimate for telemetry degradation

When telemetry data is transmitted through a communication link, some degradation in telemetry performance occurs as a result of the imperfect frequency response of the channel. The term telemetry degradation as used here is the increase in received signal power required to offset this filtering. The usual approach to assessing this degradation is to assume that it is equal to the signal power loss in the filtering, which is easily calculated. However, this approach neglects the effects of the nonlinear phase response of the filter, the effect of any reduction of the receiving system noise due to the filter, and intersymbol interference. Here, an exact calculation of the telemetry degradation, which includes all of the above effects, is compared with the signal power loss calculation for RF filtering of NRZ data on a carrier. The signal power loss calculation is found to be a reasonable approximation when the filter follows the point at which the receiving system noise is introduced, especially if the signal power loss is less than 0.5 dB. The signal power loss approximation is less valid when the receiving system noise is not filtered.

Koerner, M. A.↗

NASA NTRS2010

Numerical Simulation of Vitiation Effects on a Hydrogen-Fueled Dual-Mode Scramjet

The Wind-US computational fluid dynamics (CFD) flow solver was used to simulate dual-mode direct-connect ramjet/scramjet engine flowpath tests conducted in the University of Virginia (UVa) Supersonic Combustion Facility (SCF). The objective was to develop a computational capability within Wind-US to aid current hypersonic research and provide insight to flow as well as chemistry details that are not resolved by instruments available. Computational results are compared with experimental data to validate the accuracy of the numerical modeling. These results include two fuel-off non-reacting and eight fuel-on reacting cases with different equivalence ratios, split between one set with a clean (non-vitiated) air supply and the other set with a vitiated air supply (12 percent H2O vapor). The Peters and Rogg hydrogen-air chemical kinetics model was selected for the scramjet simulations. A limited sensitivity study was done to investigate the choice of turbulence model and inviscid flux scheme and led to the selection of the k-epsilon model and Harten, Lax and van Leer (for contact waves) (HLLC) scheme for general use. Simulation results show reasonably good agreement with experimental data and the overall vitiation effects were captured.

Vyas, Manan A.↗

DOE OSTI2022

Hybrid large-eddy simulation with adaptive downstream anisotropic eddy viscosity model

Turbulence pervades most flows of engineering interest, and its prediction remains a challenge on both accuracy and cost. One promising predictive approach that reduces cost combines large-eddy simulation (LES) with simulation based on Reynolds averaged Navier–Stokes equations (RANS). This study presents a method to overcome stability and accuracy issues associated with these hybrid LES-RANS methods. The method developed involves extracting the Reynolds stresses from the upstream LES solution and then using that information to convert the downstream RANS model from a scalar eddy-viscosity model to an anisotropic nonlinear eddy-viscosity model. The method developed differs from the downstream tensorial eddy-viscosity model by being independent of the coordinate system. The method developed was evaluated by computing film cooling of a flat plate with the coolant injected through one row of circular holes. Results obtained show instabilities at the LES-to-RANS interface to be eliminated. Also, the method developed yielded solutions that compare reasonably well with those from LES, even though a significant portion of the flow is computed by the adapted anisotropic RANS model instead of LES, which significantly reduced the number of grid points and computational time needed. Since modification of the downstream RANS model is based on information extracted from the upstream LES solution, the method developed is adaptive to the problem being studied.

42 ENGINEERING↗

DOE OSTI2020

Viscous Dissipation and Apparent Permeability of Gas Flow in Nanoporous Media

The nanopore confinement effect can increase the mobility for gas flow in shale rocks, and the apparent permeability is usually adopted to account for this effect. However, most of the apparent permeability calculation models are derived based on straight channels (capillaries) and neglect the influence of end effect. The end effect is caused by the bending of streamlines which yields more viscous dissipation and additional flow resistance for gas flow in nanoporous media. Here, for the first time the viscous dissipation for gas flow in nanoporous media is analyzed. In addition, an improved apparent permeability calculation model is proposed based on the Beskok-Karniadakis (B-K) model by introducing the influence of end effect, which is characterized by the ratio of the length to the width of a short channel (L/H). The accuracy of this model is validated by lattice Boltzmann simulations of gas flow in three different geometries: an infinite-length straight channel, a finite-length straight channel, and nanoporous media. For the infinite-length straight channel, there is no end effect therefore both the B-K model and the improved model can well predict the apparent permeability. However, for the finite-length straight channel and nanoporous media, the original B-K model overestimates the gas apparent permeability as it neglects the influence of end effect, while the improved model can still well predict the gas apparent permeability. In summary, the improved apparent permeability calculation model is more reasonable in predicting gas mobility in nanoporous media.

58 GEOSCIENCES↗

NASA NTRS1993

Constraints on the rate of discharge and duration of the Mangala Valles flood

Interest in Mangala Valles remains high within the planetary science community. This is justified because the survey mission images provide us with nearly complete coverage of the system at high resolution. Upcoming high resolution topography from the Mars Observer Laser Altimeter will enable the application of flood discharge models with an unprecedented level of detail. Previous work on the Mangala Valles problem has relied on the available low resolution topography. This has limited the ability of workers to constrain discharge calculations to only within several-order-of-magnitude estimates. Local determinations of channel depths via shadow length measurements and photoclinometric profiling are much more accurate, but can only be applied to steep slopes (and/or low sun elevations) in the case of shadow measurements, or across relatively short distances (to avoid changes in albedo along asymmetric photoclinometric profiles). We are taking stereo parallax measurements from medium resolution Viking Orbiter images, which provide a valuable intermediate check of the topography between those measurements made thus far and the upcoming Mars Observer data. The images used are from orbits 034A and 637A, and cover Mangala Valles from the source graben in Memnonia Fossae to the beginning of the bifurcated reach (at 9.5 deg lat., 151.5 deg lon.). These images are about 300 m/pixel and 250 m/pixel, respectively. Both sets of images were orthographically projected to 250 m/pixel. The separation angle between left-right pairs is approximately 52 degrees. This results in a vertical accuracy on the order of plus or minus 260 m. Though this is still somewhat coarse, the channel relief is clearly resolved. Preliminary profiles across Mangala Valles and the large topographic ridge to the east are shown on the following page. An east-west regional tilt that resulted from slight scaling differences in the digital data has been 'removed' by visually estimating a regional datum on the stereo pairs. North-south variations in scale have not yet been adjusted for, so the zero datum for each profile was simply taken to be the high water line of the channel itself. Our depth measurement for the source breach of about 750 m agrees reasonably well with values of 700-1000 m determined from shadow measurements.

Parker, T. J.↗

NASA NTRS1984

Positive-feedback photometric drift in the PDS

Digitizing flatfield images produces conditions in the Photometric Data System PDS which cause the measured density to drift by as much as .1 DN during a 10 minute interval. The drift occurs when the PDS, set up in equilibrium at fog level, subsequently scans a reasonably dense region for periods of longer than a few minutes. The drift is manifested primarily as a positive shift in density that is approximately the same for all densities. If the fog level is assumed to be in fact constant and is monitored during scans of flat fields, the PDS drift may be removed by subtracting the difference between the observed fog level and its assumed constant value for each pixel. This function is then smoothed and subtracted, as a function of scan line, from the measured density. The fog level is then adjusted to a standard value by adding a constant. The result is a flattened scan with PDS drift removed to the accuracy within which the fog level drift matches the drift at other levels.

Cornett, R. H.↗

DOE OSTI2025

Replace Human Intelligence with Fast and Smart Geometric Reasoning and Graph Neural Network to Accelerate Next Gen ModSim Workflows

We present an agent-guided approach to CAD geometry decomposition that automates hex/hybrid meshing with graph neural networks (GNNs) to accelerate next-generation ModSim workflows. Our end-to-end pipeline (i) reduces 3D boundary-representation (B-Rep) models to a 2D chordal axis skeleton (CAT) and then to a 1D bipartite graph of surface and curve nodes, (ii) assigns per node labels as Cubit® WebCut actions, (iii) trains a multi-action GNN under supervised learning, and (iv) predicts five surface-node and three curve-node actions on out-of-distribution test geometries. Each graph node carries geometric, topological, and meshing attributes drawn from the B-Rep “skin” and CAT “skeleton,” with two-way mappings across 3D↔2D↔1D representations to maintain traceability back to 3D CAD. The supervised learning model exhibits stable convergence of the binary cross-entropy loss and achieves 98.7% accuracy on unseen lattice models. To operationalize decision-making, we rank predicted commands by geometric significance and prototyped the agent-guided workflow through the Cubit® Meshing PowerTool GUI. As a stretch goal, we explore reinforcement learning (RL) to reduce or remove label requirements and to learn policies for action sequences that maximize total reward (e.g., size of hex-meshable regions and resulting hex mesh quality). When all-hex meshing is not feasible, the agent assists in producing hybrid meshes—prioritizing hex in critical regions and transitioning to tetrahedral elements (tets) elsewhere—maintaining fidelity while ensuring robustness. The overarching objective is to replace manual, heuristics-based decomposition with data-driven, reproducible automation, cutting meshing turnaround time by orders of magnitude. We anticipate direct impact on simulation workflows through intelligent, scalable decomposition of complex CAD models into hex-meshable subdomains.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Demonstrating Computational Equivalence Between Continuous and Discrete Adjoint Methods by Calculating Time-Dependent Adjoint Solutions with Neutron Diffusion Models

The continuous adjoint method and the discrete adjoint method are two alternative approaches used to calculate adjoint solutions for adjoint systems. The continuous adjoint method derives adjoint equations analytically from continuous forward equations and then solves the adjoint equations either analytically or numerically in a discretized form whereas the discrete adjoint method calculates the adjoint solutions directly from the discretized forward equations. With regard to the methodology development and calculation procedure, distinct differences are well recognized between the two methods. For certain reasons, both methods are exclusively preferred and commonly used by different computational communities, but limited studies clarify the connections between the two adjoint methods from either of the communities. Herein, this paper demonstrates the computational equivalence between the continuous and discrete adjoint methods by investigating time-dependent adjoint solutions to the two-group neutron diffusion model in nuclear reactor analysis problems using both methods. Adjoint solutions can be used to estimate system parameters for reactor safety analysis. Appropriate final state conditions for the adjoint systems are specified in both of the methods, and the conditions are clarified with proper physical explanations. With the help of an event-based case study on neutron diffusion models, the accuracy of the time-dependent adjoint fluxes obtained from both methods is verified, and the pros and cons of both adjoint methods are examined. More importantly, the computational equivalence of both methods is demonstrated when they are applied to multigroup neutron diffusion systems. The advantage of calculating time-dependent adjoint fluxes by directly solving time-dependent adjoint systems rather than taking steady-state approximations as in common practice is also demonstrated.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS1973

Oxygen fugacity values of Apollo 12, 14, and 15 lunar samples and reduced state of lunar magmas

The oxygen fugacity values of lunar samples were measured directly with an improved solid-electrolyte oxygen cell between 1000 and 1200 C with an accuracy mostly better than 0.2 log f(02) unit. The bulk rock f(02) values of basaltic igneous rocks 12009, 12053, 15058, and 15595 ranged from 10 to -15.4 power to 10 to the -15.7 power at 1000 C and from 10 to the -12.3 power to 10 to the -12.8 power at 1200 C. Those of microbreccia 14321 also fell in this range, but the change of oxygen fugacity with temperature was irregular in the first heating cycle in comparison to the smooth changes observed with the basaltic igneous rocks. Two different samples of rock 14310 also showed similar f(02) values below 1170 C, but exhibited irreversible sudden rise in f(02) at this temperature for reasons yet to be determined. The data on the phenocryst olivine and the groundmass of basalt 12009 do not conclusively indicate progressive reduction of the lunar magma during cooling.

Sato, M.↗

NASA NTRS2011

Assessment of Some Atomization Models Used in Spray Calculations

The paper presents the results from a validation study undertaken as a part of the NASA s fundamental aeronautics initiative on high altitude emissions in order to assess the accuracy of several atomization models used in both non-superheat and superheat spray calculations. As a part of this investigation we have undertaken the validation based on four different cases to investigate the spray characteristics of (1) a flashing jet generated by the sudden release of pressurized R134A from cylindrical nozzle, (2) a liquid jet atomizing in a subsonic cross flow, (3) a Parker-Hannifin pressure swirl atomizer, and (4) a single-element Lean Direct Injector (LDI) combustor experiment. These cases were chosen because of their importance in some aerospace applications. The validation is based on some 3D and axisymmetric calculations involving both reacting and non-reacting sprays. In general, the predicted results provide reasonable agreement for both mean droplet sizes (D32) and average droplet velocities but mostly underestimate the droplets sizes in the inner radial region of a cylindrical jet.

Raju, M. S.↗

NASA NTRS2008

Efficient Probabilistic Diagnostics for Electrical Power Systems

We consider in this work the probabilistic approach to model-based diagnosis when applied to electrical power systems (EPSs). Our probabilistic approach is formally well-founded, as it based on Bayesian networks and arithmetic circuits. We investigate the diagnostic task known as fault isolation, and pay special attention to meeting two of the main challenges . model development and real-time reasoning . often associated with real-world application of model-based diagnosis technologies. To address the challenge of model development, we develop a systematic approach to representing electrical power systems as Bayesian networks, supported by an easy-to-use speci.cation language. To address the real-time reasoning challenge, we compile Bayesian networks into arithmetic circuits. Arithmetic circuit evaluation supports real-time diagnosis by being predictable and fast. In essence, we introduce a high-level EPS speci.cation language from which Bayesian networks that can diagnose multiple simultaneous failures are auto-generated, and we illustrate the feasibility of using arithmetic circuits, compiled from Bayesian networks, for real-time diagnosis on real-world EPSs of interest to NASA. The experimental system is a real-world EPS, namely the Advanced Diagnostic and Prognostic Testbed (ADAPT) located at the NASA Ames Research Center. In experiments with the ADAPT Bayesian network, which currently contains 503 discrete nodes and 579 edges, we .nd high diagnostic accuracy in scenarios where one to three faults, both in components and sensors, were inserted. The time taken to compute the most probable explanation using arithmetic circuits has a small mean of 0.2625 milliseconds and standard deviation of 0.2028 milliseconds. In experiments with data from ADAPT we also show that arithmetic circuit evaluation substantially outperforms joint tree propagation and variable elimination, two alternative algorithms for diagnosis using Bayesian network inference.

Mengshoel, Ole J.↗

DOE OSTI2024

A numerical study of process complexity in permafrost dominated regions

Numerical modeling of permafrost dynamics requires adequate representation of atmospheric and surface processes, a reasonable parameter estimation strategy, and site-specific model development. The three main research objectives of the study are: (i) to propose a novel methodology that determines the required level of surface process complexity of permafrost models by conducting parameter sensitivity and calibration, (ii) to design and compare three numerical models of increasing surface process complexity, and (iii) to calibrate and validate the numerical models at the Yakou catchment on the Qinghai-Tibet Plateau as an exemplary study site. The calibration was carried out by coupling the Advanced Terrestrial Simulator (numerical model) and PEST (calibration tool). Simulation results showed that (i) A simple numerical model that considers only subsurface processes can simulate active layer development with the same accuracy as other more complex models that include surface processes. (ii) Peat and mineral soil layer permeability, Van Genuchten alpha, and porosity are highly sensitive. (iii) Liquid precipitation aids in increasing the rate of permafrost degradation. (iv) Deposition of snow insulated the subsurface during the thaw initiation period. We have developed and released an integrated code that couples the numerical software ATS to the calibration software PEST. The numerical model can be further used to determine the impacts of climate change on permafrost degradation.

Calibration↗

NASA NTRS1991

Station keeping strategy for multiple spacecraft interferometry

The feasibility of multiple spacecraft stationkeeping for submillimeter and optical interferometry is examined. A condition for interferometry is that two or more spacecraft must control their relative positions with better than 1 mn accuracy indefinitely in both radial and transverse directions although separated by as much as 1 Km in LEO and 100 Km in GEO. They must also maneuver through a useful area of the U-V plane of an arbitrary astronomical source. The problem is first outlined and a solution which utilizes gravity gradient forces to do most of the work and ion thrusters for additional maneuvering is proposed. All the perturbing forces are shown to be small compared to the ion thruster requirements. An inertial position and attitude control strategy is suggested which utilizes existing or soon to be available sensors and actuators. Finally, the fuel and power system mass requirements are estimated and found to be within reason for a 10 year mission.

Decou, Anthony B.↗

NASA NTRS2001

Twenty Years of Radiation Measurements in Low-Earth Orbit - What Have We Learned Space Radiation Environment?

The advent of the Space Shuttle program has made possible space radiation environment measurements spanning a wide range of altitudes and orbital inclinations over multiple solar cycles. These measurements range from routine integral dose measurements with thermoluminescent dosimeters to particle energy spectra measurements made with a charged particle telescope. This paper will review the new understanding about the space radiation environment gained from this diverse data set. Major findings from these measurements include: estimations of the westward drift rate of the South Atlantic Anomaly (SAA) of 0.28-0.49/y; evidence for a northward component to the SAA drift of 0.08-0.12/y; observation of the formation and decay of the pseudo-stable additional radiation belt following the Mar 1991 SPE and geomagnetic storm with an estimated decay e-folding time of 9-10 months; observation of a local geomagnetic east-west trapped proton exposure anisotropy with an estimated magnitude of 1.6-3.3; demonstration that the trapped proton exposure in low-Earth orbit (LEO) can be reasonably modeled as a power law function of atmospheric density in the SAA region, with best correlations obtained when the exospheric temperature saturates at 938-975 K; the actual solar cycle modulation of trapped proton exposure in LEO is less than predicted by the AP8 model; and the testing and validation of GCR flux models, radiation transport codes, and dynamic geomagnetic cutoff models. Long-term, time-resolved proportional counter measurements made aboard the Mir during the same period provides further demonstration of the solar cycle modulation of the trapped protons at low altitudes - the observed modulation is also well described as power law function of atmospheric density. These data and findings have helped to improve the overall accuracy of pre-mission crew exposure projections using various semi-empirical space environment models, radiation transport codes, and spacecraft radiation shielding models. During the rise phase of solar cycle 22 (1987-1991), the RMS error between preflight exposure projections and measured crew exposure was 73%. For the rise phase of cycle 23 (1997-2001), the preflight exposure projection RMS error has decreased to 23%. The launch and assembly of the Space Station has begun a new era of long-term LEO space environment monitoring. The radiation environment at the Space Station will be monitored with three external charged particle telescopes oriented in the velocity vector, anti-velocity vector, and zenith directions. Data from the telescopes will provide charge, mass, energy, and arrival direction for incident particles with energy to mass ratios of 13- 450 MeV/amu and Z of 1-24. The external environment data will be complimented by measurements from a portable charged particle telescope and proportional counter located inside the vehicle.

Golightly, Michael J.↗

NASA NTRS2013

Potassium Isotopic Compositions of NIST Potassium Standards and 40Ar/39Ar Mineral Standards

Knowledge of the isotopic ratios of standards, spikes, and reference materials is fundamental to the accuracy of many geochronological methods. For example, the 238U/235U ratio relevant to U-Pb geochronology was recently re-determined [1] and shown to differ significantly from the previously accepted value employed during age determinations. These underlying values are fundamental to accurate age calculations in many isotopic systems, and uncertainty in these values can represent a significant (and often unrecognized) portion of the uncertainty budget for determined ages. The potassium isotopic composition of mineral standards, or neutron flux monitors, is a critical, but often overlooked component in the calculation of K-Ar and 40Ar/39Ar ages. It is currently assumed that all terrestrial materials have abundances indistinguishable from that of NIST SRM 985 [2]; this is apparently a reasonable assumption at the 0.25per mille level (1σ) [3]. The 40Ar/39Ar method further relies on the assumption that standards and samples (including primary and secondary standards) have indistinguishable 40K/39K values. We will present data establishing the potassium isotopic compositions of NIST isotopic K SRM 985, elemental K SRM 999b, and 40Ar/39Ar biotite mineral standard GA1550 (sample MD-2). Stable isotopic compositions (41K/39K) were measured by the peak shoulder method with high resolution MC-ICP-MS (Thermo Scientific NEPTUNE Plus), using the accepted value of NIST isotopic SRM 985 [2] for fractionation [4] corrections [5]. 40K abundances were measured by TIMS (Thermo Scientific TRITON), using 41K/39K values from ICP-MS measurements (or, for SRM 985, values from [2]) for internal fractionation corrections. Collectively these data represent an important step towards a metrologically traceable calibration of 40K concentrations in primary 40Ar/39Ar mineral standards and improve uncertainties by ca. an order of magnitude in the potassium isotopic compositions of standards.

Morgan, Leah↗

DOE OSTI2024

Low-Cost Heliostat for High-Flux Small-Area Receivers (Final Technical Report)

This project analyzed a two-stage heliostat concept consisting of a tracking stage and a concentrating stage. The tracking stage uses mirrors mounted on a common drive that move to track the sun. The concentrating stage consists of stationary mirrors that each have a unique angle to direct rays towards a small-area, high-flux, point-focused receiver. By splitting the collection and concentrating process into two stages, multiple small, inexpensive mirrors can share a structure and be controlled by a single drive in the tracking stage. The project effort developed modeling techniques that were specifically relevant to this two-stage heliostat concept. Both field-level and unit-level models were developed. The field-level model does not explicitly consider unit-level losses which are predicted by the unit-level model and then integrated into the field-level model through a correlation referred to as an efficiency modifier. This approach is referred to as the two-model approach; the development and demonstration of this two-model approach for a multi-stage heliostat technology is a key outcome of this work. The field-level model is used to design a field that hits a specific design day power given a set of heliostat design parameters. An oversized field is simulated and then heliostat units are removed based on their annual energy production in order to generate the highest performing field. The field reduction procedure fits a smooth curve fit to annual energy production as a function of position in the field which has the effect of reducing the noise that is otherwise caused by the Monte Carlo ray tracing technique. This approach is referred to as the annual energy fit method and substantially reduces computational run time for a given field level modeling accuracy. The annual energy fit approach enables the selection of a properly sized, high-performing field using orders of magnitude fewer rays than would otherwise be possible and the development of this approach is a second key outcome of this work. These models are used within a genetic optimization algorithm in order to optimize the geometric parameters associated with a heliostat in order to achieve the lowest cost per unit of collected design day power. The cost modeling that underlies the optimization is a simple, scaling type analysis backed up by a much more detailed Design for Manufacture and Assembly (DFMA) analysis. Although the figure of merit used for optimization was not cost per mirror area, this metric is reasonable to use as a means of comparison. The optimally designed 500 kW design has a tracking mirror specific cost of $181.85/m 2 , which is significantly larger than the target value and also larger than the current state of the art. The cost of the torque-tube type linkages contributed substantially to the overall cost. Based on this observation, potentially attractive alternative design configuration utilizing a capstan type actuation system should be investigated. Finally, NREL compared the performance of the two-stage heliostat to the performance of a focused and different sized flat conventional heliostats and showed that, as expected, additional losses versus the convention heliostat caused by a worse cosine efficiency, two stages of reflection, and interstage interactions. The two-stage heliostat requires around 75% more reflective area than a flat 1x1 meter conventional heliostat (similar to a focused heliostat) and 40% more than a flat 2x2 meter conventional heliostat.

14 SOLAR ENERGY↗

DOE OSTI2022

A Logarithmic Bayesian Approach to Quantum Error Detection

We consider the problem of continuous quantum error correction from a Bayesian perspective, proposing a pair of digital filters using logarithmic probabilities that are able to achieve near-optimal performance on a three-qubit bit-flip code, while still being reasonable to implement on low-latency hardware. These practical filters are approximations of an optimal filter that we derive explicitly for finite time steps, in contrast with previous work that has relied on stochastic differential equations such as the Wonham filter. By utilizing logarithmic probabilities, we are able to eliminate the need for explicit normalization and can reduce the Gaussian noise distribution to a simple quadratic expression. The state transitions induced by the bit-flip errors are modeled using a Markov chain, which for log-probabilities must be evaluated using a LogSumExp function. We develop the two versions of our filter by constraining this LogSumExp to have either one or two inputs, which favors either simplicity or accuracy, respectively. Using simulated data, we demonstrate that the single-term and two-term filters are able to significantly outperform both a double threshold scheme and a linearized version of the Wonham filter in tests of error detection under a wide variety of error rates and time steps.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS2009

CFD Analysis of Nozzle Jet Plume Effects on Sonic Boom Signature

An axisymmetric full Navier-Stokes computational fluid dynamics study is conducted to examine nozzle exhaust jet plume effects on the sonic boom signature of a supersonic aircraft. A simplified axisymmetric nozzle geometry, representative of the nozzle on the NASA Dryden NF-15B Lift and Nozzle Change Effects on Tail Shock research airplane, is considered. The computational fluid dynamics code is validated using available wind-tunnel sonic boom experimental data. The effects of grid size, spatial order of accuracy, grid type, and flow viscosity on the accuracy of the predicted sonic boom pressure signature are quantified. Grid lines parallel to the Mach wave direction are found to give the best results. Second-order accurate upwind methods are required as a minimum for accurate sonic boom simulations. The highly underexpanded nozzle flow is found to provide significantly more reduction in the tail shock strength in the sonic boom N-wave pressure signature than perfectly expanded and overexpanded nozzle flows. A tail shock train in the sonic boom signature is observed for the highly underexpanded nozzle flow. Axisymmetric computational fluid dynamics simulations show the flow physics inside the F-15 nozzle to be nonisentropic and complex. Although the one-dimensional isentropic nozzle plume results look reasonable, they fail to capture the sonic boom shock train in the highly underexpanded nozzle flow.

Bui, Trong T.↗