Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “delta method”

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 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
NASA NTRS2017

Carbon-Carbon Nozzle Extension Development in Support of In-Space and Upper-Stage Liquid Rocket Engines

Upper stage and in-space liquid rocket engines are optimized for performance through the use of high area ratio nozzles to fully expand combustion gases to low exit pressures, increasing exhaust velocities. Due to the large size of such nozzles, and the related engine performance requirements, carbon-carbon (C-C) composite nozzle extensions are being considered to reduce weight impacts. Currently, the state-of-the-art is represented by the metallic and foreign composite nozzle extensions limited to approximately 2000 degrees F. used on the Atlas V, Delta IV, Falcon 9, and Ariane 5 launch vehicles. NASA and industry partners are working towards advancing the domestic supply chain for C-C composite nozzle extensions. These development efforts are primarily being conducted through the NASA Small Business Innovation Research (SBIR) program in addition to other low level internal research efforts. This has allowed for the initial material development and characterization, subscale hardware fabrication, and completion of hot-fire testing in relevant environments. NASA and industry partners have designed, fabricated and hot-fire tested several subscale domestically produced C-C extensions to advance the material and coatings fabrication technology for use with a variety of liquid rocket and scramjet engines. Testing at NASA's Marshall Space Flight Center (MSFC) evaluated heritage and state-of-the-art C-C materials and coatings, demonstrating the initial capabilities of the high temperature materials and their fabrication methods. This paper discusses the initial material development, design and fabrication of the subscale carbon-carbon nozzle extensions, provides an overview of the test campaign, presents results of the hot fire testing, and discusses potential follow-on development work. The follow on work includes the fabrication of ultra-high temperature materials, larger C-C nozzle extensions, material characterization, sub-element testing and hot-fire testing at larger scale.

Gradl, Paul R.↗
NASA NTRS2008

Measurement of Off-Body Velocity, Pressure, and Temperature in an Unseeded Supersonic Air Vortex by Stimulated Raman Scattering

A noninvasive optical method is used to make time-averaged (30 sec) off-body measurements in a supersonic airflow. Seeding of tracer particles is not required. One spatial component of velocity, static pressure, and static temperature are measured with stimulated Raman scattering. The three flow parameters are determined simultaneously from a common sample volume (0.3 by 0.3 by 15 mm) using concurrent measurements of the forward and backward scattered line shapes of a N2 vibrational Raman transition. The capability of this technique is illustrated with laboratory and large-scale wind tunnel testing that demonstrate 5-10% measurement uncertainties. Because the spatial resolution of the present work was improved to 1.5 cm (compared to 20 cm in previous work), it was possible to demonstrate a modest one-dimensional profiling of cross-flow velocity, pressure, and translational temperature through the low-density core of a stream-wise vortex (delta-wing model at Mach 2.8 in NASA Langley's Unitary Plan Wind Tunnel).

Herring, Gregory C.↗
NASA NTRS2014

A Novel Framework for Quantifying past Methane Recycling by Sphagnum-Methanotroph Symbiosis Using Carbon and Hydrogen Isotope Ratios of Leaf Wax Biomarkers

The concentration of atmospheric methane is strongly linked to variations in Earth's climate. Currently, we can directly reconstruct the total atmospheric concentration of methane, but not individual terms of the methane cycle. Northern wetlands, dominated by Sphagnum, are an important contributor of atmospheric methane, and we seek to understand the methane cycle in these systems. We present a novel method for quantifying the proportion of carbon Sphagnum assimilates from its methanotrophic symbionts using stable isotope ratios of leaf-wax biomarkers. Carbon isotope ratios of Sphagnum compounds are determined by two competing influences, water content and the isotope ratio of source carbon. We disentangled these effects using a combined hydrogen and carbon isotope approach. We constrained Sphagnum water content using the contrast between the hydrogen isotope ratios of Sphagnum and vascular plant biomarkers. We then used Sphagnum water content to calculate the carbon isotope ratio of Sphagnum's carbon pool. Using a mass balance equation, we calculated the proportion of recycled methane contributed to the Sphagnum carbon pool, 'PRM.' We quantified PRM in peat monoliths from three microhabitats in the Mer Bleue peatland complex. Modern studies have shown that water table depth and vegetation have strong influences on the peatland methane cycle on instrumental time scales. With this new approach, delta C-13 of Sphagnum compounds are now a useful tool for investigating the relationships among hydrology, vegetation, and methanotrophy in Sphagnum peatlands over the time scales of entire peatland sediment records, vital to our understanding of the global carbon cycle through the Late Glacial and Holocene.

air pollution↗
DOE OSTI2022

Intercomparison of Dynamically and Statistically Downscaled Climate Change Projections over the Midwest and Great Lakes Region

Downscaling of global climate model (GCMs) simulations is a key element of regional-to-local-scale climate change projections that can inform impact assessments, long-term planning, and resource management in different sectors. Here, we conduct an intercomparison between statistically and dynamically downscaled GCMs simulations using the hybrid delta (HD) and the Weather Research and Forecast (WRF) Model, respectively, over the Midwest and Great Lakes region to 1) validate their performance in reproducing extreme daily precipitation (P) and daily maximum temperature (T max ) for summer and winter and 2) evaluate projections of extremes in the future. Our results show the HD statistical downscaling approach, which includes large-scale bias correction of GCM inputs, can reproduce observed extreme P and T max reasonably well for both summer and winter. However, raw historical WRF simulations show significant bias in both extreme P and T max for both seasons. Interestingly, the convection-permitting WRF simulation at 4-km grid spacing does not produce better results for seasonal extremes than the WRF simulation at 12 km using a parameterized convection scheme. Despite a broad similarity for winter extreme P projections, the projected changes in the future summer storms are quite different between downscaling methods; WRF simulations show substantial increases in summer extreme precipitation, while the changes projected by the HD approach exhibit moderate decreases overall. The WRF simulations at 4 km also show a pronounced decoupling effect between seasonal totals and extreme daily P for summer, which suggests that there could be more intense summer extremes at two different time scales, with more severe individual convective storms combined with longer summer droughts at the end of the twenty-first century.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1981

Approximate factorization for incompressible flow

For computational solution of the incompressible Navier-Stokes equations, the approximate factorization (AF) algorithm is used to solve the vectorized momentum equation in delta form based on the pressure calculated in the previous time step. The newly calculated velocities are substituted into the pressure equation (obtained from a linear combination of the continuity and momentum equation), which is then solved by means of line SOR. Computational results are presented for the NACA 66 sub 3 018 airfoil at Reynolds numbers of 1000 and 40,000 and attack angles of 0 and 6 degrees. Comparison with wind tunnel data for Re = 40,000 indicates good qualitative agreement between measured and calculated pressure distributions. Quantitative agreement is only fair, however, with the calculations somewhat displaced from the measurements. Furthermore, the computed velocity profiles are unrealistically thick around the airfoil, due to the excessive amount of artificial viscosity needed for stability. Based on the performance of the algorithm with regard to stability, it is concluded that AF/SOR is suitable for calculations at Reynolds numbers less than 10,000. Speedwise, the method is faster than point SOR by at least a factor of two.

Bernard, R. S.↗
DOE OSTI2024

Core-edge modeling of gyrokinetic turbulence by coupling the delta-f and total-f models in the XGC code

Whole device gyrokinetic turbulence is modeled using a delta-f model in the core coupled to a total-f model in the edge. This new model is implemented in the gyrokinetic code X-Point Included Gyrokinetic Code (XGC) and works within a single simulation. It employs a canonical Maxwellian background in the core and the usual local Maxwellian in the edge. These different core and edge backgrounds are coupled over an overlap region. A thermal-bath operator is used to keep δf small in the core.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1993

On the large eddy simulation of turbulent flows in complex geometry

Application of the method of Large Eddy Simulation (LES) to a turbulent flow consists of three separate steps. First, a filtering operation is performed on the Navier-Stokes equations to remove the small spatial scales. The resulting equations that describe the space time evolution of the 'large eddies' contain the subgrid-scale (sgs) stress tensor that describes the effect of the unresolved small scales on the resolved scales. The second step is the replacement of the sgs stress tensor by some expression involving the large scales - this is the problem of 'subgrid-scale modeling'. The final step is the numerical simulation of the resulting 'closed' equations for the large scale fields on a grid small enough to resolve the smallest of the large eddies, but still much larger than the fine scale structures at the Kolmogorov length. In dividing a turbulent flow field into 'large' and 'small' eddies, one presumes that a cut-off length delta can be sensibly chosen such that all fluctuations on a scale larger than delta are 'large eddies' and the remainder constitute the 'small scale' fluctuations. Typically, delta would be a length scale characterizing the smallest structures of interest in the flow. In an inhomogeneous flow, the 'sensible choice' for delta may vary significantly over the flow domain. For example, in a wall bounded turbulent flow, most statistical averages of interest vary much more rapidly with position near the wall than far away from it. Further, there are dynamically important organized structures near the wall on a scale much smaller than the boundary layer thickness. Therefore, the minimum size of eddies that need to be resolved is smaller near the wall. In general, for the LES of inhomogeneous flows, the width of the filtering kernel delta must be considered to be a function of position. If a filtering operation with a nonuniform filter width is performed on the Navier-Stokes equations, one does not in general get the standard large eddy equations. The complication is caused by the fact that a filtering operation with a nonuniform filter width in general does not commute with the operation of differentiation. This is one of the issues that we have looked at in detail as it is basic to any attempt at applying LES to complex geometry flows. Our principal findings are summarized.

Ghosal, Sandip↗
NASA NTRS2018

Modified Geometric Truncation of the Scattering Phase Function

Phase function of light scattering on large atmospheric particles has very strong peak in forward direction constituting a challenge for accurate numerical calculations of radiance required in remote sensing problems. Scaling transformation replaces original phase function with a sum of the delta function and a new regular smooth phase function. Geometric truncation is one of the ways to construct such a smooth function. The replacement phase function coincides with the original one outside the forward cone and preserves the asymmetry parameter. It has discontinuity at the cone. Another simple functional form of the replacement phase function within the cone is suggested. It enables continuity and allows for a number of modifications. Three of them are considered in this study: preserving asymmetry parameter, providing continuity of the 1st derivative of the phase function, and preserving mean scattering angle. Yet another problem addressed in this study is objective selection of the width of the forward cone. That angle affects truncation fraction and values of the phase function within the cone. A heuristic approach providing unambiguous criterion of selection of the truncation angle is proposed. The approach has easy numerical implementation. Suggested modifications were tested on cloud phase function using discrete ordinates and Monte Carlo methods. It was shown that the modifications provide better accuracy of the radiance computation compare to the original geometric truncation with discrete ordinates while continuous derivative approach provides significant gain in computer time with Monte Carlo simulations.

Alexander Radkevich↗
NASA NTRS2002

System and Method for Measuring the Transfer Function of a Guided Wave Device

A method/system are provided for measuring the NxN scalar transfer function elements for an N-port guided wave device. Optical energy of a selected wavelength is generated at a source and directed along N reference optical paths having N reference path lengths. Each reference optical path terminates in one of N detectors such that N reference signals are produced at the N detectors. The reference signals are indicative of amplitude, phase and frequency of the optical energy carried along the N reference optical paths. The optical energy from the source is also directed to the N-ports of the guided wave device and then on to each of the N detectors such that N measurement optical paths are defined between the source and each of the N detectors. A portion of the optical energy is modified in terms of at least one of the amplitude and phase to produce N modified signals at each of the N detectors. At each of the N detectors, each of the N modified signals is combined with a corresponding one of the N reference signals to produce corresponding N combined signals at each of the N detectors. A total of N(sup 2) measurement signals are generated by the N detectors. Each of the N(sup 2) measurement signals is sampled at a wave number increment (Delta)k so that N(sup 2) sampled signals are produced. The NxN transfer function elements are generated using the N(sup 2) sampled signals. Reference and measurement path length constraints are defined such that the N combined signals at each of the N detectors are spatially separated from one another in the time domain.

Mark E Froggatt↗
NASA NTRS1994

Steps toward determination of the size and structure of the broad-line region in active galactic nuclei. 5: Variability of the ultraviolet continuum and emission lines of NGC 3783

We report on the results of intensive ultraviolet spectral monitoring of the Seyfert 1 galaxy NGC 3783. The nucleus of NGC 3783 was observed with the International Ultraviolet Explorer satellite on a regular basis for a total of 7 months, once every 4 days for the first 172 days and once every other day for the final 50 days. Significant variability was observed in both continuum and emission-line fluxes. The light curves for the continuum fluxes exhibited two well-defined local minima or 'dips,' the first lasting is less than or approximately 20 days and the second is less than or approximately 4 days, with additional episodes of relatively rapid flickering of approximately the same amplitude. As in the case of NGC 5548 (the only other Seyfert galaxy that has been the subject of such an intensive, sustained monitoring effort), the largest continuum variations were seen at the shortest wavelengths, so that the continuum became 'harder' when brighter. The variations in the continuum occurred simultaneously at all wavelengths (delta(t) is less than 2 days). Generally, the amplitude of variability of the emission lines was lower than (or comparable to) that of the continuum. Apart from Mg II (which varied little) and N V (which is relatively weak and badly blended with Ly(alpha), the light curves of the emission lines are very similar to the continuum light curves, in each case with a small systematic delay or 'lag.' As for NGC 5548, the highest ionization lines seem to respond with shorter lags than the lower ionization lines. The lags found for NGC 3783 are considerably shorter than those obtained for NGC 5548, with values of (formally) approximately 0 days for He II + O III), and approximately 4 days for Ly(alpha) and C IV. The data further suggest lags of approximately 4 days for Si IV + O IV) and 8-30 days for Si III + C III). Mg II lagged the 1460 A continuum by approximately 9 days, although this result depends on the method of measuring the line flux and may in fact be due to variability of the underlying Fe II lines. Correlation analysis further shows that the power density spectrum contains substantial unresolved power over timescales of is less than or approximately 2 days, and that the character of the continuum variability may change with time.

Reichert, G. A.↗
NASA NTRS1995

Kinetics of cycle length dependence of ventricular repolarization: effect of autonomic blockade

INTRODUCTION: Beat-to-beat adaptation of ventricular repolarization duration to cardiac cycle length and autonomic activity has not been previously characterized in the spontaneously beating human heart. METHODS AND RESULTS: The ECG of 14 healthy subjects was recorded from the supine and upright positions. Autonomic blockade was accomplished by atropine and propranolol. RR and RT intervals were measured by a computer algorithm, and the impulse response (h) from RR to RT computed. In the supine position the maximal adjustment of the RT interval occurred in the first beat following a change in cycle length (hpeak = 17.8 +/- 1.6 msec/sec), but continued to be detectable for 3.8 seconds (2.9-4.7 sec). Propranolol attenuated the peak impulse response to 15.8 +/- 4.0 msec/sec (P = NS). In the standing position the peak impulse response was increased to 25.2 +/- 5.0 msec/sec (P = 0.004 vs supine), and the impulse response duration (hdur) shortened to 1.4 seconds (1.3-1.6). This was reversed by beta blockade (hpeak = 10.7 +/- 3.6 [P = 0.005 vs standing]; hdur = 5.5 sec [4.8-6.1]). Parasympathetic and combined autonomic blockade resulted in too little residual heart rate variability to estimate the impulse response accurately. The slope of the regression of delta RT and delta RR in the supine position was 0.0177 +/- 0.0016, which was closely correlated with the peak impulse response (r = 0.91). CONCLUSIONS: System identification techniques can assist in characterizing the cycle dependence of ventricular repolarization and may provide new insights into conditions associated with abnormal repolarization.

Clinical Trial↗
NASA NTRS2011

Technique for Solving Electrically Small to Large Structures for Broadband Applications

Fast iterative algorithms are often used for solving Method of Moments (MoM) systems, having a large number of unknowns, to determine current distribution and other parameters. The most commonly used fast methods include the fast multipole method (FMM), the precorrected fast Fourier transform (PFFT), and low-rank QR compression methods. These methods reduce the O(N) memory and time requirements to O(N log N) by compressing the dense MoM system so as to exploit the physics of Green s Function interactions. FFT-based techniques for solving such problems are efficient for spacefilling and uniform structures, but their performance substantially degrades for non-uniformly distributed structures due to the inherent need to employ a uniform global grid. FMM or QR techniques are better suited than FFT techniques; however, neither the FMM nor the QR technique can be used at all frequencies. This method has been developed to efficiently solve for a desired parameter of a system or device that can include both electrically large FMM elements, and electrically small QR elements. The system or device is set up as an oct-tree structure that can include regions of both the FMM type and the QR type. The system is enclosed with a cube at a 0- th level, splitting the cube at the 0-th level into eight child cubes. This forms cubes at a 1st level, recursively repeating the splitting process for cubes at successive levels until a desired number of levels is created. For each cube that is thus formed, neighbor lists and interaction lists are maintained. An iterative solver is then used to determine a first matrix vector product for any electrically large elements as well as a second matrix vector product for any electrically small elements that are included in the structure. These matrix vector products for the electrically large and small elements are combined, and a net delta for a combination of the matrix vector products is determined. The iteration continues until a net delta is obtained that is within the predefined limits. The matrix vector products that were last obtained are used to solve for the desired parameter. The solution for the desired parameter is then presented to a user in a tangible form; for example, on a display.

Jandhyala, Vikram↗
NASA NTRS2022

Validation of the Mars 2020 Fault Protection Design: Navigating the Infinity of the Off-Nominal

On July 30th 2020, the Mars 2020 mission successfully launched out of Cape Canaveral, Florida, passed through the Earth’s shadow, and began its short cruise to Mars. Less than seven months later, the Perseverance rover touched down safely in Jezero Crater to begin its ambitious mission that includes looking for signs of ancient life and collecting samples for future return to Earth. Getting to the successful landing, or “Tango Delta Nominal,” could not have been achieved without also considering the off-nominal. One of the teams supporting this ambitious mission is the fault protection (FP) team. This team is tasked with assessing the various failures, or faults, that could prevent mission success and with ensuring that the autonomous behaviors built into the software and hardware can detect faults and recover the vehicle to a safe state. As part of its charter, the FP team designed a test campaign to provide confidence in the system’s robustness to off-nominal scenarios across all of Mars 2020’s mission phases. The greatest challenge associated with designing such a validation campaign was reducing the infinite number of anomalous scenarios into a finite test suite. In addition, the tests needed to be executed efficiently in order to utilize the team’s limited test venue access, but still needed to maintain a level of rigor that guaranteed confidence in the test outcomes. Given that each test scenario generated massive amounts of data, the team also developed methods for quickly ascertaining whether the autonomous fault protection behaviors maintained vehicle safety in the presence of an anomaly. This paper summarizes the processes that the Mars 2020 fault protection team employed to execute its off-nominal validation campaign. It captures both the methods of generating a suite of off-nominal tests, as well as reducing it to a subset that can be realistically executed within schedule and resource constraints. It also describes the various processes and philosophies that the team utilized to execute the tests efficiently, including creating a standardized procedure template, keeping the test cases modular so that they could be easily interchanged, and capturing common fault injections in a change-controlled database. Finally, it will describe the tools and processes for assessing the test data, focusing in particular on a tool that evaluated vehicle state using “secondary” sources of data to validate that the software had truly configured the spacecraft to the expected safe state.

Morantz, Chaz↗
NASA NTRS1995

A Coupled-Cluster Study of the Molecular Structure, Vibrational Spectrum, and Heats of Formation of XONO2 (X=H, F, Cl)

The equilibrium structures, harmonic vibrational frequencies, dipole moments, and IR intensities of nitric acid, fluorine nitrate, and chlorine nitrate have been investigated by using the singles and doubles coupled cluster method that also includes a perturbational estimate of the effects of connected triple excitations, CCSD(T). A standard triple-zeta double-polarized basis set was utilized. The equilibrium geometries and vibrational spectra of HONO2 and ClONO2 are shown to be in excellent agreement with the available experimental data. The ab initio vibrational spectrum of FONO2 is also shown to be in excellent agreement with experiment. Unlike the FOOF and FON molecules, but similar to the cis- and trans-FONO molecules, FONO2 is shown to possess normal bond distances. The bonding in FONO2 is shown to be more similar to that in ClONO2 than that in HONO2, although there are still significant differences, especially in the partial atomic charges as deduced from Mulliken populations. This causes FONO2 to possess almost no dipole moment, which is very different to both HONO2 and ClONO2. By using large atomic natural orbital basis sets, CCSD(T) energies are computed for four isodesmic reactions in order to determine an accurate heat of formation for FONO2. Our best estimate for DELTA H(raised circle) (sub f,298) is 3.1 plus or minus 2.0 kcal/mol, indicating that the F-ONO2 bond energy is 31.3 kcal/mol.

Lee, Timothy J.↗
NASA NTRS1997

On Efficient Parallel Implementation of Moving Body Overset Grid Methods

An investigation into the parallel performance of moving-body overset grid methods will be presented. Parallel versions of the OVERFLOW flow solver, DCF3D domain connectivity software, and SIXDO six-degree-of-freedom routine are coupled with an automatic load balance routine and tested for 3D Navier-Stokes calculations on the IBM SP2. The primary source of parallel inefficiency in moving and problems are the domain connectivity costs with DCF 3D. Although this algorithm constitutes a relatively low fraction of the total solution cost (e.g. 10-20%) in calculations on serial machines, the consequently cause a significant degradation in the overall parallel performance. The paper will highlight some approaches for improving the scalability of DCF3D. The paper will present results of a proposed new load balancing scheme that seeks more equal distribution of the inter-grid boundary points in order to more evenly load balance the donor search costs associated with DCF3D. Some preliminary results will also be given from a new solution-adaption algorithm coupled with OVERFLOW which incorporates overset cartesian grids with various levels of refinement. The measured parallel performance from a descending delta-wing configuration and a generic store-separation from a wing/pylon case will be presented.

Wissink, Andrew M.↗
NASA NTRS2016

Emissions Relationships Among Western Forest Fire Plumes: II. Plume Typing for Sources of Ozone and Aerosol Absorption

Previous studies of emission factors from biomass burning are prone to largeerrors since they ignore the interplay of mixing and varying pre-fire backgroundCO2 levels. Such complications severely affected our studies of 446 forest fireplume samples measured in the Western US by the science teams of NASAs SEAC4RS and ARCTAS airborne missions. Consequently we propose a MixedEffects Regression Emission Technique (MERET) to check techniques like theNormalized Emission Ratio Method (NERM), where use of sequentialobservations cannot disentangle emissions and mixing. We also evaluate asimpler consensus technique. All techniques relate emissions to fuel burnedusing Cburn Ctot added to the fire plume, where Ctot (CO2 + CO). Mixed-effectsregression can estimate pre-fire background values of Ctot (indexed byobservation j) simultaneously with emissions factors indexed by individualspecies i, xi (Cburn )i,j., MERET and consensus require more than twoemissions indicators. Our studies excluded samples where exogenous CO orCH4 might have been fed into a fire plume, mimicking emission.We sought to let the data on 13 gases and particulate properties suggest clustersof variables and plume types, using non-negative matrix factorization (NMF).While samples were mixtures, the NMF unmixing suggested purer burn types.Particulate properties (bscat, babs, SSA, AE) and gas-phase emissions were interrelated.Finally, we sought a simple categorization useful for modeling ozone productionin plumes. Two kinds of fires produced high ozone: those with large fuel nitrogenas evidenced by remnant CH3CN in the plumes, and also those from veryintense large burns. Fire types with optimal ratios of delta-NOydelta-HCHO associate with the highest additional ozone per unit Cburn, Perhaps theseplumes exhibit limited NOx binding to reactive organics. Perhaps these plumesexhibit limited NOx binding to reactive organics.

ozone↗
NASA NTRS1994

Ultrasound as a Noninvasive Method to Assess Changes of Intracranial Volume and Pressure During Simulated Microgravity

Headaches are commonly experienced by astronauts in microgravity and by subjects undergoing head-down tilt (simulated microgravity on Earth). Exposure to microgravity probably elevates blood pressure and flow in the head which may increase intracranial volume (ICV) and pressure (ICP) and in turn cause headache. Due to the slightly compliant nature of the cranial vault and the encasement of brain and its vasculature within this vault, any increase of ICV will increase ICP and slightly distend the cranium. Previous studies document perivascular edema and increased ICP in rhesus monkeys during head-down tilt. Elevated ICP has also been reported in humans during head-down tilt. ICP measurements in healthy humans are rare because of the invasiveness of currently-available measurement techniques. Therefore, we proposed a noninvasive ultrasound technique to assess changes of ICV and JCP. The ultrasound principle is based on compliance of the cranial vault. A 450 kHz ultrasound stimulus is transmitted through the cranium by a transducer every 7.5-10 msec. The ultrasound wave enters the brain tissue, reflects off the opposite side of the cranium and is received by the same transducer. The detected wave is compared for phase quadrature (90 deg.to transmitted wave). Because the electronic circuitry of the device maintains a 90 deg. phase (phi), any alterations in the detected wave caused by an increase of ICV and ICP will be reflected as a change in the wave frequency. Phase shift is directly proportional to path length of the wave, DELTA x, which is expressed as DELTA x = phi lambda/2 pi where lambda is wavelength. Elevated ICV and ICP expand the cranial vault and increase path length of the wave (a measure of intracranial distance). Increased path length equals reduced frequency of the detected wave. Reduced frequency is then related to elevated ICP. This technique has potential uses for ICP studies of astronauts in space and head trauma patients on Earth.

Murthy, G.↗
DOE OSTI2020

Application and Certification of Comparative Vacuum Monitoring Sensors for Structural Health Monitoring of 737 Wing Box Fittings

Multi-site fatigue damage, hidden cracks in hard-to-reach locations, disbonded joints, erosion, impact, and corrosion are among the major flaws encountered in today's extensive fleet of aging aircraft and space vehicles. The use of in-situ sensors for real-time health monitoring of aircraft structures are a viable option to overcome inspection impediments stemming from accessibility limitations, complex geometries, and the location and depth of hidden damage. Reliable, structural health monitoring systems can automatically process data, assess structural condition, and signal the need for human intervention. Prevention of unexpected flaw growth and structural failure can be improved if on-board health monitoring systems are used to continuously assess structural integrity. Such systems are able to detect incipient damage before catastrophic failures occurs. Condition-based maintenance practices could be substituted for the current time-based maintenance approach. Other advantages of on-board distributed sensor systems are that they can eliminate costly, and potentially damaging, disassembly, improve sensitivity by producing optimum placement of sensors and decrease maintenance costs by eliminating more time- consuming manual inspections. This report presents a Sandia Labs-aviation industry effort to move SHM into routine use for aircraft maintenance. This program addressed formal SHM technology validation and certification issues so that the full spectrum of concerns, including design, deployment, performance and certification were appropriately considered. The Airworthiness Assurance NDI Validation Center (AANC) at Sandia Labs, in conjunction with Boeing, Delta Air Lines, Structural Monitoring Systems Ltd., Anodyne Electronics Manufacturing Corp. and the Federal Aviation Administration (FAA) carried out a certification program to formally introduce Comparative Vacuum Monitoring (CVM) as a structural health monitoring solution to a specific aircraft wing box application. Validation tasks were designed to address the SHM equipment, the health monitoring task, the resolution required, the sensor interrogation procedures, the conditions under which the monitoring will occur, the potential inspector population, adoption of CVM into an airline maintenance program and the document revisions necessary to allow for routine use of CVM as an alternate means of performing periodic structural inspects. To carry out the validation process, knowledge of aircraft maintenance practices was coupled with an unbiased, independent evaluation. Sandia Labs designed, implemented, and analyzed the results from a focused and statistically-relevant experimental effort to quantify the reliability of the CVM system applied to the Boeing 737 Wing Box fitting application. All factors that affect SHM sensitivity were included in this program: flaw size, shape, orientation and location relative to the sensors, as well as operational and environmental variables. Statistical methods were applied to performance data to derive Probability of Detection (POD) values for CVM sensors in a manner that agrees with current nondestructive inspection (NDI) validation requirements and also is acceptable to both the aviation industry and regulatory bodies. This report presents the use of several different statistical methods, some of them adapted from NDI performance assessments and some proposed to address the unique nature of damage detection via SHM systems, and discusses how they can converge to produce a confident quantification of SHM performance An important element in developing SHM validation processes is a clear understanding of the regulatory measures needed to adopt SHM solutions along with the knowledge of the structural and maintenance characteristics that may impact the operational performance of an SHM system. This report describes the major elements of an SHM validation approach and differentiates the SHM elements from those found in NDI validation. The activities conducted in this program demonstrated the feasibility of routine SHM usage in general and CVM in particular for the application selected. They also helped establish an optimum OEM-airline-regulator process and determined how to safely adopt SHM solutions. This formal SHM validation will allow aircraft manufacturers and airlines to confidently make informed decisions about the proper utilization of CVM technology. It will also streamline the regulatory actions and formal certification measures needed to assure the safe application of SHM solutions.

42 ENGINEERING↗