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.

SourceResearch connectionFrom yearThrough year

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

At least 685 records · Page 38

All sources NASA NTRS DOE OSTI

NASA NTRS2023

Distorted-Toroidal Flux Rope Model

The 3D characterization of magnetic flux ropes observed in the heliosphere has been a challenging task for decades. This is mainly due to the limitations on inferring the 3D global topology and physical properties from the 1D time series from any spacecraft. To advance our understanding of magnetic flux ropes whose configuration departs from the typical stiff geometries, here we present an analytical solution for a 3D flux rope model with an arbitrary cross section and a toroidal global shape. This constitutes the next level of complexity following the elliptic-cylindrical (EC) geometry. The mathematical framework was established by Nieves-Chinchilla et al. with the EC flux rope model, which describes a magnetic topology with an elliptical cross section as a first approach to changes in the cross section. In the distorted-toroidal flux rope model, the cross section is described by a general function. The model is completely described by a nonorthogonal geometry and the Maxwell equations can be consistently solved to obtain the magnetic field and relevant physical quantities. As a proof of concept, this model is generalized in terms of the radial dependence of current density components. The last part of this paper is dedicated to a specific function, F (φ) = δ (1 - λ cosφ), to illustrate possibilities of the model. This model paves the way toward the investigation of complex distortions of magnetic structures in the solar wind. Future investigations will explore these distortions in depth by analyzing specific events; studying implications for physical quantities, such as magnetic fluxes, helicity, or energy; and evaluating the force balance with the ambient solar wind that allows such distortions.

Solar wind↗

NASA NTRS1976

Proton-tissue dose buildup factors

When an object is exposed to external radiation, the dose field within the object is a complicated function of the character of the external radiation, the shape of the object (including orientation), and the object's material composition. In this note, the dose conversion factors for protons in tissue are represented using buildup factors. A parametric form for the buildup factors is obtained. The values of the parameters are derived from Monte Carlo calculations of various authors. All the necessary information to estimate nuclear reaction effects in proton irration of convex objects of arbitrary shape is included.

Wilson, J. W.↗

NASA NTRS1987

Space station integrated wall design and penetration damage control

A methodology was developed to allow a designer to optimize the pressure wall, insulation, and meteoroid/debris shield system of a manned spacecraft for a given spacecraft configuration and threat environment. The threat environment consists of meteoroids and orbital debris, as specified for an arbitrary orbit and expected lifetime. An overall probability of no penetration is calculated, as well as contours of equal threat that take into account spacecraft geometry and orientation. Techniques, tools, and procedures for repairing an impacted and penetrated pressure wall were developed and tested. These techniques are applied from the spacecraft interior and account for the possibility of performing the repair in a vacuum. Hypervelocity impact testing was conducted to: (1) develop and refine appropriate penetration functions, and (2) determine the internal effects of a penetration on personnel and equipment.

Coronado, A. R.↗

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.↗

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.↗

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.↗

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↗

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.↗

NASA NTRS2011

Model Checking a Self-Stabilizing Distributed Clock Synchronization Protocol for Arbitrary Digraphs

This report presents the mechanical verification of a self-stabilizing distributed clock synchronization protocol for arbitrary digraphs in the absence of faults. This protocol does not rely on assumptions about the initial state of the system, other than the presence of at least one node, and no central clock or a centrally generated signal, pulse, or message is used. The system under study is an arbitrary, non-partitioned digraph ranging from fully connected to 1-connected networks of nodes while allowing for differences in the network elements. Nodes are anonymous, i.e., they do not have unique identities. There is no theoretical limit on the maximum number of participating nodes. The only constraint on the behavior of the node is that the interactions with other nodes are restricted to defined links and interfaces. This protocol deterministically converges within a time bound that is a linear function of the self-stabilization period.

Malekpour, Mahyar R.↗

NASA NTRS2012

Model Checking A Self-Stabilizing Synchronization Protocol for Arbitrary Digraphs

Malekpour, Mahyar R.↗