Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mathematical physics”

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 865 records · Page 48

Mathematical modeling of ice accretion on airfoils

The progress toward development of a computer model suitable for predicting icing behavior on airfoils over a wide range of environmental conditions and airfoils shapes is reported. The LEWICE program was formulated to solve a set of equations which describe the physical processes which occur during accretion of ice on an airfoil, including heat transfer in a time dependent mode, with the restriction that the flow must be describable by a two-dimensional flow code. Input data comprises the cloud liquid water content, mean droplet diameter, ambient air temperature, air velocity, and relative humidity. A potential flowfield around the airfoil is calculated, along with the droplet trajectories within the flowfield, followed by local values of water droplet collection efficiency at the impact points. Both glaze and rime ice conditions are reproduced, and comparisons with test results on icing of circular cylinders showed good agreement with the physical situation.

Macarthur, C. D.↗

Physical mechanisms of heat, momentum and turbulence fluxes

In a qualitative way, the physical mechanisms which generate fluxes of heat, momentum, and turbulence in the atmosphere are discussed. This material is presented to acquaint people with the Earth science aspects of turbulence as important processes in the atmosphere. To attempt to describe turbulent fluxes of heat, momentum, and moisture in precise mathematical detail becomes an intractable problem. It is burdened by an eighth order set of equations involving more variables than equations. It is a closure problem which requires complicated assumptions that are not necessarily always satisfied, variable boundary conditions, and sparse observational data. Therefore, we must approach the problem in a simplified manner to obtain any kind of solution involving the variables of shear, stress, and heat, moisture, and momentum fluxes. There are other problems, of course, in which the inclusion of the planetary boundary layer is extremely important. Air pollution studies, air-sea exchanges, mesoscale models, and so on, must account for the planetary layer in very specific terms. Some of the physical mechanisms that are involved in generating fluxes are described.

Theon, John S.↗

Collision-induced squeezing in a harmonic oscillator

The concept of squeezing has so far been applied mainly to light, as is evidenced by the number of research works on the subject of squeezed light. Since, in quantum mechanics, both light and the simple harmonic oscillator are described within the same mathematical framework, there is of course no difficulty in applying the concept to the simple harmonic oscillator as well. In fact, the theoretical development of squeezed states and squeezed light owes much to the physical insights that one obtains as the analogy between light and the harmonic oscillator is exploited. The example presented shows clearly that two states with different phases in general have different degrees of squeezing, even if they have the same state distribution. This means that, even if one considers collision processes that produce the same state distribution, the degree of squeezing obtained during and after the collisions can be quite different, depending on how the phases phi(sub n) of the probability amplitudes develop in time as the collisions proceed. It is therefore evident that, for a detailed study of collision-induced squeezing, further study on the time development of the phases in collisions and its relation to collision parameters such as potential energy surfaces and collision energy is needed.

Lee, Hai-Woong↗

Modeling of outgassing and matrix decomposition in carbon-phenolic composites

A new release rate equation to model the phase change of water to steam in composite materials was derived from the theory of molecular diffusion and equilibrium moisture concentration. The new model is dependent on internal pressure, the microstructure of the voids and channels in the composite materials, and the diffusion properties of the matrix material. Hence, it is more fundamental and accurate than the empirical Arrhenius rate equation currently in use. The model was mathematically formalized and integrated into the thermostructural analysis code CHAR. Parametric studies on variation of several parameters have been done. Comparisons to Arrhenius and straight-line models show that the new model produces physically realistic results under all conditions.

Mcmanus, Hugh L.↗

Annual research briefs, 1993

The 1993 annual progress reports of the Research Fellow and students of the Center for Turbulence Research are included. The first group of reports are directed towards the theory and application of active control in turbulent flows including the development of a systematic mathematical procedure based on the Navier Stokes equations for flow control. The second group of reports are concerned with the prediction of turbulent flows. The remaining articles are devoted to turbulent reacting flows, turbulence physics, experiments, and simulations.

Source record↗

Uncertainty Modeling for Structural Control Analysis and Synthesis

The development of an accurate model of uncertainties for the control of structures that undergo a change in operational environment, based solely on modeling and experimentation in the original environment is studied. The application used throughout this work is the development of an on-orbit uncertainty model based on ground modeling and experimentation. A ground based uncertainty model consisting of mean errors and bounds on critical structural parameters is developed. The uncertainty model is created using multiple data sets to observe all relevant uncertainties in the system. The Discrete Extended Kalman Filter is used as an identification/parameter estimation method for each data set, in addition to providing a covariance matrix which aids in the development of the uncertainty model. Once ground based modal uncertainties have been developed, they are localized to specific degrees of freedom in the form of mass and stiffness uncertainties. Two techniques are presented: a matrix method which develops the mass and stiffness uncertainties in a mathematical manner; and a sensitivity method which assumes a form for the mass and stiffness uncertainties in macroelements and scaling factors. This form allows the derivation of mass and stiffness uncertainties in a more physical manner. The mass and stiffness uncertainties of the ground based system are then mapped onto the on-orbit system, and projected to create an analogous on-orbit uncertainty model in the form of mean errors and bounds on critical parameters. The Middeck Active Control Experiment is introduced as experimental verification for the localization and projection methods developed. In addition, closed loop results from on-orbit operations of the experiment verify the use of the uncertainty model for control analysis and synthesis in space.

Campbell, Mark E.↗

Numerical Simulation of Protoplanetary Vortices

The fluid dynamics within a protoplanetary disk has been attracting the attention of many researchers for a few decades. Previous works include, to list only a few among many others, the well-known prescription of Shakura & Sunyaev, the convective and instability study of Stone & Balbus and Hawley et al., the Rossby wave approach of Lovelace et al., as well as a recent work by Klahr & Bodenheimer, which attempted to identify turbulent flow within the disk. The disk is commonly understood to be a thin gas disk rotating around a central star with differential rotation (the Keplerian velocity), and the central quest remains as how the flow behavior deviates (albeit by a small amount) from a strong balance established between gravitational and centrifugal forces, transfers mass and momentum inward, and eventually forms planetesimals and planets. In earlier works we have briefly described the possible physical processes involved in the disk; we have proposed the existence of long-lasting, coherent vortices as an efficient agent for mass and momentum transport. In particular, Barranco et al. provided a general mathematical framework that is suitable for the asymptotic regime of the disk; Barranco & Marcus (2000) addressed a proposed vortex-dust interaction mechanism which might lead to planetesimal formation; and Lin et al. (2002), as inspired by general geophysical vortex dynamics, proposed basic mechanisms by which vortices can transport mass and angular momentum. The current work follows up on our previous effort. We shall focus on the detailed numerical implementation of our problem. We have developed a parallel, pseudo-spectral code to simulate the full three-dimensional vortex dynamics in a stably-stratified, differentially rotating frame, which represents the environment of the disk. Our simulation is validated with full diagnostics and comparisons, and we present our results on a family of three-dimensional, coherent equilibrium vortices.

Lin, H.↗

Bidirectional Classical Stochastic Processes with Measurements and Feedback

A measurement on a quantum system is said to cause the "collapse" of the quantum state vector or density matrix. An analogous collapse occurs with measurements on a classical stochastic process. This paper addresses the question of describing the response of a classical stochastic process when there is feedback from the output of a measurement to the input, and is intended to give a model for quantum-mechanical processes that occur along a space-like reaction coordinate. The classical system can be thought of in physical terms as two counterflowing probability streams, which stochastically exchange probability currents in a way that the net probability current, and hence the overall probability, suitably interpreted, is conserved. The proposed formalism extends the . mathematics of those stochastic processes describable with linear, single-step, unidirectional transition probabilities, known as Markov chains and stochastic matrices. It is shown that a certain rearrangement and combination of the input and output of two stochastic matrices of the same order yields another matrix of the same type. Each measurement causes the partial collapse of the probability current distribution in the midst of such a process, giving rise to calculable, but non-Markov, values for the ensuing modification of the system's output probability distribution. The paper concludes with an analysis of a classical probabilistic version of the so-called grandfather paradox.

Hahne, G. E.↗

Evolutionary Computing Methods for Spectral Retrieval

A methodology for processing spectral images to retrieve information on underlying physical, chemical, and/or biological phenomena is based on evolutionary and related computational methods implemented in software. In a typical case, the solution (the information that one seeks to retrieve) consists of parameters of a mathematical model that represents one or more of the phenomena of interest. The methodology was developed for the initial purpose of retrieving the desired information from spectral image data acquired by remote-sensing instruments aimed at planets (including the Earth). Examples of information desired in such applications include trace gas concentrations, temperature profiles, surface types, day/night fractions, cloud/aerosol fractions, seasons, and viewing angles. The methodology is also potentially useful for retrieving information on chemical and/or biological hazards in terrestrial settings. In this methodology, one utilizes an iterative process that minimizes a fitness function indicative of the degree of dissimilarity between observed and synthetic spectral and angular data. The evolutionary computing methods that lie at the heart of this process yield a population of solutions (sets of the desired parameters) within an accuracy represented by a fitness-function value specified by the user. The evolutionary computing methods (ECM) used in this methodology are Genetic Algorithms and Simulated Annealing, both of which are well-established optimization techniques and have also been described in previous NASA Tech Briefs articles. These are embedded in a conceptual framework, represented in the architecture of the implementing software, that enables automatic retrieval of spectral and angular data and analysis of the retrieved solutions for uniqueness.

Terrile, Richard↗

Shape Validation and RF Performance of Inflatable Antennas

Inflatable aperture antennas are an emerging technology that is being investigated for potential use in science and exploration missions. In particular, for missions to Mars and beyond, large deployable aperture antennas can provide the antenna gain required for high data rate communications, where the necessary antenna diameter exceeds the available volume of typical launch vehicle platforms. As inflatable aperture antennas have not been proven fully qualified for space missions, the author's Master's Thesis assessed the Ruze equation in characterizing this antenna technology. Inflatable aperture antennas do not follow a parabolic shape, and so the Ruze equation is not applicable due to the macroscopic shape errors of this technology. Therefore, geometric evaluations of the surface profile cannot simply correlate antenna gain degradation with the root-mean square shape error with a parabolic surface. Consequently, the focus of this work was to derive an accurate mathematical model of an inflatable aperture antenna to in order to characterize its Radio Frequency (RF) performance. Calculus of Variations methodologies were used to derive the surface profile shape of the inflatable aperture antenna. Physical Optics techniques were used to generate the antenna pattern profile. Validation testing of the predicted inflatable antenna shape model was performed through use of Laser Radar metrology measurements on an inflatable test article. Assessments of the RF performance of the inflatable aperture antenna, compared with nominally shaped paraboloidal antennas, were obtained through simulations of both technologies using a common diameter, depth, and arc length. Assessments of the RF performance of the inflatable aperture antenna was also performed against itself for changes in distance of the antenna feed location in the axial direction. Whereas the Ruze equation is limited to assessing gain reduction, this effort will also assess beam spreading and first side lobe angle and magnitude. The ability to characterize the RF response of this antenna will provide for an improved understanding of this technology. The accurate representation of the shape of this type of antenna technology will help to identify the most appropriate ways in which this technology could be utilized in planning future communication architectures for NASA missions to Mars and beyond.

Welch, Bryan↗

Initial Ada components evaluation

The SAIC has the responsibility for independent test and validation of the SSE. They have been using a mathematical functions library package implemented in Ada to test the SSE IV and V process. The library package consists of elementary mathematical functions and is both machine and accuracy independent. The SSE Ada components evaluation includes code complexity metrics based on Halstead's software science metrics and McCabe's measure of cyclomatic complexity. Halstead's metrics are based on the number of operators and operands on a logical unit of code and are compiled from the number of distinct operators, distinct operands, and total number of occurrences of operators and operands. These metrics give an indication of the physical size of a program in terms of operators and operands and are used diagnostically to point to potential problems. McCabe's Cyclomatic Complexity Metrics (CCM) are compiled from flow charts transformed to equivalent directed graphs. The CCM is a measure of the total number of linearly independent paths through the code's control structure. These metrics were computed for the Ada mathematical functions library using Software Automated Verification and Validation (SAVVAS), the SSE IV and V tool. A table with selected results was shown, indicating that most of these routines are of good quality. Thresholds for the Halstead measures indicate poor quality if the length metric exceeds 260 or difficulty is greater than 190. The McCabe CCM indicated a high quality of software products.

Moebes, Travis↗

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↗

Robots and Manipulator Systems. Part 2

This book represents a conference on robots and manipulator systems. Topics discussed include: mathematical modeling systems integrators, the use of TV monitors for remote operation of machinery on earth and in space, an evaluation of various control modes, remotely piloted aircraft, undersea operations, and prosthetic devices for the physically handicapped.

Heer, E.↗

Microwave remote sensing: Active and passive. Volume 2 - Radar remote sensing and surface scattering and emission theory

The fundamental principles of radar backscattering measurements are presented, including measurement statistics, Doppler and pulse discrimination techniques, and associated ambiguity functions. The operation of real and synthetic aperture sidelooking airborne radar systems is described, along with the internal and external calibration techniques employed in scattering measurements. Attention is given to the physical mechanisms responsible for the scattering emission behavior of homogeneous and inhomogeneous media, through a discussion of surface roughness, dielectric properties and inhomogeneity, and penetration depth. Simple semiempirical models are presented. Theoretical models involving greater mathematical sophistication are also given for extended ocean and bare soil surfaces, and the more general case of a vegetation canopy over a rough surface.

Ulaby, F. T.↗

Measurement of Residual Flexibility for Substructures Having Prominent Flexible Interfaces

Verification of a dynamic model of a constrained structure requires a modal survey test of the physical structure and subsequent modification of the model to obtain the best agreement possible with test data. Constrained-boundary or fixed-base testing has historically been the most common approach for verifying constrained mathematical models, since the boundary conditions of the test article are designed to match the actual constraints in service. However, there are difficulties involved with fixed-base testing, in some cases making the approach impractical. It is not possible to conduct a truly fixed-base test due to coupling between the test article and the fixture. In addition, it is often difficult to accurately simulate the actual boundary constraints, and the cost of designing and constructing the fixture may be prohibitive. For use when fixed-base testing proves impractical or undesirable, alternate free-boundary test methods have been investigated, including the residual flexibility technique. The residual flexibility approach has been treated analytically in considerable detail and has had limited frequency response measurements for the method. This concern is well-justified for a number of reasons. First, residual flexibilities are very small numbers, typically on the order of 1.0E-6 in/lb for translational diagonal terms, and orders of magnitude smaller for off-diagonal values. This poses difficulty in obtaining accurate and noise-free measurements, especially for points removed from the excitation source. A second difficulty encountered in residual measurements lies in obtaining a clean residual function in the process of subtracting synthesized modal data from a measured response function. Inaccuracies occur since modes are not subtracted exactly, but only to the accuracy of the curve fits for each mode; these errors are compounded with increasing distance from the excitation point. In this paper, the residual flexibility method is applied to a simple structure in both test and analysis. Measured and predicted residual functions are compared, and regions of poor data in the measured curves are described. It is found that for accurate residual measurements, frequency response functions having prominent stiffness lines in the acceleration/force format are needed. The lack of such stiffness lines increases measurement errors. Interface drive point frequency respose functions for shuttle orbiter payloads exhibit dominant stiffness lines, making the residual test approach a good candidate for payload modal tests when constrained tests are inappropriate. Difficulties in extracting a residual flexibility value from noisy test data are discussed. It is shown that use of a weighted second order least-squares curve fit of the measured residual function allows identification of residual flexibility that compares very well with predictions for the simple structure. This approach also provides an estimate of second order residual mass effects.

Tinker, Michael L.↗

Gravitational Physics Research

Gravitational physics research at ISPAE is connected with NASA's Relativity Mission (Gravity Probe B (GP-B)) which will perform a test of Einstein's General Relativity Theory. GP-B will measure the geodetic and motional effect predicted by General Relativity Theory with extremely stable and sensitive gyroscopes in an earth orbiting satellite. Both effects cause a very small precession of the gyroscope spin axis. The goal of the GP-B experiment is the measurement of the gyroscope precession with very high precision. GP-B is being developed by a team at Stanford University and is scheduled for launch in the year 2001. The related UAH research is a collaboration with Stanford University and MSFC. This research is focussed primarily on the error analysis and data reduction methods of the experiment but includes other topics concerned with experiment systems and their performance affecting the science measurements. The hydrogen maser is the most accurate and stable clock available. It will be used in future gravitational physics missions to measure relativistic effects such as the second order Doppler effect. The HMC experiment, currently under development at the Smithsonian Astrophysical Observatory (SAO), will test the performance and capability of the hydrogen maser clock for gravitational physics measurements. UAH in collaboration with the SAO science team will study methods to evaluate the behavior and performance of the HMC. The GP-B data analysis developed by the Stanford group involves complicated mathematical operations. This situation led to the idea to investigate alternate and possibly simpler mathematical procedures to extract the GP-B measurements form the data stream. Comparison of different methods would increase the confidence in the selected scheme.

Wu, S. T.↗

Quantitative Technique for Comparing Simulant Materials through Figures of Merit

The 1989 workshop report entitled Workshop on Production and Uses of Simulated Lunar Materials and the Lunar Regolith Simulant Materials: Recommendations for Standardization, Production, and Usage, NASA Technical Publication both identified and reinforced a need for a set of standards and requirements for the production and usage of the Lunar simulant materials. As NASA prepares to return to the Moon, and set out to Mars, a set of early requirements have been developed for simulant materials and the initial methods to produce and measure those simulants have been defined. Addressed in the requirements document are: 1) a method for evaluating the quality of any simulant of a regolith, 2) the minimum characteristics for simulants of Lunar regolith, and 3) a method to produce simulants needed for NASA's Exploration mission. As an extension of the requirements document a method to evaluate new and current simulants has been rigorously defined through the mathematics of Figures of Merit (FoM). Requirements and techniques have been developed that allow the simulant provider to compare their product to a standard reference material through Figures of Merit. Standard reference material may be physical material such as the Apollo core samples or material properties predicted for any landing site. The simulant provider is not restricted to providing a single "high fidelity" simulant, which may be costly to produce. The provider can now develop "lower fidelity" simulants for engineering applications such as drilling and mobility applications.

Rickman, Doug↗

Stereo Imaging Velocimetry Technique Using Standard Off-the-Shelf CCD Cameras

Stereo imaging velocimetry is a fluid physics technique for measuring three-dimensional (3D) velocities at a plurality of points. This technique provides full-field 3D analysis of any optically clear fluid or gas experiment seeded with tracer particles. Unlike current 3D particle imaging velocimetry systems that rely primarily on laser-based systems, stereo imaging velocimetry uses standard off-the-shelf charge-coupled device (CCD) cameras to provide accurate and reproducible 3D velocity profiles for experiments that require 3D analysis. Using two cameras aligned orthogonally, we present a closed mathematical solution resulting in an accurate 3D approximation of the observation volume. The stereo imaging velocimetry technique is divided into four phases: 3D camera calibration, particle overlap decomposition, particle tracking, and stereo matching. Each phase is explained in detail. In addition to being utilized for space shuttle experiments, stereo imaging velocimetry has been applied to the fields of fluid physics, bioscience, and colloidal microscopy.

McDowell, Mark↗