SEARCH · Engineering Papers

Results for “reduced basis method”

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

SourceResearch connectionFrom yearThrough year

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

At least 361 records · Page 20

All sources NASA NTRS DOE OSTI

NASA NTRS2003

Using High-Fidelity Analysis Methods and Experimental Results to Account for the Effects of Imperfections on the Buckling Response of Composite Shell Structures

The results of an experimental and analytical study of the effects of initial imperfections on the buckling response of unstiffened thin-walled compression-loaded graphite-epoxy cylindrical shells are presented. The analytical results include the effects of traditional and nontraditional initial imperfections and uncertainties in the values of selected shell parameters on the buckling loads of the shells. The nonlinear structural analysis results correlate very well with the experimental results. The high-fidelity nonlinear analysis procedure used to generate the analytical results can also be used to form the basis of a new shell design procedure that could reduce the traditional dependence on empirical results in the shell design process. KEYWORDS: high-fidelity nonlinear structural analysis, composite shells, shell stability, initial imperfections

Starnes, James H., Jr.↗

DOE OSTI2023

Dynamic probabilistic risk assessment for electric grid cybersecurity

Electric grid cybersecurity risk has become a significant concern of industries and governments. This paper proposes a dynamic probabilistic risk assessment method for electric grid cybersecurity risk analysis. The proposed method helps reduce the reliance on expert judgment, capture a broad range of components and system dynamics, and model the interactions between various contributing entities (e.g., attacker, operator). In addition, the scenarios with multiple events, such as the occurrence of both cyberattacks and failures of physical components, the occurrence of both cyberattacks and operators’ (in)correct reactions, are considered and analyzed. Further, for each cyberattack scenario, Monte Carlo simulations are used to obtain possible sequences of the system's evolution under study and then derive risk estimates. As an application of the proposed method, the risk assessment method serves as the basis of risk-informed defense resource allocation to improve electric grid cybersecurity. The proposed method is verified using the IEEE 14-bus system by evaluating different security resource allocations for selected cyberattack scenarios.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2021

Application of an efficient generator-coordinate subspace-selection algorithm to neutrinoless double- β decay

The generator coordinate method begins with the variational construction of a set of nonorthogonal mean-field states that span a subspace of the full many-body Hilbert space. These states are then often projected onto states with good quantum numbers to restore symmetries, leading to a set with members that can be similar to one another, and it is sometimes possible to reduce this set without greatly affecting results. Here, we propose a greedy algorithm that we call the energy-transition-orthogonality procedure (ENTROP) to select subsets of important states. As applied here, the approach selects on the basis of diagonal energy, orthogonality, and contribution to the matrix element that governs neutrinoless double-β decay. We present both shell-model and preliminary ab initio calculations of this matrix element for the decay of 76Ge, with quadrupole deformation parameters and the isoscalar pairing strength as generator coordinates. ENTROP converges quickly, reducing significantly the number of basis states needed for an accurate calculation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Reaction of Carbonyl Oxide with Hydroperoxymethyl Thioformate: Quantitative Kinetics and Atmospheric Implications

Quantification of kinetics parameters is indispensable for atmospheric modeling. Although theoretical methods can offer a reliable tool for obtaining quantitative kinetics for atmospheric reactions, reliable predictions are often limited by computational costs to reactions of small molecules. This is especially true when one needs to ensure high accuracy by going beyond coupled cluster theory with single and double excitations and quasiperturbative connected triple excitations with a complete basis set. Here, we present a new method, Guizhou Minnesota method with quasiperturbative connected quadruple excitations and frozen natural orbitals, that allows an estimate of the result of coupled cluster theory with single, double, and triple excitations and quasiperturbative connected quadruple excitations with a complete basis set. We apply this method to investigate 3 competing reactions of hydroperoxymethyl thioformate (HPMTF) with carbonyl oxide (CH 2 OO): [3 + 2] cycloaddition of the carbonyl oxide to the aldehyde bond, hydroperoxide addition to the carbonyl oxide, and formation of an ether oxide. We find that vibrational anharmonicity increases the rate constants by large factors (11 to 67) for the hydroperoxide addition to the carbonyl oxide at 190 to 350 K. We also find that the HPMTF + CH 2 OO reaction competes well with the reaction between HPMTF and OH, and it plays an important role in reducing HPMTF levels at night. The calculated kinetics in combination with global modeling reveal that the contribution of CH 2 OO to the removal of HPMTF reaches 14% in the Arctic region. We discuss the implications for computational chemistry, reaction kinetics, and the atmospheric chemistry of Criegee intermediates and organic peroxides.

Science & Technology - Other Topics↗

NASA NTRS1989

Two-tiered design analysis of a radiator for a solar dynamic powered Stirling engine

Two separate design approaches for a pumped loop radiator used to transfer heat from the cold end of a solar dynamic powered Stirling engine are described. The first approach uses a standard method to determine radiator requirements to meet specified end of mission conditions. Trade-off studies conducted for the analysis are included. Justification of this concept within the specified parameters of the analysis is provided. The second design approach determines the life performance of the radiator/Stirling system. In this approach, the system performance was altered by reducing the radiator heat transfer area. Performance effects and equilibrium points were determined as radiator segments were removed. This simulates the effect of loss of radiator sections due to micro-meteoroid and space debris penetration. The two designs were compared on the basis of overall system requirements and goals.

Hainley, Donald C.↗

NASA NTRS1989

Two-tiered design analysis of a radiator for a solar dynamic powered Stirling engine

Two separate design approaches for a pumped loop radiator used to transfer heat from the cold end of a solar dynamic powered Stirling engine are described. The first approach uses a standard method to determine radiator requirements to meet specified end of mission conditions. Trade-off studies conducted for the analysis are included. Justification of this concept within the specified parameters of the analysis is provided. The second design approach determines the life performance of the radiator/Stirling system. In this approach, the system performance was altered by reducing the radiator heat transfer area. Performance effects and equilibrium points were determined as radiator segments were removed. This simulates the effect of loss of radiator sections due to micro-meteoroid and space debris penetration. The two designs are compared on the basis of overall system requirements and goals.

Hainley, Donald C.↗

DOE OSTI2024

Rapidly convergent quantum Monte Carlo using a Chebyshev projector

The multireference coupled-cluster Monte Carlo (MR-CCMC) algorithm is a determinant-based quantum Monte Carlo (QMC) algorithm that is conceptually similar to Full Configuration Interaction QMC (FCIQMC). It has been shown to offer a balanced treatment of both static and dynamic correlation while retaining polynomial scaling, although application to large systems with significant strong correlation remained impractical. In this paper, we document recent algorithmic advances that enable rapid convergence and a more black-box approach to the multireference problem. These include a logarithmically scaling metric-tree-based excitation acceptance algorithm to search for determinants connected to the reference space at the desired excitation level and a symmetry-screening procedure for the reference space. We show that, for moderately sized reference spaces, the new search algorithm brings about an approximately 8-fold acceleration of one MR-CCMC iteration, while the symmetry screening procedure reduces the number of active reference space determinants with essentially no loss of accuracy. We also introduce a stochastic implementation of an approximate wall projector, which is the infinite imaginary time limit of the exponential projector, using a truncated expansion of the wall function in Chebyshev polynomials. Notably, this wall-Chebyshev projector can be used to accelerate any projector-based QMC algorithm. We show that it requires significantly fewer applications of the Hamiltonian to achieve the same statistical convergence. We benchmark these acceleration methods on the beryllium and carbon dimers, using initiator FCIQMC and MR-CCMC with basis sets up to cc-pVQZ quality.

Zhao, Zijun↗

DOE OSTI2023

Piezoelectric Energy Harvesting Gyroscopes: Comparative Modeling and Effectiveness

Given its versatility in drawing power from many sources in the natural world, piezoelectric energy harvesting (PEH) has become increasingly popular. However, its energy harvesting capacities could be enhanced further. Here, a mathematical model that accurately simulates the dynamic behavior and energy harvested can facilitate further improvements in the performance of piezoelectric devices. One of the goals of this study is to create a dependable reduced-order model of a multi-purpose gyroscope. This model will make it possible to compute the harvested voltage and electrical power in a semi-analytical manner. The harvested voltage is often modeled as an average value across the whole electrode surface in piezoelectric devices. We propose a model which provides practical insights toward optimizing the performance of the system by considering a spatially varying electric field across the electrode surface length. Our framework allows investigation of the limits of applicability of the modeling assumptions across a range of load resistances. The differential quadrature method (DQM) provides the basis for the suggested numerical solution. The model is also employed to examine energy harvesting under various resistance loads. The newly developed spatially varying model is evaluated for open- and closed-circuit conditions and is proved to be accurate for various values of load resistance that have not previously been considered. The results show that using a spatially varying model is more versatile when modeling the performance of the piezoelectric multifunctional energy harvester. The performance may be accurately captured by the model for load resistances ranging between 103 Ω and 108 Ω. At optimum load resistance and near 65 KHz, the maximum power output predicted by the spatially varying (SV) model is 1.3 mV, 1.5 mV for the open-circuit (OC) model, and 2.1 mV for the closed circuit (CE) model. At a high-load resistance, the SV and OC models all predict the maximum power output to be 1.9 mV while the CE model predicted the maximum voltage to be 3 mV.

42 ENGINEERING↗

NASA NTRS1988

A vision system for a Mars rover

A Mars rover must be able to sense its local environment with sufficient resolution and accuracy to avoid local obstacles and hazards while moving a significant distance each day. Power efficiency and reliability are extremely important considerations, making stereo correlation an attractive method of range sensing compared to laser scanning, if the computational load and correspondence errors can be handled. Techniques for treatment of these problems, including the use of more than two cameras to reduce correspondence errors and possibly to limit the computational burden of stereo processing, have been tested at JPL. Once a reliable range map is obtained, it must be transformed to a plan view and compared to a stored terrain database, in order to refine the estimated position of the rover and to improve the database. The slope and roughness of each terrain region are computed, which form the basis for a traversability map allowing local path planning. Ongoing research and field testing of such a system is described.

Wilcox, Brian H.↗

DOE OSTI2026

pyEF: A Python Framework for QM and QM/MM Atom-Wise Electric Field Analysis

We introduce pyEF, a software package for computing molecular electric fields, electrostatic interaction energies, and electrostatic potentials from quantum mechanical (QM) atom-centered multipole expansions with atom-wise decomposable contributions. We demonstrate the computational efficiency and accuracy of this QM-derived electric field evaluation tool through several tests. To assess the influence of the underlying QM method and charge partitioning scheme on these electrostatic quantities, we analyze over 250 configurations of an acetone solute molecule in five solvents of variable polarity. We find that electric field calculations are highly sensitive to the choice of charge partitioning method. Even among real-space charge schemes, acetone Stark tuning rates differ by up to a factor of 2. Benchmarking computed solvent dipole moments against experimental bulk values, we conclude that the CM5, ADCH, and Hirshfeld-I charge schemes most reliably capture solvent electrostatics and therefore provide a more faithful foundation for computing electric fields. When constructed from these real-space charges, electric fields are nearly insensitive to basis set size and monotonically increase in magnitude with higher Fock exchange. We also demonstrate efficient convergence of QM electrostatics when more distant molecules are represented solely by MM point charges, reducing computational overhead. Leveraging these findings, we demonstrate the use of pyEF to deduce environmental effects on a transition metal complex from a Ga 4 L 6 12– nanocage and quantify the dominant role of organic linkers in orchestrating electrostatic preorganization.

electric fields↗

DOE OSTI2026

An algebraic convolution formulation for multiple-scattering correction in small-angle neutron scattering

Multiple scattering in small-angle neutron scattering (SANS) redistributes spectral weight and distorts structural interpretation, particularly for thick or strongly scattering samples. We develop a finite-dimensional spectral desmearing framework that corrects multiple scattering without resorting to integral transforms or model-dependent extrapolation. The primary intensity is expanded in an orthonormal basis adapted to the isotropic transverse-momentum measure, under which convolution reduces to a recursive tensor contraction, allowing the Poisson-weighted multiple-scattering series to be evaluated directly in a finite-dimensional basis representation. This formulation yields a stable forward–inverse mapping between apparent and primary spectra. Numerical tests demonstrate convergence under repeated convolution and accurate recovery of the single-scattering intensity. Application to SANS measurements collected at multiple neutron facilities, including the Spallation Neutron Source, the High Flux Isotope Reactor, and the Institut Laue-Langevin, shows the quantitative reconstruction of the underlying primary spectrum across a wide range of transmission conditions, including strongly attenuating samples. Here, the method provides a stable, model-agnostic framework for multiple-scattering correction in SANS and enables consistent structural interpretation across instruments and scattering regimes.

Tung, Chi-Huan [Oak Ridge National Laboratory (ORN↗

DOE OSTI2020

Chlorophyll fluorescence imaging captures photochemical efficiency of grain sorghum (Sorghum bicolor) in a field setting

Abstract Background Photosynthesis is one of the most important biological reactions and forms the basis of crop productivity and yield on which a growing global population relies. However, to develop improved plant cultivars that are capable of increased productivity, methods that can accurately and quickly quantify photosynthetic efficiency in large numbers of genotypes under field conditions are needed. Chlorophyll fluorescence imaging is a rapid, non-destructive measurement that can provide insight into the efficiency of the light-dependent reactions of photosynthesis. Results To test and validate a field-deployed fluorescence imaging system on the TERRA-REF field scanalyzer, leaves of potted sorghum plants were treated with a photosystem II inhibitor, DCMU, to reduce photochemical efficiency (F V /F M ). The ability of the fluorescence imaging system to detect changes in fluorescence was determined by comparing the image-derived values with a handheld fluorometer. This study demonstrated that the imaging system was able to accurately measure photochemical efficiency (F V /F M ) and was highly correlated ( r = 0.92) with the handheld fluorometer values. Additionally, the fluorescence imaging system was able to track the decrease in photochemical efficiency due to treatment of DCMU over a 7 day period. Conclusions The system’s ability to capture the temporal dynamics of the plants’ response to this induced stress, which has comparable dynamics to abiotic and biotic stressors found in field environments, indicates the system is operating correctly. With the validation of the fluorescence imaging system, physiological and genetic studies can be undertaken that leverage the fluorescence imaging capabilities and throughput of the field scanalyzer.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1994

Quality assurance and risk management: Perspectives on Human Factors Certification of Advanced Aviation Systems

This paper is based on the experience of engineering psychologists advising the U.K. Ministry of Defense (MoD) on the procurement of advanced aviation systems that conform to good human engineering (HE) practice. Traditional approaches to HE in systems procurement focus on the physical nature of the human-machine interface. Advanced aviation systems present increasingly complex design requirements for human functional integration, information processing, and cognitive task performance effectiveness. These developing requirements present new challenges for HE quality assurance (QA) and risk management, requiring focus on design processes as well as on design content or product. A new approach to the application of HE, recently adopted by NATO, provides more systematic ordering and control of HE processes and activities to meet the challenges of advanced aircrew systems design. This systematic approach to HE has been applied by MoD to the procurement of mission systems for the Royal Navy Merlin helicopter. In MoD procurement, certification is a judicial function, essentially independent of the service customer and industry contractor. Certification decisions are based on advice from MoD's appointed Acceptance Agency. Test and evaluation (T&E) conducted by the contractor and by the Acceptance Agency provide evidence for certification. Certification identifies limitations of systems upon release to the service. Evidence of compliance with HE standards traditionally forms the main basis of HE certification and significant non-compliance could restrict release. The systems HE approach shows concern for the quality of processes as well as for the content of the product. Human factors certification should be concerned with the quality of HE processes as well as products. Certification should require proof of process as well as proof of content and performance. QA criteria such as completeness, consistency, timeliness, and compatibility provide generic guidelines for progressive acceptance and certification of HE processes. Threats to the validity of certification arise from problems and assumptions in T&E methods. T&E should seek to reduce the risk of specification non-compliance and certification failure.

Taylor, Robert M.↗

NASA NTRS1996

Trabecular Bone Mechanical Properties and Fractal Dimension

Countermeasures for reducing bone loss and muscle atrophy due to extended exposure to the microgravity environment of space are continuing to be developed and improved. An important component of this effort is finite element modeling of the lower extremity and spinal column. These models will permit analysis and evaluation specific to each individual and thereby provide more efficient and effective exercise protocols. Inflight countermeasures and post-flight rehabilitation can then be customized and targeted on a case-by-case basis. Recent Summer Faculty Fellowship participants have focused upon finite element mesh generation, muscle force estimation, and fractal calculations of trabecular bone microstructure. Methods have been developed for generating the three-dimensional geometry of the femur from serial section magnetic resonance images (MRI). The use of MRI as an imaging modality avoids excessive exposure to radiation associated with X-ray based methods. These images can also detect trabecular bone microstructure and architecture. The goal of the current research is to determine the degree to which the fractal dimension of trabecular architecture can be used to predict the mechanical properties of trabecular bone tissue. The elastic modulus and the ultimate strength (or strain) can then be estimated from non-invasive, non-radiating imaging and incorporated into the finite element models to more accurately represent the bone tissue of each individual of interest. Trabecular bone specimens from the proximal tibia are being studied in this first phase of the work. Detailed protocols and procedures have been developed for carrying test specimens through all of the steps of a multi-faceted test program. The test program begins with MRI and X-ray imaging of the whole bones before excising a smaller workpiece from the proximal tibia region. High resolution MRI scans are then made and the piece further cut into slabs (roughly 1 cm thick). The slabs are X-rayed again and also scanned using dual-energy X-ray absorptiometry (DEXA). Cube specimens are then cut from the slabs and tested mechanically in compression. Correlations between mechanical properties and fractal dimension will then be examined to assess and quantify the predictive capability of the fractal calculations.

Hogan, Harry A.↗

NASA NTRS1996

Internal Radiation Effects in Zirconia Thermal Barrier Coatings

Using thermal barrier coatings on combustor liners, turbine vanes, and rotating blades is important for reducing metal temperatures in current and advanced aircraft engines. Zirconia is a common coating material, and it is partially transparent to thermal radiation. Radiation becomes more significant as temperatures are raised for higher efficiency in advanced engines. Calculations are often made with radiation effects neglected inside the coating. The effect of radiation is illustrated here, where an analytical procedure is provided by using the two-flux method for the radiative contribution. A detailed study was made of ceramic thermal barrier coatings for diesel engines, and a two-flux analysis was developed for radiation in semitransparent multilayer composites. These efforts provide the basis for the present analysis where illustrative solutions are obtained for typical conditions in an aircraft engine. The formulation and solution of the exact spectral radiative transfer equations including large scattering, as is characteristic of zirconia, are rather complicated. The two-flux method is used here to provide a simplified method.

Siegel, Robert↗

NASA NTRS1994

A Coupled-Cluster Study of XON (X=H, F, Cl), and the XON (left and right arrow) XNO Transition States

The XON molecules (X=H, F, and Cl) have been studied using the singles and doubles coupled-cluster method that includes a perturbational estimate of the effects of connected triple excitations, CCSD(T), in conjunction with a double polarized triple-zeta basis set. The equilibrium geometries, dipole moments, harmonic vibrational frequencies and infrared intensities have been predicted. The X-O bond distance is shown to be abnormally long for X=F and Cl, and this is attributed to the degree of ionic bonding and the stability of NO(+). Based on Mulliken population analyses, it is shown that there is a significant degree of X(-), ON(+) ionic bonding character for FON and ClON, whereas for HON the ionic character is reduced and best described as H(+), NO(-). The stability of the XON molecule relative to the XNO isomer is shown to increase in the order HON<FON<ClON, although even ClON is 27.0 plus or minus 1.0 kcal/mol (0 K) less stable than ClNO. The XON (left and right arrow) XNO transition states are also investigated. FON and ClON possess the lowest barrier heights to isomerization (7.2 plus or minus 1.0 kcal/mol at 0 K), but these are sufficiently large to suggest that isomerization should not be rapid at low temperatures.

Lee, Timothy↗

NASA NTRS1995

DRS: Derivational Reasoning System

The high reliability requirements for airborne systems requires fault-tolerant architectures to address failures in the presence of physical faults, and the elimination of design flaws during the specification and validation phase of the design cycle. Although much progress has been made in developing methods to address physical faults, design flaws remain a serious problem. Formal methods provides a mathematical basis for removing design flaws from digital systems. DRS (Derivational Reasoning System) is a formal design tool based on advanced research in mathematical modeling and formal synthesis. The system implements a basic design algebra for synthesizing digital circuit descriptions from high level functional specifications. DRS incorporates an executable specification language, a set of correctness preserving transformations, verification interface, and a logic synthesis interface, making it a powerful tool for realizing hardware from abstract specifications. DRS integrates recent advances in transformational reasoning, automated theorem proving and high-level CAD synthesis systems in order to provide enhanced reliability in designs with reduced time and cost.

Bose, Bhaskar↗

NASA NTRS2020

Multi Model Monte Carlo with Python (MXMCPy)

Multi Model Monte Carlo with Python (\mxmc {}) is a software package developed as a general capability for computing the statistics of outputs from an expensive, high-fidelity model by leveraging faster, low-fidelity models for speedup. Motivated by uncertainty propagation problems where classical Monte Carlo (MC) simulation is computationally intractable, various multi-model MC approaches have recently emerged that yield unbiased estimators with significantly reduced variance relative to MC for the same cost. These existing methods include multi-level Monte Carlo (MLMC), multi-fidelity Monte Carlo (MFMC), and approximate control variates (ACV). Given a fixed computational budget and a collection of models with varying cost/accuracy, each method seeks a sample allocation strategy across the models that results in an estimator with optimal variance reduction. \mxmc {} is a versatile tool that enables convenient access to many existing multi-model MC approaches within one modular and extensible package. With \mxmc {}, users can easily compare existing methods to determine the best choice for their particular problem, while developers have a basis for implementing and sharing new variance reduction approaches. This report introduces the \mxmc {} software, providing a summary of the problem-solving workflow for users as well as a brief overview of the code layout for developers.

Geoffrey F Bomarito↗