SEARCH · Engineering Papers
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.
Explore the source record for details and available documents.
A procedure is described for the analysis of the stable carbon isotopic composition of dissolved organic carbon (DOC) in natural waters from marine and higher-salinity environments. Rapid (less than 5 min) and complete oxidation of DOC is achieved using a modification of previous photochemical oxidation techniques. The CO2 evolved from DOC oxidation can be collected in less than 10 min for isotopic analysis. The procedure is at present suitable for oxidation and collection of 1-5 micromoles of carbon and has an associated blank of 0.1-0.2 micromole of carbon. Complete photochemical oxidation of DOC standards was demonstrated by quantitative recovery of CO2 as measured manometrically. Isotopic analyses of standards by photochemical and high-temperature sealed-tube combustion methods agreed to within 0.3%. Photochemical oxidation of DOC in a representative sediment pore-water sample was also quantitative, as shown by the excellent agreement between the photochemical and sealed-tube methods. The delta 13C values obtained for pore-water DOC using the two methods of oxidation were identical, suggesting that the modified photochemical method is adequate for the isotopically non-fractionated oxidation of pore-water DOC. The procedure was evaluated through an analysis of DOC in pond and pore waters from a hypersaline microbial mat environment. Concentrations of DOC in the water column over the mat displayed a diel pattern, but the isotopic composition of this DOC remained relatively constant (average delta 13C = -12.4%). Pore-water DOC exhibited a distinct concentration maximum in the mat surface layer, and delta 13C of pore-water DOC was nearly 8% lighter at 1.5-2.0-cm depth than in the mat surface layer (0-0.5-cm depth). These results demonstrate the effectiveness of the method in elucidating differences in DOC concentration and delta 13C over biogeochemically relevant spatial and temporal scales. Carbon isotopic analysis of DOC in natural waters, especially pore waters, should be a useful probe of biogeochemical processes in recent environments.
A discrete-vortex method is developed to investigate the wing rock problem associated with highly swept wings. The method uses two logarithmic vortices placed above the wing to represent the vortex flow field and uses boundary conditions based on conical flow, vortex rate of change of momentum, and other considerations to position the vortices and determine their strengths. A relationship based on the time analogy and conical-flow assumptions is used to determine the hysteretic positions of the vortices during roll oscillations. Static and dynamic vortex positions and wing rock amplitudes and frequencies calculated by using the method are generally in good agreement with available experimental data. The results verify that wing rock is caused by hysteretic deflections of the vortices and indicate that the stabilizing moments that limit wing rock amplitudes are the result of the one primary vortex moving outboard of the wing where it has little influence on the wing.
The objective is to characterize the fracture behavior and to define the fracture mechanisms for new Al-Li-Cu alloys, with emphasis on the role of indium additions and cryogenic temperatures. Three alloys were investigated in rolled product form: 2090 baseline and 2090 + indium produced by Reynolds Metals, and commercial AA 2090-T81 produced by Alcoa. The experimental 2090 + In alloy exhibited increases in hardness and ultimate strength, but no change in tensile yield strength, compared to the baseline 2090 composition in the unstretched T6 condition. The reason for this behavior is not understood. Based on hardness and preliminary Kahn Tear fracture experiments, a nominally peak-aged condition was employed for detailed fracture studies. Crack initiation and growth fracture toughness were examined as a function of stress state and microstructure using J(delta a) methods applied to precracked compact tension specimens in the LT orientation. To date, J(delta a) experiments have been limited to 23 C. Alcoa 2090-T81 exhibited the highest toughness regardless of stress state. Fracture was accompanied by extensive delamination associated with high angle grain boundaries normal to the fatigue precrack surface and progressed microscopically by a transgranular shear mechanism. In contrast the two peak-aged Reynolds alloys had lower toughness and fracture was intersubgranular without substantial delamination. The influences of cryogenic temperature, microstructure, boundary precipitate structure, and deformation mode in governing the competing fracture mechanisms will be determined in future experiments. Results contribute to the development of predictive micromechanical models for fracture modes in Al-Li alloys, and to fracture resistant materials.
Reliable predictions of climate change impacts on water use, irrigation requirements and yields of irrigated sugarcane in South Africa (a water-scarce country) are necessary to plan adaptation strategies. Although previous work has been done in this regard, methodologies and results vary considerably. The objectives were (1) to estimate likely impacts of climate change on sugarcane yields, water use and irrigation demand at three irrigated sugarcane production sites in South Africa (Malelane, Pongola and La Mercy) for current (1980-2010) and future (2070-2100) climate scenarios, using an approach based on the Agricultural Model Inter-comparison and Improvement Project (AgMIP) protocols; and (2) to assess the suitability of this methodology for investigating climate change impacts on sugarcane production. Future climate datasets were generated using the Delta downscaling method and three Global Circulation Models (GCMs) assuming atmospheric CO2 concentration [CO2] of 734 ppm(A2 emissions scenario). Yield and water use were simulated using the DSSAT-Canegro v4.5 model. Irrigated cane yields are expected to increase at all three sites (between 11 and 14%), primarily due to increased interception of radiation as a result of accelerated canopy development. Evapotranspiration and irrigation requirements increased by 11% due to increased canopy cover and evaporative demand. Sucrose yields are expected to decline because of increased consumption of photo-assimilate for structural growth and maintenance respiration. Crop responses in canopy development and yield formation differed markedly between the crop cycles investigated. Possible agronomic implications of these results include reduced weed control costs due to shortened periods of partial canopy, a need for improved efficiency of irrigation to counter increased demands, and adjustments to ripening and harvest practices to counter decreased cane quality and optimize productivity. Although the Delta climate data downscaling method is considered robust, accurate and easily-understood, it does not change the future number of rain-days per month. The impacts of this and other climate data simplifications ought to be explored in future work. Shortcomings of the DSSAT-Canegro model include the simulated responses of phenological development, photosynthesis and respiration processes to high temperatures, and the disconnect between simulated biomass accumulation and expansive growth. Proposed methodology refinements should improve the reliability of predicted climate change impacts on sugarcane yield.
A system and method for optimizing a fiber optic sensor by properly clocking or rotationally orienting a window with a fiber optic focuser. This method may include shining light through the focuser into the window and receiving with a reflectometer a first reflection from a first surface of the window and a second reflection from a second surface of the window. Next, the method may include calculating a delta of intensities between the first and second reflections, and then rotating an orientation of the window to a next rotational orientation. The steps of rotating the orientation of the window, receiving reflections, and calculating their delta may be repeated any plurality of times to solve for a plurality of deltas. Then, the method may include the steps of determining which one of the deltas is the largest and fixing the window to the focuser at a rotational orientation associated therewith.
Modern advanced Phasor Measurement Units (PMUs) are developed with ultra-high reporting rates to meet the demand for monitoring the power systems dynamics in detail. Due to the large volume of data, the communication and storage systems are seriously challenged with the presence of Ultra-High-Density (UHD) synchrophasor and Point on Wave (POW) data. Therefore, it is an urgent task to compress the UHD data for more efficient communication and data storage. This paper proposes several methods to compress the synchrophasor and POW data in a lossless manner. First, an Improved-Time-Series-Special Compression (ITSSC) method is proposed to compress the UHD frequency data. Second, a Delta-difference Huffman method is combined with the TSSC algorithm to compress the UHD phase angle data. Finally, a cyclical high-order delta modulation method is proposed to compress the UHD POW data. The proposed models are extensively tested and compared with different existing lossless compression algorithms using the field-collected synchrophasor and POW data at different reporting rates. The results indicate that the proposed algorithms are efficient in performing lossless compression for the UHD synchrophasor and POW data in real time.
The coming massive parallelism of exascale computing presents a pressing challenge for the many DOE simulations of time-dependent partial differential equations (PDEs), which typically use traditional sequential time stepping methods. Since this traditional approach is inherently serial, it presents a sequential bottleneck when moving to exascale computing, because future performance gains will come through greater concurrency, not faster clock speeds. Thus, the goal of this work is to research parallelism in time, i.e., methods that compute multiple time values simultaneously, not sequentially. The focus will be on hyperbolic and chaotic problems of interest to DOE, with the goal of enabling scalable simulations of time-dependent hyperbolic and chaotic problems on future architectures. The chosen methodology for solving these problems parallel-in-time is multigrid, because multigrid (when it works) is a powerful, optimal, and scalable solver for discretized PDEs. Multigrid is already commonly used in many DOE simulations for scalably and optimally solving space-only PDE problems. The areas of hyperbolic and chaotic problems are chosen because of their relevance to problems of programmatic interest to DOE. However, these problems are also well-known to be difficult for parallelin-time methods, with the most common method, parareal, diverging in many cases. The current state-of-the-art for parallel-in-time at LLNL is the multigrid reduction in time (MGRIT) XBraid package, which also struggles for such problems, while still showing some improvement over parareal. In summary, new methods are needed for an efficient parallel-in-time scheme for hyperbolic and chaotic problems, and this work shall research promising new multigrid methods in this area. In particular, this work shall continue researching the directions from the current collaboration with Dr. Falgout, which are laid out in the work Toward Parallel in Time for Chaotic Dynamical Systems and showed the first known results of a parallel-in-time speedup for a chaotic problem. This work outlines two key improvements to XBraid for chaotic problems, the so-called “theta” and “delta-correction” methods. Here, these two improvements will be implemented in a high-performance but general way in XBraid and explored for more complicated problems. We will additionally research, as time allows, improvements to these techniques, as well as multigrid relaxation techniques based on Least Squares Shadowing (LSS by Wang) and a nonintrusive block tridiagonal solver based on MGRIT, called TriMGRIT.
The coming massive parallelism of exascale computing presents a pressing challenge for the many DOE simulations of time-dependent partial differential equations (PDEs), which typically use traditional sequential time stepping methods. Since this traditional approach is inherently serial, it presents a sequential bottleneck when moving to exascale computing, because future performance gains will come through greater concurrency, not faster clock speeds. Thus, the goal of this work is to research parallelism in time, i.e., methods that compute multiple time values simultaneously, not sequentially. The focus will be on hyperbolic and chaotic problems of interest to DOE, with the goal of enabling scalable simulations of time-dependent hyperbolic and chaotic problems on future architectures. The chosen methodology for solving these problems parallel-in-time is multigrid, because multigrid (when it works) is a powerful, optimal, and scalable solver for discretized PDEs. Multigrid is already commonly used in many DOE simulations for scalably and optimally solving space-only PDE problems. The areas of hyperbolic and chaotic problems are chosen because of their relevance to problems of programmatic interest to DOE. However, these problems are also well-known to be difficult for parallel-in-time methods, with the most common method, parareal, diverging in many cases. The current state of-the-art for parallel-in-time at LLNL is the multigrid reduction in time (MGRIT) XBraid package, which also struggles for such problems, while still showing some improvement over parareal. In summary, new methods are needed for an efficient parallel-in-time scheme for hyperbolic and chaotic problems, and this work shall research promising new multigrid methods in this area. In particular, this work shall continue researching the directions from the current collaboration with Dr. Falgout, which are laid out in the work Toward Parallel in Time for Chaotic Dynamical Systems and showed the first known results of a parallel-in-time speedup for a chaotic problem. This work outlines two key improvements to XBraid for chaotic problems, the so-called “theta” and “delta-correction” methods. Here, these two improvements will be further researched and improved (including with a new relaxation method inspired by on Least Squares Shadowing (LSS)) and explored for more complicated problems.
Counts of a specific characteristic were obtained within regions defined on an object that was manufactured in a proprietary setting. The count regions were altered during production and resulted in misaligned or overlapping count data. A closed-formula maximum likelihood estimator (MLE) of the new region means is derived using all of the available count data and an independent Poisson model. The MLE is shown to be preferable to estimators constructed using generalized linear models for the overlapping data setting. This closed-form estimator extends to over-dispersed overlapping count data as the quasi-MLE and also performs well with correlated overlapping count data. Standard errors for the estimator are approximated and are validated with a simulation study. Additionally, the methods are extended to overlapping multinomial data. Illustrative examples of the methods are provided throughout the paper and are reproducible with the supplemental R code. Additionally, proofs of the paper’s results are also included in the supplemental material.
A nonsteady vortex-lattice method is introduced for predicting the dynamic stability derivatives of a delta wing undergoing an oscillatory motion. The analysis is applied to several types of small oscillations in pitch. The angle of attack varied between + or - 1 deg, with the mean held at 0 deg when the flow was assumed to be attached and between + or - 1 deg and the mean held at 15 deg when both leading-edge separation and wake roll-up were included. The computed results for damping in pitch are compared with several other methods and with experiments, and are found to be consistent and in good agreement.
Fatigue crack growth rate (FCGR) tests were conducted on IN 718 in the solution annealed and aged condition at room temperature in accordance with E647-87. As part of each test, the FCGR threshold was measured using the decreasing Delta K method. A new heat of material was being tested and some of this material was sent to a different laboratory which wanted to use a specimen with a 127 mm width. Threshold data previously had been established on specimens with a width of 50.8 mm. As a check of the laboratory, tests were conducted at room temperature and R equal to 0.1 for comparison with the earlier data. The results were a threshold significantly higher than previously observed. Interchanging of specimen sizes and laboratories showed that the results were not due to a heat-to-heat or lab-to-lab variation. The results to be presented here are those obtained at the original laboratory. Growth rates were measured using the electric potential drop technique at R values of 0.1, 0.7, and 0.9. Compact tension specimen sizes with planer dimensions of 25.4 mm, 50.8 mm, and 127 mm were used. Crack growth rates at threshold were generally below 2.5 X 10(exp -8) mm / cycle. Closure measurements were made on some of the specimens by a manual procedure using a clip gage. When the crack growth rate data for the specimens tested at R equal to 0.1 were plotted as a function of applied Delta K, the thresholds varied with specimen width. The larger the width, the higher the threshold. The thresholds varied from 6.5 MPa-m(exp 1/2) for the 25.4 mm specimen to 15.4 MPa-m(exp 1/2) for the 127 mm specimen. At R equal to 0.7, the 25.4 mm and 50.8 mm specimens had essentially the same threshold, about 2.9 MPa-m(exp 1/2)while the 127 mm specimen had a threshold of 4.5 MPa-m(exp 1/2). When plotted as a function of effective Delta K, the R equal to 0.1 data are essentially normalized. Various aspects of the test procedure will be discussed as well as the results of analysis of the data using some different closure models.
The LINAC at the Los Alamos Neutron Science Center (LANSCE) has been utilizing the Delta-t method to match the RF cavities to the design acceleration parameters since its commissioning in 1972. The differences in time-of-flight between two subsequent Beam Position and Phase Monitors (BPPMs) are measured with both accelerated and drifting beams, depending on the whether the module is set to on or off. The algorithm optimizes the module amplitude and phase via iterative measurements if the initial phase is in the vicinity of the design value. With an upgrade to a faster readout system, a scan over the whole RF cavity phase range requires relatively less time than the classical optimization procedure. The Phase Scan Signature Matching (PSSM) method provides a time-efficient method that ensures the phase selection lands on the bunching side and empowers future analyses to build module-specific models. The PSSM also utilizes a direct model to determine the correct amplitude to sub-percent level instead of using linearized matrices. Furthermore, lacking a reliable energy measurement method in the LINAC, we measure the beam phases at two downstream locations to increase the precision of energy measurements. In this letter, we also discuss the sensitivities of PSSM, error propagation, and the implementation results for the 2019 and 2020 beam cycles.
Let $\mathcal{H}$ be an $r$-uniform hypergraph. The minimum positive co-degree of $\mathcal{H}$, denoted by $\delta_{r-1}^+(\mathcal{H})$, is the minimum $k$ such that if $S$ is an $(r-1)$-set contained in a hyperedge of $\mathcal{H}$, then $S$ is contained in at least $k$ hyperedges of $\mathcal{H}$. For $r\geq k$ fixed and $n$ sufficiently large, we determine the maximum possible size of an intersecting $r$-uniform $n$-vertex hypergraph with minimum positive co-degree $\delta_{r-1}^+(\mathcal{H}) \geq k$ and characterize the unique hypergraph attaining this maximum. This generalizes the Erd\Hos--Ko--Rado theorem which corresponds to the case $k=1$. Furthermore, our proof is based on the delta-system method.
Radiographic, penetrant, ultrasonic shear wave and ultrasonic delta inspection methods for critical defect determination
As part of a continuing program of numerical tests of convergence accelerators, the iterated Aitken's Delta-squared method, Wynn's epsilon algorithm, Brezinski's theta algorithm, and Levin's u transform are compared on a broad range of test problems: linearly convergence alternating, monotone, and irregular-sign series, logarithmically convergent series, power method and Bernoulli method sequences, alternating and monotone asymptotic series, and some perturbation series arising in applications. In each category either the epsilon algorithm or the u transform gives the best results of the four methods tested. In some cases differences among methods are slight, and in others they are quite striking.
This paper first examines the extension of two-frequency, or Delta-k near-nadir remote sensing techniques to off-nadir radar altimetry. A different approach to sea state sensing is investigated which appears to offer much higher accuracy. The second section examines the Delta-k method of sensing rough surface area correlation or ocean wavenumber spectrum. The technique shows promise of directly sensing this parameter, in contrast with synthetic-aperture radar methods which are based on complex wave interaction mechanisms and mathematical transformations requiring the acquisition of voluminous data.
Evolutionary sequences for low-metallicity stars (Z ranging from 0.001 to 0.0001) to study the effects of internal stellar rotation on the evolutionary time scales in the pre-main sequence, the main sequence (MS), and around the MS turnoff. Although a substantial amount of angular momentum remains in the interior, rotation is only a minor perturbation on the structure and ages of globular cluster stars. Even models with large initial angular momenta have MS lifetimes that are within 1 percent of those of standard models of the same mass and composition. Therefore, rotation does not affect age estimates of globular clusters from isochrone fitting. Furthermore, the models suggest that because rotation is not likely to affect horizontal-branch (HB) morphology, it does not affect significantly age estimates from the Delta-V method. Nevertheless, the internal angular momentum in the models is consistent with observations of surface rotational velocities on the HB, which require the preservation of a large reservoir of internal angular momentum.