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 559 records · Page 31

All sources NASA NTRS DOE OSTI

NASA NTRS1987

Covariance and the hierarchy of frame bundles

This is an essay on the general concept of covariance, and its connection with the structure of the nested set of higher frame bundles over a differentiable manifold. Examples of covariant geometric objects include not only linear tensor fields, densities and forms, but affinity fields, sectors and sector forms, higher order frame fields, etc., often having nonlinear transformation rules and Lie derivatives. The intrinsic, or invariant, sets of forms that arise on frame bundles satisfy the graded Cartan-Maurer structure equations of an infinite Lie algebra. Reduction of these gives invariant structure equations for Lie pseudogroups, and for G-structures of various orders. Some new results are introduced for prolongation of structure equations, and for treatment of Riemannian geometry with higher-order moving frames. The use of invariant form equations for nonlinear field physics is implicitly advocated.

Estabrook, Frank B.↗

NASA NTRS1989

Relative tracking of multiple spacecraft by interferometry

An analysis is conducted of earth-based radiometric tracking of one spacecraft relative to an angularly nearby second spacecraft. Two cases are studied: relative positioning between a lander and a rover on the surface of Mars and relative tracking between a Mars lander and a Mars orbiter. All spacecraft signals are simultaneously received in the same beamwidth of an earth tracking antenna. Differential interferometric measurement errors are predicted. Errors which scale with angular separation between sources and errors which scale with temporal separation between measurement epochs are reduced virtually to zero. System thermal noise and systematic phase shifts introduced by receiver electronics typically dominate the error budget. Solar plasma delays become dominant for signal paths which pass close to the sun. Precise line-of-sight range measurements, differenced between stations, are also considered. Meter-level accuracy is obtained for lander/rover relative position by combining interferometric and precise range measurements. Either data type alone, for geometries where earth is not near zero declination as seen from Mars and Mars is not near zero declination as seen from earth, can provide accuracy at the 10-100-m level.

Border, James S.↗

DOE OSTI2026

Magnus force induced magnetic diode effect in skyrmion systems

We show that skyrmions can exhibit a “magnetic diode effect,” where there is a nonreciprocal response in the transport when the magnetic field is reversed. This effect can be achieved for skyrmions moving in channels with a sawtooth potential on one side and a reversed sawtooth potential on the other side. We consider the cases of both spin-transfer torque (STT) and spin–orbit torque (SOT). When the magnetic field is held fixed, the velocity response of the skyrmion is the same for current applied in either direction for both torques, so there is no current diode effect. When the magnetic field is reversed, under STT driving the velocity of the skyrmion reverses and its absolute value changes. Under SOT driving, the velocity remains in the same direction but drops to a much lower value, resulting in negative differential conductivity. For a fixed current, we find a nonreciprocal skyrmion velocity as a function of the applied field’s sign, in analogy to the velocity–current curves observed in the usual diode effect. The nonreciprocity is generated by the Magnus force, which causes skyrmions to interact preferentially with one side of the channel. Since the channel sides have opposite asymmetry, a positive magnetic field can cause the skyrmion to interact with the “hard” asymmetry side of the channel, while a negative magnetic field causes the skyrmion to interact with the “easy” asymmetry side. This geometry could be used to create new kinds of magnetic-field-induced diode effects that can be harnessed in new types of skyrmion-based devices.

36 MATERIALS SCIENCE↗

DOE OSTI2021

The PENELOPE Physics Models and Transport Mechanics. Implementation into Geant4

A translation of the PENELOPE physics subroutines to C++, designed as an extension of the GEANT4 toolkit, is presented. The Fortran code system PENELOPE performs Monte Carlo simulation of coupled electron-photon transport in arbitrary materials for a wide energy range, nominally from 50 eV up to 1 GeV. PENELOPE implements the most reliable interaction models that are currently available, limited only by the required generality of the code. In addition, the transport of electrons and positrons is simulated by means of an elaborate class II scheme in which hard interactions (involving deflection angles or energy transfers larger than pre-defined cutoffs) are simulated from the associated restricted differential cross sections. After a brief description of the interaction models adopted for photons and electrons/positrons, we describe the details of the class-II algorithm used for tracking electrons and positrons. The C++ classes are adapted to the specific code structure of GEANT4. They provide a complete description of the interactions and transport mechanics of electrons/positrons and photons in arbitrary materials, which can be activated from the G4ProcessManager to produce simulation results equivalent to those from the original PENELOPE programs. The combined code, named PENG4, benefits from the multi-threading capabilities and advanced geometry and statistical tools of GEANT4.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1983

The prediction of the noise of supersonic propellers in time domain - New theoretical results

In this paper, a new formula for the prediction of the noise of supersonic propellers is derived in the time domain which is superior to the previous formulations in several respects. The governing equation is based on the Ffowcs Williams-Hawkings (FW-H) equation with the thickness source term replaced by an equivalent loading source term derived by Isom (1975). Using some results of generalized function theory and simple four-dimensional space-time geometry, the formal solution of the governing equation is manipulated to a form requiring only the knowledge of blade surface pressure data and geometry. The final form of the main result of this paper consists of some surface and line integrals. The surface integrals depend on the surface pressure, time rate of change of surface pressure, and surface pressure gradient. These integrals also involve blade surface curvatures. The line integrals which depend on local surface pressure are along the trailing edge, the shock traces on the blade, and the perimeter of the airfoil section at the inner radius of the blade. The new formulation is for the full blade surface and does not involve any numerical observer time differentiation. The method of implementation on a computer for numerical work is also discussed.

Farassat, F.↗

NASA NTRS1983

Finite Element Analysis of a Deployable Space Structure

The dynamic characteristics of the Solar Array Flight Experiment (SAFE) structure during deployment and retraction are investigated. The SAFE structure consists of a deployable mast with an attached solar blanket designed with accordion type folds to permit packaging in a small volume. The planar form of the blanket geometry during deployment is maintained by a blanket tension/guidewire system. Structurally, the mast is modeled as an Euler beam column with inplane and out of plane bending and finite torsional stiffness. For out of plane motion, the blanket is modeled as a distributed mass uniformly supported by the three guidewires. For inplane motion the blanket displacements are assumed to vary linearly from the mast base to the mast tip. The mathematical model uses a virtual work formulation, required because the axial loading on the mast is nonconservative, combined with assumed beam modes to derive the differential equations of motion. Consideration of the time dependent boundary conditions results in an infinite set of ODE with time dependent coefficients. Finally, correlation of mast tip accelerations to mast base bending moments for specified modal motions are indicated.

Weeks, G. E.↗

NASA NTRS1991

Computational attributes of the integral form of the equation of transfer

Difficulties can arise in radiative and neutron transport calculations when a highly anisotropic scattering phase function is present. In the presence of anisotropy, currently used numerical solutions are based on the integro-differential form of the linearized Boltzmann transport equation. This paper, departs from classical thought and presents an alternative numerical approach based on application of the integral form of the transport equation. Use of the integral formalism facilitates the following steps: a reduction in dimensionality of the system prior to discretization, the use of symbolic manipulation to augment the computational procedure, and the direct determination of key physical quantities which are derivable through the various Legendre moments of the intensity. The approach is developed in the context of radiative heat transfer in a plane-parallel geometry, and results are presented and compared with existing benchmark solutions. Encouraging results are presented to illustrate the potential of the integral formalism for computation. The integral formalism appears to possess several computational attributes which are well-suited to radiative and neutron transport calculations.

Frankel, J. I.↗

NASA NTRSDate not supplied

Progress on Inverse Estimation Technique of Non-Linear Pitch Damping Coefficient Curves Using Free-Flight CFD Generated Trajectories

Characterization of entry vehicle pitch damping coefficient curves is crucial to ensure appropriate re-entry and overall mission success. The pitch damping coefficient (C_(m_q )+C_(m_α ̇ )) is used to encapsulate the oscillatory growth or decay of a body during a trajectory. The inverse estimation technique utilizes an existing Free-Flight CFD (FF-CFD) dataset and wraps a reconstruction algorithm in an optimizer. The reconstruction integrates the planar equations of motion derived by Schoenenberger, Queen [1] using Python’s scipy.integrate.solve_ivp. The optimizer’s objective function is the normalized 𝐿2 residual of the angle of attack peaks between the reconstructed trajectory and the original data produced with FF-CFD. Inclusion of the peak times in this residual calculation allows for simultaneous optimization of the pitch moment coefficient, C_(m_α ). This residual equation is shown below in Eq. 1. The optimizer scipy.optimize.minimize was used with the gradient-based Powell method for the analysis presented, however the differential evolution method was investigated as means of comparison, and was found to produce marginally lower residual values with prohibitively longer run times. Further, the pitch damping curve is found by fitting a cubic interpolation function to a set of (α, (C_(m_q )+C_(m_α ̇ ))) control points, where the α points are held constant and the (C_(m_q )+C_(m_α ̇ )) values are the optimized parameters. The pitch moment curve uses a linear interpolation between the minimum and maximum α in the dataset. FF-CFD generated trajectories using the Dragonfly capsule geometry with the Genesis ballistic range model parameters were simulated and used for this analysis. These FF-CFD trajectories simulate planar motion, as restricted by the reconstructing the equations of motion, of three different cases: 1-DoF (free-to-pitch), 2-DoF (free-to-pitch and heave), and 3-DoF (free-to-pitch, heave, and decelerate). Pitch damping coefficient curves generated using this inverse estimation curve technique with FF-CFD 1-DoF Dragonfly data are found in Fig. 1. Preliminary results reconstructing ballistic range shots using these FF-CFD derived predictions of the pitch damping curve (Fig. 1) are shown in Fig. 2. It should be noted that the ballistic range shot used a Genesis model whereas the FF-CFD data used a Dragonfly geometry, however these geometries are similar.

entry↗

NASA NTRS1997

Aerodynamic Shape Sensitivity Analysis and Design Optimization of Complex Configurations Using Unstructured Grids

A three-dimensional unstructured grid approach to aerodynamic shape sensitivity analysis and design optimization has been developed and is extended to model geometrically complex configurations. The advantage of unstructured grids (when compared with a structured-grid approach) is their inherent ability to discretize irregularly shaped domains with greater efficiency and less effort. Hence, this approach is ideally suited for geometrically complex configurations of practical interest. In this work the nonlinear Euler equations are solved using an upwind, cell-centered, finite-volume scheme. The discrete, linearized systems which result from this scheme are solved iteratively by a preconditioned conjugate-gradient-like algorithm known as GMRES for the two-dimensional geometry and a Gauss-Seidel algorithm for the three-dimensional; similar procedures are used to solve the accompanying linear aerodynamic sensitivity equations in incremental iterative form. As shown, this particular form of the sensitivity equation makes large-scale gradient-based aerodynamic optimization possible by taking advantage of memory efficient methods to construct exact Jacobian matrix-vector products. Simple parameterization techniques are utilized for demonstrative purposes. Once the surface has been deformed, the unstructured grid is adapted by considering the mesh as a system of interconnected springs. Grid sensitivities are obtained by differentiating the surface parameterization and the grid adaptation algorithms with ADIFOR (which is an advanced automatic-differentiation software tool). To demonstrate the ability of this procedure to analyze and design complex configurations of practical interest, the sensitivity analysis and shape optimization has been performed for a two-dimensional high-lift multielement airfoil and for a three-dimensional Boeing 747-200 aircraft.

Taylor, Arthur C., III↗

DOE OSTI2022

Slow manifold reduction as a systematic tool for revealing the geometry of phase space

Many non-dissipative reduced plasma models can be derived from more fundamental non-dissipative models by restricting to an approximate invariant manifold. I present a general systematic procedure for finding the Hamiltonian formulation of a plasma model that can be derived in this manner. Several illustrative examples are considered in detail.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Tip-based proximity ferroelectric switching and piezoelectric response in wurtzite multilayers

Proximity ferroelectricity is a paradigm for inducing ferroelectricity when a nonferroelectric polar material (such as Al⁢ N), which is unswitchable with an external field below the dielectric breakdown field, becomes a practically switchable ferroelectric in direct contact with a thin switchable ferroelectric layer (such as Al 1−𝑥 ⁢Sc 𝑥 ⁢N). Here, we develop a Landau-Ginzburg-Devonshire approach to study the proximity effect of local piezoelectric response and polarization reversal in wurtzite ferroelectric multilayers under a sharp electrically biased tip. Using finite-element modeling, we analyze the probe-induced nucleation of nanodomains, the features of local polarization hysteresis loops and coercive fields in the Al 1−𝑥 ⁢Sc 𝑥 ⁢N/Al⁢ N bilayers and three-layers. Similar to the wurtzite multilayers sandwiched between two parallel electrodes, the regimes of “proximity switching” (when all layers collectively switch) and the regime of “proximity suppression” (when they collectively do not switch) are the only two possible regimes in the probe-electrode geometry. However, the parameters and asymmetry of the local piezoresponse and polarization hysteresis loops depend significantly on the sequence of the layers with respect to the probe. The physical mechanism of proximity ferroelectricity in the local probe geometry is a depolarizing electric field determined by the polarization of the layers and their relative thickness. The field, whose direction is opposite to the polarization vector in the layer(s) with the larger spontaneous polarization (such as Al⁢ N), renormalizes the double-well ferroelectric potential to lower the steepness of the switching barrier in the “otherwise unswitchable” polar layers. Tip-based control of domains in otherwise nonferroelectric layers using proximity ferroelectricity can provide nanoscale control of domain reversal in memory, actuation, sensing, and optical applications. The ability of the tip-induced proximity switching to differentially switch multilayers, based on the order of the layers, provides a powerful tool for selective domain engineering.

36 MATERIALS SCIENCE↗

NASA NTRS1987

Expansion of epicyclic gear dynamic analysis program

The multiple mesh/single stage dynamics program is a gear tooth analysis program which determines detailed geometry, dynamic loads, stresses, and surface damage factors. The program can analyze a variety of both epicyclic and single mesh systems with spur or helical gear teeth including internal, external, and buttress tooth forms. The modifications refine the options for the flexible carrier and flexible ring gear rim and adds three options: a floating Sun gear option; a natural frequency option; and a finite element compliance formulation for helical gear teeth. The option for a floating Sun incorporates two additional degrees of freedom at the Sun center. The natural frequency option evaluates the frequencies of planetary, star, or differential systems as well as the effect of additional springs at the Sun center and those due to a flexible carrier and/or ring gear rim. The helical tooth pair finite element calculated compliance is obtained from an automated element breakup of the helical teeth and then is used with the basic gear dynamic solution and stress postprocessing routines. The flexible carrier or ring gear rim option for planetary and star spur gear systems allows the output torque per carrier and ring gear rim segment to vary based on the dynamic response of the entire system, while the total output torque remains constant.

Boyd, Linda Smith↗

NASA NTRS1992

Is the Oort A-value a universal growth rate limit for accretion disk shear instabilities?

A weak-field local MHD instability that is of importance to accretion disks is examined. The maximum growth rate of the instability is found to be not only independent of the magnetic field strength but independent of field geometry as well. In particular, all Keplerian disks are unstable in the presence of any weak poloidal field, with the ratio of the maximum growth rate to disk angular velocity given by 3/4. The maximum growth rate of any weak field configuration that is not purely toroidal is given by the local Oort A-value of the disk. The behavior is studied by using a form of the dynamical Hill equations. It is conjectured that the Oort A-value is an upper bound to the growth rate of any instability feeding upon the free energy of differential rotation.

Balbus, Steven A.↗

NASA NTRS1993

Performance evaluation of the JPL TOPSAR system: An across-track interferometric SAR system for topographic mapping

An across-track interferometric (XTI) SAR has the capability of determining the three dimensional target location from the recorded radar echos thus providing high spatial resolution both horizontally and vertically. The output of an XTI SAR will typically include both a radar backscatter map and a digital elevation model (DEM). The accuracy of the derived DEM is determined by factors such as data acquisition geometry, signal-to-noise ratio, knowledge of platform position and attitude as well as the accuracy of the processing system. In the summer of 1992 the authors acquired TOPSAR data over the Ft. Irwin area in California, a desert area with significant relief (height standard deviation 150 m). Very accurate DEM's have been derived for this area by the Topographic Engineering Center (TEC) using digital correlation methods on 1:20,000 scale digitized photographs. Corner reflectors were deployed in the area, and their locations were determined to cm accuracies by the Defense Mapping Agency (DMA) using differential GPS techniques. DEM's generated from the acquired radar data have been rotated and translated to overlay the reference DEM's provided by TEC. A detailed description of the errors and their characteristics will be given. The standard deviations of the radar data sets we acquired measured over a 5.6 by 7 km area were 1.9 and 2.3 m. The corresponding numbers for flat areas were 1.1 and 2.0 m and for mountain areas 3.3 and 2.2 m.

Madsen, Soren N.↗

DOE OSTI2025

Computational Modeling of Graphite Degradation due to Molten Salt Infiltration and Wear

Molten-salt reactors (MSRs) represent a promising next-generation reactor design, with graphite serving as a moderator and/or reflector in several designs. However, due to limited experimental data and operational experience, a technical understanding of the structural integrity of graphite in molten salt environments remains incomplete. This report presents a modeling-based evaluation of graphite degradation in MSR environments, focusing on the effects of salt infiltration in fuel salt-based designs and surface wear in pebble bed reactor designs. The objective of this study is to enhance understanding of the structural integrity challenges posed by these degradation mechanisms and to provide a framework for assessing graphite behavior in MSRs. The first part of the report investigates the phenomenon of molten salt infiltration into graphite. This infiltration occurs when molten salt permeates the interconnected pore structure of the graphite moderator, driven by factors such as pressure differentials and the physical properties of both the salt and graphite. The infiltration process is influenced by characteristics of the pore structure, viscosity of the molten salt, and the interfacial energies between the graphite, salt, and the atmosphere within the graphite pore. Utilizing a coupled multiphysics modeling approach with Grizzly software, the study evaluates the stress induced by internal heat sources due to infiltration, which can lead to structural concerns. This evaluation is crucial for understanding how infiltration affects the mechanical integrity of graphite components in MSRs. The study considers the Molten-Salt Reactor Experiment (MSRE) graphite stringer geometry due to the availability of relevant data. Through detailed finite element analysis, the study examines stress distributions at varying infiltration percentages, revealing that stress levels increase with higher amounts of infiltration. Rare-event simulations, using the parallel subset simulation (PSS) framework, further quantify the failure probabilities under input uncertainties, with a user-specified failure metric. The PSS framework also identifies critical input parameters that significantly affect the stress values, including infiltration amount, thermal conductivity, and power density. Additionally, considering realistic reactor scenarios, the analysis was performed to account for the combined effects of radiation and infiltration, and modeling strategies on how to analyze new reactor designs or new graphite grades are discussed. The second part of the report focuses on wear mechanisms in pebble bed-based MSRs. As graphite fuel pebbles interact with the graphite reflector block, wear can result in material loss and the formation of surface defects, which may act as stress concentrators. A similar multiphysics modeling framework is employed to assess the impact of wear on the structural integrity of graphite components. This study considers a generic fluoride-cooled high-temperature reactor (gFHR) design due to the availability of comprehensive data. Worst-case scenario dimensions of the reflector blocks were analyzed under thermal and radiation conditions. Subsequently, wear in the form of idealized pits and grooves is modeled on the inner surface of the graphite block, with the maximum stress from previous simulations. The simulations show that groove-type defects are more detrimental than pits, leading to higher stress concentrations. Considering worst-case simulation scenarios and experimental wear rates, it was determined that the formation of a surface defect critical enough to affect the stress may not be possible in a gFHR design. Overall, the findings of this research contribute to the development of robust modeling tools for predicting graphite behavior under various operational conditions in MSRs.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS1992

Evaluation and application of the Baldwin-Lomax turbulence model in two-dimensional, unsteady, compressible boundary layers with and without separation in engine inlets

There is a practical need to model high speed flows that exist in jet engine inlets. The boundary layers that form in these inlets may be turbulent or laminar and either separated or attached. Also, unsteady supersonic inlets may be subject to frequent changes in operating conditions. Some changes in the operating conditions of the inlets may include varying the inlet geometry, bleeds and bypasses, and rotating or translating the centerbody. In addition, the inlet may be either started or unstarted. Therefore, a CFD code, used to model these inlets, may have to run for several different cases. Also, since the flow conditions through an unsteady inlet may be continually fluctuating, the CFD code which models these flows may have to be run over many time steps. Therefore, it would be beneficial that the code run quickly. Many turbulence models, however, are cumbersome to implement and require a lot of computer time to run, since they add to the number of differential equations to be solved to model a flow. The Baldwin-Lomax turbulence model is a popular model. It is an algebraic, eddy viscosity model. The Baldwin-Lomax model is used in many CFD codes because it is quick and easy to implement. In this paper, we will discuss implementing the Baldwin-Lomax turbulence model for both steady and unsteady compressible flows. In addition, these flows may be either separated or attached. In order to apply this turbulence model to flows which may be subjected to these conditions, certain modifications should be made to the original Baldwin-Lomax model. We will discuss these modifications and determine whether the Baldwin-Lomax model is a viable turbulence model that produces reasonably accurate results for high speed flows that can be found in engine inlets.

Sakowski, Barbara↗

NASA NTRS1993

Two wavelength satellite laser ranging using SPAD

When ranging to satellites with lasers, there are several principal contributions to the error budget: from the laser ranging system on the ground, from the satellite retroarray geometry, and from the atmosphere. Using a single wavelength, we have routinely achieved a ranging precision of 8 millimeters when ranging to the ERS-1 and Starlette satellites. The systematic error of the atmosphere, assuming the existing dispersion models, is expected to be of the order of 1 cm. Multiple wavelengths ranging might contribute to the refinement of the existing models. Taking into account the energy balance, the existing picosecond lasers and the existing receiver and detection technology, several pairs or multiple wavelengths may be considered. To be able to improve the atmospheric models to the subcentimeter accuracy level, the differential time interval (DTI) has to be determined within a few picoseconds depending on the selected wavelength pair. There exist several projects based on picosecond lasers as transmitters and on two types of detection techniques: one is based on photodetectors, like photomultipliers or photodiodes connected to the time interval meters. Another technique is based on the use of a streak camera as an echo signal detector, temporal analyzer, and time interval vernier. The temporal analysis at a single wavelength using the streak camera showed the complexity of the problem.

Prochazka, Ivan↗

NASA NTRS1997

Constrained Multipoint Aerodynamic Shape Optimization Using an Adjoint Formulation and Parallel Computers

An aerodynamic shape optimization method that treats the design of complex aircraft configurations subject to high fidelity computational fluid dynamics (CFD), geometric constraints and multiple design points is described. The design process will be greatly accelerated through the use of both control theory and distributed memory computer architectures. Control theory is employed to derive the adjoint differential equations whose solution allows for the evaluation of design gradient information at a fraction of the computational cost required by previous design methods. The resulting problem is implemented on parallel distributed memory architectures using a domain decomposition approach, an optimized communication schedule, and the MPI (Message Passing Interface) standard for portability and efficiency. The final result achieves very rapid aerodynamic design based on a higher order CFD method. In order to facilitate the integration of these high fidelity CFD approaches into future multi-disciplinary optimization (NW) applications, new methods must be developed which are capable of simultaneously addressing complex geometries, multiple objective functions, and geometric design constraints. In our earlier studies, we coupled the adjoint based design formulations with unconstrained optimization algorithms and showed that the approach was effective for the aerodynamic design of airfoils, wings, wing-bodies, and complex aircraft configurations. In many of the results presented in these earlier works, geometric constraints were satisfied either by a projection into feasible space or by posing the design space parameterization such that it automatically satisfied constraints. Furthermore, with the exception of reference 9 where the second author initially explored the use of multipoint design in conjunction with adjoint formulations, our earlier works have focused on single point design efforts. Here we demonstrate that the same methodology may be extended to treat complete configuration designs subject to multiple design points and geometric constraints. Examples are presented for both transonic and supersonic configurations ranging from wing alone designs to complex configuration designs involving wing, fuselage, nacelles and pylons.

Reuther, James↗