Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “generalized algorithm”

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 973 records · Page 54

Lattice Green’s Functions for High-Order Finite Difference Stencils

Lattice Green's Functions (LGFs) are fundamental solutions to discretized linear operators, and as such they are a useful tool for solving discretized elliptic PDEs on domains that are unbounded in one or more directions. The majority of existing numerical solvers that make use of LGFs rely on a second-order discretization and operate on domains with free-space boundary conditions in all directions. Under these conditions, fast expansion methods are available that enable precomputation of 2D or 3D LGFs in linear time, avoiding the need for brute-force multi-dimensional quadrature of numerically unstable integrals. Here we focus on higher-order discretizations of the Laplace operator on domains with more general boundary conditions, by (1) providing an algorithm for fast and accurate evaluation of the LGFs associated with high-order dimension-split centered finite differences on unbounded domains, and (2) deriving closed-form expressions for the LGFs associated with both dimension-split and Mehrstellen discretizations on domains with one unbounded dimension. Through numerical experiments we demonstrate that these techniques provide LGF evaluations with near machine-precision accuracy, and that the resulting LGFs allow for numerically consistent solutions to high-order discretizations of the Poisson's equation on fully or partially unbounded 3D domains.

97 MATHEMATICS AND COMPUTING↗

Improve Learning from Crowds via Generative Augmentation

Crowdsourcing provides an efficient label collection schema for supervised machine learning. However, to control annotation cost, each instance in the crowdsourced data is typically annotated by a small number of annotators. This creates a sparsity issue and limits the quality of machine learning models trained on such data. In this paper, we study how to handle sparsity in crowdsourced data using data augmentation. Specifically, we propose to directly learn a classifier by augmenting the raw sparse annotations. We implement two principles of high-quality augmentation using Generative Adversarial Networks: 1) the generated annotations should follow the distribution of authentic ones, which is measured by a discriminator; 2) the generated annotations should have high mutual information with the ground-truth labels, which is measured by an auxiliary network. Extensive experiments and comparisons against an array of state-of-the-art learning from crowds methods on three real-world datasets proved the effectiveness of our data augmentation framework. It shows the potential of our algorithm for low-budget crowdsourcing in general.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗

Clustering Algorithm for AM Parts using GSH and EDT with Autoencoder

SAND2025-10103O The Clustering Algorithm for AM Parts Using GSH and (EDT With Autoencoder is a software tool. It uses a clustering algorithm for additive manufacturing (AM) parts using generalized spherical harmonics (GSH) and Euclidean distance transform (EDT) with an autoencoder to quantify material microstructure. The tool offers improved sensitivity to microstructural changes compared to traditional approaches. The tool integrates multiple microstructural properties, such as grain morphology, crystallographic orientation, and material phase information, to provide a comprehensive analysis of material microstructures. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Rodgers, Theron [Sandia National Lab. (SNL-CA), Li↗

Parallel Time Integration for Constrained Optimization

The number of transistors in an average processor continues to increase, but individual clock speeds have plateaued. Those transistors are instead going into additional cores, increasing the number of different things that a processor can do at once and placing an emphasis on parallel computation. Many problems in scientific computing follow a time-evolution model, and it can be difficult to solve such problems in parallel across the temporal domain. The Multi-Grid Reduction In Time (MGRIT) algorithm, developed at Lawrence Livermore National Laboratory (LLNL), solves differential equations with a method designed specifically to take advantage of extreme numbers of processors by parallelizing across time. The Tri-diagonal MGRIT (TriMGRIT) algorithm, also developed at LLNL, is a generalization of MGRIT which enables parallel-in-time solving of a greater number of problems. Constrained optimization problems, in particular, may be solved in parallel using TriMGRIT. These consist of choosing a control function such that an objective functional is minimized, constrained by a differential-equation. We consider two such problems: applying torque to a pendulum to bring it to a gentle stop and moving a crowd of people from one distribution into another. We also perform some miscellaneous theoretical and practical research, including investigating the use of a line-search subroutine to refine intermediate TriMGRIT results and preliminary work on strategies for choosing operators for TriMGRIT to use.

97 MATHEMATICS AND COMPUTING↗

Toward Accelerated Nuclear-physics Parameter Estimation from Binary Neutron Star Mergers: Emulators for the Tolman–Oppenheimer–Volkoff Equations

Abstract Gravitational-wave observations of binary neutron-star (BNS) mergers have the potential to revolutionize our understanding of the nuclear equation of state (EOS) and the fundamental interactions that determine its properties. However, Bayesian parameter estimation frameworks do not typically sample over microscopic nuclear-physics parameters that determine the EOS. One of the major hurdles in doing so is the computational cost involved in solving the neutron-star structure equations, known as the Tolman–Oppenheimer–Volkoff (TOV) equations. In this paper, we explore approaches to emulating solutions for the TOV equations: multilayer perceptrons (MLPs), Gaussian processes, and a data-driven variant of the reduced basis method (RBM). We implement these emulators for three different parameterizations of the nuclear EOS, each with a different degree of complexity represented by the number of model parameters. We find that our MLP-based emulators are generally more accurate than the other two algorithms, whereas the RBM results in the largest speedup with respect to the full high-fidelity TOV solver. We employ these emulators for a simple parameter inference using a potentially loud BNS observation and show that the posteriors predicted by our emulators are in excellent agreement with those obtained from the full TOV solver.

79 ASTRONOMY AND ASTROPHYSICS↗

The MSFC UNIVAC 1108 EXEC 8 simulation model

A model is presented which simulates the MSFC Univac 1108 multiprocessor system. The hardware/operating system is described to enable a good statistical measurement of the system behavior. The performance of the 1108 is evaluated by performing twenty-four different experiments designed to locate system bottlenecks and also to test the sensitivity of system throughput with respect to perturbation of the various Exec 8 scheduling algorithms. The model is implemented in the general purpose system simulation language and the techniques described can be used to assist in the design, development, and evaluation of multiprocessor systems.

Williams, T. G.↗

Design of structures for optimum geometry

A general method is proposed for geometric optimization of finite element structures on the assumption that a reasonable initial geometry is specified. The number of finite elements, number of joints, and the element-joint relationships are specified and are not changed during the optimization process. The structure may be statically indeterminate and may support multiple loading conditions. Design variables include geometric and member sizing parameters. Design objectives may be minimum weight or cost, and constraints include strength and stiffness limitations. The optimization procedure is a general mathematical programming approach. An iterative design algorithm is outlined, and the method is specialized for truss structures. Design examples are given.

Vanderplaats, G. N.↗

Some aspects of the inversion problem in remote sensing

Several commonly used methods for inversion--constrained linear inversion, synthesis (Backus-Gilbert) methods and nonlinear iterative techniques for the Chahine type--are discussed. It is demonstrated that a very close connection exists between Backus-Gilbert solutions and those given by constrained linear inversion. A number of examples of the application of such methods are presented, showing that resolution is not greatly different for quite different algorithms, a result quite in accord with general theoretical considerations: more resolution can be achieved at the expense of introducing greater a priori bias in the procedure.

Twomey, S.↗

Mathematical description and program documentation for CLASSY, an adaptive maximum likelihood clustering method

Discussed in this report is the clustering algorithm CLASSY, including detailed descriptions of its general structure and mathematical background and of the various major subroutines. The report provides a development of the logic and equations used with specific reference to program variables. Some comments on timing and proposed optimization techniques are included.

Lennington, R. K.↗

An algorithm for finite element analysis of partly wrinkled membranes

A generalization and numerical implementation is presented of the Stein-Hedgepeth (1961) continuum theory for the analysis of partly wrinkled membranes, such as the tensioned membrane surfaces which are found in spacecraft structural components. The approach has its basis in experimental observations showing that when wrinkles develop in a membrane parallel to the x direction, the associated overall contraction in the y direction exceeds that predicted by the Poisson's ratio effect.

Miller, R. K.↗

Validation of a transonic analysis code for use in preliminary design of advanced transport configurations

The WBPPW code has the capability of analyzing flow-field effects about configurations which include wing pylons and engine nacelles or pods in addition to the basic wing/fuselage combination. Using the concept of grid embedding, the code solves the extended small disturbance transonic flow equation for complex flow interactions of the various configuration components. A general description of the code and solution algorithm is included. Results are presented and compared with experiment for various configurations which encompass the code capabilities. These include wing planform and wing contour modifications and variations in nacelle position beneath a high-aspect-ratio wing. Results are analyzed in the light of preliminary design, where the capability to accurately compute flow-field effects resulting from various configuration perturbations is important. The comparisons show that the computational results are sensitive to subtle design modifications and that the code could be used as an effective guide during the design process for transport configurations.

Waggoner, E. G.↗

Theoretical and software considerations for general dynamic analysis using multilevel substructured models

The dynamic analysis of complex structural systems using the finite element method and multilevel substructured models is presented. The fixed-interface method is selected for substructure reduction because of its efficiency, accuracy, and adaptability to restart and reanalysis. This method is extended to reduction of substructures which are themselves composed of reduced substructures. The implementation and performance of the method in a general purpose software system is emphasized. Solution algorithms consistent with the chosen data structures are presented. It is demonstrated that successful finite element software requires the use of software executives to supplement the algorithmic language. The complexity of the implementation of restart and reanalysis porcedures illustrates the need for executive systems to support the noncomputational aspects of the software. It is shown that significant computational efficiencies can be achieved through proper use of substructuring and reduction technbiques without sacrificing solution accuracy. The restart and reanalysis capabilities and the flexible procedures for multilevel substructured modeling gives economical yet accurate analyses of complex structural systems.

Schmidt, R. J.↗

Constitutive modeling of cyclic plasticity and creep, using an internal time concept

Using the concept of an internal time as related to plastic strains, a differential stress-strain relation for elastoplasticity is rederived, such that (1) the concept of a yield-surface is retained; (2) the definitions of elastic and plastic processes are analogous to those in classical plasticity theory; and (3) its computational implementation, via a 'tangent-stiffness' finite element method and a 'generalized-midpoint-radial-return' stress-integration algorithm, is simple and efficient. Also, using the concept of an internal time, as related to both the inelastic strains as well as the Newtonian time, a constitutive model for creep-plasticity interaction, is discussed. The problem of modeling experimental data for plasticity and creep, by the present analytical relations, as accurately as desired, is discussed. Numerical examples which illustrate the validity of the present relations are presented for the cases of cyclic plasticity and creep.

Watanabe, O.↗

Embedding of attitude determination in n-dimensional spaces

The problem of attitude determination in n-dimensional spaces is addressed. The proper parameters are found, and it is shown that not all three-dimensional methods have useful extensions to higher dimensions. It is demonstrated that Rodriguez parameters are conveniently extendable to other dimensions. An algorithm for using these parameters in the general n-dimensional case is developed and tested with a four-dimensional example. The correct mathematical description of angular velocities is addressed, showing that angular velocity in n dimensions cannot be represented by a vector but rather by a tensor of the second rank. Only in three dimensions can the angular velocity be described by a vector.

Bar-Itzhack, Itzhack Y.↗

Hypercube matrix computation task

A major objective of the Hypercube Matrix Computation effort at the Jet Propulsion Laboratory (JPL) is to investigate the applicability of a parallel computing architecture to the solution of large-scale electromagnetic scattering problems. Three scattering analysis codes are being implemented and assessed on a JPL/California Institute of Technology (Caltech) Mark 3 Hypercube. The codes, which utilize different underlying algorithms, give a means of evaluating the general applicability of this parallel architecture. The three analysis codes being implemented are a frequency domain method of moments code, a time domain finite difference code, and a frequency domain finite elements code. These analysis capabilities are being integrated into an electromagnetics interactive analysis workstation which can serve as a design tool for the construction of antennas and other radiating or scattering structures. The first two years of work on the Hypercube Matrix Computation effort is summarized. It includes both new developments and results as well as work previously reported in the Hypercube Matrix Computation Task: Final Report for 1986 to 1987 (JPL Publication 87-18).

Calalo, Ruel H.↗

Rocket noise filtering system using digital filters

A set of digital filters is designed to filter rocket noise to various bandwidths. The filters are designed to have constant group delay and are implemented in software on a general purpose computer. The Parks-McClellan algorithm is used. Preliminary tests are performed to verify the design and implementation. An analog filter which was previously employed is also simulated.

Mauritzen, David↗

Numerical Arc-Segmentation Algorithm For A Radio Conference

NASARC computer program developed from general planning principles and decisions of both sessions of World Administrative Radio Conference on Use of Geostationary Satellite Orbit and on Planning of Space Services Utilizing It (WARC-85 and WARC-88). Written to help countries satisfy requirements for nationwide fixed-satellite services from at least one orbital position within predetermined arc. Written in ANSI standard FORTRAN 77.

Whyte, W. A.↗