Search NASA⌕ Search

SEARCH · Search NASA

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.

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

At least 361 records · Page 20

Spectral element methods: Algorithms and architectures

Spectral element methods are high-order weighted residual techniques for partial differential equations that combine the geometric flexibility of finite element methods with the rapid convergence of spectral techniques. Spectral element methods are described for the simulation of incompressible fluid flows, with special emphasis on implementation of spectral element techniques on medium-grained parallel processors. Two parallel architectures are considered: the first, a commercially available message-passing hypercube system; the second, a developmental reconfigurable architecture based on Geometry-Defining Processors. High parallel efficiency is obtained in hypercube spectral element computations, indicating that load balancing and communication issues can be successfully addressed by a high-order technique/medium-grained processor algorithm-architecture coupling.

Fischer, Paul↗

Spectral element methods - Algorithms and architectures

Spectral element methods are high-order weighted residual techniques for partial differential equations that combine the geometric flexibility of finite element methods with the rapid convergence of spectral techniques. Spectral element methods are described for the simulation of incompressible fluid flows, with special emphasis on implementation of spectral element techniques on medium-grained parallel processors. Two parallel architectures are considered; the first, a commercially available message-passing hypercube system; the second, a developmental reconfigurable architecture based on Geometry-Defining Processors. High parallel efficiency is obtained in hypercube spectral element computations, indicating that load balancing and communication issues can be successfully addressed by a high-order technique/medium-grained processor algorithm-architecture coupling.

Fischer, Paul↗

Adhesive-Bonded Composite Joint Analysis with Delaminated Surface Ply Using Strain-Energy Release Rate

This paper presents an analytical model to determine the strain energy release rate due to an interlaminar crack of the surface ply in adhesively bonded composite joints subjected to axial tension. Single-lap shear-joint standard test specimen geometry with thick bondline is followed for model development. The field equations are formulated by using the first-order shear-deformation theory in laminated plates together with kinematics relations and force equilibrium conditions. The stress distributions for the adherends and adhesive are determined after the appropriate boundary and loading conditions are applied and the equations for the field displacements are solved. The system of second-order differential equations is solved to using the symbolic computation tool Maple 9.52 to provide displacements fields. The equivalent forces at the tip of the prescribed interlaminar crack are obtained based on interlaminar stress distributions. The strain energy release rate of the crack is then determined by using the crack closure method. Finite element analyses using the J integral as well as the crack closure method are performed to verify the developed analytical model. It has been shown that the results using the analytical method correlate well with the results from the finite element analyses. An attempt is made to predict the failure loads of the joints based on limited test data from the literature. The effectiveness of the inclusion of bondline thickness is justified when compared with the results obtained from the previous model in which a thin bondline and uniform adhesive stresses through the bondline thickness are assumed.

Chadegani, Alireza↗

Machine Learning to Select Experiments Driven by Fundamental Science and Applications for Targeted Nuclear Data Improvement

This work describes a blueprint for a process that accelerates progress in science by quantitatively answering the following question: What is the optimal combination of fundamental-science and application-driven experiments to maximally reduce pertinent data uncertainties? Answering this question entails solving a high-dimensional and complex optimization problem that is best solved with advanced statistic techniques often classified as machine learning. We apply this process within the framework of nuclear data with the aim to select an experiment combination that will reduce uncertainties in 239 Pu nuclear data for neutron energies between 1 and 600 keV. In this field, fundamental-physics driven data, called differential, look at one nuclear physics observable at a time. They are contrasted to application-driven, integral, data where one or few resulting values inform a broad set of nuclear data across several nuclides and energies. The candidates for integral experiments are criticality measurements that were refined by a genetic algorithm to be maximally sensitive to 239 Pu fission cross sections in the desired energy range. Twenty-three candidate differential experiments were investigated and span multiple nuclear physics observables (e.g., total, capture cross sections) for isotopes appearing in the integral experiments. The optimal combination among these candidate experiments was investigated via generalized least squares fitting, augmented with Gaussian processes to ameliorate statistical irregularities in data, and the D-optimality criterion. The latter evaluates for each pair of candidates the joint reduction in uncertainties of all 12200 nuclear data appearing in the integral experiments compared to the knowledge we have from 168 past experiments, theory, and nuclear data. We chose as differential measurements those that investigate 63 Cu and 239 Pu total cross sections, based on D-optimality rank and feasibility constraints. Two integral (criticality) experiments were selected: An experiment with Al 2 ⁢O 3 and graphite interleaved with Pu and a thick Cu reflector explores 1–30 keV, while we target the 30–600 keV range with an experiment that swaps boron in place of graphite with a different geometry.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Compatibility of molten plutonium with wrought and additively manufactured metal crucibles

Understanding plutonium’s interaction with metals is crucial for optimizing pyrochemical operations, nuclear fuel containment, and various actinide processing techniques. Traditionally, tantalum crucibles are employed for plutonium processing due to their high durability, excellent temperature stability, and low solubility in plutonium. However, tantalum faces challenges such as plutonium wetting and diffusion, making surface coatings particularly important for crucibles in pyrochemical applications to enhance corrosion resistance against plutonium. Tantalum is also expensive and difficult to machine, prompting the need for advanced manufacturing techniques to address these challenges. Here, in this work, we investigate the interaction of Pu with tantalum and titanium crucibles fabricated using both traditional machining methods and laser powder bed fusion (LPBF) additive manufacturing (AM). LPBF-AM is an advanced technique that allows for the creation of complex geometries from traditionally difficult-to-machine metals by using a high-powered laser to build parts. Previous studies of conventional manufactured tantalum have utilized oxidation and carburization of the surface to mitigate plutonium wetting; however, no studies of surface modified LPBF-AM material have been undertaken. These studies are crucial, given the typical differences in the grain structure between conventional and LPBF-AM materials. All crucibles underwent differential scanning calorimetry to confirm the melting of plutonium. Subsequently, the crucibles were sectioned and mounted in epoxy for microstructural analysis using optical microscopy and scanning electron microscopy. This investigation, comparing the performance of wrought vs AM metal crucibles, provides a basis for future tooling applications in actinide processing techniques and can address the challenges associated with traditional machining, particularly in pyrochemical applications.

Actinides↗

Coulomb Branch Amplitudes from a Deformed Amplituhedron Geometry

The amplituhedron provides, via geometric means, the all-loop integrand of scattering amplitudes in maximally supersymmetric Yang-Mills theory. Unfortunately, dimensional regularization, used conventionally for integration, breaks the beautiful geometric picture. This motivates us to propose a “deformed” amplituhedron. Focusing on the four-particle amplitude, we introduce two deformation parameters, which can be interpreted as particle masses. We provide evidence that the mass pattern corresponds to a specific choice of vacuum expectation values on the Coulomb branch. The deformed amplitude is infrared finite, making the answer well defined in four dimensions. Leveraging four-dimensional integration techniques based on differential equations, we compute the amplitude up to two loops. In the limit where the deformation parameters are taken to zero, we recover the known Bern-Dixon-Smirnov amplitude. In the limit where only one deformation parameter is taken to zero, we find a connection to the angle-dependent cusp anomalous dimension. Published by the American Physical Society 2024

Physics↗

Nature of convection-stabilized dc arcs in dual-flow nozzle geometry. I - The cold flow field and dc arc characteristics. II - Optical diagnostics and theory

Steady-state low-current air arcs in a dual-flow nozzle system are studied experimentally. The cold flow field with no arc is investigated using a 12.7-mm diameter dual-flow nozzle in a steady-flow facility. Mach number and mass flux distributions are determined for various nozzle-pressure ratios and nozzle-gap spacing. It is found that the shock waves in the converging-diverging nozzles result in a decrease in overal resistance by about 15 percent. Also, Schlieren and differential interferometry techniques are used to visualize the density gradients within the arc plasma and thermal mantle. Both optical techniques reveal a laminar arc structure for a reservoir pressure of 1 atm at various current levels. Experimentally determined axial static pressure and cold-flow mass flux rate distributions and a channel-flow model with constant arc temperatre are used to solve the energy integral for the arc radius as a function of axial distance. The arc electric field strength, voltage, resistance, and power are determined with Ohm's law and the total heat transfer is related to arc power.

Serbetci, Ilter↗

Hypercomplex Automatic Differentiation in the Eulerian Hydrocode PAGOSA

Enabling the computation of partial derivatives or sensitivities in production hydrocodes is beneficial for design, optimization, sensitivity analysis, and uncertainty quantification. Traditional finite difference approximations of these sensitivities are inefficient since convergence studies of the step size is required for each parameter of interest. For these reasons, HYPercomplex Automatic Differentiation (HYPAD) was implemented in the Eulerian hydrocode PAGOSA. HYPAD is analogous to forward-mode automatic differentiation except hypercomplex numbers (numbers with multiple imaginary parts) are used instead of dual numbers. Accurate partial derivatives can be computed of all state variables with respect to multiple input variables in a single run. The method was implemented using operator overloading to handle hypercomplex algebra. HYPAD was demonstrated and verified on Sod’s shock tube problem to compute derivatives of the state variables with respect to a material parameter, initial conditions, and geometry.

97 MATHEMATICS AND COMPUTING↗

Current-State Constrained Filter Bank for Wald Testing of Spacecraft Conjunctions

We propose a filter bank consisting of an ordinary current-state extended Kalman filter, and two similar but constrained filters: one is constrained by a null hypothesis that the miss distance between two conjuncting spacecraft is inside their combined hard body radius at the predicted time of closest approach, and one is constrained by an alternative complementary hypothesis. The unconstrained filter is the basis of an initial screening for close approaches of interest. Once the initial screening detects a possibly risky conjunction, the unconstrained filter also governs measurement editing for all three filters, and predicts the time of closest approach. The constrained filters operate only when conjunctions of interest occur. The computed likelihoods of the innovations of the two constrained filters form a ratio for a Wald sequential probability ratio test. The Wald test guides risk mitigation maneuver decisions based on explicit false alarm and missed detection criteria. Since only current-state Kalman filtering is required to compute the innovations for the likelihood ratio, the present approach does not require the mapping of probability density forward to the time of closest approach. Instead, the hard-body constraint manifold is mapped to the filter update time by applying a sigma-point transformation to a projection function. Although many projectors are available, we choose one based on Lambert-style differential correction of the current-state velocity. We have tested our method using a scenario based on the Magnetospheric Multi-Scale mission, scheduled for launch in late 2014. This mission involves formation flight in highly elliptical orbits of four spinning spacecraft equipped with antennas extending 120 meters tip-to-tip. Eccentricities range from 0.82 to 0.91, and close approaches generally occur in the vicinity of perigee, where rapid changes in geometry may occur. Testing the method using two 12,000-case Monte Carlo simulations, we found the method achieved a missed detection rate of 0.1%, and a false alarm rate of 2%.

Carpenter, J. Russell↗

Analysis of Triangular Rollable and Collapsible Composite Booms Under the Effects of Gravity and Twist

Finite element analysis is used to investigate Triangular Rollable and Collapsible (TRAC) composite booms. The TRAC booms are subjected to a suite of loading conditions, including in-plane bending, out-of-plane bending, and axial compression. Simulations containing as-designed boom geometry are performed, as are simulations with geometric deviations associated with the as-built conditions of a sample TRAC boom. Namely, a cross-sectional twist along the length of the boom is incorporated into the analysis to study the effect of the twist on the deformation characteristics of the boom. Twisting in the sample boom was observed after short-term storage in a rolled configuration and may be attributed to differential creep deformation between the inner and outer flanges of the TRAC boom. Simulations are conducted with and without the effects of gravity to understand the significance of gravity on testing conditions. Booms of three lengths are evaluated: 3 meters, 7.5 meters and 30 meters. The mechanical responses of the TRAC booms are discussed with reference to the effects of gravity, boom length and cross-sectional twisting. Gravity is found to have a more significant effect on boom mechanical response as boom length is increased. Cross-sectional twisting is found to have a more significant effect on boom mechanical response in the shorter booms.

High-strain composites↗

Analysis of Triangular Rollable and Collapsible (TRAC) Composite Booms under the Effects of Gravity and Twist

Finite element analysis is used to investigate Triangular Rollable and Collapsible (TRAC) composite booms. The TRAC booms are subjected to a suite of loading conditions, including in-plane bending, out-of-plane bending, and axial compression. Simulations containing as-designed boom geometry are performed, as are simulations with geometric deviations associated with the as-built conditions of a sample TRAC boom. Namely, a cross-sectional twist along the length of the boom is incorporated into the analysis to study the effect of the twist on the deformation characteristics of the boom. Twisting in the sample boom was observed after short-term storage in a rolled configuration and may be attributed to differential creep deformation between the inner and outer flanges of the TRAC boom. Simulations are conducted with and without the effects of gravity to understand the significance of gravity on testing conditions. Booms of three lengths are evaluated: 3 meters, 7.5 meters and 30 meters. The mechanical responses of the TRAC booms are discussed with reference to the effects of gravity, boom length and cross-sectional twisting. Gravity is found to have a more significant effect on boom mechanical response as boom length is increased. Cross-sectional twisting is found to have a more significant effect on boom mechanical response in the shorter booms.

High-strain composites↗

Direct, simple and efficient computation of all components of the virtual-casing magnetic field in axisymmetric geometries with Kapur–Rokhlin quadrature

In a recent publication (Toler et al., J. Plasma Phys., vol. 89, issue 2, 2023, p. 905890210), we demonstrated that for axisymmetric geometries, the Kapur–Rokhlin quadrature rule provided an efficient and high-order accurate method for computing the normal component, on the plasma surface, of the magnetic field due to the toroidal current flowing in the plasma, via the virtual-casing principle. The calculation was indirect, as it required the prior computation of the magnetic vector potential from the virtual-casing principle, followed by the computation of its tangential derivative by Fourier differentiation, to obtain the normal component of the magnetic field. Our approach did not provide the other components of the virtual-casing magnetic field. In this letter, we show that a more direct and more general approach is available for the computation of the virtual-casing magnetic field. The Kapur–Rokhlin quadrature rule accurately calculates the principal value integrals in the expression for all the components of the magnetic field on the plasma boundary, and the numerical error converges at a rate nearly as high as the indirect method we presented previously.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Recognition of 3-D symmetric objects from range images in automated assembly tasks

A new technique is presented for the three dimensional recognition of symmetric objects from range images. Beginning from the implicit representation of quadrics, a set of ten coefficients is determined for symmetric objects like spheres, cones, cylinders, ellipsoids, and parallelepipeds. Instead of using these ten coefficients trying to fit them to smooth surfaces (patches) based on the traditional way of determining curvatures, a new approach based on two dimensional geometry is used. For each symmetric object, a unique set of two dimensional curves is obtained from the various angles at which the object is intersected with a plane. Using the same ten coefficients obtained earlier and based on the discriminant method, each of these curves is classified as a parabola, circle, ellipse, or hyperbola. Each symmetric object is found to possess a unique set of these two dimensional curves whereby it can be differentiated from the others. It is shown that instead of using the three dimensional discriminant which involves evaluation of the rank of its matrix, it is sufficient to use the two dimensional discriminant which only requires three arithmetic operations.

Alvertos, Nicolas↗

Three-dimensional aspects of interplanetary shock waves

The three-dimensional geometries of two classes of interplanetary shock waves are analyzed: (1) the more or less spherical shock wave originating from some short-lived solar event, such as a solar flare, and (2) the corotating shock pair formed by the interaction of long-lived solar wind streams. The systematic distortion of a spherical wave due to persistent latitudinal solar wind structure should amount to 25 deg for an equator-to-pole solar wind speed differential of 200 km/sec. A geometrical argument is given for determining the heliocentric distance of the onset of corotating shocks formed by contiguous long-lived solar streams. Estimates are obtained for the heliocentric distance to the inner edges of shock waves as a function of speed difference between the streams and as a function of heliographic latitude. Means by which the three-dimensional structures studied can be observed by out-of-the-ecliptic space missions are discussed.

Siscoe, G. L.↗

Energy transfer in volume-reflecting heat shields

An approximate analysis of radiative transfer in highly scattering materials was developed based on the Kubelka-Munk differential equations--a set of two differential equations representing the spatial rate of change of radiative half-fluxes within the scattering media. These approximate solutions of the Kubelka-Munk equations together with analytic solutions for the steady state temperature distribution for two types of boundary conditions are given. These solutions show the influence of back surface reflectance, scattering power, incident radiative flux parameter and boundary conductive flux parameter on overall reflectance and temperature distributions. This radiation field analysis, adapted to spherical geometry, was applied to the evaluation of the thermal performance of teflon and fritted quartz as heat protection materials for entry into the atmosphere of Jupiter.

Weston, K. C.↗

Dynamics of early planetary gear trains

A method to analyze the static and dynamic loads in a planetary gear train was developed. A variable-variable mesh stiffness (VVMS) model was used to simulate the external and internal spur gear mesh behavior, and an equivalent conventional gear train concept was adapted for the dynamic studies. The analysis can be applied either involute or noninvolute spur gearing. By utilizing the equivalent gear train concept, the developed method may be extended for use for all types of epicyclic gearing. The method is incorporated into a computer program so that the static and dynamic behavior of individual components can be examined. Items considered in the analysis are: (1) static and dynamic load sharing among the planets; (2) floating or fixed Sun gear; (3) actual tooth geometry, including errors and modifications; (4) positioning errors of the planet gears; (5) torque variations due to noninvolute gear action. A mathematical model comprised of power source, load, and planetary transmission is used to determine the instantaneous loads to which the components are subjected. It considers fluctuating output torque, elastic behavior in the system, and loss of contact between gear teeth. The dynamic model has nine degrees of freedom resulting in a set of simultaneous second order differential equations with time varying coefficients, which are solved numerically. The computer program was used to determine the effect of manufacturing errors, damping and component stiffness, and transmitted load on dynamic behavior. It is indicated that this methodology offers the designer/analyst a comprehensive tool with which planetary drives may be quickly and effectively evaluated.

August, R.↗

Initial On-Orbit Spatial Resolution Characterization of OrbView-3 Panchromatic Images

Characterization was conducted under the Memorandum of Understanding among Orbital Sciences Corp., ORBIMAGE, Inc., and NASA Applied Sciences Directorate. Acquired five OrbView-3 panchromatic images of the permanent Stennis Space Center edge targets painted on a concrete surface. Each image is available at two processing levels: Georaw and Basic. Georaw is an intermediate image in which individual pixels are aligned by a nominal shift in the along-scan direction to adjust for the staggered layout of the panchromatic detectors along the focal plane array. Georaw images are engineering data and are not delivered to customers. The Basic product includes a cubic interpolation to align the pixels better along the focal plane and to correct for sensor artifacts, such as smile and attitude smoothing. This product retains satellite geometry - no rectification is performed. Processing of the characterized images did not include image sharpening, which is applied by default to OrbView-3 image products delivered by ORBIMAGE to customers. Edge responses were extracted from images of tilted edges in two directions: along-scan and cross-scan. Each edge response was approximated with a superposition of three sigmoidal functions through a nonlinear least-squares curve-fitting. Line Spread Functions (LSF) were derived by differentiation of the analytical approximation. Modulation Transfer Functions (MTF) were obtained after applying the discrete Fourier transform to the LSF.

Blonski, Slawomir↗

Streamlined Convergence Acceleration for CFD Codes

Enigma, a simplified interface to the PETSc library, is shown to enable the rapid solution of discrete partial differential equations. Two CFD codes, LAURA and HyperSolve, use Enigma to compute steady solutions of the Navier-Stokes equations. Using PETSc, Enigma is shown to provide a Jacobian-Free Newton-Krylov method (JFNK), globalized with pseudotransient continuation, that improves efficiency over the point-implicit relaxation method traditionally used by LAURA. It is shown that iterative error has a large impact on surface heat transfer predicted by LAURA on an axisymmetric sphere-cone geometry. Also, the convergence rate of HyperSolve simulating subsonic flow over a delta wing geometry with the JFNK method is shown to be more efficient than employing a defect correction method as the nonlinear solver.

Thompson, Kyle B.↗