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 217 records · Page 12

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Assessment of MOOSE-Based Tools for Calculating Radial Core Expansion

Radial core expansion in liquid-metal cooled fast reactor systems is a well-known phenomenon that produces strong reactivity feedback effects. An inherently safe reactor design takes advantage of negative reactivity feedback in accident scenarios by utilizing a core restraint system which produces a bowed shape that allows for radial expansion of the fuel regions. Detailed modeling and simulation of core radial expansion itself as well as subsequent reactivity feedback is a challenging task involving contact of many fuel assembly elements and physics feedback from neutronics, thermal hydraulics, and thermal mechanical response. A variety of physics codes have been developed to model aspects of radial core expansion but in general invoke geometrical or physics approximations. No code system currently exists which tightly and robustly couples these physics with enough detail to fully resolve the complex core radial expansion reactivity feedback effects. The future availability of such a code system is of vital importance to fully understanding the reactivity feedback effects that occur due to radial expansion, and consequently to optimizing the design of the core restraint system. A high-fidelity code will also be used to benchmark existing lower fidelity, faster running models to understand their benefits, limitations, and range of applications. A code development path using MOOSE-based tools is proposed in order to leverage the detailed geometry capabilities and natural tight coupling and robustness of MOOSE-based applications for modeling this complex phenomena. While simulation of the full phenomenon involves several physics, an assessment has been initiated on the capabilities and readiness of the currently available Tensor Mechanics module within MOOSE for calculation of the structural mechanical responses which occur within the reactor core. This report focuses on modeling the force-deformation response which mimics the physics of a duct contact deformation, as well as differential thermal expansion which produces a thermal bowed shaped for the fuel assemblies in a core. Simple examples were initially performed such as simple supported beam bending under load. The complexity of examples was progressively increased to better mimic the duct behavior by including a differential thermal example and inclusion of hexagonal cross-sections in the geometry. Further assessment of the structural mechanical response simulation capability is still required for modeling duct contact interactions and irradiation creep and swelling. Companion thermal hydraulic and neutronics assessments will also be required; these activities are planned for future years. Finally, integration of the multiple physics components through MOOSE is required to predict the core radial expansion and subsequent feedback effects.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2024

MatCal Users Guide: Release 1.3.0

Any continuum mechanics model will require three components: (1) a discretized geometry of the boundary value problem being studied, (2) the partial differential equations to be solved, and (3) the initial conditions and boundary conditions for the problem. To describe material behavior in these computational models, material models contribute to (2) the underlying equations and, occasionally, to (3) the initial conditions for the simulation. These material models can exhibit a mathematical form that is empirically based, based on first principles, or developed from both empirical observations and known physics. In general, these models are meant to represent a class of materials with well understood behavior. As a result, material models have parameters that must be tuned or calibrated so that the model response matches characterization data available for the specific material it is intended to represent when used to simulate a specific system. For simple models, such as isotropic, linear elastic materials in solid mechanics, this calibration process can be a simple analytical calculation directly extracting the parameters from experimental measurements. For complex models that have many inputs and require many characterization datasets to adequately identify the material behavior, the model calibration process can require an inverse problem approach where an optimization is performed to tune the model parameters to the available data.

36 MATERIALS SCIENCE↗

NASA NTRS1992

Three-dimensional co-culture process

The present invention relates to a 3-dimensional co-culture process, more particularly to methods or co-culturing at least two types of cells in a culture environment, either in space or in unit gravity, with minimum shear stress, freedom for 3-dimensional spatial orientation of the suspended particles and localization of particles with differing or similar sedimentation properties in a similar spatial region to form 3-dimensional tissue-like structures. Several examples of multicellular 3-dimensional experiences are included. The protocol and procedure are also set forth. The process allows simultaneous culture of multiple cell types and supporting substrates in a manner which does not disrupt the 3-dimensional spatial orientation of these components. The co-cultured cells cause a mutual induction effect which mimics the natural hormonal signals and cell interactions found in the intact organism. This causes the tissues to differentiate and form higher 3-dimensional structures such as glands, junctional complexes polypoid geometries, and microvilli which represent the corresponding in-vitro structures to a greater degree than when the cell types are cultured individually or by conventional processes. This process was clearly demonstrated for the case of two epithelial derived colon cancer lines, each co-cultured with normal human fibroblasts and with microcarrier bead substrates. The results clearly demonstrate increased 3-dimensional tissue-like structure and biochemical evidence of an increased differentiation state. With the present invention a variety of cells may be co-cultured to produce tissue which has 3-dimensionality and has some of the characteristics of in-vitro tissue. The process provides enhanced 3-dimensional tissue which create a multicellular organoid differentiation model.

Wolf, David A.↗

DOE OSTI2026

Ground and excited state gradients with end-to-end differentiable semiempirical quantum chemistry

Accurate and efficient gradients of molecular energy with respect to nuclear degrees of freedom are essential for geometry optimization and molecular dynamics, including simulations that go beyond the Born–Oppenheimer regime. A common approach involves deriving analytical formulas for new electronic structure methods, which is often conceptually difficult and requires tedious coding. Here, we implement analytical, semi-numerical, and automatic differentiation (AD)-based gradient pathways for semiempirical Hamiltonian models in the PYSEQM software package, leveraging both graphics processing unit (GPU) and central processing unit (CPU) architectures. We further extend these capabilities to excited states calculated using the configuration interaction singles and time-dependent Hartree–Fock ansätze. We benchmark wall time, peak memory usage, and accuracy across three molecular families of varying chemical complexity, including systems of up to a thousand atoms. For ground-state simulations, analytical and AD gradients achieve near-identical GPU runtimes, while semi-numerical gradients are slower on GPU but remain competitive on CPU. For excited states, both analytical and custom AD approaches using implicit differentiation show similar performance and low memory requirements, whereas gradients with full AD are memory-limited. AD gradients match analytical ones in accuracy across all tested systems, aided by a quaternion-based diatomic frame rotation for two-center quantities that ensures smooth energy surfaces. Overall, automatic differentiation emerges as a practical alternative to analytical gradients in semiempirical quantum chemistry, offering high accuracy while allowing seamless integration in AI-driven workflows and popular packages, such as PyTorch and JAX. Our results provide actionable guidance for selecting optimal gradient strategies in large-scale ground- and excited-state molecular dynamics simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1983

Angular distribution for electron excitation of the 4(2)S yields 4(2)P transition in Zn II - Comparison of experiment and theory

Differential electron-scattering cross sections for inelastic excitation of an ion have been measured for the first time. Experiments were carried out in a cross electron-ion beam geometry for the 4(2)S yields 4(2)P transition in Zn II at 75 eV. In addition, differential cross sections were calculated at energies between 15 and 100 eV in a five-state close-coupling approximation in which 4s, 4p, 3d(9)4s(2), 5s, and 4 d states were included. Agreement in shape between theory and experiment at 75 eV is excellent.

Chutjian, A.↗

NASA NTRS1990

Time-optimal Aircraft Pursuit-evasion with a Weapon Envelope Constraint

The optimal pursuit-evasion problem between two aircraft including a realistic weapon envelope is analyzed using differential game theory. Six order nonlinear point mass vehicle models are employed and the inclusion of an arbitrary weapon envelope geometry is allowed. The performance index is a linear combination of flight time and the square of the vehicle acceleration. Closed form solution to this high-order differential game is then obtained using feedback linearization. The solution is in the form of a feedback guidance law together with a quartic polynomial for time-to-go. Due to its modest computational requirements, this nonlinear guidance law is useful for on-board real-time implementation.

Menon, P. K. A.↗

DOE OSTI2024

Improving River Routing Using a Differentiable Muskingum‐Cunge Model and Physics‐Informed Machine Learning

Recently, rainfall-runoff simulations in small headwater basins have been improved by methodological advances such as deep neural networks (NNs) and hybrid physics-NN models—particularly, a genre called differentiable modeling that intermingles NNs with physics to learn relationships between variables. However, hydrologic routing simulations, necessary for simulating floods in stem rivers downstream of large heterogeneous basins, had not yet benefited from these advances and it was unclear if the routing process could be improved via coupled NNs. We present a novel differentiable routing method (δMC-Juniata-hydroDL2) that mimics the classical Muskingum-Cunge routing model over a river network but embeds an NN to infer parameterizations for Manning's roughness (n) and channel geometries from raw reach-scale attributes like catchment areas and sinuosity. The NN was trained solely on downstream hydrographs. Synthetic experiments show that while the channel geometry parameter was unidentifiable, n can be identified with moderate precision. With real-world data, the trained differentiable routing model produced more accurate long-term routing results for both the training gage and untrained inner gages for larger subbasins (>2,000 km2) than either a machine learning model assuming homogeneity, or simply using the sum of runoff from subbasins. The n parameterization trained on short periods gave high performance in other periods, despite significant errors in runoff inputs. The learned n pattern was consistent with literature expectations, demonstrating the framework's potential for knowledge discovery, but the absolute values can vary depending on training periods. The trained n parameterization can be coupled with traditional models to improve national-scale hydrologic flood simulations.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1993

A simulation of GPS and differential GPS sensors

The Global Positioning System (GPS) is a revolutionary advance in navigation. Users can determine latitude, longitude, and altitude by receiving range information from at least four satellites. The statistical accuracy of the user's position is directly proportional to the statistical accuracy of the range measurement. Range errors are caused by clock errors, ephemeris errors, atmospheric delays, multipath errors, and receiver noise. Selective Availability, which the military uses to intentionally degrade accuracy for non-authorized users, is a major error source. The proportionality constant relating position errors to range errors is the Dilution of Precision (DOP) which is a function of the satellite geometry. Receivers separated by relatively short distances have the same satellite and atmospheric errors. Differential GPS (DGPS) removes these errors by transmitting pseudorange corrections from a fixed receiver to a mobile receiver. The corrected pseudorange at the moving receiver is now corrupted only by errors from the receiver clock, multipath, and measurement noise. This paper describes a software package that models position errors for various GPS and DGPS systems. The error model is used in the Real-Time Simulator and Cockpit Technology workstation simulations at NASA-LaRC. The GPS/DGPS sensor can simulate enroute navigation, instrument approaches, or on-airport navigation.

Rankin, James M.↗

DOE OSTI2023

Multislice forward modeling of coherent surface scattering imaging on surface and interfacial structures

To study nanostructures on substrates, surface-sensitive reflection-geometry scattering techniques such as grazing incident small angle X-ray scattering are commonly used to yield an averaged statistical structural information of the surface sample. Grazing incidence geometry can probe the absolute three-dimensional structural morphology of the sample if a highly coherent beam is used. Coherent surface scattering imaging (CSSI) is a powerful yet non-invasive technique similar to coherent X-ray diffractive imaging (CDI) but performed at small angles and grazing-incidence reflection geometry. A challenge with CSSI is that conventional CDI reconstruction techniques cannot be directly applied to CSSI because the Fourier-transform-based forward models cannot reproduce the dynamical scattering phenomenon near the critical angle of total external reflection of the substrate-supported samples. To overcome this challenge, we have developed a multislice forward model which can successfully simulate the dynamical or multi-beam scattering generated from surface structures and the underlying substrate. The forward model is also demonstrated to be able to reconstruct an elongated 3D pattern from a single shot scattering image in the CSSI geometry through fast-performing CUDA-assisted PyTorch optimization with automatic differentiation.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

DOE OSTI2021

Deformation and strength of mantle relevant garnets: Implications for the subduction of basaltic-rich crust

Garnet is an important mineral phase in the upper mantle as it is both a key component in bulk mantle rocks, and a primary phase at high pressure within subducted basalt. In this work, we focus on the strength of garnet and the texture that develops within garnet during accommodation of differential deformational strain. We use X-ray diffraction in a radial geometry to analyze texture development in situ in three garnet compositions under pressure at 300 K: a natural garnet (Prp 60 Alm 37 ) to 30 GPa, and two synthetic majorite-bearing compositions (Prp 59 Maj 41 and Prp 42 Maj 58 ) to 44 GPa. All three garnets develop a modest (100) texture at elevated pressure under axial compression. Elasto-viscoplastic self-consistent (EVPSC) modeling suggests that two slip systems are active in the three garnet compositions at all pressures studied: {110}<$1\bar{1}1$> and {001}<110>. We determine a flow strength of ~5 GPa at pressures between 10 to 15 GPa for all three garnets; these values are higher than previously reported yield strengths measured on natural and majoritic garnets. Strengths calculated using the experimental lattice strain differ from the strength generated from those calculated using EVPSC. Prp 67 Alm 33 , Prp 59 Maj 41 , and Prp 42 Maj 58 are of comparable strength to each other at room temperature, which indicates that majorite substitution does not greatly affect the strength of garnets. Additionally, all three garnets are of similar strength as lower mantle phases such as bridgmanite and ferropericlase, suggesting that garnet may not be notably stronger than the surrounding lower mantle/deep upper mantle phases at the base of the upper mantle.

58 GEOSCIENCES↗

DOE OSTI2022

Taylor wave solution for a general equation of state

This document describes a solution procedure for calculating the Taylor wave behind an unsupported Chapman–Jouguet (CJ) detonation in planar, cylindrical, and spherical geometries given a general equation of state. The resulting semi-analytic solution can be utilized to examine new equation of state models for detonation products and during the verification of hydrodynamic codes. The governing partial differential equations are reduced to ordinary differential equations in both characteristic and self-similar forms. The first-order systems corresponding to each geometry are amenable to solution numerically using commonly available methods. A difficulty arises at the CJ point in radial coordinates where the similarity equations become singular. Two separate strategies are proposed to integrate the first-order system. The first one uses an asymptotic approximation near the CJ point that can be used to perturb the boundary conditions. The second one applies a change of variables which removes the singularity at the expense of an additional equation to be integrated. A test problem is provided for the Davis products equation of state to illustrate the qualitative features of the Taylor wave in each geometric configuration and compared with a Lagrangian hydrodynamics research code. A Python code listing gives an implementation using the SciPy library to assists users in generating the results.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1999

The Magnetospheric Multiscale Mission...Resolving Fundamental Processes in Space Plasmas

The Magnetospheric Multiscale (MMS) mission is a multiple-spacecraft Solar-Terrestrial Probe designed to study the microphysics of magnetic reconnection, charged particle acceleration, and turbulence in key boundary regions of Earth's magnetosphere. These three processes, which control the flow of energy, mass, and momentum within and across plasma boundaries, occur throughout the universe and are fundamental to our understanding of astrophysical and solar system plasmas. Only in Earth's magnetosphere, however, are they readily accessible for sustained study through in-situ measurement. MMS will employ five co-orbiting spacecraft identically instrumented to measure electric and magnetic fields, plasmas, and energetic particles. The initial parameters of the individual spacecraft orbits will be designed so that the spacecraft formation will evolve into a three-dimensional configuration near apogee, allowing MMS to differentiate between spatial and temporal effects and to determine the three dimensional geometry of plasma, field, and current structures. In order to sample all of the magnetospheric boundary regions, MMS will employ a unique four-phase orbital strategy involving carefully sequenced changes in the local time and radial distance of apogee and, in the third phase, a change in orbit inclination from 10 degrees to 90 degrees. The nominal mission operational lifetime is two years. Launch is currently scheduled for 2006.

Curtis, S.↗

NASA NTRSDate not supplied

Formal Verification of the Interaction Between Semi-Algebraic Sets and Real Analytic Functions

Semi-algebraic sets and real analytic functions are fundamental concepts in Real Algebraic Geometry and Real Analysis, respectively. These concepts interact in the study of Differential Equations, where the real analytic solution to a differential equation is known to enter or exit a semi-algebraic set in a predicable way. Motivated to enhance the capability to reason about differential equations in the Prototype Verification System (PVS), a formalization of multivariate polynomials, semi-algebraic sets, and real analytic functions is developed. The favorable way that a real analytic function enters and exits a semi-algebraic set is proven. It is further shown that if the function is assumed to be smooth, a slightly weaker assumption than real analytic, these favorable interactions with semi-algebraic sets may fail.

Real analytic functions↗

NASA NTRS1984

Development of MCAERO wing design panel method with interactive graphics module

A reliable and efficient iterative method has been developed for designing wing section contours corresponding to a prescribed subcritical pressure distribution. The design process is initialized by using MCAERO (MCAIR 3-D Subsonic Potential Flow Analysis Code) to analyze a baseline configuration. A second program DMCAERO is then used to calculate a matrix containing the partial derivative of potential at each control point with respect to each unknown geometry parameter by applying a first-order expansion to the baseline equations in MCAERO. This matrix is calculated only once but is used in each iteration cycle to calculate the geometry perturbation and to analyze the perturbed geometry. The potential on the new geometry is calculated by linear extrapolation from the baseline solution. This extrapolated potential is converted to velocity by numerical differentiation, and velocity is converted to pressure by using Bernoulli's equation. There is an interactive graphics option which allows the user to graphically display the results of the design process and to interactively change either the geometry or the prescribed pressure distribution.

Hawk, J. D.↗

NASA NTRS1974

Three-dimensional elastic stress and displacement analysis of finite geometry solids containing cracks

The line method of analysis is applied to the Navier-Cauchy equations of elastic equilibrium to calculate the displacement distributions in various bodies containing cracks. The application of this method to these equations leads to coupled sets of simultaneous ordinary differential equations whose solutions are obtained along sets of lines in a discretized region. When decoupling the equations and their boundary conditions is not possible, the use of a successive approximation procedure permits the analytical solution of the resulting ordinary differential equations. The results obtained show a considerable potential for using this method in the three-dimensional analysis of finite geometry solids and suggest a possible extension of this technique to nonlinear material behavior.

Gyekenyesi, J. P.↗

NASA NTRS1975

Three-dimensional elastic stress and displacement analysis of finite geometry solids containing cracks

Gyekenyesi, J. P.↗

NASA NTRS1990

The method of lines in analyzing solids containing cracks

A semi-numerical method is reviewed for solving a set of coupled partial differential equations subject to mixed and possibly coupled boundary conditions. The line method of analysis is applied to the Navier-Cauchy equations of elastic and elastoplastic equilibrium to calculate the displacement distributions in various, simple geometry bodies containing cracks. The application of this method to the appropriate field equations leads to coupled sets of simultaneous ordinary differential equations whose solutions are obtained along sets of lines in a discretized region. When decoupling of the equations and their boundary conditions is not possible, the use of a successive approximation procedure permits the analytical solution of the resulting ordinary differential equations. The use of this method is illustrated by reviewing and presenting selected solutions of mixed boundary value problems in three dimensional fracture mechanics. These solutions are of great importance in fracture toughness testing, where accurate stress and displacement distributions are required for the calculation of certain fracture parameters. Computations obtained for typical flawed specimens include that for elastic as well as elastoplastic response. Problems in both Cartesian and cylindrical coordinate systems are included. Results are summarized for a finite geometry rectangular bar with a central through-the-thickness or rectangular surface crack under remote uniaxial tension. In addition, stress and displacement distributions are reviewed for finite circular bars with embedded penny-shaped cracks, and rods with external annular or ring cracks under opening mode tension. The results obtained show that the method of lines presents a systematic approach to the solution of some three-dimensional mechanics problems with arbitrary boundary conditions. The advantage of this method over other numerical solutions is that good results are obtained even from the use of a relatively coarse grid.

Gyekenyesi, John P.↗

NASA NTRS1996

Geometric Model for a Parametric Study of the Blended-Wing-Body Airplane

A parametric model is presented for the blended-wing-body airplane, one concept being proposed for the next generation of large subsonic transports. The model is defined in terms of a small set of parameters which facilitates analysis and optimization during the conceptual design process. The model is generated from a preliminary CAD geometry. From this geometry, airfoil cross sections are cut at selected locations and fitted with analytic curves. The airfoils are then used as boundaries for surfaces defined as the solution of partial differential equations. Both the airfoil curves and the surfaces are generated with free parameters selected to give a good representation of the original geometry. The original surface is compared with the parametric model, and solutions of the Euler equations for compressible flow are computed for both geometries. The parametric model is a good approximation of the CAD model and the computed solutions are qualitatively similar. An optimal NURBS approximation is constructed and can be used by a CAD model for further refinement or modification of the original geometry.

Mastin, C. Wayne↗