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.

SourceResearch connectionFrom yearThrough year

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

At least 847 records · Page 47

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Using Homogenized Models to Explore the Effect of Fracture Densities on Weathering

Despite its importance, only a few researchers have incorporated the effects of fracturing into models of reactive transport for rock weathering. Here we explore 2D simulations that describe weathering under conditions of diffusive and advective transport within heterogeneous media consisting of rocky blocks and fractures. In our simulations, the Darcy velocities vary in space and time and depend on weathering processes within the rock matrix. We explore simulations with saturated and unsaturated flow for weathering bedrock that consists of blocks separated by inert or weathered material. The simulations show that a simplified homogenized model can approximate exact solutions for some of the simulated columns and hills and can allow exploration of coupling between flow and reaction in fractured rock. These hillslope simulations document that, even in the presence of 2D water flow, i) an increase in fracture density results in faster weathering advance rates; and ii) the water table locates deeper for a rock system that is weathered and fractured rather than weathered and unfractured. Some of these patterns have also been observed for natural systems. But these simulations also highlight how simplified models that do not use appropriate averaging of heterogeneities can be inaccurate in predicting weathering rate for natural systems. For example, if water flows both vertically and laterally through the vadose zone of a hill, then a prediction of the depth of regolith that is based on modeling strictly unidirectional downward infiltration will be unrealistically large. Likewise, if the fracture density observed near the land surface is used in a model to predict depth of weathering for a system where the fracture density decreases downward, the model will overestimate regolith depth. Learning how to develop accurately homogenized models could thus enable better conceptual models and predictions of weathering advance in natural systems.

58 GEOSCIENCES↗

NASA NTRS1974

Evaluation of atomic constants for optical radiation, volume 1

Atomic constants for optical radiation are discussed which include transition probabilities, line strengths, and oscillator strengths for both dipole and quadrupole transitions, as well as the associated matrix elements needed for line broadening calculations. Atomic constants were computed for a wide selection of elements and lines. An existing computer program was used, with modifications to include, in an approximate manner, the effect of equivalent electrons, and to enable reordering and restructuring of the output for publication. This program is suitable for fast, low cost computation of the optical constants, using the Coulomb approximation formalism for LS coupling.

Kylstra, C. D.↗

NASA NTRS1983

CO 1-0 band isotopic lines as intensity standards

Fundamental band vibration-rotational line intensities are reported for the (C-12)(O-16), (C-13)(O-16), and (C-12)(O-18) molecules. The intensities were computed using numerically obtained matrix elements of a recently, determined electric dipole moment function. New intensity measurements for some of these lines obtained via high-resolution Fourier spectra are also reported. The overall consistency (approximately 2.2 percent) between calculations, which are based on recent experimental determinations of CO 1-0 band intensities, not including our own, those of the authors, and independent experimental results of the authors suggests that the lines of the fundamental band of CO can be used as absolute intensity standards.

Chackerian, C., Jr.↗

NASA NTRS1985

The substitution of nickel for cobalt in hot isostatically pressed powder metallurgy UDIMET 700 alloys

Nickel was substituted in various proportions for cobalt in a series of five hot-isostatically-pressed powder metallurgy alloys based on the UDIMET 700 composition. These alloys were given 5-step heat treatments appropriate for use in turbine engine disks. The resultant microstructures displayed three distinct sizes of gamma-prime particles in a gamma matrix. The higher cobalt-content alloys contained larger amounts of the finest gamma-prime particles, and had the lowest gamma-gamma-prime lattice mismatch. While all alloys had approximately the same tensile properties at 25 and 650 gamma C, the rupture lives at 650 and 760 C peaked in the alloys with cobalt contents between 12.7 and 4.3 pct. Minimum creep rates increased as cobalt contents were lowered, suggesting their correlation with the gamma-prime particle size distribution and the gamma-gamma-prime mismatch. It was also found that, on overaging at temperatures higher than suitable for turbine disk use, the high cobalt-content alloys were prone to sigma phase formation.

Harf, F. H.↗

NASA NTRS1986

Stability of semidiscrete approximations for hyperbolic initial-boundary-value problems: An eigenvalue analysis

A hyperbolic initial-boundary-value problem can be approximated by a system of ordinary differential equations (ODEs) by replacing the spatial derivatives by finite-difference approximations. The resulting system of ODEs is called a semidiscrete approximation. A complication is the fact that more boundary conditions are required for the spatially discrete approximation than are specified for the partial differential equation. Consequently, additional numerical boundary conditions are required and improper treatment of these additional conditions can lead to instability. For a linear initial-boundary-value problem (IBVP) with homogeneous analytical boundary conditions, the semidiscrete approximation results in a system of ODEs of the form du/dt = Au whose solution can be written as u(t) = exp(At)u(O). Lax-Richtmyer stability requires that the matrix norm of exp(At) be uniformly bounded for O less than or = t less than or = T independent of the spatial mesh size. Although the classical Lax-Richtmyer stability definition involves a conventional vector norm, there is no known algebraic test for the uniform boundedness of the matrix norm of exp(At) for hyperbolic IBVPs. An alternative but more complicated stability definition is used in the theory developed by Gustafsson, Kreiss, and Sundstrom (GKS). The two methods are compared.

Warming, Robert F.↗

NASA NTRS1989

1200 to 1400 K slow strain rate compressive properties of NiAl/Ni2AlTi-base materials

An attempt to apply the Martin Marietta Corporation's XD technology to the fabrication of NiAl-Ni2AlTi materials with improved creep properties is presented. Composite materials, containing from 0 to 30 vol pct of nominally 1-micron-diameter TiB2 particles in the intermetallic matrix have been produced by the XD process and compacted by hot pressing. Such composites demonstrated significant strength increases, approaching 3-fold for the 20 vol pct materials, in comparison to the unreinforced aluminide. This behavior was accomplished without deleterious side effects as the grain boundaries and particle-matrix interfaces were intact after compressive deformation to 10 percent or more strain. Typical true compressive stress-strain diagrams for materials tested in air between 1200 and 1400 K at approximate strain rates of 1.7 x 10 to the -6th/sec are presented.

Whittenberger, J. Daniel↗

NASA NTRS1990

Fresh look at floating shock fitting

A fast implicit upwind procedure for the two-dimensional Euler equations is described that allows accurate computations of shocked flows on nonadapted meshes. Away from shocks, the second-order accurate upwinding is based on the split-coefficient-matrix (SCM) method. In the presence of shocks, the difference stencils are modified using a floating shock fitting technique. Rapid convergence to steady-state solutions is attained with a diagonalized approximate factorization (AF) algorithm. Results are presented for Riemann's problem, for a regular shock reflection at an inviscid wall, for supersonic flow past a cylinder, and for a transonic airfoil. All computed shocks are ideally sharp and in excellent agreement with other numerical results or 'exact' solutions. Most importantly, this has been accomplished on unusually crude meshes without any attempt to align grid lines with shock fronts or to cluster grid lines around shocks.

Hartwich, PETER-M.↗

DOE OSTI2024

Luther-Emery liquid and dominant singlet superconductivity in the hole-doped Haldane spin-1 chain

We investigate the pairing tendencies in the hole-doped Haldane spin-1 chain. To allow for doping, we extend the original spin chain Hamiltonian into a fermionic model involving a two-orbital Hubbard chain at intermediate or strong repulsive interaction strengths U and for degenerate orbitals. At half filling and large U, the ferromagnetic Hund's coupling, J H , generates effective spin-1 moments, with antiferromagnetic correlations between sites. Using large-scale density matrix renormalization group calculations, we accurately study the system's behavior under light hole-doping. For U = 1.6 in units of the noninteracting bandwidth and for J H /U ≳ 0.275, we find that singlet pairing dominates the long-distance physics, establishing this system as a promising platform for repulsively mediated superconductivity. We provide concrete examples of materials that could realize the physics described here. We also provide evidence that the system approaches a Luther-Emery liquid state at large system sizes, reminiscent of the behavior of doped one-orbital two-leg ladders at weak coupling, which also have superconducting tendencies. The numerically calculated central charge approaches one in the thermodynamic limit, indicating a single gapless mode as is expected for the Luther-Emery state. Exponents characterizing the power-law decays of singlet pair-pair and charge density-density correlations are determined, and found to approximately satisfy the Luther-Emery identity.

1-dimensional systems↗

NASA NTRS1992

An approximately factored incremental strategy for calculating consistent discrete aerodynamic sensitivity derivatives

An incremental strategy is presented for iteratively solving very large systems of linear equations, which are associated with aerodynamic sensitivity derivatives for advanced CFD codes. It is shown that the left-hand side matrix operator and the well-known factorization algorithm used to solve the nonlinear flow equations can also be used to efficiently solve the linear sensitivity equations. Two airfoil problems are considered as an example: subsonic low Reynolds number laminar flow and transonic high Reynolds number turbulent flow.

Korivi, V. M.↗

NASA NTRS2014

Optics of Water Cloud Droplets Mixed with Black-Carbon Aerosols

We use the recently extended superposition T-matrix method to calculate scattering and absorption properties of micrometer-sized water droplets contaminated by black carbon. Our numerically exact results reveal that, depending on the mode of soot-water mixing, the soot specific absorption can vary by a factor exceeding 6.5. The specific absorption is maximized when the soot material is quasi-uniformly distributed throughout the droplet interior in the form of numerous small monomers. The range of mixing scenarios captured by our computations implies a wide range of remote sensing and radiation budget implications of the presence of black carbon in liquid-water clouds. We show that the popular Maxwell-Garnett effective-medium approximation can be used to calculate the optical cross sections, single-scattering albedo, and asymmetry parameter for the quasi-uniform mixing scenario, but is likely to fail in application to other mixing scenarios and in computations of the elements of the scattering matrix.

soot↗

DOE OSTI2025

Sensitivity of magnetic islands in permanent magnet stellarators using the gradient and Hessian methods

Stellarator plasmas are known to be very sensitive to perturbations in the magnetic field. The permanent magnet stellarator was in part developed as a solution to high machining tolerances placed on the shape properties of electromagnetic coils in traditional stellarators. However, as a consequence of this high sensitivity to the field structure, sensitivities of permanent magnet stellarator plasmas to perturbations of permanent magnet properties must necessarily be well-understood. The gradient and Hessian matrix methods have been previously demonstrated to be useful sensitivity analysis methods for modular coils. We apply these two methods to the study of island width sensitivities in both the MUSE and PM4STELL permanent magnet stellarator projects. These sensitivity methods were used to determine the relative impacts of permanent magnet parameter perturbations on island widths in the vacuum field approximation of both stellarator equilibria. The square of resonant magnetic field perturbation is used here as a proxy for island width. In particular, gradients of magnetizations of individual magnets were examined in MUSE, as well as gradients of magnet group displacements informed by device design. Three different forms of permanent magnet magnetization perturbations are investigated for MUSE, and the flux surface response to perturbations is demonstrated. The Hessian matrix method is applied to PM4STELL, illustrating the sensitivity of dominant island widths to displacements of toroidal wedge structures. These methods allow for selective direction of experimental resources toward regions of heightened sensitivity, while constraints on less impactful permanent magnet parameters can be relaxed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Machine Learning Surrogates of a Fuel Matrix Degradation Process Model for Performance Assessment of a Nuclear Waste Repository

Spent nuclear fuel repository simulations are currently not able to incorporate detailed fuel matrix degradation (FMD) process models due to their computational cost, especially when large numbers of waste packages breach. The current paper uses machine learning to develop artificial neural network and k-nearest neighbor regression surrogate models that approximate the detailed FMD process model while being computationally much faster to evaluate. Further, using fuel cask temperature, dose rate, and the environmental concentrations of CO 3 2- , O 2 , Fe 2+ , and H 2 as inputs, these surrogates show good agreement with the FMD process model predictions of the UO 2 degradation rate for conditions within the range of the training data. A demonstration in a full-scale shale repository reference case simulation shows that the incorporation of the surrogate models captures local and temporal environmental effects on fuel degradation rates while retaining good computational efficiency.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2023

Polyester-based epoxy vitrimer integrating spent coffee ground as a natural filler

Coffee is a widely consumed beverage, but approximately 80% mass of coffee cherry is discarded as waste. Here, in this study, spent coffee grounds (SCG) were recycled and used as a natural filler in a polyester-based epoxy matrix made from reactions of diglycidyl ethers of Bisphenol A and polyethylene glycol with methyl nadic anhydride. The properties of SCG/epoxy composites with various SCG contents up to 40% w/w were evaluated using various techniques such as FT-IR spectroscopy, tensile testing, optical microscopy, differential scanning calorimetry, and rheology. The SCG filled polyester-based epoxy matrix demonstrated excellent dimensional stability and retained a permanently crosslinked network, but re-shaping was possible at ~150 °C through transesterification, which occurred due to a dynamic chain exchange reaction between the ester group of the epoxy matrix and the hydroxyl group of the SCG fillers. The transesterification also improved the interfacial adhesion between the matrix and SCG fillers, resulting in tensile strength of over 20 MPa and a modulus of around 2.5 GPa, even with high SCG contents up to 40% w/w. While neat, cured epoxy shows thermoset behavior, filling the matrix with SCG enables reforming and origami-inspired spontaneous change in shape.

42 ENGINEERING↗

DOE OSTI2026

Data-Driven Performance Optimization of Gamma Spectrometers With Many Channels

In gamma spectrometers with variable spectroscopic performance across many channels (e.g., many pixels or voxels), a tradeoff exists between including data from successively worse-performing readout channels and increasing efficiency. Brute-force calculation of the optimal set of included channels is exponentially infeasible as the number of channels grows, and approximate methods are required. In this work, we present a data-driven framework for attempting to find near-optimal sets of included detector channels. The framework leverages non-negative matrix factorization (NMF) to learn the behavior of gamma spectra across the detector and clusters similarly-performing detector channels together. Performance comparisons are then made between spectra with channel clusters removed, which is more feasible than brute force. The framework is general and can be applied to arbitrary, user-defined performance metrics depending on the application. We apply this framework to optimizing gamma spectra measured by H3D M400 CdZnTe (CZT) spectrometers, which exhibit variable performance across their crystal volumes. In particular, we show several examples optimizing various performance metrics for uranium and plutonium gamma spectra in non-destructive assay (NDA) for nuclear safeguards, and explore trends in performance versus parameters such as clustering algorithm type. We also compare the NMF + clustering pipeline to several non-machine-learning (ML) algorithms, including several greedy algorithms. Although, we find that the NMF + clustering pipeline tends to find the best-performing set of detector voxels, significantly improving over the unoptimized spectra, but that a greedy accumulation of spectra segmented by detector depth can, in some cases, give similar performance improvements in much less computation time.

Energy resolution↗

NASA NTRS1979

Gust alleviation - Criteria and control laws

The relationships between criteria specified for aircraft gust alleviation and the form of the control laws that result from the criteria are considered. Open-loop gust alleviation based on the linearized, small perturbation equations of aircraft motion is discussed, and an approximate solution of the open-loop control law is presented for the case in which the number of degrees of freedom of the aircraft exceeds the rank of the control effectiveness matrix. Excessive actuator lag is compensated for by taking into account actuator dynamics in the equations of motion, resulting in the specification of a general load network. Criteria for gust alleviation when output motions are gust alleviated and the closed-loop control law derived from them are examined and linear optimal control law is derived. Comparisons of the control laws reveal that the effectiveness of an open-loop control law is greatest at low aircraft frequencies but deteriorates as the natural frequency of the actuators is approached, while closed-loop methods are found to be more effective at higher frequencies.

Rynaski, E. G.↗

NASA NTRS1980

Possible isotopic fractionation effects in sputtered minerals

A model which makes definite predictions for the fractionation of isotopes in sputtered material is discussed. The fractionation patterns are nonlinear, and the pattern for a particular set of isotopes depends on the chemical matrix within which those isotopes are contained. Calculations are presented for all nonmonoisotopic elements contained in the minerals perovskite, anorthite, ackermanite, enstatite, and troilite. All isotopes are fractionated at the level of approximately 4-6 deg/o per atomic mass unit. Oxygen is always positively fractionated (heavier isotopes sputtered preferentially), and heavier elements are generally negatively fractioned (light isotopes sputtered preferentially). The value of Delta (O-18:O-16) is always less by about 1.8 deg/o than a linear extrapolation based upon the calculated delta (O-17:O-16) value would suggest. The phenomenon of both negative and positive fractionation patterns from a single target mineral are used to make an experimental test of the proposed model.

Haff, P. K.↗

NASA NTRS1985

Deformation and failure of transversely loaded composite plates

Quasi-isotropic graphite/epoxy clamped circular plates were loaded under control point loading. Deflections and strains at various locations were monitored up to failure. Various damage mechanisms, such as intralaminar matrix cracking, delaminations and fiber breaks, are correlated with the strain records. The damage at various stages of deformation was characterized by means of X-radiography. Results are compared with an approximate isotropic plate theory.

Daniel, I. M.↗

NASA NTRS1990

A superlinear interior points algorithm for engineering design optimization

We present a quasi-Newton interior points algorithm for nonlinear constrained optimization. It is based on a general approach consisting of the iterative solution in the primal and dual spaces of the equalities in Karush-Kuhn-Tucker optimality conditions. This is done in such a way to have primal and dual feasibility at each iteration, which ensures satisfaction of those optimality conditions at the limit points. This approach is very strong and efficient, since at each iteration it only requires the solution of two linear systems with the same matrix, instead of quadratic programming subproblems. It is also particularly appropriate for engineering design optimization inasmuch at each iteration a feasible design is obtained. The present algorithm uses a quasi-Newton approximation of the second derivative of the Lagrangian function in order to have superlinear asymptotic convergence. We discuss theoretical aspects of the algorithm and its computer implementation.

Herskovits, J.↗