Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical differentiation”

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 1,297 records · Page 72
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Plastic work constrained elastoplastic topology optimization

An elastoplastic topology optimization framework for limiting plastic work generation while maximizing stiffness is presented. The kinematics and constitutive model are based on finite strain linear isotropic hardening plasticity, and the balance laws are solved using a total Lagrangian finite element formulation. Aggregation of the specific plastic work combined with an adaptive normalization scheme efficiently constrains the maximum specific plastic work. The optimization problem is regularized using an augmented partial differential equation filter, and is solved by the method of moving asymptotes where path-dependent sensitivities are derived using the adjoint method. The numerical examples show a clear dependence on the optimized maximum stiffness structures for different levels of constrained specific plastic work. It is also shown that due to the history dependency of the plasticity, the load path significantly influences the structural performance and optimized topology.

42 ENGINEERING↗
DOE OSTI2025

Fast and Scalable FFT-Based GPU-Accelerated Algorithms for Block-Triangular Toeplitz Matrices with Application to Linear Inverse Problems Governed by Autonomous Dynamical Systems

In this work, we present an efficient and scalable algorithm for performing matrix-vector multiplications (matvecs) for block Toeplitz matrices. Such matrices, which are shift-invariant with respect to their blocks, arise in the context of solving inverse problems governed by autonomous systems, and time-invariant systems in particular. In this article, we consider inverse problems that infer unknown parameters from observational data of a linear time-invariant dynamical system given in the form of partial differential equations (PDEs). Matrix-free Newton-conjugate-gradient methods are often the gold standard for solving these inverse problems, but they require numerous actions of the Hessian on a vector. Matrix-free adjoint-based Hessian matvecs require solution of a pair of linearized forward/adjoint PDE solves per Hessian action, which may be prohibitive for large-scale inverse problems. Time invariance of the forward PDE problem leads to a block Toeplitz structure of the discretized parameter-to-observable (p2o) map defining the mapping from inputs (parameters) to outputs (observables) of the PDEs. This block Toeplitz structure enables us to exploit two key properties: (1) compact storage of the p2o map and its adjoint, and (2) efficient fast Fourier transform–based Hessian matvecs. The proposed algorithm is mapped onto large multi-GPU clusters and achieves more than 80% of peak bandwidth on NVIDIA A100 GPUs. Excellent weak scaling is shown for up to 48 A100 GPUs. For the targeted problems, the implementation executes Hessian matvecs within fractions of a second, which is orders of magnitude faster than can be achieved by conventional matrix-free Hessian matvecs via forward/adjoint PDE solves.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

A User Guide to PARET/ANL

PARET was originally created in 1969 at what is now Idaho National Laboratory (INL), to analyze reactivity insertion events in research and test reactor cores cooled by light or heavy water, with fuel composed of either plates or pins. The use of PARET is also appropriate for fuel assemblies with curved fuel plates when their radii of curvatures are large with respect to the fuel plate thickness. The PARET/ANL version of the code has been developed at Argonne National Laboratory (ANL) under the sponsorship of the U.S. Department of Energy/NNSA since the inception of the Reactor Conversion Program. Since Reduced Enrichment for Research and Test Reactors (RERTR) began in 1978, PARET/ANL has been benchmarked to experimental data including SPERT testing, and used to determine the expected transient behavior of a large number of reactors regardless of enrichment both inside and outside the Reactor Conversion Program. This document provides the pertinent information for the use of PARET/ANL Version 7.6. PARET/ANL models the various fueled regions of a reactor core as channels. Each of these channels consists of a single flat fuel plate/pin (including cladding and, optionally, a gap) with water coolant on each side. In slab geometry the coolant channels for a given fuel plate are of identical dimensions (mirror symmetry), but they can be of different thickness in each channel. There can be many channels, but each channel is independent and coupled only through reactivity feedback effects to the whole core. The time-dependent differential equations that represent the system are replaced by an equivalent set of finite-difference equations in space and time, which are integrated numerically. PARET/ANL uses fundamentally the same numerical scheme as RELAP5 for the time-integration of the point-kinetics equations. The one-dimensional thermal-hydraulic model includes temperature-dependent thermal properties of the solid materials, such as heat capacity and thermal conductivity, as well as the transient heat production and heat transfer from the fuel meat to the coolant. Temperature- and pressure-dependent thermal properties of the coolant such as enthalpy, density, thermal conductivity, and viscosity are also used in determining parameters such as friction factors and heat transfer coefficients. The code first determines the steady-state solution for the initial state. Then the solution of the transient is obtained by integration in time and space. Multiple heat transfer, DNB and flow instability correlations are available. The code was originally developed to model reactors cooled by an open loop, which was adequate for rapid transients in pool-type cores. An external loop model appropriate for Miniature Neutron Source Reactors (MNSR’s) was also added to PARET/ANL to model natural circulation within the vessel, heat transfer from the vessel to pool and heat loss by evaporation from the pool. PARET/ANL also contains models for decay heat after shutdown, control rod reactivity versus time or position, time-dependent pump flow, and loss-of-flow event with flow reversal as well as logic for trips on period, power, and flow. Feedback reactivity effects from coolant density changes and temperature changes are represented by tables. Feedback reactivity from fuel heat-up (Doppler Effect) is represented by a four-term polynomial in powers of fuel temperature. Photo-neutrons produced in beryllium or in heavy water may be included in the point-kinetics equations by using additional delayed neutron groups.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1971

Nutational stability of a dual-spin satellite under the influence of applied reaction torques

It is assumed that since the solar paddle attachments to the hub of the spacecraft are not rigidly locked, the effect of the solar panels can be replaced by a constant reaction torque acting on the hub of the spacecraft. This could result in the satellite having an equilibrium motion about an axis displaced from the nominal axis of symmetry. The variational equations of motion are developed about such an equilibrium position using the SAS-A spacecraft as a model. Energy dissipation on the rotor as well as the main body is included. This nonautonomous set of differential equations are linearized and transformed to an autonomous set using the Liapunov reducibility theorem. The stability of the kinematically similar system is examined numerically using representative SAS-A parameters for the case when either pair of solar panels is assumed to be loosely attached. Stability is verified for small system nutation angles (0.1 degree).

Bainum, P. M.↗
NASA NTRS1972

Nutational stability of a dual-spin satellite under the influence of applied reaction torques.

The nutational stability of a dual-spin satellite under the influence of an applied reaction torque is analyzed. It is assumed that since the solar paddle attachments to the hub of the spacecraft are not rigidly locked, the effect of the solar panels can be replaced by a constant reaction torque acting on the hub of the spacecraft. This could result in the satellite having an equilibrium motion about an axis displaced from the nominal axis of symmetry. The variational equations of motion are developed about such an equilibrium position using the SAS-A spacecraft as a model. Energy dissipation on the rotor as well as the main body is included. This nonautonomous set of differential equations are linearized and transformed to an autonomous set using the Lyapunov Reducibility Theorem. The stability of the kinematically similar system is examined numerically using representative SAS-A parameters for the case when either pair of solar panels is assumed to be loosely attached. Stability is verified for small system nutation angles (0.1 degree) although the time constant associated with the least damped mode is approximately one order of magnitude larger than that for the SAS-A system without the influence of the reaction torques.

Bainum, P. M.↗
NASA NTRS1977

Development of prediction techniques for multi-jet thermal ground flow fields and fountain formation

Avoiding detrimental ground interaction is important for practical V/STOL aircraft. This paper reports recent developments in a numerical method for estimating thermal ground footprints. Upwash and fountain formation for arbitrarily oriented jet arrangements is predicted. Flow asymmetry due to roll, pitch, differential thrust or ground inclination is included. The prediction methodology uses simple inviscid relations for energy and momentum conservation along with an empirical entrainment law, applied in independent sectors of the wall jet and upwash. Asymmetrical stagnation line prediction is compared with experiment. Detailed flow measurements for a three-jet interaction are also presented.

Siclari, M. J.↗
NASA NTRS1980

Computation of three-dimensional viscous supersonic flow in inlets

A new approach has been developed for the computation of the three-dimensional viscous supersonic flow with embedded subsonic regions adjacent to solid boundaries and is applied to a mixed-compression supersonic inlet typical of current designs. The approach uses a reduced form of the three-dimensional Navier-Stokes equations so that the resultant equations can be treated as an initial boundary value problem and thus be solved by non-iterative forward marching in space. The numerical procedure utilizes an efficient consistently-split linearized block implicit technique to solve the finite difference analogues to the set of governing partial differential equations.

Buggeln, R. C.↗
NASA NTRS1988

Applications of an exponential finite difference technique

An exponential finite difference scheme first presented by Bhattacharya for one dimensional unsteady heat conduction problems in Cartesian coordinates was extended. The finite difference algorithm developed was used to solve the unsteady diffusion equation in one dimensional cylindrical coordinates and was applied to two and three dimensional conduction problems in Cartesian coordinates. Heat conduction involving variable thermal conductivity was also investigated. The method was used to solve nonlinear partial differential equations in one and two dimensional Cartesian coordinates. Predicted results are compared to exact solutions where available or to results obtained by other numerical methods.

Handschuh, Robert F.↗
NASA NTRS1990

Fourier analysis of finite element preconditioned collocation schemes

The spectrum of the iteration operator of some finite element preconditioned Fourier collocation schemes is investigated. The first part of the paper analyses one-dimensional elliptic and hyperbolic model problems and the advection-diffusion equation. Analytical expressions of the eigenvalues are obtained with use of symbolic computation. The second part of the paper considers the set of one-dimensional differential equations resulting from Fourier analysis (in the tranverse direction) of the 2-D Stokes problem. All results agree with previous conclusions on the numerical efficiency of finite element preconditioning schemes.

Deville, Michel O.↗
NASA NTRS1993

Geology of Galileo Regio quadrangle, Ganymede

Galileo Regio quadrangle, including most of Galileo Regio, contains two basic geologic units: older dark, furrowed, and heavily cratered material and younger light, grooved, relatively less-cratered material. Dark material dominates, the light material occurring to the W and SW in about five percent of the quadrangle. Dominant structural features, of uncertain origin, are the NW-trending arcuate furrows and associated orthogonal and oblique furrows in the dark material; grooves in the light material are not as well-ordered nor as intensely developed as elsewhere on Ganymede. Fifteen palimpsests, located partly or totally in the quadrangle and representing the period when the crust became sufficiently strong to record impact events, are subdivided into three relative-age classes, as are the numerous impact craters, some of which are dome, moat, and rampart-like. Galileo Regio is believed to have formed when, because of cooling, differentiation, and expansion, the thin, weak, icy lithosphere broke into blocks and light material filled the space between, or overspread those blocks that subsided, or both. Grooves formed in most of the light material and both light and dark material were cratered by impact, the most recent of which created bright-ray craters, the youngest features visible in the quadrangle.

Underwood, James R., Jr.↗
NASA NTRS1993

Approximate optimal guidance for the advanced launch system

A real-time guidance scheme for the problem of maximizing the payload into orbit subject to the equations of motion for a rocket over a spherical, non-rotating earth is presented. An approximate optimal launch guidance law is developed based upon an asymptotic expansion of the Hamilton - Jacobi - Bellman or dynamic programming equation. The expansion is performed in terms of a small parameter, which is used to separate the dynamics of the problem into primary and perturbation dynamics. For the zeroth-order problem the small parameter is set to zero and a closed-form solution to the zeroth-order expansion term of Hamilton - Jacobi - Bellman equation is obtained. Higher-order terms of the expansion include the effects of the neglected perturbation dynamics. These higher-order terms are determined from the solution of first-order linear partial differential equations requiring only the evaluation of quadratures. This technique is preferred as a real-time, on-line guidance scheme to alternative numerical iterative optimization schemes because of the unreliable convergence properties of these iterative guidance schemes and because the quadratures needed for the approximate optimal guidance law can be performed rapidly and by parallel processing. Even if the approximate solution is not nearly optimal, when using this technique the zeroth-order solution always provides a path which satisfies the terminal constraints. Results for two-degree-of-freedom simulations are presented for the simplified problem of flight in the equatorial plane and compared to the guidance scheme generated by the shooting method which is an iterative second-order technique.

Feeley, T. S.↗
NASA NTRS1995

Effects of pressure gradients on turbulent premixed flames

The influence of a constant acceleration on a turbulent premixed flame is studied by direct numerical simulation. This acceleration induces a mean pressure gradient across the flame brush, leading to a modification of the turbulent flame structure due to differential buoyancy mechanisms between heavy cold fresh and light hot burnt gases. Such a pressure gradient may be encountered in practical applications in ducted flames. A favorable pressure gradient, i.e. the pressure decreases from unburnt to burnt gases, is found to decrease the flame wrinkling, the flame brush thickness, and the turbulent flame speed. A favorable pressure gradient also promotes counter-gradient turbulent transport. On the other hand, adverse pressure gradients tend to increase the flame brush thickness and turbulent flame speed, and promote classical gradient turbulent transport. The balance equation for the turbulent flux of the Favre averaged progress variable is also analyzed. The first results show that the fluctuating pressure term, cannot be neglected as generally assumed in models. Simple models assuming that a high mean pressure gradient may only be balanced by the cross-dissipation term seem too approximate. This analysis has to be continued to compare simulation data and closure schemes proposed for the transport equation. The analysis developed by Veynante et al.(1995) has been extended to imposed acceleration and mean pressure gradients. A simple model for the turbulent flux is proposed and validated from simulation data. Then, a modified criterion is derived to delineate between counter-gradient and gradient turbulent diffusion. In fact, counter-gradient diffusion may occur in most practical applications, especially for ducted flames.

Veynante, D.↗
NASA NTRS1996

Approach to Synchronization Control of Magnetic Bearings Using Fuzzy Logic

This paper presents a fuzzy-logic approach to the synthesis of synchronization control for magnetically suspended rotor system. The synchronization control enables a whirling rotor to undergo synchronous motion along the magnetic bearing axes; thereby avoiding the gyroscopic effect that degrade the stability of rotor systems when spinning at high speed. The control system features a fuzzy controller acting on the magnetic bearing device, in which the fuzzy inference system trained through fuzzy rules to minimize the differential errors between four bearing axes so that an error along one bearing axis can affect the overall control loop for the motion synchronization. Numerical simulations of synchronization control for the magnetically suspended rotor system are presented to show the effectiveness of the present approach.

Yang, Li-Farn↗
NASA NTRS2000

A Dream of a Mission: Stellar Imager and Seismic Probe

The Stellar Imager and Seismic Probe (SISP) is a mission to understand the various effects of magnetic fields of stars, the dynamos that generate them, and the internal structure and dynamics of the stars in which they exist. The ultimate goal is to achieve the best-possible forecasting of solar activity on times scales ranging up to decades, and an understanding of the impact of stellar magnetic activity on astrobiology and life in the Universe. The road to that goal will revolutionize our understanding of stars and stellar systems, the building blocks of the Universe. SISP will zoom in on what today - with few exceptions - we only know as point sources, revealing processes never before seen, thus providing a tool to astrophysics as fundamental as the microscope is to the study of life on Earth. SISP is an ultraviolet aperture-synthesis imager with 8-10 telescopes with meter-class apertures, and a central hub with focal-plane instrumentation that allows spectrophotometry in passbands as narrow as a few Angstroms up to hundreds of Angstroms. SISP will image stars and binaries with one hundred to one thousand resolution elements on their surface, and sound their interiors through asteroseismology to image internal structure, differential rotation, and large-scale circulations; this will provide accurate knowledge of stellar structure and evolution and complex transport processes, and will impact numerous branches of (astro)physics ranging from the Big Bang to the future of the Universe. Fitting naturally within the NASA long-term time line, SISP complements defined missions, and with them will show us entire other solar systems, from the central star to their orbiting planets.

Carpenter, Kenneth G.↗
NASA NTRS1956

On the Theory of Anisotropic Shallow Shells

A numerical analysis of thin-walled shallow shells is presented. Equations of equilibrium and relations between deformations and stresses are included along with fundamental differential equations.

Ambartsumyan, S. A.↗
NASA NTRS2001

Auxins and tropisms

Differential growth of plants in response to the changes in the light and gravity vectors requires a complex signal transduction cascade. Although many of the details of the mechanisms by which these differential growth responses are induced are as yet unknown, auxin has been implicated in both gravitropism and phototropism. Specifically, the redistribution of auxin across gravity or light-stimulated tissues has been detected and shown to be required for this process. The approaches by which auxin has been implicated in tropisms include isolation of mutants altered in auxin transport or response with altered gravitropic or phototropic response, identification of auxin gradients with radiolabeled auxin and auxin-inducible gene reporter systems, and by use of inhibitors of auxin transport that block gravitropism and phototropism. Proteins that transport auxin have been identified and the mechanisms which determine auxin transport polarity have been explored. In addition, recent evidence that reversible protein phosphorylation controls this process is summarized. Finally, the data in support of several hypotheses for mechanisms by which auxin transport could be differentially regulated during gravitropism are examined. Although many details of the mechanisms by which plants respond to gravity and light are not yet clear, numerous recent studies demonstrate the role of auxin in these processes.

NASA Discipline Plant Biology↗
NASA NTRS2009

Developing Lightning Prediction Tools for the CCAFS Dual-Polarimetric Radar

NASA Marshall Space Flight Center and the University of Alabama Huntsville are collaborating with the 45th Weather Squadron (45WS) to develop improved lightning prediction capabilities for the new C-band dual-polarimetric weather radar being acquired for use by 45WS and launch weather forecasters at Cape Canaveral Air Force Station (CCAFS). In particular, these algorithms will focus on lightning onset, cessation and combined lightning-radar applications for convective winds assessment. Research using radar reflectivity (Z) data for prediction of lightning onset has been extensively discussed in the literature and subsequently applied by launch weather forecasters as it pertains to lightning nowcasting. Currently the forecasters apply a relatively straight forward but effective temperature-Z threshold algorithm for assessing the likelihood of lightning onset in a given storm. In addition, a layered VIL above the freezing level product is used as automated guidance for the onset of lightning. Only limited research and field work has been conducted on lightning cessation using Z and vertically-integrated Z for determining cessation. Though not used operationally vertically-integrated Z (basis for VIL) has recently shown promise as a tool for use in nowcasting lightning cessation. The work discussed herein leverages and expands upon these and similar reflectivity-threshold approaches via the application/addition of over two decades of polarimetric radar research focused on distinct multi-parameter radar signatures of ice/mixed-phase initiation and ice-crystal orientation in highly electrified convective clouds. Specifically, our approach is based on numerous previous studies that have observed repeatable patterns in the behavior of the vertical hydrometeor column as it relates to the temporal evolution of differential reflectivity and depolarization (manifested in either LDR or p(sub hv)), development of in-situ mixed and ice phase microphysics, electric fields, and ensuing lightning in the sub-tropical/tropical convection typical of the southeastern U.S., Maritime Continent, and southwestern Amazon. The polarimetric signatures detected in this setting provide a basis for automated 3-D detection of hydrometeor types in fuzzy logic hydrometeor identification algorithms (HID). Our working hypothesis is that improvement in lightning onset warning lead time and specificity for a given storm, relative to application of a Z-threshold algorithm, should arise as a consequence of the ability of dual-polarimetric radar to unambiguously detect and identify (through HID algorithms) the updraft elevation of rain-water cores above the freezing level and subsequent onset of drop freezing, riming, and robust mixed phase processes leading to significant charge separation and lightning. This type of algorithm, though dependent on the quality of the polarimetric data should be less susceptible to variable Z-calibration that can impact a given Z-threshold approach. To facilitate development of the algorithm while the 45WS dual-pol radar is in its current test stages and to evaluate the impact of polarimetric data quality (e.g., modified scan parameters and sampling) on the ensuing algorithms, we are using the ARMOR C-band dual-pol radar in Huntsville combined with N. Alabama LMA data and ARMOR HID algorithms [NCAR algorithm modified for application at C-band] in a testbed fashion. For lightning cessation we are revisiting the application of differential propagation phase variables for the monitoring of ice crystal alignment driven by in-cloud electric fields combined with metrics of ice water path (i.e., vertically integrated reflectivity). Importantly it should be noted that this approach is still very much a research topic and as such, we will explore operational applications that involve radar frequencies other than C-Band by using the UAH MAX X-band dual-pol radar in slow staring modes.

Petersen, W. A.↗