SEARCH · Engineering Papers

Results for “solution regularity”

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 325 records · Page 18

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Paired autoencoders for likelihood-free estimation in inverse problems

Abstract We consider the solution of nonlinear inverse problems where the forward problem is a discretization of a partial differential equation. Such problems are notoriously difficult to solve in practice and require minimizing a combination of a data-fit term and a regularization term. The main computational bottleneck of typical algorithms is the direct estimation of the data misfit. Therefore, likelihood-free approaches have become appealing alternatives. Nonetheless, difficulties in generalization and limitations in accuracy have hindered their broader utility and applicability. In this work, we use a paired autoencoder framework as a likelihood-free estimator (LFE) for inverse problems. We show that the use of such an architecture allows us to construct a solution efficiently and to overcome some known open problems when using LFEs. In particular, our framework can assess the quality of the solution and improve on it if needed. We demonstrate the viability of our approach using examples from full waveform inversion and inverse electromagnetic imaging.

Chung, Matthias (ORCID:0000000178224539)↗

DOE OSTI2024

Dynamic and relativistic effects on Pt–Pt indirect spin–spin coupling in aqueous solution studied by ab initio molecular dynamics and two- vs four-component density functional NMR calculations

Treating 195Pt nuclear magnetic resonance parameters in solution remains a considerable challenge from a quantum chemistry point of view, requiring a high level of theory that simultaneously takes into account the relativistic effects, the dynamic treatment of the solvent–solute system, and the dynamic electron correlation. A combination of Car-Parrinello molecular dynamics (CPMD) and relativistic calculations based on two-component zeroth order regular approximation spin–orbit Kohn–Sham (2c-ZKS) and four-component Dirac–Kohn–Sham (4c-DKS) Hamiltonians is performed to address the solvent effect (water) on the conformational changes and JPtPt1 coupling. A series of bridged PtIII dinuclear complexes [L1–Pt2(NH3)4(Am)2–L2]n+ (Am = α–pyrrolidonate and pivalamidate; L = H2O, Cl−, and Br−) are studied. The computed Pt–Pt coupling is strongly dependent on the conformational dynamics of the complexes, which, in turn, is correlated with the trans influence among axial ligands and with the angle N–C–O from the bridging ligands. The J-coupling is decomposed in terms of dynamic contributions. The decomposition reveals that the vibrational and explicit solvation contributions reduce JPtPt1 of diaquo complexes (L1 = L2 = H2O) in comparison to the static gas-phase magnitude, whereas the implicit solvation and bulk contributions correspond to an increase in JPtPt1 in dihalo (L1 = L2 = X−) and aquahalo (L1 = H2O; L2 = X−) complexes. Relativistic treatment combined with CPMD shows that the 2c-ZKS Hamiltonian performs as well as 4c-DKS for the JPtPt1 coupling.

Chemistry↗

DOE OSTI2021

Integrating Electromagnetic Acoustic Transducers in a Modular Robotic Gripper for Inspecting Tubular Components

Tubular structures are critical components in infrastructure such as power plants. Throughout their life, they are subjected to extreme conditions or suffer from defects such as corrosion and cracks. Although regular inspection of these components is necessary, such inspection is limited by safety-related risks and limited access for human inspection. Robots can provide a solution for automatic inspection. The main challenge, however, lies in integrating sensors for nondestructive evaluation with robotic platforms. As part of developing a versatile lizard-inspired tube inspector robot, in this study the authors propose to integrate electromagnetic acoustic transducers into a modular robotic gripper for use in automated ultrasonic inspection. In particular, spiral coils with cylindrical magnets are integrated into a novel friction-based gripper to excite Lamb waves in thin cylindrical structures. To evaluate the performance of the integrated sensors, the gripper was attached to a robotic arm manipulator and tested on pipes of different outer diameters. Two sets of tests were carried out on both defect-free pipes and pipes with simulated defects, including surface partial cracking and corrosion. The inspection results indicated that transmitted and received signals could be acquired with an acceptable signal-to-noise ratio in the time domain. Moreover, the simulated defects could be successfully detected using the integrated robotic sensing system.

Materials Science↗

DOE OSTI2021

In Vivo Study of Osteochondral Defect Regeneration Using Innovative Composite Calcium Phosphate Biocement in a Sheep Model

This study aimed to clarify the therapeutic effect and regenerative potential of the novel, amino acids-enriched acellular biocement (CAL) based on calcium phosphate on osteochondral defects in sheep. Eighteen sheep were divided into three groups, the treated group (osteochondral defects filled with a CAL biomaterial), the treated group with a biocement without amino acids (C cement), and the untreated group (spontaneous healing). Cartilages of all three groups were compared with natural cartilage (negative control). After six months, sheep were evaluated by gross appearance, histological staining, immunohistochemical staining, histological scores, X-ray, micro-CT, and MRI. Treatment of osteochondral defects by CAL resulted in efficient articular cartilage regeneration, with a predominant structural and histological characteristic of hyaline cartilage, contrary to fibrocartilage, fibrous tissue or disordered mixed tissue on untreated defect (p < 0.001, modified O’Driscoll score). MRI results of treated defects showed well-integrated and regenerated cartilage with similar signal intensity, regularity of the articular surface, and cartilage thickness with respect to adjacent native cartilage. We have demonstrated that the use of new biocement represents an effective solution for the successful treatment of osteochondral defects in a sheep animal model, can induce an endogenous regeneration of cartilage in situ, and provides several benefits for the design of future therapies supporting osteochondral defect healing.

Kresakova, Lenka↗

DOE OSTI2025

Molecular aspect ratio effect on axial thermal transport in solution-spun carbon nanotube fibers

Neat, densely packed, and highly aligned carbon nanotube fibers (CNTFs) have appealing room-temperature axial thermal conductivity (⁠k⁠) and thermal diffusivity (⁠α⁠) for applications in lightweight heat spreading, flexible thermal connections, and thermoelectric active cooling. Although CNTFs are regularly produced from different input carbon nanotubes (CNTs), prior work has not quantified how the CNT molecular aspect ratio r (i.e., molecular length-to-diameter ratio) influences k and in well-aligned, packed CNTFs. Here, we perform self-heated steady-state and three-omega thermal measurements at room temperature on CNTF suspended in vacuum. Our results show that k increases from 150 to 380 W / mK for viscosity-averaged molecular aspect ratios increasing from r = 960 to 5600 and nanotube diameters of ~2 nm, which we attribute to the effects of thermal resistances between CNT bundles. CNTFs made with varying volume fraction Φ of constituent high-r and low-r CNT have properties that fall within or below the typical macroscopic rule-of-mixtures bounds. The thermal diffusivity α scales with k, leading to a sample-averaged volumetric heat capacity of 1.5 ± 0.3 MJ / m 3 K⁠. This work's findings that fibers made from longer CNT have larger k and α at room temperature motivate further investigation into thermal transport in solution-spun CNTF.

42 ENGINEERING↗

NASA NTRS1997

Regularization for Atmospheric Temperature Retrieval Problems

Passive remote sensing of the atmosphere is used to determine the atmospheric state. A radiometer measures microwave emissions from earth's atmosphere and surface. The radiance measured by the radiometer is proportional to the brightness temperature. This brightness temperature can be used to estimate atmospheric parameters such as temperature and water vapor content. These quantities are of primary importance for different applications in meteorology, oceanography, and geophysical sciences. Depending on the range in the electromagnetic spectrum being measured by the radiometer and the atmospheric quantities to be estimated, the retrieval or inverse problem of determining atmospheric parameters from brightness temperature might be linear or nonlinear. In most applications, the retrieval problem requires the inversion of a Fredholm integral equation of the first kind making this an ill-posed problem. The numerical solution of the retrieval problem requires the transformation of the continuous problem into a discrete problem. The ill-posedness of the continuous problem translates into ill-conditioning or ill-posedness of the discrete problem. Regularization methods are used to convert the ill-posed problem into a well-posed one. In this paper, we present some results of our work in applying different regularization techniques to atmospheric temperature retrievals using brightness temperatures measured with the SSM/T-1 sensor. Simulation results are presented which show the potential of these techniques to improve temperature retrievals. In particular, no statistical assumptions are needed and the algorithms were capable of correctly estimating the temperature profile corner at the tropopause independent of the initial guess.

Velez-Reyes, Miguel↗

DOE OSTI2021

Extracting Material Property Measurement Data from Scientific Articles

Machine learning-based prediction of material properties is often hampered by the lack of sufficiently large training datasets. The majority of such measurement data is embedded in scientific literature and the ability to automatically extract these data is essential to support the development of reliable property prediction methods. In this work, we describe a methodology for an automatic property extraction framework using material solubility as the target property. We create an annotated dataset containing tags for solubility-related entities using a combination of regular expressions and manual tagging. We then compare five entity recognition models leveraging both token-level and span-level architectures on the task of classifying solute names, solubility values, and solubility units. Additionally, we explore a novel pretraining approach that leverages automated chemical name and quantity extraction tools to generate large datasets that do not rely on intensive manual effort. Finally, we perform an analysis to identify the causes of classification errors.

Panapitiya, Gihan U.↗

NASA NTRS1983

Transformations of the perturbed two-body problem to unperturbed harmonic oscillators

Singular, nonlinear, and Liapunov unstable equations are made regular and linear through transformations that change the perturbed planar problem of two bodies into unperturbed and undamped harmonic oscillators with constant coefficients, so that the stable solution may be immediately written in terms of the new variables. The use of arbitrary and special functions for the transformations allows the systematic discussion of previously introduced and novel anomalies. For the case of the unperturbed two-body problem, it is proved that if transformations are power functions of the radial variable, only the eccentric and the true anomalies (with the corresponding transformations of the radial variable) will result in harmonic oscillators. The present method significantly reduces computation requirements in autonomous space operations.

Szebehely, V.↗

DOE OSTI2024

A $C^1$-Conforming Arbitrary-Order Two-Dimensional Virtual Element Method for the Fourth-Order Phase-Field Equation

We present a two-dimensional conforming virtual element method for the fourth-order phase-field equation. Our proposed numerical approach to the solution of this high-order phase-field (HOPF) equation relies on the design of an arbitrary-order accurate, virtual element space with $C^1$ global regularity. Such regularity is guaranteed by taking the values of the virtual element functions and their full gradient at the mesh vertices as degrees of freedom. Attaining high-order accuracy requires also edge polynomial moments of the trace of the virtual element functions and their normal derivatives. In this work, we detail the scheme construction, and prove its convergence by deriving error estimates in different norms. A set of representative test cases allows us to assess the behavior of the method.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Sampling frequency thresholds for the quantum advantage of the quantum approximate optimization algorithm

We compare the performance of the Quantum Approximate Optimization Algorithm (QAOA) with state-of-the-art classical solvers Gurobi and MQLib to solve the MaxCut problem on 3-regular graphs. We identify the minimum noiseless sampling frequency and depth p required for a quantum device to outperform classical algorithms. There is potential for quantum advantage on hundreds of qubits and moderate depth with a sampling frequency of 10 kHz. We observe, however, that classical heuristic solvers are capable of producing high-quality approximate solutions in linear time complexity. In order to match this quality for large graph sizes N, a quantum device must support depth p > 11. Additionally, multi-shot QAOA is not efficient on large graphs, indicating that QAOA p ≤ 11 does not scale with N. These results limit achieving quantum advantage for QAOA MaxCut on 3-regular graphs. Other problems, such as different graphs, weighted MaxCut, and 3-SAT, may be better suited for achieving quantum advantage on near-term quantum devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Space-Split Algorithm for Sensitivity Analysis of Discrete Chaotic Systems With Multidimensional Unstable Manifolds

Accurate approximations of the change of a system's output and its statistics with respect to the input are highly desired in computational dynamics. Ruelle's linear response theory provides breakthrough mathematical machinery for computing the linear response of chaotic dynamical systems. In this paper, we propose an algorithm for sensitivity analysis of discrete chaos with an arbitrary number of positive Lyapunov exponents. We combine the concept of perturbation space-splitting, which regularizes Ruelle's original expression, together with measure-based parameterization of the expanding subspace. We use these tools to rigorously derive trajectory-following recursive relations that converge exponentially fast, and construct a memory-efficient Monte Carlo scheme for derivatives of the output statistics. Thanks to the regularization and lack of simplifying assumptions on the system's behavior, our method is immune to the common problems of other popular methods such as the exploding tangent solutions and unphysical shadowing directions. Here, we provide a ready-to-use algorithm, analyze its complexity, and demonstrate several numerical examples of sensitivity computation using physically-inspired low-dimensional systems.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1988

Higher order sensitivity of solutions to convex programming problems without strict complementarity

Consideration is given to a family of convex programming problems which depend on a vector parameter. It is shown that the solutions of the problems and the associated Lagrange multipliers are arbitrarily many times directionally differentiable functions of the parameter, provided that the data of the problems are sufficiently regular. The characterizations of the respective derivatives are given.

Malanowski, Kazimierz↗

NASA NTRS1994

Algorithms for bilevel optimization

General multilevel nonlinear optimization problems arise in design of complex systems and can be used as a means of regularization for multi-criteria optimization problems. Here, for clarity in displaying our ideas, we restrict ourselves to general bi-level optimization problems, and we present two solution approaches. Both approaches use a trust-region globalization strategy, and they can be easily extended to handle the general multilevel problem. We make no convexity assumptions, but we do assume that the problem has a nondegenerate feasible set. We consider necessary optimality conditions for the bi-level problem formulations and discuss results that can be extended to obtain multilevel optimization formulations with constraints at each level.

Alexandrov, Natalia↗

NASA NTRS1990

The shallow structure of the Martian lithosphere in the vicinity of the ridged plains

On Mars, the most prominent assemblage of regularly spaced plain ridges occurs in the plains of the Coprates and Lunae Planum regions. In this paper, a number of tectonic models is developed and quantitatively evaluated, that relate the regular spacing of these ridges to the shallow internal structure of Mars. It is shown that models with both rigid and deformable megaregolith-basement interface conditions can provide solutions which can explain the ridge spacing within the constraint provided by the estimated thickness of the smooth plains materials. Implications of the models for the structure and state of stress in the Tharsis region at the time of ridge formation are discussed.

Zuber, M. T.↗

NASA NTRS1995

Spacecraft Availability Enhancement by In-flight Testing of Spare Parts

This paper describes an approach to improve availability by testing redundant parts at pre-determined intervals. The purpose of the testing is to detect nonfunctional back-up equipment and develop work-around measures or replacement spacecraft before a failure of the primary equipment reduces availability. The work reported here is an outgrowth of the NASA Space Network's program for the maintenance and replenishment of the Tracking and Data Relay Satellite System. The approach is based on a standby factor and a cyclic stress factor. The standby factor accounts for the effects of adverse storage conditions encountered as part of the in-flight environment. The stress factor accounts for the effects of physical or thermal cycling due to the application of force or power that characterizes the operation or use of a component. By the quantitative consideration of standby and cycling risks, a regular testing interval can be calculated - an interval for testing that furnishes information for availability planning but does not subject the spacecraft to undue risk. This paper includes quantitative solutions for appropriate testing intervals for equipment configurations and failure rates that are representative of a Tracking and Relay Data Satellite. The effect of a single test on the availability of certain equipment is also illustrated.

Cutler, Robert M.↗

NASA NTRS2003

Using CFD as Rocket Injector Design Tool: Recent Progress at Marshall Space Flight Center

The choice of tools used for injector design is in a transitional phase between exclusive reliance on the empirically based correlations and extensive use of computational fluid dynamics (CFD). The Next Generation Launch Technology (NGLT) Program goals emphasizing lower costs and increased reliability have produced a need to enable CFD as an injector design tool in a shorter time frame. This is the primary objective of the Staged Combustor Injector Technology Task currently under way at Marshall Space Flight Center (MSFC). The documentation of this effort begins with a very brief status of current injector design tools. MSFC's vision for use of CFD as a tool for combustion devices design is stated and discussed with emphasis on the injector. The concept of the Simulation Readiness Level (SRL), comprised of solution fidelity, robustness and accuracy, is introduced and discussed. This quantitative measurement is used to establish the gap between the current state of demonstrated capability and that necessary for regular use in the design process. MSFC's view of the validation process is presented and issues associated with obtaining the necessary data are noted and discussed. Three current experimental efforts aimed at generating validation data are presented. The importance of uncertainty analysis to understand the data quality is also demonstrated. First, a brief status of current injector design tools is provided as context for the current effort. Next, the MSFC vision for using CFD as an injector design tool is stated. A generic CFD-based injector design methodology is also outlined and briefly discussed. Three areas where MSFC is using injector CFD analyses for program support will be discussed. These include the Integrated Powerhead Development (IPD) engine which uses hydrogen and oxygen propellants in a full flow staged combustion (FFSC) cycle and the TR-107 and the RS84 engine both of which use RP-1 and oxygen in an ORSC cycle. Finally, an attempt is made to objectively summarize what progress has been made at MSFC in enabling CFD as an injector design tool.

Tucker, Kevin↗

NASA NTRS2009

Comparison of Node-Centered and Cell-Centered Unstructured Finite-Volume Discretizations: Viscous Fluxes - Part 1

Discretization of the viscous terms in current finite-volume unstructured-grid schemes are compared using node-centered and cell-centered approaches in two dimensions. Accuracy and efficiency are studied for six nominally second-order accurate schemes: a node-centered scheme, cell-centered node-averaging schemes with and without clipping, and cell-centered schemes with unweighted, weighted, and approximately mapped least-square face gradient reconstruction. The grids considered range from structured (regular) grids to irregular grids composed of arbitrary mixtures of triangles and quadrilaterals, including random perturbations of the grid points to bring out the worst possible behavior of the solution. Two classes of tests are considered. The first class of tests involves smooth manufactured solutions on both isotropic and highly anisotropic grids with discontinuous metrics, typical of those encountered in grid adaptation. The second class concerns solutions and grids varying strongly anisotropically over a curved body, typical of those encountered in high-Reynolds number turbulent flow simulations. Results from the first class indicate the face least-square methods, the node-averaging method without clipping, and the node-centered method demonstrate second-order convergence of discretization errors with very similar accuracies per degree of freedom. The second class of tests are more discriminating. The node-centered scheme is always second order with an accuracy and complexity in linearization comparable to the best of the cell-centered schemes. In comparison, the cell-centered node-averaging schemes are less accurate, have a higher complexity in linearization, and can fail to converge to the exact solution when clipping of the node-averaged values is used. The cell-centered schemes using least-square face gradient reconstruction have more compact stencils with a complexity similar to the complexity of the node-centered scheme. For simulations on highly anisotropic curved grids, the least-square methods have to be amended either by introducing a local mapping of the surface anisotropy or modifying the scheme stencil to reflect the direction of strong coupling.

Diskin, Boris↗

NASA NTRS2010

Comparison of Node-Centered and Cell-Centered Unstructured Finite-Volume Discretizations: Viscous Fluxes

Discretization of the viscous terms in current finite-volume unstructured-grid schemes are compared using node-centered and cell-centered approaches in two dimensions. Accuracy and complexity are studied for four nominally second-order accurate schemes: a node-centered scheme and three cell-centered schemes - a node-averaging scheme and two schemes with nearest-neighbor and adaptive compact stencils for least-square face gradient reconstruction. The grids considered range from structured (regular) grids to irregular grids composed of arbitrary mixtures of triangles and quadrilaterals, including random perturbations of the grid points to bring out the worst possible behavior of the solution. Two classes of tests are considered. The first class of tests involves smooth manufactured solutions on both isotropic and highly anisotropic grids with discontinuous metrics, typical of those encountered in grid adaptation. The second class concerns solutions and grids varying strongly anisotropically over a curved body, typical of those encountered in high-Reynolds number turbulent flow simulations. Tests from the first class indicate the face least-square methods, the node-averaging method without clipping, and the node-centered method demonstrate second-order convergence of discretization errors with very similar accuracies per degree of freedom. The tests of the second class are more discriminating. The node-centered scheme is always second order with an accuracy and complexity in linearization comparable to the best of the cell-centered schemes. In comparison, the cell-centered node-averaging schemes may degenerate on mixed grids, have a higher complexity in linearization, and can fail to converge to the exact solution when clipping of the node-averaged values is used. The cell-centered schemes using least-square face gradient reconstruction have more compact stencils with a complexity similar to that of the node-centered scheme. For simulations on highly anisotropic curved grids, the least-square methods have to be amended either by introducing a local mapping based on a distance function commonly available in practical schemes or modifying the scheme stencil to reflect the direction of strong coupling. The major conclusion is that accuracies of the node centered and the best cell-centered schemes are comparable at equivalent number of degrees of freedom.

Diskin, Boris↗