SEARCH · Engineering Papers

Results for “factorization methods”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 1,171 records · Page 65

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Improvements to a class of hybrid methods for radiation transport: Nyström reconstruction and defect correction methods

In this study, two modifications are introduced for improving the accuracy, versatility, and robustness of a class of hybrid methods for radiation transport. In general, such methods are constructed by splitting the radiative flux into collided and uncollided components to which low- and high-resolution angular approximations are applied, respectively. In this work we focus on discrete ordinates discretizations of high and low order. The first modification we introduce changes the way in which the collided component is mapped into the uncollided component at the end of each time step in a simulation. The new mapping is a Nyström-type reconstruction that is applicable to arbitrary discrete ordinates quadratures, is guaranteed to preserve positivity of the solution provided that all ordinate weights are positive, is significantly more accurate than previous methods, and can be readily extended to other discretizations such as moment methods, finite element methods, and diffusion approximations. The second modification leverages integral deferred correction (IDC) to iteratively correct for the splitting error introduced by the inconsistency in angular discretization between the collided and uncollided components, in addition to improving the accuracy of the low-order temporal error that is treated by traditional IDC methods. Numerical tests in one- and two-dimensional geometries are used to demonstrate the increased accuracy and efficiency of the proposed modifications. It is found that the two techniques combined yield methods with solution accuracy and memory requirements comparable to that of monolithic discrete ordinates methods while reducing runtime by as much as a factor of between two and ten, depending on the problem.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS1978

Strainrange partitioning behavior of the nickel-base superalloys, Rene 80 and 100

A study has been made to assess the ability of the method of Strainrange Partitioning (SRP) to both correlate and predict high-temperature low-cycle fatigue lives of nickel-base superalloys for gas turbine applications. SRP is shown to correlate the cyclic lives of the baseline tests to within factors of nearly two. The partitioned strainrage versus life relationships for uncoated Rene 80 and cast IN 100 have also been determined from the ductility normalized-strainrange partitioning equations. These were used to predict the cyclic lives of the baseline tests. Predicted and observed cyclic lives agreed to within factors of nearly three. The life predictability of the method is also verified for cast IN 100 by applying the baseline results to the cyclic-life prediction of a series of complex strain-cycling tests with multiple hold periods at constant strain. It is concluded that the method of SRP can correlate and predict the cyclic lives of laboratory specimens of the nickel-base superalloys evaluated in this program.

Halford, G. R.↗

NASA NTRS2014

The AgMIP Coordinated Climate-Crop Modeling Project (C3MP): Methods and Protocols

Climate change is expected to alter a multitude of factors important to agricultural systems, including pests, diseases, weeds, extreme climate events, water resources, soil degradation, and socio-economic pressures. Changes to carbon dioxide concentration ([CO2]), temperature, and water (CTW) will be the primary drivers of change in crop growth and agricultural systems. Therefore, establishing the CTW-change sensitivity of crop yields is an urgent research need and warrants diverse methods of investigation. Crop models provide a biophysical, process-based tool to investigate crop responses across varying environmental conditions and farm management techniques, and have been applied in climate impact assessment by using a variety of methods (White et al., 2011, and references therein). However, there is a significant amount of divergence between various crop models' responses to CTW changes (Rotter et al., 2011). While the application of a site-based crop model is relatively simple, the coordination of such agricultural impact assessments on larger scales requires consistent and timely contributions from a large number of crop modelers, each time a new global climate model (GCM) scenario or downscaling technique is created. A coordinated, global effort to rapidly examine CTW sensitivity across multiple crops, crop models, and sites is needed to aid model development and enhance the assessment of climate impacts (Deser et al., 2012). To fulfill this need, the Coordinated Climate-Crop Modeling Project (C3MP) (Ruane et al., 2014) was initiated within the Agricultural Model Intercomparison and Improvement Project (AgMIP; Rosenzweig et al., 2013). The submitted results from C3MP Phase 1 (February 15, 2013-December 31, 2013) are currently being analyzed. This chapter serves to present and update the C3MP protocols, discuss the initial participation and general findings, comment on needed adjustments, and describe continued and future development. AgMIP aims to improve substantially the climate, crop, and economic simulation tools that are used to characterize the agricultural sector, to assess future world food security under changing climate conditions, and to enhance adaptation capacity both globally and regionally. To understand better and improve the modeled crop responses, AgMIP has conducted detailed crop model intercomparisons at closely observed field sites for wheat (Asseng et al., 2013), rice (Li et al., in review), maize (Bassu et al., 2014), and sugarcane (Singels et al., 2013). A coordinated modeling exercise was one of the original motivations for AgMIP, and C3MP provides rapid estimation of crop responses to CO2, water, and temperature (CTW) changes, adding dimension and insight into the crop model intercomparisons, while facilitating interactions within the global community of modelers. C3MP also contributes a fast-track, multi-model climate sensitivity assessment for the AgMIP climate and crop modeling teams on Research Track 2 (Fig. 1), which seeks to understand the impact of projected climatic changes on crop production and food security (Rosenzweig et al., 2013; Ruane et al., 2014).

agriculture↗

DOE OSTI2023

Combining variational autoencoders and physical bias for improved microscopy data analysis *

Electron and scanning probe microscopy produce vast amounts of data in the form of images or hyperspectral data, such as electron energy loss spectroscopy or 4D scanning transmission electron microscope, that contain information on a wide range of structural, physical, and chemical properties of materials. To extract valuable insights from these data, it is crucial to identify physically separate regions in the data, such as phases, ferroic variants, and boundaries between them. In order to derive an easily interpretable feature analysis, combining with well-defined boundaries in a principled and unsupervised manner, here we present a physics augmented machine learning method which combines the capability of variational autoencoders to disentangle factors of variability within the data and the physics driven loss function that seeks to minimize the total length of the discontinuities in images corresponding to latent representations. Our method is applied to various materials, including NiO-LSMO, BiFeO 3 , and graphene. The results demonstrate the effectiveness of our approach in extracting meaningful information from large volumes of imaging data. The customized codes of the required functions and classes to develop phyVAE is available at https://github.com/arpanbiswas52/phy-VAE.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS2020

Progress of Coating Stress Compensation of Silicon Mirrors for Lynx X-ray Telescope Mission Concept Using Thermal Oxide Patterning Method

A thermal oxide patterning method has proven to be effective for correcting coating-stress-induced distortion on flat silicon wafers. We report progress on developing this method for correcting curved silicon mirrors distorted by front-side iridium coatings. Owing to the difference in geometry, a finite element model has been established to calculate the appropriate duty cycle maps in thermal oxide hexagon patterns used for compensation. In addition, a photolithographic process, along with three-dimensional printed equipment, has been developed for creating patterns precisely on the back side of curved mirrors. The developed method has been used to recover the original surface shape of two silicon mirrors which are 100-mm long, 0.5-mm thick, having 312-mm radius of curvature, and 30 deg in azimuthal span (Wolter-I geometry). These mirrors’ front sides are sputter-coated by 20-nm iridium layers with ∼-70 N∕m integrated stress. Measurement results show that the developed method can mitigate coating-induced distortion by a factor of ∼5 in RMS height and ∼4 in RMS slope error, corresponding to ∼0.5 arc sec RMS slope error. Residual errors after correction are dominated by midfrequency ripples created during the annealing process, which will be resolved in the future. The presented method is precise and inexpensive and a potential candidate for resolving the coating stress issue for Lynx optics in the future.

x-ray↗

NASA NTRS1989

Parallel linear equation solvers for finite element computations

The overall objective of this research is to develop efficient methods for the solution of linear and nonlinear systems of equations on parallel and supercomputers, and to apply these methods to the solution of problems in structural analysis. Attention has been given so far only to linear equations. The methods considered for the solution of the stiffness equation Kx=f have been Choleski factorization and the conjugate gradient iteration with SSOR and Incomplete Choleski preconditioning. More detail on these methods will be given on subsequent slides. These methods have been used to solve for the static displacements for the mast and panel focus problems in conjunction with the CSM testbed system based on NICE/SPAR.

Ortega, James M.↗

DOE OSTI2022

Zeolitic Imidazolate Framework Membranes: Novel Synthesis Methods and Progress Toward Industrial Use

In the last decade, zeolitic imidazolate frameworks (ZIFs) have been studied extensively for their potential as selective separation membranes. In this review, we highlight unique structural properties of ZIFs that allow them to achieve certain important separations, like that of propylene from propane, and summarize the state of the art in ZIF thin-film deposition on porous substrates and their modification by postsynthesis treatments. We also review the reported membrane performance for representative membrane synthesis approaches and attempt to rank the synthesis methods with respect to potential for scalability. To compare the dependence of membrane performance on membrane synthesis methods and operating conditions, we map out fluxes and separation factors of selected ZIF-8 membranes for propylene/propane separation. Finally, we provide future directions considering the importance of further improvements in scalability, cost effectiveness, and stable performance under industrially relevant conditions.

Chemistry↗

NASA NTRS1998

Prediction of Transitional Flows in the Low Pressure Turbine

Current turbulence models tend to give too early and too short a length of flow transition to turbulence, and hence fail to predict flow separation induced by the adverse pressure gradients and streamline flow curvatures. Our discussion will focus on the development and validation of transition models. The baseline data for model comparisons are the T3 series, which include a range of free-stream turbulence intensity and cover zero-pressure gradient to aft-loaded turbine pressure gradient flows. The method will be based on the conditioned N-S equations and a transport equation for the intermittency factor. First, several of the most popular 2-equation models in predicting flow transition are examined: k-e [Launder-Sharina], k-w [Wilcox], Lien-Leschiziner and SST [Menter] models. All models fail to predict the onset and the length of transition, even for the simplest flat plate with zero-pressure gradient(T3A). Although the predicted onset position of transition can be varied by providing different inlet turbulent energy dissipation rates, the appropriate inlet conditions for turbulence quantities should be adjusted to match the decay of the free-stream turbulence. Arguably, one may adjust the low-Reynolds-number part of the model to predict transition. This approach has so far not been very successful. However, we have found that the low-Reynolds-number model of Launder and Sharma [1974], which is an improved version of Jones and Launder [1972] gave the best overall performance. The Launder and Sharma model was designed to capture flow re-laminarization (a reverse of flow transition), but tends to give rise to a too early and too fast transition in comparison with the physical transition. The three test cases were for flows with zero pressure gradient but with different free-stream turbulent intensities. The same can be said about the model when considering flows subject to pressure gradient(T3C1). To capture the effects of transition using existing turbulence models, one approach is to make use of the concept of the intermittency to predict the flow transition. It was originally based on the intermittency distribution of Narasimha [1957], and then gradually evolved into a transport equation for the intermittency factor. Gostelow and associates [1994,1995] have made some improvements to Narasimha's method in an attempt to account for both favorable and adverse pressure gradients. Their approach is based on a linear, explicit combination of laminar and turbulent solutions. This approach fails to predict the overshoot of the skin friction on a flat plate near the end of transition zone, even though the length of transition is well predicted. The major flaw of Gostelow's approach is that it assumes the non-turbulent part being the laminar solution and the turbulent part being the turbulent solution and they do not interact across the transitional region. The technique in condition averaging the flow equations in intermittent flows was first introduced by Libby [1975] and Dopazo [1977] and further refined by Dick and associates [1988, 1996]. This approach employs two sets of transport equations for the non-turbulent part and the other for the turbulent part. The advantage of this approach is that it allows the interaction of non-turbulent and turbulent velocities through the introduction of additional source terms in the continuity and momentum equations for the non-turbulent and turbulent velocities. However, the strong coupling of the two sets of equations has caused some numerical difficulties, which requires special attention. The prediction of the skin friction can be improved by this approach via the implicit coupling of non-turbulent and turbulent velocity flelds. Another improvement of the interrmittency model can be further made by allowing the intermittency to vary in the cross-stream direction. This is one step prior to testing any proposal for the transport equation for the intermittency factor. Instead of solving the transport equation for the intermittency factor, the distribution for the intermittency factor is prescribed by Klebanoff's empirical formula [1955]. The skin friction is very well predicted by this new modification, including the overshoot of the profile near the end of the transition zone. The outcome of this study is very encouraging since it indicates that the proper description of the intermittency distribution is the key to the success of the model prediction. This study will be used to guide us on the modelling of the intermittency transport equation.

Huang, George↗

DOE OSTI2022

Gapped-filtering for efficient Chebyshev expansion of the density projection operator

In this paper, we develop the gapped-filtering method, whereby a short Chebyshev expansion accurately represents the density-matrix operator. The method optimizes the Chebyshev coefficients to give the correct density matrix at all energies except within the gapped region where there are no eigenstates. Gapped filtering reduces the number of required terms in the Chebyshev expansion compared to traditional expansion methods, as long as one knows or can determine efficiently the HOMO and LUMO positions. The reduction is especially noticeable (factors of 2-3) when high accuracy is sought. To exemplify the method, we use gapped-filtering to increase the efficiency of stochastic-GW calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1993

Effects of free-stream turbulence on boundary-layer transition

Traditional methods for the prediction of the transition location rest either on lump parameters like the momentum thickness or on N factors obtained from the approximate linear stability characteristics of the flow. These methods account for the effects of the disturbance environment by empirical correlation with observed transition locations. We have developed a highly efficient and accurate approach to stability analysis and transition simulation based on parabolized stability equations (PSE). Here we present a first application of this approach to realistic problems of practical interest. We study the effects of the turbulence level on transition and heat transfer in flows over flat plates, curved plates, and the stator blade of a gas-turbine model. We discuss the formulation of an input model for the PSE analysis in the light of the present knowledge of receptivity mechanisms. The input model is used to compute transition locations for various configurations and flow conditions. The results compare favorably with experimental data.

Herbert, TH.↗

DOE OSTI2022

Extreme events, energy security and equality through micro- and macro-levels: Concepts, challenges and methods

Low-income households face long-standing challenges of energy insecurity and inequality (EII). During extreme events (e.g., disasters and pandemics) these challenges are especially severe for vulnerable populations reliant on energy for health, education, and well-being. However, many EII studies rarely incorporate the micro- and macro-perspectives of resilience and reliability of energy and internet infrastructure and social-psychological factors. To remedy this gap, we first address the impacts of extreme events on EII among vulnerable populations. Second, we evaluate the driving factors of EII and how they change during disasters. Third, we situate these inequalities within broader energy systems and pinpoint the importance of equitable infrastructure systems by examining infrastructure reliability and resilience and the role of renewable technologies. Then, we consider the factors influencing energy consumption, such as energy practices, socio-psychological factors, and internet access. Finally, we propose interdisciplinary research methods to study these issues during extreme events and provide recommendations.

Chen, CF↗

NASA NTRS1993

A transient FETI methodology for large-scale parallel implicit computations in structural mechanics, part 2

Explicit codes are often used to simulate the nonlinear dynamics of large-scale structural systems, even for low frequency response, because the storage and CPU requirements entailed by the repeated factorizations traditionally found in implicit codes rapidly overwhelm the available computing resources. With the advent of parallel processing, this trend is accelerating because explicit schemes are also easier to parallellize than implicit ones. However, the time step restriction imposed by the Courant stability condition on all explicit schemes cannot yet and perhaps will never be offset by the speed of parallel hardware. Therefore, it is essential to develop efficient and robust alternatives to direct methods that are also amenable to massively parallel processing because implicit codes using unconditionally stable time-integration algorithms are computationally more efficient than explicit codes when simulating low-frequency dynamics. Here we present a domain decomposition method for implicit schemes that requires significantly less storage than factorization algorithms, that is several times faster than other popular direct and iterative methods, that can be easily implemented on both shared and local memory parallel processors, and that is both computationally and communication-wise efficient. The proposed transient domain decomposition method is an extension of the method of Finite Element Tearing and Interconnecting (FETI) developed by Farhat and Roux for the solution of static problems. Serial and parallel performance results on the CRAY Y-MP/8 and the iPSC-860/128 systems are reported and analyzed for realistic structural dynamics problems. These results establish the superiority of the FETI method over both the serial/parallel conjugate gradient algorithm with diagonal scaling and the serial/parallel direct method, and contrast the computational power of the iPSC-860/128 parallel processor with that of the CRAY Y-MP/8 system.

Farhat, Charbel↗

NASA NTRS1992

A transient FETI methodology for large-scale parallel implicit computations in structural mechanics

Explicit codes are often used to simulate the nonlinear dynamics of large-scale structural systems, even for low frequency response, because the storage and CPU requirements entailed by the repeated factorizations traditionally found in implicit codes rapidly overwhelm the available computing resources. With the advent of parallel processing, this trend is accelerating because explicit schemes are also easier to parallelize than implicit ones. However, the time step restriction imposed by the Courant stability condition on all explicit schemes cannot yet -- and perhaps will never -- be offset by the speed of parallel hardware. Therefore, it is essential to develop efficient and robust alternatives to direct methods that are also amenable to massively parallel processing because implicit codes using unconditionally stable time-integration algorithms are computationally more efficient when simulating low-frequency dynamics. Here we present a domain decomposition method for implicit schemes that requires significantly less storage than factorization algorithms, that is several times faster than other popular direct and iterative methods, that can be easily implemented on both shared and local memory parallel processors, and that is both computationally and communication-wise efficient. The proposed transient domain decomposition method is an extension of the method of Finite Element Tearing and Interconnecting (FETI) developed by Farhat and Roux for the solution of static problems. Serial and parallel performance results on the CRAY Y-MP/8 and the iPSC-860/128 systems are reported and analyzed for realistic structural dynamics problems. These results establish the superiority of the FETI method over both the serial/parallel conjugate gradient algorithm with diagonal scaling and the serial/parallel direct method, and contrast the computational power of the iPSC-860/128 parallel processor with that of the CRAY Y-MP/8 system.

Farhat, Charbel↗

DOE OSTI2021

A Targeted Sequencing Assay for Serotyping Escherichia coli Using AgriSeq Technology

The gold standard method for serotyping Escherichia coli has relied on antisera-based typing of the O- and H-antigens, which is labor intensive and often unreliable. In the post-genomic era, sequence-based assays are potentially faster to provide results, could combine O-serogrouping and H-typing in a single test, and could simultaneously screen for the presence of other genetic markers of interest such as virulence factors. Whole genome sequencing is one approach; however, this method has limited multiplexing capabilities, and only a small fraction of the sequence is informative for subtyping or identifying virulence potential. A targeted, sequence-based assay and accompanying software for data analysis would be a great improvement over the currently available methods for serotyping. The purpose of this study was to develop a high-throughput, molecular method for serotyping E. coli by sequencing the genes that are required for production of O- and H-antigens, as well as to develop software for data analysis and serotype identification. To expand the utility of the assay, targets for the virulence factors, Shiga toxins ( stx 1 , and stx 2 ) and intimin ( eae ) were included. To validate the assay, genomic DNA was extracted from O-serogroup and H-type standard strains and from Shiga toxin-producing E. coli , the targeted regions were amplified, and then sequencing libraries were prepared from the amplified products followed by sequencing of the libraries on the Ion S5™ sequencer. The resulting sequence files were analyzed via the SeroType Caller™ software for identification of O-serogroup, H-type, and presence of stx 1 , stx 2 , and eae . We successfully identified 169 O-serogroups and 41 H-types. The assay also routinely detected the presence of stx 1a,c,d (3 of 3 strains), stx 2c−e,g (8 of 8 strains), stx 2f (1 strain), and eae (6 of 6 strains). Taken together, the high-throughput, sequence-based method presented here is a reliable alternative to antisera-based serotyping methods for E. coli .

Elder, Jacob R.↗

DOE OSTI2024

Helium-4 gravitational form factors: Exchange currents

We evaluate the leading exchange corrections to the helium-4 gravitational form factors (GFFs) to momenta of the order of the nucleon mass. We use both the K-harmonic method with simple pair nucleon potential, and a Jastrow trial function using the Argonne 𝑣 14 potential, to evaluate the helium-4 GFFs. The exchange current contributions include the pair interaction, plus the seagull and the pion exchange interactions, modulo the recoil corrections. To estimate the off-shellness of the pion nucleon coupling in this momenta range, we discuss the results using either the pseudoscalar (PS) or pseudovector (PV) pion-nucleon couplings. When the PV coupling is used, the pair diagram contribution is higher order in the nonrelativistic expansion. The results for the helium-4 A-GFF are comparable to those given by the impulse approximation, especially for the PS coupling using both the K-harmonic method and variational method. The exchange current contributions with the PS coupling for the charge form factor of helium-4, yield better agreement with the existing data over a broad range of momenta, especially when the Argonne 𝑣 14 potential including the D-wave admixture is used.

A ≤ 5↗

NASA NTRS1983

Stability analysis of intermediate boundary conditions in approximate factorization schemes

In many cases, approximate factorization schemes have provided a significant increase in efficiency over previously used solution methods in certain problems. The present investigation is concerned with the importance of intermediate boundary conditions in approximate factorization schemes, taking into account a specific example regarding a boundary-induced stability restriction in a scheme for the transonic full-potential equation. The considered scheme has been discussed by Holst (1979). Holst's scheme is a variation of the AF2 schemes described by Ballhaus and Steger (1975). The application of the AF2 scheme to the two-dimensional Laplace equation in a rectangle is studied giving attention to the stability of the scheme in connection with various boundary conditions for the intermediate variable.

South, J. C., Jr.↗

NASA NTRS1989

GOES I/M image navigation and registration

Image Navigation and Registration (INR) is the system that will be used on future Geostationary Operational Environmental Satellite (GOES) missions to locate and register radiometric imagery data. It consists of a semiclosed loop system with a ground-based segment that generates coefficients to perform image motion compensation (IMC). The IMC coefficients are uplinked to the satellite-based segment, where they are used to adjust the displacement of the imagery data due to movement of the imaging instrument line-of-sight. The flight dynamics aspects of the INR system is discussed in terms of the attitude and orbit determination, attitude pointing, and attitude and orbit control needed to perform INR. The modeling used in the determination of orbit and attitude is discussed, along with the method of on-orbit control used in the INR system, and various factors that affect stability. Also discussed are potential error sources inherent in the INR system and the operational methods of compensating for these errors.

Fiorello, J. L., Jr.↗

DOE OSTI2023

Laser Driven Hydrothermal Processing (LDHP) for Rapid Dissolution of Glassy Nuclear Debris

Discovered in 2009 on the SLAB laser system at LLNL, Laser-Driven Hydrothermal Processing (LDHP) both spalls and transiently dissolves a wetted surface (e.g., rock, concrete, paint) when that surface is struck by a laser beam pulse under specific conditions. Previous projects have explored the application of LDHP to expediting sample dissolution for post-detonation nuclear forensics. Early results were very promising, indicating that this method could decrease material dissolution time by a factor of 10 or more. However, the large and complicated lasers previously used are impractical for this application. Under the R^3 venture we aimed to benchmark the LDHP capabilities of turn-key commercial off-the-shelf (COTS) lasers, and to establish how this process could be integrated into an analytical workflow for nuclear forensic analysis. We demonstrated that geological rock standards could be quickly processed into fine powder by LDHP followed by complete acid digestion in 15 minutes, without elemental fractionation that would skew analytical results. We determined that the COTS laser can be used for LDHP that could potentially reduce dissolution time for solid nuclear debris samples in post detonation nuclear forensics. One downside to this approach, however, is that the laser interactions are different with every material, so extensive method validation would be required for all likely material compositions. The LDHP processing rates of several materials were tested and determined to be favorable in the beginning with significant slowdown due to particulate scattering. This could be mitigated either by stopping and centrifuging the sample at specific intervals, or through design of a continuous solvent flow through mechanism.

36 MATERIALS SCIENCE↗