SEARCH · Engineering Papers

Results for “differential geometry”

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.

641 records · Page 36

All sources NASA NTRS DOE OSTI

DOE OSTI2026

From Detector Layout to Signal Analysis: Geometry Optimization and Neutron-Gamma Tagging in Plastic Scintillator Detectors

Neutrinos are elementary particles with many properties still unknown. Their masses so far have only upper and lower limits. Still, due to neutrino oscillations, it is clear that they are not massless, as stated by the Standard Model of Elementary Particles. Neutrinos are also present in the Universe in vast amounts, but they rarely interact with the surrounding matter. Their abundance makes them very interesting for many theories beyond the Standard Model, e.g., dark matter searches and CP symmetry violation in the leptonic sector, which could be (partially) responsible for the observed matter-antimatter asymmetry in today's Universe.\\ The Deep Underground Neutrino Experiment (DUNE) is a next-generation accelerator-based neutrino oscillation experiment that will study neutrinos with unprecedented precision and may answer many open questions. DUNE uses a powerful neutrino beam from Fermilab. It consists of a Near Detector (ND) complex to measure neutrinos before oscillat ion, and a Far Detector (FD) complex $1300\ \mathrm{km}$ away to measure them after oscillation. As part of the ND complex, measurements are also possible with different angles to the neutrino beam. This enables excellent control of systematic uncertainties of e.g., neutrino cross section measurements.\\ One detector in the near detector complex is The Muon Spectrometer (TMS), an extension of a Liquid Argon (LAr) detector that measures the charge and momentum of muons produced in neutrino interactions within the LAr. The design of this detector, which consists of alternating layers of steel and plastic scintillator bars, must be optimized for the expected muon energies.\\ In this thesis, a study of the optimal module orientation plan is presented, which is necessary for the physics performance of the near detector complex and, by extension, DUNE. As part of this study, the event reconstruction was also developed and improved. Simulated muons are then reconstructed, and the performan ce of different module orientation plans is tested.\\ As a second part, a study of neutron and gamma tagging using a Pulse Shape Discriminating (PSD) plastic scintillator is presented. The properties of this material allow particle differentiation based on the temporal distribution of emitted light. A novel approach to using the individual light signals was successfully tested using data from a small, local test setup.

Nehm, Asa [Mainz U.]↗

DOE OSTI2025

Inverse Trans Influence and Uranium-Arene σ-Bonding Drive Molecular Geometry: Ligand Modification from Hard to Soft Flips the Oxide from Axial to Equatorial

A rare example of an equatorially bound terminal uranium(V) oxo complex in a chelating sulfur-based ligand environment, namely [(mes( Me,Ad ArS) 3 )U V (O eq )(THF)] (2), is presented. Octahedrally coordinated 2 is obtained by reaction of the mesitylene-anchored tris-thiophenolate-coordinated uranium(III) complex [U III ((SAr Ad,Me ) 3 mes)] (1) with the oxygen-atom transfer reagent N 2 O. The observed, equatorially bound oxo ligand in 2 is in stark contrast to its known tris-aryloxide analog, [(mes( Me,Ad ArO) 3 )U V (O ax )(THF)] (A), where the oxo ligand occupies the typically observed axial coordination site. Complexes 1 and 2 are characterized by single-crystal X-ray diffraction analyses and spectroscopic and magnetochemical methods, including 1 H NMR, UV/vis/NIR electronic absorption, as well as EPR spectroscopy and SQUID magnetometry, thus confirming the C S symmetry and the pentavalent oxidation state of 2. Encompassing quantum chemical calculations (DFT and CASPT2) on 2 and its tris-phenolate analog A, support and rationalize the structural and electronic differences. The molecular orbital pictures show that a stabilizing σ-bonding interaction arising from the U–O eq inverse trans influence (ITI) is present in 2 but missing in A. In 2, the sulfur 3p orbitals are closer in energy to the uranium 5f manifold than the arene π-system, leading to an ITI, while U–arene σ- or δ-bonding is not observed. Although the arene orbitals remain separated from the uranium 5f orbitals in A, the absence of an ITI allows the arene a 2u orbital to engage in a σ-type interaction with the metal. Thus, incorporating a tris-thiophenolate to an arene anchor introduces a new design concept in molecular f-element chemistry. This approach stabilizes an equatorially bound U(V) oxo center, contrasting with its tris-phenolate counterpart, where oxo coordination is axial. The observed geometric divergence, driven by competing ITI and U–arene interactions, not only tunes electronic structure but also leads to differentiated reactivity: only the phenolate analogs activate H 2 O, while the thiolates do not.

Hydrocarbons↗

NASA NTRS1998

Trellises and Trellis-Based Decoding Algorithms for Linear Block Codes

A code trellis is a graphical representation of a code, block or convolutional, in which every path represents a codeword (or a code sequence for a convolutional code). This representation makes it possible to implement Maximum Likelihood Decoding (MLD) of a code with reduced decoding complexity. The most well known trellis-based MLD algorithm is the Viterbi algorithm. The trellis representation was first introduced and used for convolutional codes [23]. This representation, together with the Viterbi decoding algorithm, has resulted in a wide range of applications of convolutional codes for error control in digital communications over the last two decades. There are two major reasons for this inactive period of research in this area. First, most coding theorists at that time believed that block codes did not have simple trellis structure like convolutional codes and maximum likelihood decoding of linear block codes using the Viterbi algorithm was practically impossible, except for very short block codes. Second, since almost all of the linear block codes are constructed algebraically or based on finite geometries, it was the belief of many coding theorists that algebraic decoding was the only way to decode these codes. These two reasons seriously hindered the development of efficient soft-decision decoding methods for linear block codes and their applications to error control in digital communications. This led to a general belief that block codes are inferior to convolutional codes and hence, that they were not useful. Chapter 2 gives a brief review of linear block codes. The goal is to provide the essential background material for the development of trellis structure and trellis-based decoding algorithms for linear block codes in the later chapters. Chapters 3 through 6 present the fundamental concepts, finite-state machine model, state space formulation, basic structural properties, state labeling, construction procedures, complexity, minimality, and sectionalization of trellises. Chapter 7 discusses trellis decomposition and subtrellises for low-weight codewords. Chapter 8 first presents well known methods for constructing long powerful codes from short component codes or component codes of smaller dimensions, and then provides methods for constructing their trellises which include Shannon and Cartesian product techniques. Chapter 9 deals with convolutional codes, puncturing, zero-tail termination and tail-biting.Chapters 10 through 13 present various trellis-based decoding algorithms, old and new. Chapter 10 first discusses the application of the well known Viterbi decoding algorithm to linear block codes, optimum sectionalization of a code trellis to minimize computation complexity, and design issues for IC (integrated circuit) implementation of a Viterbi decoder. Then it presents a new decoding algorithm for convolutional codes, named Differential Trellis Decoding (DTD) algorithm. Chapter 12 presents a suboptimum reliability-based iterative decoding algorithm with a low-weight trellis search for the most likely codeword. This decoding algorithm provides a good trade-off between error performance and decoding complexity. All the decoding algorithms presented in Chapters 10 through 12 are devised to minimize word error probability. Chapter 13 presents decoding algorithms that minimize bit error probability and provide the corresponding soft (reliability) information at the output of the decoder. Decoding algorithms presented are the MAP (maximum a posteriori probability) decoding algorithm and the Soft-Output Viterbi Algorithm (SOVA) algorithm. Finally, the minimization of bit error probability in trellis-based MLD is discussed.

Lin, Shu↗

DOE OSTI2022

Comprehensive uncertainty quantification (UQ) for full engineering models by solving probability density function (PDF) equation

This report details a new method for propagating parameter uncertainty (forward uncertainty quantification) in partial differential equations (PDE) based computational mechanics applications. The method provides full-field quantities of interest by solving for the joint probability density function (PDF) equations which are implied by the PDEs with uncertain parameters. Full-field uncertainty quantification enables the design of complex systems where quantities of interest, such as failure points, are not known apriori. The method, motivated by the well-known probability density function (PDF) propagation method of turbulence modeling, uses an ensemble of solutions to provide the joint PDF of desired quantities at every point in the domain. A small subset of the ensemble is computed exactly, and the remainder of the samples are computed with approximation of the driving (dynamics) term of the PDEs based on those exact solutions. Although the proposed method has commonalities with traditional interpolatory stochastic collocation methods applied directly to quantities of interest, it is distinct and exploits the parameter dependence and smoothness of the dynamics term of the governing PDEs. The efficacy of the method is demonstrated by applying it to two target problems: solid mechanics explicit dynamics with uncertain material model parameters, and reacting hypersonic fluid mechanics with uncertain chemical kinetic rate parameters. A minimally invasive implementation of the method for representative codes SPARC (reacting hypersonics) and NimbleSM (finite- element solid mechanics) and associated software details are described. For solid mechanics demonstration problems the method shows order of magnitudes improvement in accuracy over traditional stochastic collocation. For the reacting hypersonics problem, the method is implemented as a streamline integration and results show very good accuracy for the approximate sample solutions of re-entry flow past the Apollo capsule geometry at Mach 30.

42 ENGINEERING↗

DOE OSTI2025

Coupled geomechanical investigation of depletion-induced fault reactivation

Fault reactivation during subsurface fluid production pose significant challenges to safe and sustainable resource extraction. Here, this study presents a three-dimensional coupled geomechanical framework to investigate the processes driving fault reactivation, capturing the interactions between reservoir dynamics and geomechanical responses. Verification against theoretical estimations based on linear poroelasticity confirms the model's capacity in representing reservoir background stress responses. However, the study reveals that relying solely on background stress states can underestimate or overestimate fault reactivation potential, emphasizing the importance of including localized stress perturbations such as differential compaction and stress redistribution. Applied to a fault (M1) inspired by the geological characteristics of the Groningen field, the model shows slip initiation at 2965 m depth with 16.0 MPa depletion, aligning with field observations where seismicity occurred at approximately 3 km depth after 15.8 MPa depletion. Parametric studies reveal: (1) inelastic reservoir compaction delays fault reactivation and mitigates fault slip by reducing stress concentration, (2) higher intermediate in-situ stress magnitudes decrease the Coulomb Failure Stress (CFS) increase rate and reduce fault slip, (3) larger fault offsets amplify shear stress near the offset zone, promoting earlier reactivation and longer rupture propagation, and (4) fault permeability significantly influences pressure diffusion, with low-permeability faults leading to sharper stress changes and earlier fault destabilization. These insights highlight the critical role of geological and mechanical parameters in fault reactivation and provide a predictive framework for mitigating induced seismicity risks.

Coupled geomechanical modeling↗

NASA NTRS1993

Structured adaptive grid generation using algebraic methods

The accuracy of the numerical algorithm depends not only on the formal order of approximation but also on the distribution of grid points in the computational domain. Grid adaptation is a procedure which allows optimal grid redistribution as the solution progresses. It offers the prospect of accurate flow field simulations without the use of an excessively timely, computationally expensive, grid. Grid adaptive schemes are divided into two basic categories: differential and algebraic. The differential method is based on a variational approach where a function which contains a measure of grid smoothness, orthogonality and volume variation is minimized by using a variational principle. This approach provided a solid mathematical basis for the adaptive method, but the Euler-Lagrange equations must be solved in addition to the original governing equations. On the other hand, the algebraic method requires much less computational effort, but the grid may not be smooth. The algebraic techniques are based on devising an algorithm where the grid movement is governed by estimates of the local error in the numerical solution. This is achieved by requiring the points in the large error regions to attract other points and points in the low error region to repel other points. The development of a fast, efficient, and robust algebraic adaptive algorithm for structured flow simulation applications is presented. This development is accomplished in a three step process. The first step is to define an adaptive weighting mesh (distribution mesh) on the basis of the equidistribution law applied to the flow field solution. The second, and probably the most crucial step, is to redistribute grid points in the computational domain according to the aforementioned weighting mesh. The third and the last step is to reevaluate the flow property by an appropriate search/interpolate scheme at the new grid locations. The adaptive weighting mesh provides the information on the desired concentration of points to the grid redistribution scheme. The evaluation of the weighting mesh is accomplished by utilizing the weight function representing the solution variation and the equidistribution law. The selection of the weight function plays a key role in grid adaptation. A new weight function utilizing a properly weighted boolean sum of various flowfield characteristics is defined. The redistribution scheme is developed utilizing Non-Uniform Rational B-Splines (NURBS) representation. The application of NURBS representation results in a well distributed smooth grid by maintaining the fidelity of the geometry associated with boundary curves. Several algebraic methods are applied to smooth and/or nearly orthogonalize the grid lines. An elliptic solver is utilized to smooth the grid lines if there are grid crossings. Various computational examples of practical interest are presented to demonstrate the success of these methods.

Yang, Jiann-Cherng↗

DOE OSTI2025

An implementation of a high-order generalized finite difference method for solving the time-harmonic cold plasma wave equation in toroidal geometry

A high-order physics-informed meshless finite difference numerical technique is introduced for solving the time-harmonic cold plasma wave equation in toroidal geometries, presenting a novel application of the generalized finite difference (GFD) method to plasma wave simulations. The algorithm employs an irregular distribution of computational points, with local point density informed by the shortest wavelength derived from the cold plasma dispersion relation. Numerical stability and robustness are addressed using regularization techniques. The algorithm, implemented for two spatial dimensions, solves for the wave electric field and is demonstrated to achieve convergence rates of $\mathcal{O}$($\mathcal{h}$ $\mathcal{P}$ )⁠. Verification tests reproduce plane wave solutions, and example simulations of ion cyclotron resonance heating and electron cyclotron resonance heating demonstrate its capability, approaching realistic tokamak plasma scenarios. This work contributes to laying a foundation for the GFD method to be used in more sophisticated, optimized, and physically realistic full-wave simulations in time-harmonic plasma wave research.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Performance Evaluation of AG-1 FC HEPA Filters and Medium in Nuclear Complex Facilities - 20434

High Efficiency Particulate Air (HEPA) filters are credited as the final barrier against the release of radioactive aerosol contamination in nearly every operating U. S. Department of Energy (DOE) and National Nuclear Security Agency (NNSA) nuclear facility. The Institute for Clean Energy Technology (ICET) at Mississippi State University maintains a research program that studies various aspects of these components of containment systems and seeks to answer key questions from across the industry. This programmatic overview will include results from recent and ongoing studies, including a study on filter design and performance envelope, degradation due to aging or fatigue, fire event impact on HEPA filters, and accelerated aging of medium and component parts. ICET has completed a study helping to define the threshold for combined elevated temperature and relative humidity resistance of separator style and separator-less style HEPA filters that are used in DOE complexes. The study provides experimental data that considers elevated temperature, elevated relative humidity (RH), and target differential pressure (dP) in an attempt to gain a better understanding or more comprehensive insight into the operating envelope of different configurations of ASME (American Society for Mechanical Engineers) AG-1 Section FC filters. Test variables including air temperature ranges of 48.88 deg. C, 54.44 deg. C, or 60 deg. C (120, 130, or 140 deg. F); air RH ranges including 60-70%, 80%, or 90+% with initial filter dP of either 497.68 Pa, 746.52 Pa, or 995.36 Pa (2 inches water column (in. w. c.), 3 in. w. c., or 4 in. w. c.) included for full bracketing of each set of conditions. A study examining the effects of aging and fatigue on nuclear grade HEPA filters and medium elucidates the physical properties of media along with testing of new and aged ASME AG-1 HEPA filters. In addition to performing autopsies of the tested, aged filters to help better understand the service life of their performance, newly designed accelerated aging chambers allow new medium to undergo accelerated aging and exposure treatments. Evaluation of degradation in physical properties and functionality requires a high population of aged filters of different ages, designs, manufacturers, and operational histories. The limited availability of this population and the impropriety of looking solely at formulations and components from past years requires experimental design that will provide foresight into future filter performance. This undertaking therefore involves not only properties analysis of aged medium from different operational histories, ages, designs, and manufacturers, but also prescience in analytical goals via accelerated aging studies of newly manufactured media. A more complete understanding of the mechanisms of degradation from past formulations as well as current and future formulations is possible. Accelerated Aging of newly manufactured media utilizes exposure treatments based on the Arrhenius equation to artificially age medium and therefore glimpse into the future. Evaluation of current media properties This testing allows for the comparison of performance and durability of new filters under upset or design basis conditions with aged filters that were in service under ambient conditions and other aged filters retained in storage. Susceptibility of HEPA filters to the effects of combustion byproducts, heat, and water is also being studied in an attempt to gain understanding on the prevention of filter failures during fire events in nuclear facilities. ICET seeks to determine the effect of filter performance due to smoke loading using characterized smoke from variable fuel compositions, geometries, and loads; the smoke capacity of various filters with different burning conditions, heat release rates and transport, as well as mass transport from the air stream will be studied in low flow containment systems. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

NASA NTRS1999

Effects of Melt Processing on Evolution of Structure in PEEK

We report on the effects of melt processing temperature on structure formation in Poly(ether-ether-ketone), PEEK. Real time Small Angle X-ray Scattering, SAXS, and thermal analysis are used to follow the melting behavior after various stages of processing. Assignment of peaks to structural entities within the material, the relative perfection of the crystals, and the possibility of their reorganization, are all influenced by the melt processing history. With the advent of high intensity synchrotron sources of X-radiation, polymer scientists gain a research tool which, when used along with thermal analysis, provides additional structural information about the crystals during growth and subsequent melting. PEEK is an engineering thermoplastic polymer with a very high glass transition temperature (145 C) and crystal melting point (337 C). PEEK has been the subject of recent studies by X-ray scattering in which melt and cold crystallization were followed in real-time. X-ray scattering and thermal studies have been used to address the formation of dual endothermic response which has been variously ascribed to lamellar insertion, dual crystal populations, or melting followed by re-crystallization. Another important issue is whether all of the amorphous phase is located in interlamellar regions, or alternatively whether some is located in "pockets" away from the crystalline lamellar stacks. The interpretation of scattering from lamellar stacks varies depending upon whether such amorphous pockets are formed. Some groups believe all of the amorphous phase is interlamellar. This leads to selection of a smaller thickness for the crystals. Other groups suggest that most amorphous phase is not interlamellar, and this leads to the suggestion that the crystal thickness is larger than the amorphous layer within the stacks. To investigate these ideas, we used SAXS and Differential Scanning Calorimetry to compare results of single and dual stage melt crystallization of PEEK using a treatment scheme involving annealing/crystallization at T(sub a1) followed by annealing at T(sub a2) where either T(sub a1) < T(sub a2) or T(sub a1) > T(sub a2). We proposed a model to explain multiple melting endotherms in PPS, treated according to one or two-stage melt or cold crystallization. Key features of this model are that multiple endotherms: (1) are due to reorganization/recrystallization after cold crystallization; and, (2) are dominated by crystal morphology after melt crystallization at high T. In other words, multiple distinct crystal populations are formed by the latter treatment, leading to observation of multiple melting. PEEK 45OG pellets (ICI Americas) were the starting material for this study. Films were compression molded at 400 C, then quenched to ice water. Samples were heated to 375 C in a Mettler FP80 hot stage and held for three min. to erase crystal seeds before cooling them to T(sub a1) = 280 C . Samples were held at T(sub a2) for a period of time, then immediately heated to 360 C. In the second treatment samples were held at T(sub a1) = 31 C for different crystallization times t(sub c) then cooled to 295 C and held 15 min. In situ (SAXS) experiments were performed at the Brookhaven National Synchrotron Light Source with the sample located inside the Mettler hot stage. The system was equipped with a two-dimensional position sensitive detector. The sample to detector distance was 172.7 cm and the X-ray wavelength was 1.54 Angstroms. SAXS data were taken continuously during the isothermal periods and during the heating to 360 C at 5 C/min. Each SAXS scan was collected for 30 sec. Since the samples were isotropic, circular integration was used to increase the signal to noise ratio. After dual stage melt crystallization with T(sub a1) < T(sub a2) the lower melting endotherm arises from holding at T(sub a1). During cooling a broad distribution of crystals forms, and the low-melting tail is perfected during T(sub al). Heating to T(sub a2) melts these imperfect crystals and allows others with greater average long spacing to form in their place. After dual stage crystallization with T(sub a1) > T(sub a2), the amount of space remaining for additional growth at T(sub a2) depends upon the holding time at T(sub a1). The long period of crystals formed at T(sub a2) is smaller than that formed at T(sub a1) due to growth in a now-restricted geometry. Perfection of crystals is seen as an increase of the intensity of the population scattering at higher s, while the intensity of the population scattering at lower s stays constant. During heating from below to above the minor endotherm, we see rapid decrease of the intensity of the X-ray scattering corresponding to the population of crystals scattering in the shoulder. Another important observation is that after the sample is annealed at 295 C, the shoulder intensity is restored once again. The population scattering at higher s remains longer before it disappears in the sample treated to the second stage of melt crystallization, compared to the sample crystallized with a single stage. This could be interpreted as an effect of continued perfection of the less perfect population, which is also reflected in the increased melting temperature of the smaller endotherm as the holding time at 295 C increases. In the corresponding DSC scans we see a shift in the area and the peak temperature of the lower melting endotherm with an increase of the annealing time.

Georgiev, Georgi↗

NASA NTRS2008

Spiking Neurons for Analysis of Patterns

Artificial neural networks comprising spiking neurons of a novel type have been conceived as improved pattern-analysis and pattern-recognition computational systems. These neurons are represented by a mathematical model denoted the state-variable model (SVM), which among other things, exploits a computational parallelism inherent in spiking-neuron geometry. Networks of SVM neurons offer advantages of speed and computational efficiency, relative to traditional artificial neural networks. The SVM also overcomes some of the limitations of prior spiking-neuron models. There are numerous potential pattern-recognition, tracking, and data-reduction (data preprocessing) applications for these SVM neural networks on Earth and in exploration of remote planets. Spiking neurons imitate biological neurons more closely than do the neurons of traditional artificial neural networks. A spiking neuron includes a central cell body (soma) surrounded by a tree-like interconnection network (dendrites). Spiking neurons are so named because they generate trains of output pulses (spikes) in response to inputs received from sensors or from other neurons. They gain their speed advantage over traditional neural networks by using the timing of individual spikes for computation, whereas traditional artificial neurons use averages of activity levels over time. Moreover, spiking neurons use the delays inherent in dendritic processing in order to efficiently encode the information content of incoming signals. Because traditional artificial neurons fail to capture this encoding, they have less processing capability, and so it is necessary to use more gates when implementing traditional artificial neurons in electronic circuitry. Such higher-order functions as dynamic tasking are effected by use of pools (collections) of spiking neurons interconnected by spike-transmitting fibers. The SVM includes adaptive thresholds and submodels of transport of ions (in imitation of such transport in biological neurons). These features enable the neurons to adapt their responses to high-rate inputs from sensors, and to adapt their firing thresholds to mitigate noise or effects of potential sensor failure. The mathematical derivation of the SVM starts from a prior model, known in the art as the point soma model, which captures all of the salient properties of neuronal response while keeping the computational cost low. The point-soma latency time is modified to be an exponentially decaying function of the strength of the applied potential. Choosing computational efficiency over biological fidelity, the dendrites surrounding a neuron are represented by simplified compartmental submodels and there are no dendritic spines. Updates to the dendritic potential, calcium-ion concentrations and conductances, and potassium-ion conductances are done by use of equations similar to those of the point soma. Diffusion processes in dendrites are modeled by averaging among nearest-neighbor compartments. Inputs to each of the dendritic compartments come from sensors. Alternatively or in addition, when an affected neuron is part of a pool, inputs can come from other spiking neurons. At present, SVM neural networks are implemented by computational simulation, using algorithms that encode the SVM and its submodels. However, it should be possible to implement these neural networks in hardware: The differential equations for the dendritic and cellular processes in the SVM model of spiking neurons map to equivalent circuits that can be implemented directly in analog very-large-scale integrated (VLSI) circuits.

Huntsberger, Terrance↗

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↗