Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computational mathematics”

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,243 records · Page 69

A mathematical problem and a Spacecraft Control Laboratory Experiment (SCOLE) used to evaluate control laws for flexible spacecraft. NASA/IEEE design challenge

The problen of controlling large, flexible space systems has been evaluated using computer simulation. In several cases, ground experiments have also been used to validate system performance under more realistic conditions. There remains a need, however, to test additional control laws for flexible spacecraft and to directly compare competing design techniques. A program is discussed which has been initiated to make direct comparisons of control laws for, first, a mathematical problem, then and experimental test article being assembled under the cognizance of the Spacecraft Control Branch at the NASA Langley Research Center with the advice and counsel of the IEEE Subcommittee on Large Space Structures. The physical apparatus will consist of a softly supported dynamic model of an antenna attached to the Shuttle by a flexible beam. The control objective will include the task of directing the line-of-sight of the Shuttle antenna configuration toward a fixed target, under conditions of noisy data, control authority and random disturbances.

Taylor, Lawrence W., Jr.↗

CaTs Lab (CHAOS and Thermal Sciences Laboratory)

The CHAOS and Thermal Sciences Laboratory (CaTs) at James Madison University evolved into a noteworthy effort to increase minority representation in the sciences and mathematics. Serving ten students and faculty directly, and nearly 50 students indirectly, CaTs, through recruitment efforts, workshops, mentoring programs, tutorial services and research and computational laboratories, fulfilled its intent to initiate an academically enriched research program aimed at strengthening the academic and self-actualization skills of undergraduate students with potential to pursue doctoral study in the sciences. The stated goal of the program was to increase by 5% the number of enrolled mathematics and science students into the program. Success far exceeded the program goals by producing 100% graduation rate of all supported recipients during its tenure, with 30% of the students subsequently in pursuit of graduate degrees. Student retention in the program exceeded 90% and faculty participation exceeded the three members involved in mentoring and tutoring, gaining multi-disciplinary support. Aggressive marketing of the program resulted in several paid summer internships and commitments from NASA and an ongoing relationship with CHROME, a nationally recognized organization which focuses on developing minority students in the sciences and mathematics. Success of the program was only limited by the limited fiscal resources at NASA which resulted in phasing out of the program.

Teate, Anthony A.↗

A finite element program for postbuckling calculations (PSTBKL)

The object of the research reported herein was to develop a general mathematical model and solution methodologies for analyzing the structural response of thin, metallic shell structures under large transient, cyclic, or static thermochemical loads. This report describes the computer program resulting from the research. Among the system responses associated with these loads and conditions are thermal buckling, creep buckling, and ratcheting. Thus geometric and material nonlinearities (of high order) have been anticipated and are considered in developing the mathematical model. The methodology is demonstrated through different problems of extension, shear, and of planar curved beams. Moreover, importance of the inclusion of large strains is clearly demonstrated, through the chosen applications.

Simitses, G. T.↗

Modelling Equilibrium and Fractional Crystallization in the System MgO-FeO-CaO-Al2O3-SiO2

A mathematical modelling technique for use in petrogenesis calculations in the system MgO-FeO-CaO-Al2O3-SiO2 is reported. Semiempirical phase boundary and elemental distribution information was combined with mass balance to compute approximate equilibrium crystallization paths for arbitrary system compositions. The calculation is applicable to a range of system compositions and fractionation calculations are possible. The goal of the calculation is the computation of the composition and quantity of each phase present as a function of the degree of solidification. The degree of solidification is parameterized by the heat released by the solidifying phases. The mathematical requirement for the solution of this problem is: (1) An equation constraining the composition of the magma for each solid phase in equilibrium with the liquidus phase, and (2) an equation for each solid phase and each component giving the distribution of that element between that phase and the magma.

Herbert, F.↗

Improved Timing Scheme for Spaceborne Precipitation Radar

An improved timing scheme has been conceived for operation of a scanning satellite-borne rain-measuring radar system. The scheme allows a real-time-generated solution, which is required for auto targeting. The current timing scheme used in radar satellites involves pre-computing a solution that allows the instrument to catch all transmitted pulses without transmitting and receiving at the same time. Satellite altitude requires many pulses in flight at any time, and the timing solution to prevent transmit and receive operations from colliding is usually found iteratively. The proposed satellite has a large number of scanning beams each with a different range to target and few pulses per beam. Furthermore, the satellite will be self-targeting, so the selection of which beams are used will change from sweep to sweep. The proposed timing solution guarantees no echo collisions, can be generated using simple FPGA-based hardware in real time, and can be mathematically shown to deliver the maximum number of pulses per second, given the timing constraints. The timing solution is computed every sweep, and consists of three phases: (1) a build-up phase, (2) a feedback phase, and (3) a build-down phase. Before the build-up phase can begin, the beams to be transmitted are sorted in numerical order. The numerical order of the beams is also the order from shortest range to longest range. Sorting the list guarantees no pulse collisions. The build-up phase begins by transmitting the first pulse from the first beam on the list. Transmission of this pulse starts a delay counter, which stores the beam number and the time delay to the beginning of the receive window for that beam. The timing generator waits just long enough to complete the transmit pulse plus one receive window, then sends out the second pulse. The second pulse starts a second delay counter, which stores its beam number and time delay. This process continues until an output from the first timer indicates there is less than one transmit pulse width until the start of the next receive event. This blocks future transmit pulses in the build-up phase. The feedback phase begins with the first timer paying off and starting the first receive window. When the first receive window is complete, the timing generator transmits the next beam from the list. When the second timer pays off, the second receive event is started. Following the second receive event, the timing generator will transmit the next beam on the list and start an additional timer. The timers work in a circular buffer fashion so there only need to be enough to cover the maximum number of echoes in flight.

Berkun, Andrew↗

Three-axis attitude determination from vector observations

Two computationally efficient algorithms are presented for determining three-axis attitude from two or more vector observations. The first of these, the TRIAD algorithm, provides a deterministic (i.e., nonoptimal) solution for the attitude based on two vector observations. The second, the QUEST algorithm, is an optimal algorithm which determines the attitude that achieves the best weighted overlap of an arbitrary number of reference and observation vectors. Analytical expressions are given for the covariance matrices for the two algorithms using a fairly realistic model for the measurement errors. The mathematical relationship of the two algorithms and their relative merits are discussed and numerical examples are given. The advantage of computing the covariance matrix in the body frame rather than in the inertial frame (e.g., in terms of Euler angles) is emphasized. These results are valuable when a single-frame attitude must be computed frequently. They will also be useful to the mission analyst or spacecraft engineer for the evaluation of launch-window constraints or of attitude accuracies for different attitude sensor configurations.

Shuster, M. D.↗

Mathematical modeling of the aerodynamics of high-angle-of-attack maneuvers

This paper is a review of the current state of aerodynamic mathematical modeling for aircraft motions at high angles of attack. The mathematical model serves to define a set of characteristic motions from whose known aerodynamic responses the aerodynamic response to an arbitrary high angle-of-attack flight maneuver can be predicted. Means are explored of obtaining stability parameter information in terms of the characteristic motions, whether by wind-tunnel experiments, computational methods, or by parameter-identification methods applied to flight-test data. A rationale is presented for selecting and verifying the aerodynamic mathematical model at the lowest necessary level of complexity. Experimental results describing the wing-rock phenomenon are shown to be accommodated within the most recent mathematical model by admitting the existence of aerodynamic hysteresis in the steady-state variation of the rolling moment with roll angle. Interpretation of the experimental results in terms of bifurcation theory reveals the general conditions under which aerodynamic hysteresis must exist.

Schiff, L. B.↗

Graph Theory Roots of Spatial Operators for Kinematics and Dynamics

Spatial operators have been used to analyze the dynamics of robotic multibody systems and to develop novel computational dynamics algorithms. Mass matrix factorization, inversion, diagonalization, and linearization are among several new insights obtained using such operators. While initially developed for serial rigid body manipulators, the spatial operators and the related mathematical analysis have been shown to extend very broadly including to tree and closed topology systems, to systems with flexible joints, links, etc. This work uses concepts from graph theory to explore the mathematical foundations of spatial operators. The goal is to study and characterize the properties of the spatial operators at an abstract level so that they can be applied to a broader range of dynamics problems. The rich mathematical properties of the kinematics and dynamics of robotic multibody systems has been an area of strong research interest for several decades. These properties are important to understand the inherent physical behavior of systems, for stability and control analysis, for the development of computational algorithms, and for model development of faithful models. Recurring patterns in spatial operators leads one to ask the more abstract question about the properties and characteristics of spatial operators that make them so broadly applicable. The idea is to step back from the specific application systems, and understand more deeply the generic requirements and properties of spatial operators, so that the insights and techniques are readily available across different kinematics and dynamics problems. In this work, techniques from graph theory were used to explore the abstract basis for the spatial operators. The close relationship between the mathematical properties of adjacency matrices for graphs and those of spatial operators and their kernels were established. The connections hold across very basic requirements on the system topology, the nature of the component bodies, the indexing schemes, etc. The relationship of the underlying structure is intimately connected with efficient, recursive computational algorithms. The results provide the foundational groundwork for a much broader look at the key problems in kinematics and dynamics. The properties of general graphs and trees of nodes and edge were examined, as well as the properties of adjacency matrices that are used to describe graph connectivity. The nilpotency property of such matrices for directed trees was reviewed, and the adjacency matrices were generalized to the notion of block weighted adjacency matrices that support block matrix elements. This leads us to the development of the notion of Spatial Kernel Operator SKO kernels. These kernels provide the basis for the development of SKO resolvent operators.

Jain, Abhinandan↗

Multiwavelength pyrometry for nongray bodies

A multiwavelength technique was developed and applied to measure the temperatures of nongray surfaces. The instruments required are a spectral radiometer, a dedicated auxiliary radiation source, and a computer. In general, three radiation spectra are recorded: (1) spectrum S sub 0 of the auxiliary radiation source; (2) spectrum S sub 1 of the surface-emitted radiation; and (3) spectrum S sub 2, the sum of the radiation of S sub 1 plus the reflected radiation due to the incidence of the auxiliary radiation source on the surface. Subtracting spectrum S sub 1 from spectrum Sub 2 yields the reflection spectrum resulting from the incidence radiation. From these spectra, a quantity z(lambda) is derived and is related to the reflectivity r(lambda) by r(lambda) = z(lambda)/f, where f is a constant. Spectrum S sub 1 is represented mathematically as the product of a wavelength-dependent emissivity obtained from Kirchhoff's law and a Planck function of temperature T. Application of two-variable (lambda and z), nonlinear, least-squares curve-fitting computer software to fit spectrum S sub 1 to this mathematical expression yielded the surface temperature. This technique also measured the spectral reflectivity and emissivity of the surface. Instrumentation necessary to extend measurement to elevated temperatures and in the presence of reflective interference is discussed.

Ng, Daniel↗

A mathematical basis for the design and design optimization of adaptive trusses in precision control

A mathematical basis for the optimal design of adaptive trusses to be used in supporting precision equipment is provided. The general theory of adaptive structures is introduced, and the global optimization problem of placing a limited number, q, of actuators, so as to maximally achieve precision control and provide prestress, is stated. Two serialized optimization problems, namely, optimal actuator placement for prestress and optimal actuator placement for precision control, are addressed. In the case of prestressing, the computation of a 'desired' prestress is discussed, the interaction between actuators and redundants in conveying the prestress is shown in its mathematical form, and a methodology for arriving at the optimal placement of actuators and additional redundants is discussed. With regard to precision control, an optimal placement scheme (for q actuators) for maximum 'authority' over the precision points is suggested. The results of the two serialized optimization problems are combined to give a suboptimal solution to the global optimization problem. A method for improving this suboptimal actuator placement scheme by iteration is presented.

Das, S. K.↗

Criteria for Modeling in LES of Multicomponent Fuel Flow

A report presents a study addressing the question of which large-eddy simulation (LES) equations are appropriate for modeling the flow of evaporating drops of a multicomponent liquid in a gas (e.g., a spray of kerosene or diesel fuel in air). The LES equations are obtained from the direct numerical simulation (DNS) equations in which the solution is computed at all flow length scales, by applying a spatial low-pass filter. Thus, in LES the small scales are removed and replaced by terms that cannot be computed from the LES solution and instead must be modeled to retain the effect of the small scales into the equations. The mathematical form of these models is a subject of contemporary research. For a single-component liquid, there is only one LES formulation, but this study revealed that for a multicomponent liquid, there are two non-equivalent LES formulations for the conservation equations describing the composition of the vapor. Criteria were proposed for selecting the multicomponent LES formulation that gives the best accuracy and increased computational efficiency. These criteria were applied in examination of filtered DNS databases to compute the terms in the LES equations. The DNS databases are from mixing layers of diesel and kerosene fuels. The comparisons resulted in the selection of one of the multicomponent LES formulations as the most promising with respect to all criteria.

Bellan, Josette↗

An investigation of constraint-based component-modeling for knowledge representation in computer-aided conceptual design

Originally, computer programs for engineering design focused on detailed geometric design. Later, computer programs for algorithmically performing the preliminary design of specific well-defined classes of objects became commonplace. However, due to the need for extreme flexibility, it appears unlikely that conventional programming techniques will prove fruitful in developing computer aids for engineering conceptual design. The use of symbolic processing techniques, such as object-oriented programming and constraint propagation, facilitate such flexibility. Object-oriented programming allows programs to be organized around the objects and behavior to be simulated, rather than around fixed sequences of function- and subroutine-calls. Constraint propagation allows declarative statements to be understood as designating multi-directional mathematical relationships among all the variables of an equation, rather than as unidirectional assignments to the variable on the left-hand side of the equation, as in conventional computer programs. The research has concentrated on applying these two techniques to the development of a general-purpose computer aid for engineering conceptual design. Object-oriented programming techniques are utilized to implement a user-extensible database of design components. The mathematical relationships which model both geometry and physics of these components are managed via constraint propagation. In addition, to this component-based hierarchy, special-purpose data structures are provided for describing component interactions and supporting state-dependent parameters. In order to investigate the utility of this approach, a number of sample design problems from the field of aerospace engineering were implemented using the prototype design tool, Rubber Airplane. The additional level of organizational structure obtained by representing design knowledge in terms of components is observed to provide greater convenience to the program user, and to result in a database of engineering information which is easier both to maintain and to extend.

Kolb, Mark A.↗

Evaluation of Mathematical Turbulence Models

Simplified models for internal flow described, and their predictions compared with experimental results. Report presents account of various models used in computation of turbulent flows. Applications of these models to internal flows evaluated by analysis of predictions of various turbulence models in some important flow configurations.

Nallasamy, M.↗

GEODYN system support program, volume 4

The GEODYN Orbit Determination and Geodetic Parameter Estimation System consists of a set of computer programs designed to determine and analyze definitive satellite orbits and their associated geodetic and measurement parameters. This manual describes the Support Programs used by the GEODYN System. The mathematics and programming descriptions are detailed. The operational procedures of each program are presented. GEODYN ancillary analysis programs may be grouped into three different categories: (1) orbit comparison - DELTA (2) data analysis using reference orbits - GEORGE, and (3) pass geometry computations - GROUNDTRACK. All of the above three programs use one or more tapes written by the GEODYN program in either a data reduction or orbit generator run.

Mullins, N. E.↗

The solar radiation pressure on the Mariner 9 Mars orbiter.

The refined mathematical model of the force created by the light pressure of the sun has been used to compute the solar radiation pressure force acting on the Mariner 9 (Mariner Mars 1971) spacecraft, taking into account the reflectivity characteristics of all its components. The results have been compared with values obtained from Mariner 9 observations during the cruise phase and are found to be in agreement within 0.1% of the values.

Georgevic, R. M.↗

The solar radiation pressure on the Mariner 9 Mars orbiter

The refined mathematical model of the force created by the light pressure of the sun has been used to compute the solar radiation pressure force acting on the Mariner 9 (Mariner Mars 1971) spacecraft, taking into account the reflectivity characteristics of all its components. The results have been compared with values obtained from Mariner 9 observations during the cruise phase and are found to be in agreement within 0.1% of the values.

Georgevic, R. M.↗

Antenna evaluation study for the shuttle multispectral radar, phase 2

The results of the second phase of the Antenna Evaluation Study for the Shuttle Imaging Radar are presented. The objectives of Phase II were (1) to complete the specifications for the subarray test panels, (2) to begin a study of the effects of electrical and mechanical tolerance variations on overall SIRA performance, (3) to initiate the development of a mathematical model which adequately described the array performance and (4) to begin the development of a comprehensive computer program which will eventually simulate the performance characteristics of the antenna in a spaceborne environment. Items (2), (3), and (4) were begun in Phase I (ahead of schedule), and because of this, it has been possible to accelerate the Phase II modeling/simulation objectives to the point where simulations of expected mechanical/electrical errors have already been produced.

Coffey, E. L., III↗