Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear embedding”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

A globally well-posed finite element algorithm for aerodynamics applications

A finite element CFD algorithm is developed for Euler and Navier-Stokes aerodynamic applications. For the linear basis, the resultant approximation is at least second-order-accurate in time and space for synergistic use of three procedures: (1) a Taylor weak statement, which provides for derivation of companion conservation law systems with embedded dispersion-error control mechanisms; (2) a stiffly stable second-order-accurate implicit Rosenbrock-Runge-Kutta temporal algorithm; and (3) a matrix tensor product factorization that permits efficient numerical linear algebra handling of the terminal large-matrix statement. Thorough analyses are presented regarding well-posed boundary conditions for inviscid and viscous flow specifications. Numerical solutions are generated and compared for critical evaluation of quasi-one- and two-dimensional Euler and Navier-Stokes benchmark test problems.

Iannelli, G. S.↗
NASA NTRS1978

Residual stresses in angleplied laminates and their effects on laminate behavior

Evidence of the presence of lamination residual stresses in angleplied laminates were transply cracks and warpage of unsymmetric laminates which occur prior to application of any mechanical load. Lamination residual strains were measured using the embedded strain gage technique. These strains result from the temperature differences between cure and room temperature and vary linearly within this temperature range. Lamination residual stresses were usually present in angleplied fiber composites laminates; they were also present in unidirectional hybrids and superhybrids. For specific applications, the magnitudes of lamination residual stresses were determined and evaluated relative to the anticipated applied stresses. Particular attention was given to cyclic thermal loadings in applications where the thermal cycling takes place over a wide temperature range.

Chamis, C. C.↗
NASA NTRS1993

Resonant frequencies of irregularly shaped microstrip antennas using method of moments

This paper describes an application of the method of moments to determine resonant frequencies of irregularly shaped microstrip patches embedded in a grounded dielectric slab. For analysis, the microstrip patch is assumed to be excited by a linearly polarized plane wave that is normal to the patch. The surface-current density that is induced on the patch because of the incident field is expressed in terms of subdomain functions by dividing the patch into identical rectangular subdomains. The amplitudes of the subdomain functions, as a function of frequency, are determined using the electric-field integral equation (EFIE) approach in conjunction with the method of moments. The resonant frequencies of the patch are then obtained by selecting the frequency at which the amplitude of the surface-current density is real. The resonant frequencies of the equilateral triangular and other nonrectangular patches are computed using the present technique, and these frequencies are compared with measurements and other independent calculations.

Deshpande, Manohar D.↗
DOE OSTI2020

Track Seeding and Labelling with Embedded-space Graph Neural Networks

To address the unprecedented scale of HL-LHC data, the Exa.TrkX project is investigating a variety of machine learning approaches to particle track reconstruction. The most promising of these solutions, graph neural networks (GNN), process the event as a graph that connects track measurements (detector hits corresponding to nodes) with candidate line segments between the hits (corresponding to edges). Detector information can be associated with nodes and edges, enabling a GNN to propagate the embedded parameters around the graph and predict node-, edge- and graph-level observables. Previously, message-passing GNNs have shown success in predicting doublet likelihood, and we here report updates on the state-of-the-art architectures for this task. In addition, the Exa.TrkX project has investigated innovations in both graph construction, and embedded representations, in an effort to achieve fully learned end-to-end track finding. Hence, we present a suite of extensions to the original model, with encouraging results for hitgraph classification. In addition, we explore increased performance by constructing graphs from learned representations which contain non-linear metric structure, allowing for efficient clustering and neighborhood queries of data points. We demonstrate how this framework fits in with both traditional clustering pipelines, and GNN approaches. The embedded graphs feed into high-accuracy doublet and triplet classifiers, or can be used as an end-to-end track classifier by clustering in an embedded space. A set of post-processing methods improve performance with knowledge of the detector physics. Finally, we present numerical results on the TrackML particle tracking challenge dataset, where our framework shows favorable results in both seeding and track finding.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS1994

Langmuir turbulence in the auroral ionosphere 1: Linear theory

Intense bursts of Langmuir waves with electric fields of 50 to 500 mV / m have been frequently observed at altitudes greater than 500 km in the auroral ionosphere. These bursts are driven by 20 eV to 4 keV field-aligned electrons, which are embedded in an approximately isotropic nonthermal tail of scattered electrons. The Langmuir bursts are often observed at altitudes where the ionosphere is moderately magnetized (OMEGA (sub e) approximately equals omega (sub pe)). Both the moderate magnetization and the scattered electrons have a major influence on the linear dispersion and damping of Langmuir waves. In particular, the linear dispersion is topologically different depending on whether the magnetic field is subcritical (OMEGA (sub e) less than omega (sub pe)) or supercritical (OMEGA (sub e) greater than omega (sub pe)). The correct dispersion and damping can account for the observed polarization of the Langmuir waves, which is very nearly parallel to the geomagnetic field. Inferred properties of the linear instability driven by the field-aligned electrons are discussed. The linear dispersion and damping derived here provide the basis for a nonlinear turbulence study described in a companion paper (Newman et al., this issue).

Newman, D. L.↗
DOE OSTI2025

The revised elastic field of an edge dislocation

We modify the traditional linear-elastic field of an edge dislocation. The modifications stem from symmetry and energy requirements imposed in terms of embedded (deformed) coordinates. These requirements are satisfied if a line force is added to the dislocation field. The modification to the stress field is expressed by coefficients that are a function only of Poisson’s ratio. Qualitatively, the field of the dislocation is increased in the glide direction and decreased in the climb direction. There are small changes to the field of a screw dislocation that do not entail the addition of a line force. Effects of nonlinearity, anisotropic elasticity, and image forces are briefly considered.

Hirth, J. P.↗
NASA NTRS1985

Modeling the effects of latitudinal gradients in stellar winds, with application to the solar wind

A steady, axisymmetric, quasi-radial, global model previously developed for stellar winds with embedded magnetic fields has been extended to include latitudinal gradient effects on the azimuthal velocity and magnetic field. The linear results at large radii are presented for large-amplitude latitudinal variations in the radial magnetic field, mass loss rate, and radial velocity of the wind. The magnetohydrodynamic (MHD) equations predict meridional flows that develop naturally from internal magnetic stresses. The flows open flux tubes in the star's equatorial plane, redistributing mass and magnetic flux as a function of stellar latitude. The plasma spins up to conserve angular momentum in fields and plasma. The results are generally applicable to stellar winds (including radiatively driven winds), provided that the internal structure is not dominated by rotation. The asymptotic solutions do not explicitly depend on the form of the energy equation, although the assumed O(1) state which drives these solutions depends on the deposition of energy and momentum throughout the wind.

Nerney, S.↗
NASA NTRS1990

Preliminary thermal/structural analysis of a carbon-carbon/refractory-metal heat-pipe-cooled wing leading edge

This study presents preliminary thermal/structural analyses of a carbon-carbon/refractory-metal heat-pipe-cooled wing leading edge concept designed for an air breathing single-stage-to-orbit hypersonic vehicle. The concept features chordwise (i.e., normal to the leading edge) and spanwise (i.e., parallel to the leading edge) refractory-metal heat pipes which are completely embedded within a carbon-carbon primary structure. Studies of the leading edge were performed using nonlinear thermal and linear structural three-dimensional finite element analyses. The concept was shown to be thermally feasible within the limits of the assumptions made in the analyses when internal radiative cooling is present during ascent, and a three-dimensional carbon-carbon architecture is used. In addition, internal radiative cooling was found not to be necessary during descent. The linear stress analysis indicated excessively large thermal stresses in the rafractory metal walls of the heat pipes even though a soft layer of carbon was included between the heat pipe and the carbon-carbon structure in an attempt to reduce the thermal stresses. A nonlinear structural analysis may be necessary to properly model the response of the refractory-metal heat pipes.

Glass, David E.↗
NASA NTRS2012

Free-Free Transitions in the Presence of Laser Fields and Debye Potential at Very Low Incident Electron Energies

We study the free-free transition in electron-helium ion in the ground state and embedded in a Debye potential in the presence of an external laser field at very low incident electron energies. The laser field is treated classically while the collision dynamics is treated quantum mechanically. The laser field is chosen as monochromatic, linearly polarized and homogeneous. The incident electron is considered to be dressed by the laser field in a nonperturbative manner by choosing Volkov wave function for it. The scattering wave function for the incident electron on the target embedded in a Debye potential is solved numerically by taking into account the effect of electron exchange. We calculate the laser-assisted differential and total cross sections for free-free transition for absorption/emission of a single photon or no photon exchange. The results will be presented at the conference.

Bhatia, Anand↗
DOE OSTI2022

Star–Galaxy Image Separation with Computationally Efficient Gaussian Process Classification

Abstract We introduce a novel method for discerning optical telescope images of stars from those of galaxies using Gaussian processes (GPs). Although applications of GPs often struggle in high-dimensional data modalities such as optical image classification, we show that a low-dimensional embedding of images into a metric space defined by the principal components of the data suffices to produce high-quality predictions from real large-scale survey data. We develop a novel method of GP classification hyperparameter training that scales approximately linearly in the number of image observations, which allows for application of GP models to large-size Hyper Suprime-Cam Subaru Strategic Program data. In our experiments, we evaluate the performance of a principal component analysis embedded GP predictive model against other machine-learning algorithms, including a convolutional neural network and an image photometric morphology discriminator. Our analysis shows that our methods compare favorably with current methods in optical image classification while producing posterior distributions from the GP regression that can be used to quantify object classification uncertainty. We further describe how classification uncertainty can be used to efficiently parse large-scale survey imaging data to produce high-confidence object catalogs.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1986

Models of polarized infrared emission from bipolar nebulae

Many stars with circumstellar dust shells show a high degree of linear polarization (Sato et al. 1985). We are developing a model which assumes that the polarization arises from scattering by circumstellar dust. Our model assumes a geometry in which the star is surrounded by an optically thin spherical dust shell and embedded within an optically thick disk. This geometry is consistent with that proposed for objects with bipolar molecular outflow. This is important because many bipolar flow objects have also been observed to be highly polarized. The high degree of linear polarization is produced because the disk differentially attenuates the light from the star. The light incident from the point source is attenuated by a factor of exp(-tau/cos theta) where theta is the angle between a ray from the point source to the scatterer and a ray normal to the disk; tau is the optical depth at the wavelength of interest. Hence, the light scattered from the regions directly above and below the disk give the largest contribution to the total flux. The scattering angle for light from these regions is near 90 deg., so the light is strongly polarized and, in the Rayleigh scattering regime, is polarized parallel to the disk. The Stokes parameters for the scattered light from each particle in the shell are calculated by using the scattering matrix elements generated by a Mie scattering program. After the Stokes parameters for each particle are computed they are summed to give the Stokes parameters for the entire shell. Two graphs are presented which show the intensity and polarization spectrum generated by our model using the optical constants for astronomical silicates as defined by Draine and Lee (1984).

Burns, M. S.↗
DOE OSTI2021

Hydro-chemo-mechanical phase field formulation for corrosion induced cracking in reinforced concrete

Steel corrosion induced cracking in reinforced concrete structures is often caused by volumetric straining at the steel-concrete interface. This volumetric strain is mainly controlled by the rate of generation of corrosion product, which depends on the diffusion of Cl{sup −} and dissolved O{sub 2} as well as the chemical reaction at the steel-concrete interface. While many numerical modeling techniques have been proposed to tackle this complex phenomenon, most have been limited to using staggered solution schemes, where the diffusion, chemical reaction, mechanical strain, and crack fields are solved for separately. Unfortunately, such approaches introduce numerical errors that can lead to unrealistic predictions. Therefore, a coupled solution, is needed to remove the formulation induced errors and understand this complex phenomenon. To provide a unifying solution framework for corrosion induced cracking in reinforced concrete, a phase field formulation is presented and implemented with FEM algorithms. The governing equations are derived from the system's free energy and Fick's law of diffusion. The proposed formulation was first validated in 1D by comparing the modeling results with the semi-analytical solution approximated via Laurent series to avoid singularities in the complex domain. 2D plane strain solutions to the embedded rebar undergoing corrosion induced cracking with pre-cracks are also presented and compared with the approximated solutions from linear elastic fracture mechanics theory (LEFM). The modeling of crack propagation enabled by the proposed formulation is also presented along with parametric studies to reveal the roles of the fracture toughness and permeability of the concrete on the corrosion induced cracking. A case study was also conducted on H-pile steel with concrete jackets to demonstrate the feasibility of modeling corrosion induced fracture in reinforced concrete structures with complex steel-concrete interfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Nonlinear convergence in contact mechanics: Immersed boundary finite volume

In this report we present an immersed boundary finite volume (IBM) method for simulating quasistatic contact mechanics of linearly elastic domains at small strains. In IBM, all external boundaries and internal contacts of an object are represented by embedded surfaces inside a Cartesian mesh, which need not conform to the grid lines. The contact constraints consist of the non-penetrability condition and Coulomb’s friction law, which are discretized using special interpolation stencils and enforced via penalty parameters. The resulting nonlinear system depends on displacement unknowns only. To solve it, we use the Newton method but find that it diverges frequently. To understand the divergence pattern, we analyze a simplified 2-cell problem and show that the global convergence of Newton cannot be ensured for any choice of penalty parameters. We thus propose a modified Newton solver, which guarantees convergence for the 2-cell problem and is numerically verified to converge for all the challenging simulations considered herein. While both 1 st - and 2 nd -order variants of IBM, in displacement unknowns, are proposed, the modified Newton solver applies only to the 1 st -order variant.

42 ENGINEERING↗
NASA NTRS1985

An efficient iteration strategy for the solution of the Euler equations

A line Gauss-Seidel (LGS) relaxation algorithm in conjunction with a one-parameter family of upwind discretizations of the Euler equations in two-dimensions is described. The basic algorithm has the property that convergence to the steady-state is quadratic for fully supersonic flows and linear otherwise. This is in contrast to the block ADI methods (either central or upwind differenced) and the upwind biased relaxation schemes, all of which converge linearly, independent of the flow regime. Moreover, the algorithm presented here is easily enhanced to detect regions of subsonic flow embedded in supersonic flow. This allows marching by lines in the supersonic regions, converging each line quadratically, and iterating in the subsonic regions, thus yielding a very efficient iteration strategy. Numerical results are presented for two-dimensional supersonic and transonic flows containing both oblique and normal shock waves which confirm the efficiency of the iteration strategy.

Walters, R. W.↗
NASA NTRS1994

Discrete-Layer Piezoelectric Plate and Shell Models for Active Tip-Clearance Control

The objectives of this work were to develop computational tools for the analysis of active-sensory composite structures with added or embedded piezoelectric layers. The targeted application for this class of smart composite laminates and the analytical development is the accomplishment of active tip-clearance control in turbomachinery components. Two distinct theories and analytical models were developed and explored under this contract: (1) a discrete-layer plate theory and corresponding computational models, and (2) a three dimensional general discrete-layer element generated in curvilinear coordinates for modeling laminated composite piezoelectric shells. Both models were developed from the complete electromechanical constitutive relations of piezoelectric materials, and incorporate both displacements and potentials as state variables. This report describes the development and results of these models. The discrete-layer theories imply that the displacement field and electrostatic potential through-the-thickness of the laminate are described over an individual layer rather than as a smeared function over the thickness of the entire plate or shell thickness. This is especially crucial for composites with embedded piezoelectric layers, as the actuating and sensing elements within these layers are poorly represented by effective or smeared properties. Linear Lagrange interpolation polynomials were used to describe the through-thickness laminate behavior. Both analytic and finite element approximations were used in the plane or surface of the structure. In this context, theoretical developments are presented for the discrete-layer plate theory, the discrete-layer shell theory, and the formulation of an exact solution for simply-supported piezoelectric plates. Finally, evaluations and results from a number of separate examples are presented for the static and dynamic analysis of the plate geometry. Comparisons between the different approaches are provided when possible, and initial conclusions regarding the accuracy and limitations of these models are given.

Heyliger, P. R.↗
DOE OSTI2023

Polarization consistent dielectric screening in polarizable continuum model calculations of solvation energies

A polarization consistent framework, where dielectric screening is affected consistently in polarizable continuum model (PCM) calculations, is employed for the study of solvation energies. The computational framework combines a screened range-separated-hybrid functional (SRSH) with PCM calculations, SRSH-PCM, where dielectric screening is imposed in both PCM self-consistent reaction field (SCRF) iterations and the electronic structure Hamiltonian. We begin by demonstrating the impact of modifying the Hamiltonian to include such dielectric screening in SCRF iterations by considering the solutions of electrostatically embedded Hartree–Fock (HF) exact exchange equations. Long-range screened HF-PCM calculations are shown to capture properly the linear dependence of gap energy of frontier orbitals on the inverse of the dielectric constant, whereas unscreened HF-PCM orbital energies are fallaciously semi-constant with respect to the dielectric constant and, therefore, inconsistent with the ionization energy gaps. Similar trends affect density functional theory (DFT) calculations that aim to achieve predictive quality. Importantly, the dielectric screened calculations are shown to significantly affect DFT- and HF PCM-based solvation energies, where screened solvation energies are smaller compared to the unscreened values. Importantly, SRSH-PCM, therefore, appears to reduce the tendency of DFT-PCM to overestimate solvation energies, where we find the effect to increase with the dielectric constant and the polarity of the molecular solute, trends that enhance the quality of DFT-PCM calculations of solvation energy. Understanding the relationship of dielectric screening in the Hamiltonian and DFT-PCM calculations can ultimately benefit on-going efforts for the design of predictive and parameter free descriptions of solvation energies.

Chemistry↗
NASA NTRS1989

Integrated tools for control-system analysis

The basic functions embedded within a user friendly software package (MATRIXx) are used to provide a high level systems approach to the analysis of linear control systems. Various control system analysis configurations are assembled automatically to minimize the amount of work by the user. Interactive decision making is incorporated via menu options and at selected points, such as in the plotting section, by inputting data. There are five evaluations such as the singular value robustness test, singular value loop transfer frequency response, Bode frequency response, steady-state covariance analysis, and closed-loop eigenvalues. Another section describes time response simulations. A time response for random white noise disturbance is available. The configurations and key equations used for each type of analysis, the restrictions that apply, the type of data required, and an example problem are described. One approach for integrating the design and analysis tools is also presented.

Ostroff, Aaron J.↗
NASA NTRS2016

Probabilistic Modeling Of Ocular Biomechanics In VIIP: Risk Stratification

Visual Impairment and Intracranial Pressure (VIIP) syndrome is a major health concern for long-duration space missions. Currently, it is thought that a cephalad fluid shift in microgravity causes elevated intracranial pressure (ICP) that is transmitted along the optic nerve sheath (ONS). We hypothesize that this in turn leads to alteration and remodeling of connective tissue in the posterior eye which impacts vision. Finite element (FE) analysis is a powerful tool for examining the effects of mechanical loads in complex geometries. Our goal is to build a FE analysis framework to understand the response of the lamina cribrosa and optic nerve head to elevations in ICP in VIIP. To simulate the effects of different pressures on tissues in the posterior eye, we developed a geometric model of the posterior eye and optic nerve sheath and used a Latin hypercubepartial rank correlation coef-ficient (LHSPRCC) approach to assess the influence of uncertainty in our input parameters (i.e. pressures and material properties) on the peak strains within the retina, lamina cribrosa and optic nerve. The LHSPRCC approach was repeated for three relevant ICP ranges, corresponding to upright and supine posture on earth, and microgravity [1]. At each ICP condition we used intraocular pressure (IOP) and mean arterial pressure (MAP) measurements of in-flight astronauts provided by Lifetime Surveillance of Astronaut Health Program, NASA Johnson Space Center. The lamina cribrosa, optic nerve, retinal vessel and retina were modeled as linear-elastic materials, while other tissues were modeled as a Mooney-Rivlin solid (representing ground substance, stiffness parameter c1) with embedded collagen fibers (stiffness parameters c3, c4 and c5). Geometry creationmesh generation was done in Gmsh [2], while FEBio was used for all FE simulations [3]. The LHSPRCC approach resulted in correlation coefficients in the range of 1. To assess the relative influence of the uncertainty in an input parameter on the peak strains, we ranked and then normalized these coefficients, considering that normalized values 0.5 implied a substantial influence on the range of the peak strains in the optic nerve head (ONH). IOP and ICP were found to have a major influence on the peak strains in the ONH, as did optic nerve and LC stiffness. Interestingly, the stiffness of the sclera far from the scleral canal did not have a large influence on peak strains in ONH tissues; however, the collagen fiber stiffness in the peripapillary sclera and annular ring both influenced the peak strains within the ONH. We have created a physiologically relevant model that incorporated collagen fibers to study the effects of elevated ICP. Elevated ICP resulted in strains in the optic nerve that are not predicted to occur on earth: the upright or supine conditions. We found that IOP, ICP, lamina cribrosa stiffness and optic nerve stiffness had the highest association with these extreme strains in the ONH. These extreme strains may activate mechanosensitive cells that induce tissue remodeling and are a risk factor for the development of VIIP.

biomechanics↗