Engineering Papers⌕ Search

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.

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

1,020 records · Page 57
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Technical Note: On the impact of the incident electron beam energy on the primary dose component of flattening filter free photon beams

Purpose: For commercially available linear accelerators (Linacs), the electron energies of flattening filter free (FFF) and flattened (FF) beams are either identical or the electron energy of the FFF beam is increased to match the percentage depth dose curve (PDD) of the FF beam (in reference geometry). This study focuses on the primary dose components of FFF beams for both kinds of settings, studied on the same Linac. Methods: The measurements were conducted on a VersaHD Linac (Elekta, Crawley, UK) for both FF and FFF beams with nominal energies of 6 and 10 MV. In the clinical setting of the VersaHD, the energy of FFF{sub M} (Matched) beams is set to match the PDDs of the FF beams. In contrast the incident electron beam of the FFF{sub U} beam was set to the same energy as for the FF beam. Half value layers (HVLs) and a dual parameter beam quality specifier (DPBQS) were determined. Results: For the 6 MV FFF{sub M} beam, HVL and DPBQS values were very similar compared to those of the 6 MV FF beam, while for the 10 MV FFF{sub M} and FF beams, only %dd(10){sub x} and HVL values were comparable (differences below 1.5%). This shows that matching the PDD at one depth does not guarantee other beam quality dependent parameters to be matched. For FFF{sub U} beams, all investigated beam quality specifiers were significantly different compared to those for FF beams of the same nominal accelerator potential. The DPBQS of the 6 MV FF and FFF{sub M} beams was equal within the measurement uncertainty and was comparable to published data of a machine with similar TPR{sub 20,10} and %dd(10){sub x}. In contrast to that, the DPBQS’s two parameters of the 10 MV FFF{sub M} beam were substantially higher compared to those for the 10 MV FF beam. Conclusions: PDD-matched FF and FFF beams of both nominal accelerator potentials were observed to have similar HVL values, indicating similarity of their primary dose components. Using the DPBQS revealed that the mean attenuation coefficient was found to be the same within the uncertainty of 0.8% for 6 MV FF and 6 MV FFF{sub M} beams, while for 10 MV beams, they differed by 6.4%. This shows that the DPBQS can provide a differentiation of photon beam characteristics that would remain hidden by the use of a single beam quality specifier, such as %dd(10){sub x} or HVL.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2012

High-Order Curvilinear Finite Element Methods for Lagrangian Hydrodynamics [High Order Curvilinear Finite Elements for Lagrangian Hydrodynamics]

The numerical approximation of the Euler equations of gas dynamics in a movingLagrangian frame is at the heart of many multiphysics simulation algorithms. Here, we present a general framework for high-order Lagrangian discretization of these compressible shock hydrodynamics equations using curvilinear finite elements. This method is an extension of the approach outlined in [Dobrev et al., Internat. J. Numer. Methods Fluids, 65 (2010), pp. 1295--1310] and can be formulated for any finite dimensional approximation of the kinematic and thermodynamic fields, including generic finite elements on two- and three-dimensional meshes with triangular, quadrilateral, tetrahedral, or hexahedral zones. We discretize the kinematic variables of position and velocity using a continuous high-order basis function expansion of arbitrary polynomial degree which is obtained via a corresponding high-order parametric mapping from a standard reference element. This enables the use of curvilinear zone geometry, higher-order approximations for fields within a zone, and a pointwise definition of mass conservation which we refer to as strong mass conservation. Moreover, we discretize the internal energy using a piecewise discontinuous high-order basis function expansion which is also of arbitrary polynomial degree. This facilitates multimaterial hydrodynamics by treating material properties, such as equations of state and constitutive models, as piecewise discontinuous functions which vary within a zone. To satisfy the Rankine--Hugoniot jump conditions at a shock boundary and generate the appropriate entropy, we introduce a general tensor artificial viscosity which takes advantage of the high-order kinematic and thermodynamic information available in each zone. Finally, we apply a generic high-order time discretization process to the semidiscrete equations to develop the fully discrete numerical algorithm. Our method can be viewed as the high-order generalization of the so-called staggered-grid hydrodynamics (SGH) approach and we show that under specific low-order assumptions, we exactly recover the classical SGH method. We also present numerical results from an extensive series of verification tests that demonstrate several important practical advantages of using high-order finite elements in this context.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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 OSTI2016

EMAPS: An Efficient Multiscale Approach to Plasma Systems with Non-MHD Scale Effects

Global interactions of energetic ions with magnetoplasmas and neutral gases lie at the core of many space and laboratory plasma phenomena ranging from solar wind entry into and transport within planetary magnetospheres and exospheres to fast-ion driven instabilities in fusion devices to astrophysics-in-lab experiments. The ability of computational models to properly account for physical effects that underlie such interactions, namely ion kinetic, ion cyclotron, Hall, collisional and ionization processes is important for the success and planning of experimental research in plasma physics. Understanding the physics of energetic ions, in particular their nonlinear resonance interactions with Alfvén waves, is central to improving the heating performance of magnetically confined plasmas for future energy generation. Fluid models are not adequate for high-beta plasmas as they cannot fully capture ion kinetic and cyclotron physics (e.g., ion behavior in the presence of magnetic nulls, shock structures, plasma interpenetration, etc.). Recent results from global reconnection simulations show that even in a MHD-like regime there may be significant differences between kinetic and MHD simulations. Therefore, kinetic modeling becomes essential for meeting modern day challenges in plasma physics. The hybrid approximation is an intermediate approximation between the fluid and fully kinetic approximations. It eliminates light waves, removes the electron inertial temporal and spatial scales from the problem and enables full-orbit ion kinetics. As a result, hybrid codes have become effective tools for exploring ion-scale driven phenomena associated with ion beams, shocks, reconnection and turbulence that control the large-scale behavior of laboratory and space magnetoplasmas. A number of numerical issues, however, make three-dimensional (3D) large-scale hybrid simulations of inhomogeneous magnetized plasmas prohibitively expensive or even impossible. To resolve these difficulties we have developed a novel Event-driven Multiscale Asynchronous Parallel Simulation (EMAPS) technology that replaces time stepping with self-adaptive update events. Local calculations are carried out only on an “as needed basis”. EMAPS (i) guarantees accurate and stable processing of physical variables in time accurate simulations, and (ii) eliminates unnecessary computation. Applying EMAPS to the hybrid model has resulted in the development of a unique parallel code, dimension-independent (compile-time-configurable) HYPERS (Hybrid Parallel Event-Resolved Simulator) that scales to hundreds of thousands of parallel processors. HYPERS advances electromagnetic fields and particles asynchronously on time scales determined by local physical laws and mesh properties. To achieve high computational accuracy in complex device geometries, HYPERS employs high-fidelity Cartesian grids with masked conductive cells. The HYPERS model includes multiple ion species, energy and momentum conserving ion-ion collisions, and provides a number of approximations for plasma resistivity and vacuum regions. Both local and periodic boundary conditions are allowed. The HYPERS solver preserves zero divergence of magnetic field. The project has demonstrated HYPERS capabilities on a number of applications of interest to fusion and astrophysical plasma physics applications listed below. 1. Theta-pinch formation of FRCs The formation, spontaneous spin-up, and stability of theta-pinch formed field-reversed configurations have been studied self-consistently in 3D. The end-to-end hybrid simulations reveal poloidal profiles of implosion-driven fast toroidal plasma rotation and demonstrate three discharge regimes as a function of experimental parameters: the decaying stable configuration, the tilt unstable configuration, and the nonlinear evolution of a fast growing tearing mode. 2. FRC collisions with magnetic mirrors Interactions of fast plasma streams and objects with magnetic obstacles (dipoles, mirrors, etc) lie at the core of many space and laboratory plasma phenomena ranging from magnetoshells and solar wind interactions with planetary magnetospheres to compact fusion plasmas. HYPERS simulations are compared with data from the MSX experiment (LANL) that focuses on the physics of magnetized collisionless shocks through the acceleration and subsequent stagnation of FRC plasmoids against a strong magnetic mirrors and flux-conserving boundaries. 3. Exploding magnetoplasmas Results from hybrid simulations of two experiments at the LAPD and Nevada Terawatt Facility are discussed where short-pulse lasers are used to ablate solid targets to produce plasmas that expand across external magnetic fields. The first simulation recreates flutelike density striations observed at the leading edge of a carbon plasma and predicts an early destruction of the magnetic cavity in agreement with experimental evidence. In the second simulation a polyethylene target is ablated into a mixture of protons and carbon ions. A mechanism is demonstrated that allows protons to penetrate the magnetic field in the form of a collimated flow. The results are compared to experimental data and single-fluid MHD simulations. The EMAPS framework has the potential for wide application in many other engineering and scientific fields, such as climate models, biological systems, electronic devices, seismic events, oil reservation simulators that all involve advancing solutions of partial differential equations in time where the rate of activity can be adapted widely over the spatial domain depending on locally space/time phenomena (“events”).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Quantitative Characterization of Impacts of Coupled Geomechanics and Flow on Safe and Permanent Geological Storage of CO 2 in Fractured Aquifers

This is the final scientific report for the award DE-FE0023305, entitled “Quantitative Characterization of Impacts of Coupled Geomechanics and Flow on Safe and Permanent Geological Storage of CO 2 in Fractured Aquifers.” The work has been divided into six tasks. In Task 2, we characterized rock properties, which are important when developing a quantitative approach for understanding and predicting geomechanical effects on large-scale CO 2 injection and long-term storage in the subsurface. Rock properties of interest for this characterization include porosity, permeability, elastic constants, strength, and heat capacity. We measured rock properties for three different rock types: concrete, sandstone and shale. These properties were obtained from the acoustic measurement, permeability and porosity measurement, Brazilian test, the uniaxial compressive test, and heat capacity measurement. Then, using the Brazilian technique to fracture samples, we measured their permeability under brine injection for varying effective stresses, followed by similar measurements of permeability under sc-CO 2 injection. Permeability to brine and sc-CO 2 decreased as the effective stress increased. The apparent permeability to sc-CO 2 was an order of magnitude less than that for brine, a consequence of CO 2 being the non-wetting fluid. In Task 3, we developed understandings and correlations for CO 2 injection pressure induced fracturing. We initially fractured four eight inch cubes of concrete with brine to establish a standard test procedure. Then, additional concrete samples were fractured using brine. Confining stresses were exerted on the samples and brine was injected at a constant rate into a borehole that was previously drilled into the concrete cube. The borehole pressure was measured and fracture initiation was identified as a peak in the pressure profile. That was followed by fracture propagation, which was identified as a plateau in the pressure profile, and finally the fracture reached the sample surface. The sample surfaces were photographed before and after fracturing, which enabled determination of where the fracture reached the sample surface. Afterwards, aqueous food dye solution was injected into the fractured concrete samples to color the fracture planes. The borehole was filled with food dye solution in advance, a gas pulse was transmitted to the borehole, and the borehole was then shut in to allow the food dye to be driven into existing fractures. Then, the concrete samples were broken down by high-pressure gas injection to reveal the geometry and morphology of the hydraulic fractures. In addition, acoustic measurements were conducted on multiple locations on each face of the concrete samples before and after injection for comparison. This comparison was used to confirm changes in the internal structure of the concrete sample, namely the formation of fractures. There were six concrete samples fractured with brine. Afterwards, concrete samples were fracture using sc-CO 2 . The main change in experimental procedure was the inclusion of temperature control, since in the field, we mostly encounter temperature and pressure conditions above the supercritical point of CO 2 . The concrete samples were pre-heated in an oven to elevate their temperature, and thermal tape was used to heat the CO 2 to a temperature above its supercritical temperature. Two concrete samples were fractured using liquid CO 2 , four concrete samples were fractured using gaseous CO 2 , and twenty two concrete samples were fractured using sc-CO 2 . For the concrete samples fractured using sc-CO 2 , three had pre-existing hydraulic fractures on the surface near the wellbore, and two were water saturated (as opposed to containing only air). Two concrete samples were composite samples, with a high permeability ball in the center, in order to represent CO 2 injection into a high permeability zone surrounded by a low permeability sealing formation. Later, we fractured five shale samples, obtained from the Niobrara shale outcrop, with sc-CO 2 . The injection rate was constant, confining stresses were applied, the shale samples were pre-heated in an oven to the desired temperature, and the sample was then fractured. An analysis of the failure type was done for many of the concrete samples, based on the predicted break down pressure for tensile and shear failure, and it was found that most samples fractured due to shear failure. Also, the presence of induced fractures has little to no effect on the fracture initiation pressure, but significantly changes the fracture morphology. The CO 2 fracturing in water saturated samples behaved much similar to brine induced fracturing. There was a significant difference of breakdown pressure for injection of CO 2 and brine, where the breakdown pressure for CO 2 induced fracturing is generally around the minimum horizontal stress and that for brine is much higher. Finally, the fracture orientation for most samples with large confining stress differential was perpendicular to the minimum horizontal stress. When the stress differential, especially the difference between the two directions with the smallest stresses was small, the fracture orientation was determined not only by the minimum stress direction, but also the magnitude of the confining stress, the tensile strength of the rock, and the breakdown pressure. In Task 4, we modified our coupled flow-geomechanical models to model fracture growth and propagation in storage formations and caprocks. These flow-geomechanical models are TOUGH2-CSM and TOUGH2-FLAC. The TOUGH2-CSM fluid and heat flow formulation is based on the TOUGH2 formulation of mass and energy conservation equations that govern fluid and heat flow in general multiphase, multicomponent, multi-porosity systems. The TOUGH2-CSM simplified geomechanical formulation is based on the linear theory of elasticity applied to multi-porosity non-isothermal (thermo-multi-poroelastic) media. We previously derived, from the fundamental relations of the linear theory of elasticity, an equation relating mean stress, pore pressures, temperatures, and body force (the Mean Stress equation) that we added to the TOUGH2-CSM fluid and heat flow equations. We extended the TOUGH2-CSM simplified geomechanical formulation to calculate the entire stress tensor by deriving equations for stress tensor components from derivatives of the Cartesian thermo-multi-poroelastic Navier equation components. In addition, we derived equations for stress tensor components in rz-coordinates. We solved these geomechanical equations using the integral finite difference method. This method utilizes momentum fluxes obtained by “factoring” a divergence operator out of each geomechanical equation. In Task 5, we did a literature survey to determine suitable correlations for fracture initiation, growth, and propagation. After surveying numerous papers, we modified TOUGH2-CSM to model stress dependent fracture initiation and growth using the Mohr-Coulomb criterion for shear failure in faults and caprock and the condition of negative minimum effective stress for tensile failure. Tensile failure results in a fracture whose width depends on the difference between fracture pressure and minimum stress and whose propagation is based on the difference between the stress intensity factor at the fracture tip and the rock toughness. The TOUGH-FLAC model is based on the FLAC3D ubiquitous joint model, extended for the modeling of permeability changes induced by changes in effective normal stress across fractures as well permeability changes caused by shearing along fractures. In Task 6, we validated our coupled flow-geomechanical models using analytical solutions and problems from the literature. We validated the TOUGH2-CSM modifications using analytical solutions to the displacement from a uniform load on a semi-infinite elastic medium and the two-dimensional Mandel-Cryer effect. Those were followed by a sample problem to test our stress tensor calculations, which was for injection into a single-phase reservoir with constant properties. Comparisons of our simulator to published results were done on for the depletion of a single-phase reservoir with stress dependent porosity, two examples from the In Salah gas project, an axisymmetric baserock–reservoir–caprock system with a normal faulting stress regime that had thermally induced shear stresses, simulations of CO 2 leakage through caprock in a two-dimensional reservoir in Cartesian coordinates, and predictions of caprock failure from CO 2 injection into an axisymmetric reservoir. Finally, we simulated CO 2 injection pressure-induced fracturing from one of our laboratory studies, Sample 39, and obtained a fracture profile. We also validated the TOUGH2-FLAC model for fluid driven fracture growth against solutions based on the KGD model and a case with an inclined fracture that was loaded from the boundaries to achieve wing-crack propagation. The strain softening tensile behavior and softening of modulus considering a simple damage approach was verified by a simulation representing hydraulic fracturing stress measurement around a vertical well. Finally, the model was validated against deep fracture zone opening and surface uplift at In Salah with inverses analysis used to improve the match of simulation with field data. In Task 7, we developed a scheme based on inverse modeling that can be used to identify caprock leakage. Injection of fluid into a reservoir results in a time varying pressure profile that depends on the properties of the rock and the fluid. The presence of significant caprock leakage would effect this pressure profile and we used inverse modeling, namely the Levenberg-Marquardt method, to determine the leakage location from the pressure profile. A two-dimensional and a three-dimensional example based on published data were presented. The published simulation was run and its results were considered as “field” data. In the two-dimensional example, a simulation of CO 2 injection into a two-dimensional aquifer-caprock system, the caprock had a vertical fault through which leakage could occur. The location of this leakage was obtained from our inverse modeling scheme. In the three-dimensional example, the simulation domain had four geological layers, one of which was the injection zone that contained a horizontal injection well. We simulated fractures as high permeability gridblocks, introduced a fracture that spanned the caprock depth, and ran a simulation to obtain “field” data. Then, using our inverse modeling scheme, we obtained the location of that fracture.

58 GEOSCIENCES↗
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 OSTI2016

Parameterization retrieval of trace gas volume mixing ratios from Airborne MAX-DOAS

We present a parameterization retrieval of volume mixing ratios (VMRs) from differential slant column density (dSCD) measurements by Airborne Multi-AXis Differential Optical Absorption Spectroscopy (AMAX-DOAS). The method makes use of the fact that horizontally recorded limb spectra (elevation angle 0°) are strongly sensitive to the atmospheric layer at instrument altitude. These limb spectra are analyzed using reference spectra that largely cancel out column contributions from above and below the instrument, so that the resulting limb dSCDs, i.e., the column integrated concentration with respect to a reference spectrum, are almost exclusively sensitive to the atmospheric layers around instrument altitude. The conversion of limb dSCDs into VMRs is then realized by calculating box air mass factors (Box-AMFs) for a Rayleigh atmosphere and applying a scaling factor constrained by O 4 dSCDs to account for aerosol extinction. An iterative VMR retrieval scheme corrects for trace gas profile shape effects. Benefits of this method are (1) a fast conversion that only requires the computation of Box-AMFs in a Rayleigh atmosphere; (2) neither local aerosol extinction nor the slant column density in the DOAS reference (SCD ref ) needs to be known; and (3) VMRs can be retrieved for every measurement point along a flight track, thus increasing statistics and adding flexibility to capture concentration gradients. Sensitivity studies are performed for bromine monoxide (BrO), iodine monoxide (IO) and nitrogen dioxide (NO 2 ), using (1) simulated dSCD data for different trace gas and aerosol profiles and (2) field measurements from the Tropical Ocean tRoposphere Exchange of Reactive halogen species and Oxygenated VOC (TORERO) field experiment. For simulated data in a Rayleigh atmosphere, the agreement between the VMR from the parameterization method (VMR para ) and the true VMR (VMR true ) is excellent for all trace gases. Offsets, slopes and R 2 values for the linear fit of VMR para over VMR true are, respectively (0.008 ± 0.001) pptv, 0.988 ± 0.001, 0.987 for BrO; (-0.0066 ± 0.0001) pptv, 1.0021 ± 0.0003, 0.9979 for IO; (-0.17 ± 0.03) pptv, 1.0036 ± 0.0001, 0.9997 for NO 2 . The agreement for atmospheres with aerosol shows comparable R 2 values to the Rayleigh case, but slopes deviate a bit more from one: (0.093 ± 0.002) pptv, 0.933 ± 0.002, 0.907 for BrO; (0.0021 ± 0.0004) pptv, 0.887 ± 0.001, 0.973 for IO; (8.5 ± 0.1) pptv, 0.8302 ± 0.0006, 0.9923 for NO 2 . VMR para from field data are further compared with optimal estimation retrievals (VMR OE ). Least orthogonal distance fit of the data give the following equations: BrO para = (0.1 ± 0.2) pptv + (0.95 ± 0.14) × BrO OE ; IO para = (0.01 ± 0.02) pptv + (1.00 ± 0.12) × IO OE ; NO 2para = (3.9 ± 2.5) pptv + (0.87 ± 0.15) × NO 2OE . Overall, we conclude that the parameterization retrieval is accurate with an uncertainty of 20 % for IO, 30 % for BrO and NO 2 , but not better than 0.05 pptv IO, 0.5 pptv BrO and 10 pptv NO 2 . Finally, the retrieval is applicable over a wide range of atmospheric conditions and measurement geometries and not limited to the interpretation of vertical profile measurements in the remote troposphere.

54 ENVIRONMENTAL SCIENCES↗
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↗
DOE OSTI2011

Statistical variability and confidence intervals for planar dose QA pass rates

Purpose: The most common metric for comparing measured to calculated dose, such as for pretreatment quality assurance of intensity-modulated photon fields, is a pass rate (%) generated using percent difference (%Diff), distance-to-agreement (DTA), or some combination of the two (e.g., gamma evaluation). For many dosimeters, the grid of analyzed points corresponds to an array with a low areal density of point detectors. In these cases, the pass rates for any given comparison criteria are not absolute but exhibit statistical variability that is a function, in part, on the detector sampling geometry. In this work, the authors analyze the statistics of various methods commonly used to calculate pass rates and propose methods for establishing confidence intervals for pass rates obtained with low-density arrays. Methods: Dose planes were acquired for 25 prostate and 79 head and neck intensity-modulated fields via diode array and electronic portal imaging device (EPID), and matching calculated dose planes were created via a commercial treatment planning system. Pass rates for each dose plane pair (both centered to the beam central axis) were calculated with several common comparison methods: %Diff/DTA composite analysis and gamma evaluation, using absolute dose comparison with both local and global normalization. Specialized software was designed to selectively sample the measured EPID response (very high data density) down to discrete points to simulate low-density measurements. The software was used to realign the simulated detector grid at many simulated positions with respect to the beam central axis, thereby altering the low-density sampled grid. Simulations were repeated with 100 positional iterations using a 1 detector/cm{sup 2} uniform grid, a 2 detector/cm{sup 2} uniform grid, and similar random detector grids. For each simulation, %/DTA composite pass rates were calculated with various %Diff/DTA criteria and for both local and global %Diff normalization techniques. Results: For the prostate and head/neck cases studied, the pass rates obtained with gamma analysis of high density dose planes were 2%-5% higher than respective %/DTA composite analysis on average (ranging as high as 11%), depending on tolerances and normalization. Meanwhile, the pass rates obtained via local normalization were 2%-12% lower than with global maximum normalization on average (ranging as high as 27%), depending on tolerances and calculation method. Repositioning of simulated low-density sampled grids leads to a distribution of possible pass rates for each measured/calculated dose plane pair. These distributions can be predicted using a binomial distribution in order to establish confidence intervals that depend largely on the sampling density and the observed pass rate (i.e., the degree of difference between measured and calculated dose). These results can be extended to apply to 3D arrays of detectors, as well. Conclusions: Dose plane QA analysis can be greatly affected by choice of calculation metric and user-defined parameters, and so all pass rates should be reported with a complete description of calculation method. Pass rates for low-density arrays are subject to statistical uncertainty (vs. the high-density pass rate), but these sampling errors can be modeled using statistical confidence intervals derived from the sampled pass rate and detector density. Thus, pass rates for low-density array measurements should be accompanied by a confidence interval indicating the uncertainty of each pass rate.

61 RADIATION PROTECTION AND DOSIMETRY↗