Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “matrix approximation”

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.

At least 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Towards excitations and dynamical quantities in correlated lattices with density matrix embedding theory

Density matrix embedding theory (DMET) provides a framework to describe ground-state expectation values in strongly correlated systems, but its extension to dynamical quantities is still an open problem. We show one route to obtaining excitations and dynamical spectral functions by using the techniques of DMET to approximate the matrix elements that arise in a single-mode inspired excitation ansatz. We demonstrate this approach in the one-dimensional Hubbard model, comparing the neutral excitations, single-particle density of states, charge, and spin dynamical structure factors to benchmarks from the Bethe ansatz and density matrix renormalization group. Finally, our work highlights the potential of these ideas in building computationally efficient approaches for dynamical quantities.

1-dimensional systems↗
NASA NTRS1996

Some Experiences with Nonoverlapping Schur Complement Parallel Preconditioning for CFD Calculations

In this work we consider solving matrices which arise from the discretization of advection-diffusion field equations on arbitrary triangulated domains using stabilized numerical methods. The talk will discuss several candidate matrix preconditioning algorithms based on the 2 x 2 block factorization induced by an apriori partitioning of the triangulated domain. Application of the 2 x 2 block preconditioner requires the formation and inversion of the Schur complement submatrix. We consider several strategies for simplifying this task: incomplete Schur complement factorizations, drop tolerance element filling, Schur complement probing, and localized Schur complement inversion. Numerical results will be shown comparing performance and efficiency of these approximations. The matrix preconditioner has also been embedded into a Newton algorithm for solving the nonlinear Euler and Navier-Stokes equations governing compressible flow. The remainder of the talk will show numerous examples in CFD to demonstrate the efficiency and robustness of the techniques.

Barth, Timothy J.↗
NASA NTRS1996

Characterization of the proteins comprising the integral matrix of Strongylocentrotus purpuratus embryonic spicules

In the present study, we enumerate and characterize the proteins that comprise the integral spicule matrix of the Strongylocentrotus purpuratus embryo. Two-dimensional gel electrophoresis of [35S]methionine radiolabeled spicule matrix proteins reveals that there are 12 strongly radiolabeled spicule matrix proteins and approximately three dozen less strongly radiolabeled spicule matrix proteins. The majority of the proteins have acidic isoelectric points; however, there are several spicule matrix proteins that have more alkaline isoelectric points. Western blotting analysis indicates that SM50 is the spicule matrix protein with the most alkaline isoelectric point. In addition, two distinct SM30 proteins are identified in embryonic spicules, and they have apparent molecular masses of approximately 43 and 46 kDa. Comparisons between embryonic spicule matrix proteins and adult spine integral matrix proteins suggest that the embryonic 43-kDa SM30 protein is an embryonic isoform of SM30. An adult 49-kDa spine matrix protein is also identified as a possible adult isoform of SM30. Analysis of the SM30 amino acid sequences indicates that a portion of SM30 proteins is very similar to the carbohydrate recognition domain of C-type lectin proteins.

Non-NASA Center↗
DOE OSTI2021

Hydrodynamics of nonintegrable systems from a relaxation-time approximation

We develop a general kinetic theory framework to describe the hydrodynamics of strongly interacting, nonequilibrium quantum systems in which integrability is weakly broken, leaving a few residual conserved quantities. This framework is based on a generalized relaxation-time approximation; it gives a simple, but surprisingly accurate, prescription for computing nonequilibrium transport even in strongly interacting systems. We validate the predictions of this approximation against matrix product operator calculations on chaotic quantum spin chains, finding surprisingly good agreement. Lastly, we show that despite its simplicity, our framework can capture phenomena distinctive to strongly interacting systems, such as widely separated charge and energy diffusion constants.

1-dimensional spin chains↗
NASA NTRS2024

Shape Servoing of Deformable Objects Using Model Estimation and Barrier Lyapunov Function

An adaptive shape servoing control method is presented in this article to manipulate a deformable object into a desired shape in 3-D. A finite-point-based representation of the deformable object is used and the deformation Jacobian matrix is approximated using Fourier series basis functions. The unknown parameters of the deformation Jacobian are learned by using the velocity applied to a control point on the object and corresponding change of positions of the points describing the entire object. An integral concurrent learning (ICL)-based parameter update law is designed along with a constrained controller to satisfy the state constraints on the motion of the control point using Barrier Lyapunov function analysis. ICL-based parameter update law uses data history of velocity and corresponding positions of the points along with their current values. An efficient algorithm to update the history stack using singular value maximization is proposed based on the structure of the regressor matrix. Simulations using a physical simulator and experiments using a robot platform are performed to validate the performance of the proposed controller on two different deformable objects.

Vrithik Raj Guthikonda↗
DOE OSTI2025

Convergence Analysis of the Alternating Anderson–Picard Method for Nonlinear Fixed-Point Problems

Anderson acceleration (AA) has been widely used to solve nonlinear fixed-point problems due to its rapid convergence. This work focuses on a variant of AA in which multiple Picard iterations are performed between each AA step, referred to as the Alternating Anderson–Picard (AAP) method. Furthermore, despite introducing more “slow” Picard iterations, this method has been shown to be efficient and even more robust in both linear and nonlinear cases. However, there is a lack of theoretical analysis for AAP in the nonlinear case. In this paper, we address this gap by establishing the equivalence between AAP and a multisecant-GMRES method that uses GMRES to solve a multisecant linear system at each iteration. From this perspective, we show that AAP “converges” to the Newton-GMRES method. Specifically, as the residual approaches zero, the multisecant matrix, the approximate Jacobian inverse, the search direction, and the optimization gain of AAP converge to their counterparts in the Newton-GMRES method. These connections provide insights for analyzing the asymptotic convergence properties of AAP. Consequently, we show that AAP is locally 𝑞-linear convergent and provide an upper bound for the convergence factor of AAP. To validate the theoretical results, numerical examples are provided.

Anderson acceleration↗
DOE OSTI2025

Degradable Biocomposite Thermoplastic Polyurethanes

In this project, the team developed tough and degradable biocomposite thermoplastic polyurethanes (TPUs) by incorporating bacterial spores into TPUs as a biofunctional living filler. The team screened various bacteria and selected the Bacillus subtilis ATCC 6633 strain as the final candidate, primarily due to its genomic availability, sporulation ability and TPU assimilation activity. The heat-shock tolerance of ATCC 6633 spores was further improved through evolutionary engineering via Adaptive Laboratory Evolution (ALE), demonstrating a 17.7-fold enhanced germination efficiency post heat-shock treatment compared to the wild-type strain (WT). The team fabricated biocomposite TPUs by incorporating lyophilized powder of heat-shock tolerized (HST) spores during the hot melt extrusion (HME) of TPU at 135 °C. The baseline TPU used in this project is a commercially available soft-grade TPU (BCF45) manufactured by BASF. Colony forming unit (CFU) assays quantified that WT and HST spores in the TPU matrix retained approximately 20% and 100% survivability, respectively, after HME. Tensile testing demonstrated that the spores behaved as a polymer-reinforcing filler, positively affecting the overall tensile properties of the biocomposite TPU. For example, biocomposite TPU with WT and HST spores (BC TPU WT and BC TPU HST , respectively) exhibited up to 25% and 37% improved toughness, respectively, compared to TPU without spores. BC TPU HST showed remarkably improved disintegration in autoclaved compost (92% mass loss in 5 months), which simulated a microbially poor environment for TPU degradation. When compared to TPU without spores (44% mass loss in 5 months) the acceleration of degradation is marked. Respirometry confirmed that 72% of BC TPU HST was biomineralized into CO2 within 6 months, indicating that spores in the biocomposite TPU were germinated by utilizing nutrients in the autoclaved compost, facilitating TPU degradation at the end of the material's life. The team demonstrated the scale-up of biocomposite TPU fabrication using continuous extrusion and injection molding techniques. Processing conditions optimized in a lab-scale microcompounder were successfully transferred to a continuous extruder with a 30-fold increased throughput. Biocomposite TPUs prepared using these industry-relevant processes showed comparable toughness improvements to samples prepared in the lab-scale extruder. Excitingly, following compounding in the pilot-extruder the composite material could be injection molded, while retaining high spore viability and similar toughness improvements. The team also found that spores in biocomposite TPU served as antioxidants, preventing toughness decay during the recycled extrusion of BC TPU HST . Long-term storage tests over one year showed that the addition of spores had no negative effect on the longevity of the TPU. Furthermore, the team demonstrated the fabrication of spore-bearing biocomposite polymers with other polyesters such as PBAT, PLA, and PCL. We obtained promising preliminary data that showed overall toughness improvements for all polymers with spore addition. Finally, life cycle assessment (LCA) and techno-economic analysis (TEA) were carried out, which indicated minimal additional cost of fabrication. Overall, a tough and degradable biocomposite thermoplastic was successfully developed through this project, with all tasks completed successfully, achieving >100% of the objectives.

36 MATERIALS SCIENCE↗
NASA NTRS1980

Dynamics of rotationally periodic large space structure

The analysis of large area rotationally periodic space structures presented in the paper combines the finite element method, transfer matrix procedures, approximation methods, and periodic structure analysis to obtain computational efficiency. The computations used in the analysis indicate that additive damping mechanisms can be evaluated from the frequency response of the structure. The transient response can also be obtained from the frequency response to complete the dynamic analysis.

Mcdaniel, T. J.↗
NASA NTRS1987

Robust eigensystem assignment for flexible structures

An improved method is developed for eigenvalues and eigenvectors placement of a closed-loop control system using either state or output feedback. The method basically consists of three steps. First, the singular value of QR decomposition is used to generate an orthonormal basis that spans admissible eigenvector space corresponding to each assigned eigenvalue. Secondly, given a unitary matrix, the eigenvector set which best approximates the given matrix in the least-square sense and still satisfy eigenvalue cosntraints is determined. Thirdly, a unitary matrix is sought to minimize the error between the unitary matrix and the assignable eigenvector matrix. For use as the desired eigenvector set, two matrices, namely, the open-loop eigenvector matrix and its closest unitary matrix are proposed. The latter matrix generally encourages both minimum conditioning and control gains. In addition, the algorithm is formulated in real arithmetic for efficient implementation. To illustrate the basic concepts, numerical examples are included.

Juang, Jer-Nan↗
NASA NTRS1988

A methodology for the synthesis of robust feedback systems

A new methodology is developed for the synthesis of linear, time-variant (LTI) controllers for multivariable LTI systems. The resulting closed-loop system is nominally stable and exhibits a known level of performance. In addition, robustness of the feedback system is guaranteed, i.e., stability and performance are retained in the presence of multiple unstructured uncertainty blocks located at various points in the feedback loop. The design technique is referred to as the Causality Recovery Methodology (CRM). The CRM relies on the Youla parameterization of all stabilizing compensators to ensure nominal stability of the feedback system. A frequency-domain inequality in terms of the structured singular value mu defines the robustness specification. The optimal compensator, with respect to the mu condition, is shown to be noncausal in general. The aim of the CRM is to find a stable, causal transfer function matrix that approximates the robustness characteristics of the optimal solution. The CRM, via a series of infinite-dimensional convex programs, produces a closed-loop system whose performance robustness is at least as good as that of any initial design. The algorithm is approximated by a finite dimensional process for the purposes of implementation. Two numerical examples confirm the potential viability of the CRM concept; however, the robustness improvement comes at the expense of increased computational burden and compensator complexity.

Milich, David Albert↗
NASA NTRS1994

Sensitivity derivatives for three dimensional supersonic Euler code using incremental iterative strategy

In a recent work, an incremental strategy was proposed to iteratively solve the very large systems of linear equations that are required to obtain quasianalytical sensitivity derivatives from advanced computational fluid dynamics (CFD) codes. The technique was sucessfully demonstrated for two large two-dimensional problems: a subsonic and a transonic airfoil. The principal feature of this incremental iterative stategy is that it allows the use of the identical approximate coefficient matrix operator and algorithm to solve the nonlinear flow and the linear sensitivity equations; at convergence, the accuracy of the sensitivity derivatives is not compromised. This feature allows a comparatively straightforward extension of the methodology to three-dimensional problems; this extension is successfully demonstrated in the present study for a space-marching solution of the three-dimensional Euler equations over a Mach 2.4 blended wing-body configuration.

Korivi, Vamshi Mohan↗
DOE OSTI2025

Analytical and EZmock covariance validation for the DESI 2024 results

The estimation of uncertainties in cosmological parameters is an important challenge in Large-Scale-Structure (LSS) analyses. For standard analyses such as Baryon Acoustic Oscillations (BAO) and Full-Shape two approaches are usually considered. First: analytical estimates of the covariance matrix use Gaussian approximations and (nonlinear) clustering measurements to estimate the matrix, which allows a relatively fast and computationally cheap way to generate matrices that adapt to an arbitrary clustering measurement. On the other hand, sample covariances are an empirical estimate of the matrix based on an ensemble of clustering measurements from fast and approximate simulations. While more computationally expensive due to the large amount of simulations and volume required, these allow us to take into account systematics that are impossible to model analytically. In this work we compare these two approaches in order to enable DESI's key analyses. We find that the configuration space analytical estimate performs satisfactorily in BAO analyses and its flexibility in terms of input clustering makes it the fiducial choice for DESI's 2024 BAO analysis. On the contrary, the analytical computation of the covariance matrix in Fourier space does not reproduce the expected measurements in terms of Full-Shape analyses, which motivates the use of a corrected mock covariance for DESI's 2024 Full Shape analysis.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1995

A Reduced Dimension Static, Linearized Kalman Filter and Smoother

An approximate Kalman filter and smoother, based on approximations of the state estimation error covariance matrix, is described. Approximations include a reduction of the effective state dimension, use of a static asymptotic error limit, and a time-invariant linearization of the dynamic model for error integration. The approximations lead to dramatic computational savings in applying estimation theory to large complex systems. Examples of use come from TOPEX/POSEIDON.

Kalman filter smoother modeling ocean simulation↗
NASA NTRS1987

The matrix exponential in transient structural analysis

The primary usefulness of the presented theory is in the ability to represent the effects of high frequency linear response with accuracy, without requiring very small time steps in the analysis of dynamic response. The matrix exponential contains a series approximation to the dynamic model. However, unlike the usual analysis procedure which truncates the high frequency response, the approximation in the exponential matrix solution is in the time domain. By truncating the series solution to the matrix exponential short, the solution is made inaccurate after a certain time. Yet, up to that time the solution is extremely accurate, including all high frequency effects. By taking finite time increments, the exponential matrix solution can compute the response very accurately. Use of the exponential matrix in structural dynamics is demonstrated by simulating the free vibration response of multi degree of freedom models of cantilever beams.

Minnetyan, Levon↗
NASA NTRS2004

Polymer/Silicate Nanocomposites Used to Manufacture Gas Storage Tanks With Reduced Permeability

Over the past decade, there has been considerable research in the area of polymer-layered silicate nanocomposites. This research has shown that the dispersion of small amounts of an organically modified layered silicate improves the polymer strength, modulus, thermal stability, and barrier properties. There have been several reports on the dispersion of layered silicates in an epoxy matrix. Potential enhancements to the barrier properties of epoxy/silicate nanocomposites make this material attractive for low permeability tankage. Polymer matrix composites (PMCs) have several advantages for cryogenic storage tanks. They are lightweight, strong, and stiff; therefore, a smaller fraction of a vehicle's potential payload capacity is used for propellant storage. Unfortunately, the resins typically used to make PMC tanks have higher gas permeability than metals. This can lead to hydrogen loss through the body of the tank instead of just at welds and fittings. One approach to eliminate this problem is to build composite tanks with thin metal liners. However, although these tanks provide good permeability performance, they suffer from a substantial mismatch in the coefficient of thermal expansion, which can lead to failure of the bond between the liner and the body of the tank. Both problems could be addressed with polymersilicate nanocomposites, which exhibit reduced hydrogen permeability, making them potential candidates for linerless PMC tanks. Through collaboration with Northrop Grumman and Michigan State University, nanocomposite test tanks were manufactured for the NASA Glenn Research Center, and the helium permeability was measured. An organically modified silicate was prepared at Michigan State University and dispersed in an epoxy matrix (EPON 826/JeffamineD230). The epoxy/silicate nanocomposites contained either 0 or 5 wt% of the organically modified silicate. The tanks were made by filament winding carbon fibers with the nanocomposite resin. Helium permeability was measured by Northrop Grumman, showing that the leak rate/day of the nanocomposite matrix tank was approximately 80-percent less than that of the neat epoxy matrix tank.

Campbell, Sandi G.↗
DOE OSTI2025

The no boundary density matrix

We discuss a no-boundary proposal for a subregion of the universe. In the classical approximation, this density matrix involves finding a specific classical solution of the equations of motion with no boundary. Beyond the usual no boundary condition at early times, we also have another no boundary condition in the region we trace out. We can find the prescription by starting from the usual Hartle-Hawking proposal for the wavefunction on a full slice and tracing out the unobserved region in the classical approximation. We discuss some specific subregions and compute the corresponding solutions. These geometries lead to phenomenologically unacceptable probabilities, as expected.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗