SEARCH · Engineering Papers

Results for “differential geometry”

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 577 records · Page 32

All sources NASA NTRS DOE OSTI

DOE OSTI2020

VALVE TIMING AND AREA PROFILE SELECTION FOR HYDRAULIC PUMPS AND MOTORS

A significant amount of research has been conducted to select valve timing and area profiles that create efficient and quiet hydraulic pumps and motors. Numerous active valve architectures have been modeled and optimized, but the rationale for the final solution is often unclear. The solution is usually highly dependent on the modeled valve geometry constraints and the duty cycle of the pump or motor for which the valve was optimized. This paper first presents an ideal valve area profile and what compromises are made to design for a working pressure and speed range. The ideal valve area profile is resolved using a piston-cylinder model with a pressure dependent bulk modulus and variable valve openings. The piston-cylinder is first modeled as a pump with passive (check) valves and a specified pressure drop across the valves to determine the ideal area profile. A generalization of the ideal timing is demonstrated for positive pressure differential and positive rotation direction, also known as the first quadrant. The effects on valve timing when varying pressure are shown for the first quadrant, but the trends can be extrapolated for all quadrants of operation. The ideal valve profile is implemented as fixed valve timing, meaning the valve area is only a function of the timing angle, in a pump-motor. The pump-motor is simulated in one rotation direction through a pressure range. Performance is high in pumping operation, but when the pressure differential is reversed, cylinder pressure spikes ensue. Two strategies to modify an ideal valve area profile are presented: timing grooves and a pressure shifted valve timing. Timing grooves reduce pressure spikes and cavitation in the cylinder but generally increase throttling losses. A pressure shifted valve timing has lower throttling energy losses, making it the favored solution.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2025

Cartesian equivariant representations for learning and understanding molecular orbitals

Qualitative and quantitative orbital properties such as bonding/antibonding character, localization, and orbital energies are critical to how chemists understand reactivity, catalysis, and excited-state behavior. Despite this, representations of orbitals in deep learning models have been very underdeveloped relative to representations of molecular geometries and Hamiltonians. Here, we apply state-of-the-art equivariant deep learning architectures to the task of assigning global labels to orbitals, namely energies characterizations, given the molecular coefficients from Hartree–Fock or density functional theory. The architecture we have developed, the Cartesian Equivariant Orbital Network (CEONET), shows how molecular orbital coefficients are readily featurized as equivariant node features common to all graph-based machine-learned potentials. We find that CEONET performs well at predicting difficult quantitative labels such as the orbital energy and orbital entropy. Furthermore, we find that the CEONET representation provides an intuitive latent space for differentiating orbital character for the qualitative assignment of e.g. bonding or antibonding character. In addition to providing a useful representation for further integrating deep learning with electronic structure theory, we expect CEONET to be useful for automatizing and interpreting the results of advanced electronic structure methods such as complete active space self-consistent field theory. In particular, the ability of CEONET to infer multireference character via the orbital entropy paves the way toward the machine-learned selection of active spaces.

chemical reactions↗

DOE OSTI2024

AMReX v2024

The software framework, AMReX, supports the development of block-structured adaptive mesh refinement (AMR) algorithms for solving systems of partial differential equations. AMR reduces the computational cost and memory footprint compared to a uniform mesh while preserving the essential local descriptions of different physical processes in complex multiphysics algorithms. AMR uses a hierarchical representation of the solution at multiple levels of resolution where the solution on each level is defined on the union of data containers at that resolution. These data containers, which represent the solution over a logically rectangular subregion of the domain, can contain field data defined on a mesh, Lagrangian particles or combinations of both. In addition to these basic data types, AMReX supports a multilevel embedded boundary representation of complex geometry; linear solvers for cell-centered and nodal data; asynchronous I/O in a native format readable by ParaView, VisIt and yt; and interfaces to hypre and PETSc solvers. AMReX enables applications to run on distributed memory architectures with multicore CPUs and with GPU accelerators. AMReX uses a lightweight abstraction layer that effectively hides the details of the architecture from the application. The framework currently supports CUDA, HIP and SYCL for GPU acceleration and OpenMP for multi-core CPU architectures.

Almgren, Ann↗

NASA NTRS1983

Destabilization of a 650 km chemical boundary layer and its bearing on the evolution of the continental crust

A mechanism for the production of a chemical change in the mantle, from primordial silicate compositions above the 650 km discontinuity to differentiated compositions below, is reviewed. Some consequences of this are the stabilization of two layer convection with a temperature contrast between the anhydrous mantle solidus and the geotherm which, at 650 km depth, is lower than any other location in the mantle. With thermal contributions from the concentration of the heat producing elements U, Th and K below the 650 km mantle boundary layer and the higher geotherms in the past, widespread or catastrophic melting may have taken place at this location. An episodic breach of this boundary layer by extensive heat and mass transport may have periodically destroyed any simple two-layer convection geometry in the mantle. Such episodic injections of mass and energy into the upper mantle from below may have been the mechanism responsible for episodes of enhanced surface tectonism and thermal activity which appear to be recorded in apparent polar wandering paths and radiometric ages of continental rocks.

Herzberg, C. T.↗

DOE OSTI2021

Fracture Network Localization Preceding Catastrophic Failure in Triaxial Compression Experiments on Rocks

We quantify the spatial distribution of fracture networks throughout six in situ X-ray tomography triaxial compression experiments on crystalline rocks at confining stresses of 5–35 MPa in order to quantify how fracture development controls the final macroscopic failure of the rock, a process analogous to those that control geohazards such as earthquakes and landslides. Tracking the proportion of the cumulative volume of fractures with volumes >90th percentile to the total fracture volume, ∑ v 90 / v t o t indicates that the fracture networks tend to increase in localization toward these largest fractures for up to 80% of the applied differential stress. The evolution of this metric also matches the evolution of the Gini coefficient, which measures the deviation of a population from uniformity. These results are consistent with observations of localizing low magnitude seismicity before large earthquakes in southern California. In both this analysis and the present work, phases of delocalization interrupt the general increase in localization preceding catastrophic failure, indicating that delocalization does not necessarily indicate a reduction of seismic hazard. However, the proportion of the maximum fracture volume to the total fracture volume does not increase monotonically. Experiments with higher confining stress tend to experience greater localization. To further quantify localization, we compare the geometry of the largest fractures, with volumes >90th percentile, to the best fit plane through these fractures immediately preceding failure. The r 2 scores and the mean distance of the fractures to the plane indicate greater localization in monzonite than in granite. The smaller mean mineral diameter and lower confining stress in the granite experiments may contribute to this result. Tracking these various metrics of localization reveals a close association between macroscopic yielding and the acceleration of fracture network localization. Near yielding, ∑ v 90 / v t o t and the Gini coefficient increase while the mean distance to the final failure plane decreases. Macroscopic yielding thus occurs when the rate of fracture network localization increases.

58 GEOSCIENCES↗

NASA NTRS1988

Nonconforming mortar element methods: Application to spectral discretizations

Spectral element methods are p-type weighted residual techniques for partial differential equations that combine the generality of finite element methods with the accuracy of spectral methods. Presented here is a new nonconforming discretization which greatly improves the flexibility of the spectral element approach as regards automatic mesh generation and non-propagating local mesh refinement. The method is based on the introduction of an auxiliary mortar trace space, and constitutes a new approach to discretization-driven domain decomposition characterized by a clean decoupling of the local, structure-preserving residual evaluations and the transmission of boundary and continuity conditions. The flexibility of the mortar method is illustrated by several nonconforming adaptive Navier-Stokes calculations in complex geometry.

Maday, Yvon↗

NASA NTRS2018

Simulation of Charge Collection to Spacecraft Surfaces: Freja Satellite

As spacecraft travel through plasma, spacecraft surfaces become charged by the collection of charged particles. This process is referred to as Surface Charging. These charges can be detrimental to the vehicle's electronic subsystems as they present a threat of electrostatic discharge (ESD) to onboard circuitry. The process of Surface Charging is complex and is affected by many elements. The charging of each surface is unique. The potential of an individual surface is dependent upon many variables including but not limited to the surface's geometry, material and its location. Each surface also has unique interactions with the surrounding plasma. Other factors that play large roles in the charging process is the density and temperature of plasma ions and electrons. Using Nascap-2k, a model of the Freja satellite has been constructed, and its auroral plasma environment has been imitated to simulate surface charging characteristics. The charging process of the Freja satellite has been modeled iteratively with incremental changes in both the Maxwellian electron temperature (eV) as well as the Gaussian electron energy (eV). This study provides an analysis of the sensitivity between spacecraft surface charging and these two primary variables of electron differential flux.

Wolfley, Olivia H.↗

NASA NTRS2018

Simulation of Charge Collection to Spacecraft Surfaces: Freja Satellite

As spacecraft travel through space plasma, spacecraft surfaces become charged by the collection of charged particles. This process is referred to as Surface Charging. These charges can be detrimental to the vehicle's electronic subsystems as they present a threat of electrostatic discharge (ESD) to onboard circuitry. The process of Surface Charging is complex and is affected by many elements. The charging of each surface is unique. The potential of an individual surface is dependent upon many variables including but not limited to the surface's geometry, material and its location. Each surface also has unique interactions with the surrounding plasma. Other factors that play large roles in the charging process is the density and temperature of plasma ions and electrons. Using Nascap-2k, a model of the Freja satellite has been constructed, and its auroral plasma environment has been imitated to simulate surface charging characteristics. The charging process of the Freja satellite has been modeled iteratively with incremental changes in both the Maxwellian electron temperature (eV) as well as the Gaussian electron energy (eV). This study provides an analysis of the sensitivity between spacecraft surface charging and these two primary variables of electron differential flux.

Wolfley, Olivia Hope↗

NASA NTRS1988

Numerical simulation of axisymmetric turbulent flow in combustors and diffusors

A method for predicting turbulent flow in combustors and diffusers is developed. The Navier-Stokes equations, incorporating a turbulence kappa-epsilon model equation, were solved in a nonorthogonal curvilinear coordinate system. The solution applied the finite volume method to discretize the differential equations and utilized the SIMPLE algorithm iteratively to solve the differenced equations. A zonal grid method, wherein the flow field was divided into several subsections, was developed. This approach permitted different computational schemes to be used in the various zones. In addition, grid generation was made a more simple task. However, treatment of the zonal boundaries required special handling. Boundary overlap and interpolating techniques were used and an adjustment of the flow variables was required to assure conservation of mass, momentum and energy fluxes. The numerical accuracy was assessed using different finite differencing methods, i.e., hybrid, quadratic upwind and skew upwind, to represent the convection terms. Flows in different geometries of combustors and diffusers were simulated and results compared with experimental data and good agreement was obtained.

Yung, Chain Nan↗

NASA NTRSDate not supplied

Separation and Transition on the ROTEX-T Cone-Flare

As part of NATO STO AVT-346 “Predicting Hypersonic Boundary-Layer Transition on Complex Geometries,” coordinated experimental and computational studies were conducted on the ROTEX-T, a cone-flare geometry used in a successful flight-test experiment. At the as-flown conditions, a separation bubble existed at the compression corner. Separation, reattachment, and the multifaceted linear instability paths leading this bubble to transition to turbulence are challenging to predict, but have significant impact on surface pressure and heat flux. High-resolution background-oriented schlieren and infrared thermography measurements were made in the AFOSR–Notre Dame Large Mach-6 Quiet Tunnel at freestream unit Reynolds numbers from 5.8×10(exp 6) to 12.2×10(exp 6) m-1and nominally zero angle of attack. High-speed self-aligned focusing schlieren, infrared thermography, and focused laser differential interferometry measurements were made in the AFRL Mach-6 Ludwieg Tube from 2.2×10(exp 6) to 24.7×10(exp 6) m-1and nominally zero degrees angle of attack. The surface heat-flux and Stanton number distributions were computed. Separation and reattachment locations, as well as the flow state at each, were determined from the combination of surface and off-wall measurements. The convective and global boundary-layer instabilities of the axisymmetric laminar flow at the experimental conditions were investigated computationally. Amplification of Mack’s first and second modes were observed to have logarithmic amplification factors between 5 to 7.5 at the separation location, depending on conditions. The flow was found to be globally unstable to stationary three-dimensional disturbances concentrated in the reattachment region. Previous analysis of the ROTEX-T flight data had not assessed reattachment location or the flow state upon reattachment. Thanks to the insights gained from the coordinated, on- and off-wall ground-test measurements, these evaluations have now been made. The separation location indicated by laminar simulations is consistently numerically predicted to be upstream of the experimentally observed location for a transitional separation bubble. The cause of this difference is understood to lie in the steady-state and axisymmetric assumptions made by both solvers employed to compute the basic states analyzed, as flow topology considerations assert that unsteady two-dimensional or axisymmetric separation bubbles are structurally unstable and will become three-dimensional. Computed laminar heating rates prior to separation agreed well with experiment; transitional heating rates after reattachment were between laminar and turbulent computations.

CFD Validation↗

NASA NTRSDate not supplied

Separation and Transition on the ROTEX-T Cone-Flare

CFD Validation↗

NASA NTRSDate not supplied

Evaluating Crystallinity in Thermoplastic Composites for Aerospace Applications

Polymer matrix composites (PMCs) offer many benefits for the aerospace industry due to their potential for weight reduction when compared to metal or ceramic based materials. Most PMCs currently in flight use thermoset matrices, however, thermoplastic resins are being explored as alternatives due to their ability to be remelted, which is of particular interest due to the potential for in-situ repair and faster production. Most thermoplastic resins are semicrystalline polymers. The properties of semicrystalline thermoplastics are directly influenced by their crystallinity, which can vary due to many factors including thermal treatments, environmental conditions, and mechanical deformation. Within thermoplastic PMC parts, crystallinity gradients can arise due to variations in part geometry, across part thicknesses, and along bonded joints. Monitoring the crystallinity of thermoplastic composites is key to ensuring these materials meet the high demands required for aerospace. Several different analytical techniques exist that can be used to characterize the bulk crystallinity of thermoplastic materials. However, many existing methods lack the specificity required to identify the subtle variations in crystallinity that may play a significant role in the performance and durability of PMC parts. Because of this, a significant amount of work is still required to fully characterize and understand the crystallinity profiles of thermoplastic PMCs and the resulting impact to material properties. This talk discusses the use of multiple techniques such as Differential Scanning Calorimetry, Polarized Light Optical Microscopy, and Fourier-Transform Infrared Spectroscopy to characterize the crystallinity in carbon fiber/thermoplastic composites. Samples of different crystallinity profiles were manufactured using various cooling procedures. This work aims to provide the fundamental data necessary to understand the effects of crystallinity on thermoplastic PMCs, which is key to advancing their use in aerospace applications.

Thermoplastics↗

DOE OSTI2025

Solving sparse finite element problems on neuromorphic hardware

The finite element method (FEM) is one of the most important and ubiquitous numerical methods for solving partial differential equations (PDEs) on computers for scientific and engineering discovery. Applying the FEM to larger and more detailed scientific models has driven advances in high-performance computing for decades. Here we demonstrate that scalable spiking neuromorphic hardware can directly implement the FEM by constructing a spiking neural network that solves the large, sparse, linear systems of equations at the core of the FEM. We show that for the Poisson equation, a fundamental PDE in science and engineering, our neural circuit achieves meaningful levels of numerical accuracy and close to ideal scaling on modern, inherently parallel and energy-efficient neuromorphic hardware, specifically Intel’s Loihi 2 neuromorphic platform. We illustrate extensions to irregular mesh geometries in both two and three dimensions as well as other PDEs such as linear elasticity. Our spiking neural network is constructed from a recurrent network model of the brain’s motor cortex and, in contrast to black-box deep artificial neural network-based methods for PDEs, directly translates the well-understood and trusted mathematics of the FEM to a natively spiking neuromorphic algorithm.

Applied mathematics↗

NASA NTRS1971

Wind-driven currents in a shallow lake or sea

For shallow lakes and seas such as the great lakes (especially Lake Erie) where the depth is not much greater than the Ekman depth, the usual Ekman dynamics cannot be used to predict the wind driven currents. The necessary extension to include shallow bodies of water, given by Welander, leads to a partial differential equation for the surface displacement which in turn determines all other flow quantities. A technique for obtaining exact analytical solutions to Welander's equation for bodies of water with large class of bottom topographies which may or may not contain islands is given. It involves applying conformal mapping methods to an extension of Welander's equation into the complex plane. When the wind stress is constant (which is the usual assumption for lakes) the method leads to general solutions which hold for bodies of water of arbitrary shape (the shape appears in the solutions through a set of constants which are the coefficients in the Laurent expansion of a mapping of the particular lake geometry). The method is applied to an elliptically shaped lake and a circular lake containing an eccentrically located circular island.

Goldstein, M. E.↗

NASA NTRS1996

Chemistry by Way of Density Functional Theory

In this work we demonstrate that density functional theory (DFT) methods make an important contribution to understanding chemical systems and are an important additional method for the computational chemist. We report calibration calculations obtained with different functionals for the 55 G2 molecules to justify our selection of the B3LYP functional. We show that accurate geometries and vibrational frequencies obtained at the B3LYP level can be combined with traditional methods to simplify the calculation of accurate heats of formation. We illustrate the application of the B3LYP approach to a variety of chemical problems from the vibrational frequencies of polycyclic aromatic hydrocarbons to transition metal systems. We show that the B3LYP method typically performs better than the MP2 method at a significantly lower computational cost. Thus the B3LYP method allows us to extend our studies to much larger systems while maintaining a high degree of accuracy. We show that for transition metal systems, the B3LYP bond energies are typically of sufficient accuracy that they can be used to explain experimental trends and even differentiate between different experimental values. We show that for boron clusters the B3LYP energetics are not as good as for many of the other systems presented, but even in this case the B3LYP approach is able to help understand the experimental trends.

Bauschlicher, Charles W., Jr.↗

DOE OSTI2020

Detection, monitoring, and determination of location of changes in metallic structures using multimode acoustic signals

Methods for detection, monitoring, and determination of location of changes in rigid structures with arbitrarily complex geometries are described. Implementations include locating acoustic transducers that generate and receive acoustic signals at multiple locations along a surface of the rigid structure, wherein longitudinal spacing between the transducer locations define measurement zones. Acoustic signals with chosen amplitude-time-frequency characteristics excite multiple vibration modes in the structure within each zone. Small mechanical changes in the inspection zones lead to scattering and attenuation of broadband acoustic signals, which are detectable as changes in received signal characteristics as part of a through-transmission technique. Additional use of short, narrowband pulse acoustic signals as part of a pulse-echo technique allows determination of the relative location of the mechanical change within each zone based on the differential delay profiles.

42 ENGINEERING↗

DOE OSTI2020

Deep Minima in the Triply Differential Cross Section for Ionization of Atomic Hydrogen by Electron and Positron Impact

We investigate ionization of atomic hydrogen by electron- and positron-impact. We apply the Coulomb–Born (CB1) approximation, various modified CB1 approximations, the three body distorted wave (3DW) approximation, and the time-dependent close-coupling (TDCC) method to electron-impact ionization of hydrogen. For electron-impact ionization of hydrogen for an incident energy of approximately 76.45 eV, we obtain a deep minimum in the CB1 triply differential cross section (TDCS). However, the TDCC for 74.45 eV and the 3DW for 74.46 eV gave a dip in the TDCS. For positron-hydrogen ionization (breakup) we apply the CB1 approximation and a modified CB1 approximation. We obtain a deep minimum in the TDCS and a zero in the CB1 transition matrix element for an incident energy of 100 eV with a gun angle of 56.13 °. Corresponding to a zero in the CB1 transition matrix element, there is a vortex in the velocity field associated with this element. For both electron- and positron-impact ionization of hydrogen the velocity field rotates in the same direction, which is anticlockwise. All calculations are performed for a doubly symmetric geometry; the electron-impact ionization is in-plane and the positron-impact ionization is out-of-plane.

74 ATOMIC AND MOLECULAR PHYSICS↗

NASA NTRSDate not supplied

Windform® XT 2.0 Use as 3U CubeSat Primary Structure

CubeSats provide a platform for small-scale space research and technology demonstration at reduced complexity, cost, and development time. These advantages drove the NASA Langley Research Center (LaRC) to develop and launch the GPX2 3U CubeSat to explore the viability of using Commercial Off-The-Shelf (COTS) differential GPS (dGPS) in low earth orbit. To reduce manufacturing costs and increase design flexibility, the project chose additive manufactured Windform® XT 2.0 as the primary bus material rather than traditional subtractive manufactured (milled) metal. The bus is a two-part, Selective Laser Sintered, 3D-print structure consisting of a single-piece, five-walled chassis and single-walled cover. The bus was specially designed to allow the project to accommodate the payload electronics stack as well as antennas, receivers, and deployable mechanisms. By using an additive manufactured solution, LaRC was able to design in features unrealizable through traditional milling, with a lead-time of roughly two weeks. In comparison, traditional subtractive manufacturing limits geometry options due to toolpath reach and bus construction would have required multiple components for each wall. This would have resulted in a more costly, longer lead-time article with more joints, fasteners, and complexity with a commensurate increase in overall mass. A number of lessons-learned were captured during the design, analysis, and testing of the GPX2 CubeSat covering thermal and structural analysis, vibration modeling, and geometric tolerancing. Additionally, a variety of material testing and verification were performed before and during spacecraft design and integration to assure the suitability of Windform® XT 2.0 for the launch and mission environments. This article provides the highlights of designing and testing the GPX2 bus.

CubeSat↗