Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Arbitrary functions”

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,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
NASA NTRS2012

Discrete Adjoint-Based Design for Unsteady Turbulent Flows On Dynamic Overset Unstructured Grids

A discrete adjoint-based design methodology for unsteady turbulent flows on three-dimensional dynamic overset unstructured grids is formulated, implemented, and verified. The methodology supports both compressible and incompressible flows and is amenable to massively parallel computing environments. The approach provides a general framework for performing highly efficient and discretely consistent sensitivity analysis for problems involving arbitrary combinations of overset unstructured grids which may be static, undergoing rigid or deforming motions, or any combination thereof. General parent-child motions are also accommodated, and the accuracy of the implementation is established using an independent verification based on a complex-variable approach. The methodology is used to demonstrate aerodynamic optimizations of a wind turbine geometry, a biologically-inspired flapping wing, and a complex helicopter configuration subject to trimming constraints. The objective function for each problem is successfully reduced and all specified constraints are satisfied.

Nielsen, Eric J.↗
DOE OSTI2022

Black hole scattering and partition functions

When computing the ideal gas thermal canonical partition function for a scalar outside a black hole horizon, one encounters the divergent single-particle density of states (DOS) due to the continuous nature of the normal mode spectrum. Recasting the Lorentzian field equation into an effective 1D scattering problem, we argue that the scattering phases encode non-trivial information about the DOS and can be extracted by “renormalizing” the DOS with respect to a reference. This defines a renormalized free energy up to an arbitrary additive constant. Interestingly, we discover that the 1-loop Euclidean path integral, as computed by the Denef-Hartnoll-Sachdev formula, fixes the reference free energy to be that on a Rindler-like region, and the renormalized DOS captures the quasinormal modes for the scalar. We support these claims with the examples of scalars on static BTZ, Nariai black holes and the de Sitter static patch. For black holes in asymptotically flat space, the renormalized DOS is captured by the phase of the transmission coefficient whose magnitude squared is the greybody factor. We comment on possible connections with recent works from an algebraic point of view.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Embracing fine-root system complexity in terrestrial ecosystem modeling

Projecting the dynamics and functioning of the biosphere requires a holistic consideration of whole-ecosystem processes. However, biases toward leaf, canopy, and soil modeling since the 1970s have constantly left fine-root systems being rudimentarily treated. As accelerated empirical advances in the last two decades establish clearly functional differentiation conferred by the hierarchical structure of fine-root orders and associations with mycorrhizal fungi, a need emerges to embrace this complexity to bridge the data-model gap in still extremely uncertain models. In this work, we propose a three-pool structure comprising transport and absorptive fine roots with mycorrhizal fungi (TAM) to model vertically resolved fine-root systems across organizational and spatial–temporal scales. Emerging from a conceptual shift away from arbitrary homogenization, TAM builds upon theoretical and empirical foundations as an effective and efficient approximation that balances realism and simplicity. A proof-of-concept demonstration of TAM in a big-leaf model both conservatively and radically shows robust impacts of differentiation within fine-root systems on simulating carbon cycling in temperate forests. Theoretical and quantitative support warrants exploiting its rich potentials across ecosystems and models to confront uncertainties and challenges for a predictive understanding of the biosphere. Echoing a broad trend of embracing ecological complexity in integrative ecosystem modeling, TAM may offer a consistent framework where modelers and empiricists can work together toward this grand goal.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1974

Steady incompressible potential flow around lifting bodies immersed in a fluid

The refinement was investigated of a method for evaluating the pressure distribution on a body surface of arbitrary shape in incompressible flow. The solution was obtained in terms of the velocity potential, through numerical approximations which require the use of a high speed digital computer. The box method and the modal method are described in detail, and were applied to a very thin, rectangular wing in incompressible, steady flow. The box method is found to be more practical as it is applicable to more general geometries (the modal method requires a new set of functions for each geometry), and requires less computer time (fifty percent of that required by the modal method for the same problem).

Chiuchiolo, E. A.↗
NASA NTRS1996

Development of Finite Elements for Two-Dimensional Structural Analysis Using the Integrated Force Method

The Integrated Force Method has been developed in recent years for the analysis of structural mechanics problems. This method treats all independent internal forces as unknown variables that can be calculated by simultaneously imposing equations of equilibrium and compatibility conditions. In this paper a finite element library for analyzing two-dimensional problems by the Integrated Force Method is presented. Triangular- and quadrilateral-shaped elements capable of modeling arbitrary domain configurations are presented. The element equilibrium and flexibility matrices are derived by discretizing the expressions for potential and complementary energies, respectively. The displacement and stress fields within the finite elements are independently approximated. The displacement field is interpolated as it is in the standard displacement method, and the stress field is approximated by using complete polynomials of the correct order. A procedure that uses the definitions of stress components in terms of an Airy stress function is developed to derive the stress interpolation polynomials. Such derived stress fields identically satisfy the equations of equilibrium. Moreover, the resulting element matrices are insensitive to the orientation of local coordinate systems. A method is devised to calculate the number of rigid body modes, and the present elements are shown to be free of spurious zero-energy modes. A number of example problems are solved by using the present library, and the results are compared with corresponding analytical solutions and with results from the standard displacement finite element method. The Integrated Force Method not only gives results that agree well with analytical and displacement method results but also outperforms the displacement method in stress calculations.

Kaljevic, Igor↗
NASA NTRS1998

An Atmospheric General Circulation Model with Chemistry for the CRAY T3E: Design, Performance Optimization and Coupling to an Ocean Model

The design, implementation and performance optimization on the CRAY T3E of an atmospheric general circulation model (AGCM) which includes the transport of, and chemical reactions among, an arbitrary number of constituents is reviewed. The parallel implementation is based on a two-dimensional (longitude and latitude) data domain decomposition. Initial optimization efforts centered on minimizing the impact of substantial static and weakly-dynamic load imbalances among processors through load redistribution schemes. Recent optimization efforts have centered on single-node optimization. Strategies employed include loop unrolling, both manually and through the compiler, the use of an optimized assembler-code library for special function calls, and restructuring of parts of the code to improve data locality. Data exchanges and synchronizations involved in coupling different data-distributed models can account for a significant fraction of the running time. Therefore, the required scattering and gathering of data must be optimized. In systems such as the T3E, there is much more aggregate bandwidth in the total system than in any particular processor. This suggests a distributed design. The design and implementation of a such distributed 'Data Broker' as a means to efficiently couple the components of our climate system model is described.

Farrara, John D.↗
NASA NTRS2001

On-Orbit Calibration of Redundant Spacecraft Gyros by Optimal Reduction to Three Axes

The Aqua spacecraft will carry four single-axis gyros configured with three orthogonal axes and one skew axis. This redundancy presents a challenge for batch methods of on-orbit gyro calibration that use a spacecraft rotation model deterministically related to gyro data, in that sensor data can respond to at most three angular velocity components. When the number of gyros, N, is greater than 3, the 3xN matrix, G, that reduces the N gyro measurements to three body-frame angular-velocity components cannot be fully determined by such methods; there are many such matrices that produce essentially the same angular velocity history. In such a case, spacecraft operators require information about the Nx3 gyro linear response matrix, R, that relates gyro outputs to the body-frame angular velocities causing them. This matrix provides sufficient information to determine multiple reduced-dimension G-matrices for use in case of failure or degradation of one or more gyros, as well as to determine an optimal 3xN G for the fully-functional configuration. A viable proposal is to apply a 3xN pre-filter matrix, F, to the N gyro outputs before carrying out a conventional gyro calibration procedure. The angular-velocity history emerging from conventional calibration may then be used as input data, together with the same gyro data that generated it, to fit the alignment, scale-factor, and bias parameters of each gyro axis in turn. A difficulty of such a proposal is the arbitrariness in the choice of F. Due to gyro noise, different pre-filter matrices produce different calibrations. This paper presents a method of choosing F that is based on optimizing gyro consistency in the limit of infinite weight on gyro data, as compared to sensor data. The choice of F is independent of a priori alignment and is based on the gyro data alone. The method is applicable to any N of three or more, but reduces to conventional batch-estimation methodologies when N = 3. Results of computational comparison among calibration simulations using various choices of F will be presented for the Aqua gyro configuration with N = 4.

Radomski, M. S.↗
NASA NTRS2003

Theory of Fiber Optical Bragg Grating: Revisited

The reflected signature of an optical fiber Bragg grating is analyzed using the transfer function method. This approach is capable to cast all relevant quantities into proper places and provides a better physical understanding. The relationship between reflected signal, number of periods, index of refraction, and reflected wave phase is elucidated. The condition for which the maximum reflectivity is achieved is fully examined. We also have derived an expression to predict the reflectivity minima accurately when the reflected wave is detuned. Furthermore, using the segmented potential approach, this model can handle arbitrary index of refraction profiles and compare the strength of optical reflectivity of different profiles. The condition of a non-uniform grating is also addressed.

Tai, H.↗
DOE OSTI2020

Validation of Thermal Scattering Laws for Light Water at Elevated Temperatures with Diffusion Experiments

Light water is the most important neutron moderator in many reactor applications. Thermal neutron scattering kernels for H bound in H 2 O are represented by ENDF-format thermal scattering laws (TSLs) evaluated at discrete temperatures T . Nuclear data evaluations are conventionally validated utilizing a combination of critical benchmarks and experimental cross section data. Existing public critical benchmarks may be inadequate for testing water TSLs over the full range of T of interest in reactor applications, and experimental thermal scattering cross section data for water is sparse at elevated T . In this work, MC21 is used to simulate the decay of the equilibrium thermal neutron flux in light-water spheres of several radii. The fundamental-mode time eigenvalue $\mathcal{a}$ is calculated at 22 °C and 227 °C (at saturation pressure) for each sphere using three different H-H 2 O TSLs. Fitting $\mathcal{a}$ to a polynomial function of geometric buckling allows calculation of the thermal neutron diffusion length L . Extrapolation lengths are treated as a function of the transport mean free path and geometry. Experimental L data from 25 publications at 49 temperatures (from 10 °C to 295 °C), computed by several different time-dependent and space-dependent decay methods, is used to develop an empirical fit of L vs. T . The MC21-calculated L for the TSLs tested are within ±1% of the predicted values at 22 °C and 227 °C. This validation approach, which may be repeated at arbitrary T , constitutes an integral benchmark specific to and characterizing the detailed physics of the thermal scattering kernel applied.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1978

Advanced panel-type influence coefficient methods applied to unsteady three dimensional potential flows

A panel method for solving unsteady, subsonic wind-body-tail flow problems is formulated and partially verified. The method is applicable to general aircraft configurations consisting of arbitrary arrangements of wings, bodies, tails, and nacelles. The wake may be located arbitrarily and the unsteady, transverse component of vorticity in the wake may be assigned any covection velocity. The wake in the unsteady flow problem, therefore, can be given the location and convection velocity of the wake produced by a steady flow which is the mean flow of the unsteady flow problem. The panel method has been used as a basis for expanding the unsteady kernel function in a power series to obtain panel influence coefficients which can be integrated in closed form.

Dusto, A. R.↗
DOE OSTI2022

Conservative high-order data transfer method on generalized polygonal meshes

A conservative data transfer (remap) between two meshes is an important step of arbitrary Lagrangian-Eulerian (ALE) hydrodynamics simulations. High-order numerical methods for ALE simulations require both high-order (curvilinear) meshes and high-order remap algorithms. Here we develop a conservative and bounds-preserving method for accurate remapping of discrete fields on generalized polygonal meshes with curvilinear edges. The properties of the proposed method are studied theoretically and numerically for various (smooth and non-smooth) mesh deformations and discrete fields that represent smooth and discontinuous functions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1971

A general stagnation-point convective heating equation for arbitrary gas mixtures

The stagnation-point convective heat transfer to an axisymmetric blunt body for arbitrary gases in chemical equilibrium was investigated. The gases considered were base gases of nitrogen, oxygen, hydrogen, helium, neon, argon, carbon dioxide, ammonia, and methane and 22 gas mixtures composed of the base gases. Enthalpies ranged from 2.3 to 116.2 MJ/kg, pressures ranged from 0.001 to 100 atmospheres, and the wall temperatures were 300 and 1111 K. A general equation for the stagnation-point convective heat transfer in base gases and gas mixtures was derived and is a function of the mass fraction, the molecular weight, and a transport parameter of the base gases. The relation compares well with present boundary-layer computer results and with other analytical and experimental results. In addition, the analysis verified that the convective heat transfer in gas mixtures can be determined from a summation relation involving the heat transfer coefficients of the base gases. The basic technique developed for the prediction of stagnation-point convective heating to an axisymmetric blunt body could be applied to other heat transfer problems.

Sutton, K.↗
NASA NTRS1990

Creep and creep rupture of strongly reinforced metallic composites

A creep and creep damage theory is presented for metallic composites with strong fibers. Application is to reinforced structures in which the fiber orientation may vary throughout but a distinct fiber direction can be identified locally (local transverse isotropy). The creep deformation model follows earlier work and is based on a flow potential function that depends on invariants reflecting stress and the material symmetry. As the focus is on the interaction of creep and damage, primary creep is ignored. The creep rupture model is an extension of continuum damage mechanics and includes an isochronous damage function that depends on invariants specifying the local maximum transverse tension and the maximum longitudinal shear stress. It is posited that at high temperature and low stress, appropriate to engineering practice, these stress components damage the fiber/matrix interface through diffusion controlled void growth, eventually causing creep rupture. Experiments are outlined for characterizing a composite through creep rupture tests under transverse tension and longitudinal shear. Application is made to a thin-walled pressure vessel with reinforcing fibers at an arbitrary helical angle. The results illustrate the usefulness of the model as a means of achieving optimal designs of composite structures where creep and creep rupture are life limiting.

Robinson, D. N.↗
NASA NTRS1994

Dynamics of binary-disk interaction. 1: Resonances and disk gap sizes

We investigate the gravitational interaction of a generally eccentric binary star system with circumbinary and circumstellar gaseous disks. The disks are assumed to be coplanar with the binary, geometrically thin, and primarily governed by gas pressure and (turbulent) viscosity but not self-gravity. Both ordinary and eccentric Lindblad resonances are primarily responsible for truncating the disks in binaries with arbitrary eccentricity and nonextreme mass ratio. Starting from a smooth disk configuration, after the gravitational field of the binary truncates the disk on the dynamical timescale, a quasi-equilibrium is achieved, in which the resonant and viscous torques balance each other and any changes in the structure of the disk (e.g., due to global viscous evolution) occur slowly, preserving the average size of the gap. We analytically compute the approximate sizes of disks (or disk gaps) as a function of binary mass ratio and eccentricity in this quasi-equilibrium. Comparing the gap sizes with results of direct simulations using the smoothed particle hydrodynamics (SPH), we obtain a good agreement. As a by-product of the computations, we verify that standard SPH codes can adequately represent the dynamics of disks with moderate viscosity, Reynolds number R approximately 10(exp 3). For typical viscous disk parameters, and with a denoting the binary semimajor axis, the inner edge location of a circumbinary disk varies from 1.8a to 2.6a with binary eccentricity increasing from 0 to 0.25. For eccentricities 0 less than e less than 0.75, the minimum separation between a component star and the circumbinary disk inner edge is greater than a. Our calculations are relevant, among others, to protobinary stars and the recently discovered T Tau pre-main-sequence binaries. We briefly examine the case of a pre-main-sequence spectroscopic binary GW Ori and conclude that circumbinary disk truncation to the size required by one proposed spectroscopic model cannot be due to Linblad resonances, even if the disk is nonviscous.

Artymowicz, Pawel↗
NASA NTRS1993

Interactive method for planning constrained, fuel-optimal orbital proximity operations

An interactive graphical method for planning fuel-efficient rendezvous trajectories in the multi-spacecraft environment of the space station is presented. The method allows the operator to compose a multi-burn transfer trajectory between arbitrary initial chaser and target trajectories. The available task time of the mission is limited and the maneuver is subject to various operational constraints, such as departure, arrival, plume impingement and spatial constraints. The maneuvers are described in terms of the relataive motion experienced in a Space-Station centered coordinate system. The optimization method is based on the primer vector and its extension to non-optimal trajectories. The visual feedback of trajectory shapes, operational constraints, and optimization functions, provided by user-transparaent and continuously active background computations, allows the operator to make fast, iterative design changes which rapidly converge to fuel-efficient solutions. The optimization functions are presented. A variety of simple design examples has been presented to demonstrate the usefulness of the method. In many cases the addition of a properly positioned intermediate waypoint resulted in fuel savings of up to 30%. Furthermore, due to the counter-intuitive character of the optimization functions, most fuel-optimal solutions could not have been found without the aid of the optimization tools. Operating the system was found to be very easy, and did not require any previous in-depth knowledge of orbital dynamics or trajectory. The planning tool is an example of operator assisted optimization of nonlinear cost-functions.

Abramovitz, Adrian↗
DOE OSTI2025

Designing open quantum systems with known steady states: Davies generators and beyond

We provide a systematic framework for constructing generic models of nonequilibrium quantum dynamics with a target stationary (mixed) state. Our framework identifies (almost) all combinations of Hamiltonian and dissipative dynamics that relax to a steady state of interest, generalizing the Davies’ generator for dissipative relaxation at finite temperature to nonequilibrium dynamics targeting arbitrary stationary states. We focus on Gibbs states of stabilizer Hamiltonians, identifying local Lindbladians compatible therewith by constraining the rates of dissipative and unitary processes. Moreover, given terms in the Lindbladian not compatible with the target state, our formalism identifies the operations – including syndrome measurements and local feedback – one must apply to correct these errors. Our methods also reveal new models of quantum dynamics: for example, we provide a “measurement-induced phase transition” in which measurable two-point functions exhibit critical (power-law) scaling with distance at a critical ratio of the transverse field and rate of measurement and feedback. Time-reversal symmetry – defined naturally within our formalism – can be broken both in effectively classical and intrinsically quantum ways. Our framework provides a systematic starting point for exploring the landscape of dynamical universality classes in open quantum systems, as well as identifying new protocols for quantum error correction.

Guo, Jinkang [Department of Physics and Center for↗
NASA NTRS1980

Numerical generation of two-dimensional grids by the use of Poisson equations with grid control at boundaries

A method for generating boundary-fitted, curvilinear, two dimensional grids by the use of the Poisson equations is presented. Grids of C-type and O-type were made about airfoils and other shapes, with circular, rectangular, cascade-type, and other outer boundary shapes. Both viscous and inviscid spacings were used. In all cases, two important types of grid control can be exercised at both inner and outer boundaries. First is arbitrary control of the distances between the boundaries and the adjacent lines of the same coordinate family, i.e., stand-off distances. Second is arbitrary control of the angles with which lines of the opposite coordinate family intersect the boundaries. Thus, both grid cell size (or aspect ratio) and grid cell skewness are controlled at boundaries. Reasonable cell size and shape are ensured even in cases wherein extreme boundary shapes would tend to cause skewness or poorly controlled grid spacing. An inherent feature of the Poisson equations is that lines in the interior of the grid smoothly connect the boundary points (the grid mapping functions are second order differentiable).

Sorenson, R. L.↗
NASA NTRS1987

Theoretical studies of Resonance Enhanced Stimulated Raman Scattering (RESRS) of frequency-doubled Alexandrite laser wavelength in cesium vapor

The solutions for the imaginary susceptibility of the Raman field transition with arbitrary relaxation rates and field strengths are examined for three different sets of relaxation rates. These rates correspond to: (1) Far Infrared (FIR) Raman lasers in the diabatic collision regime without consideration of coupled population decay in a closed system, (2) Raman FIR lasers in the diabatic collision regime with coupled population conserving decay, and (3) IR Raman gain in cesium vapor. The model is further expanded to include Doppler broadening and used to predict the peak gain as a function of detuning for a frequency doubled Alexandrite laser-pumped cesium vapor gain cell.

Lawandy, Nabil M.↗