Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “rational approximations”

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 163 records · Page 9

Continuation of tailored composite structures of ordered staple thermoplastic material

The search for the cost effective composite structure has motivated the investigation of several approaches to develop composite structure from innovative material forms. Among the promising approaches is the conversion of a planar sheet to components of complex curvature through sheet forming or stretch forming. In both cases, the potential for material stretch in the fiber direction appears to offer a clear advantage in formability over continuous fiber systems. A framework was established which allows the simulation of the anisotropic mechanisms of deformation of long discontinuous fiber laminates wherein the matrix phase is a viscous fluid. Predictions for the effective viscosities of a hyper-anisotropic medium consisting of collimated, discontinuous fibers suspended in viscous matrix were extended to capture the characteristics of typical polymers including non-Newtonian behavior and temperature dependence. In addition, the influence of fiber misorientation was also modeled by compliance averaging to determine ensemble properties for a given orientation distribution. A design tool is presented for predicting the effect of material heterogeneity on the performance of curved composite beams such as those used in aircraft fuselage structures. Material heterogeneity can be induced during manufacturing processes such as sheet forming and stretch forming of thermoplastic composites. This heterogeneity can be introduced in the form of fiber realignment and spreading during the manufacturing process causing radial and tangential gradients in material properties. Two analysis procedures are used to solve the beam problems. The first method uses separate two-dimensional elasticity solutions for the stresses in the flange and web sections of the beam. The separate solutions are coupled by requiring that forces and displacements match section boundaries. The second method uses an approximate Rayleigh-Ritz technique to find the solutions for more complex beams. Analyses are performed for curved beams of various cross-sections loaded in pure bending and with a uniform distributed load. Preliminary results show that the geometry of the beam dictates the effect of heterogeneity on performance. The role of heterogeneity is larger in beams with a small average radius-to-depth ration, R/t, where R is the average radius of the beam and t is the difference between the inside and outside radii. Results of the anlysis are in the form of stresses and displacements and are compared to both mechanics of materials and numerical solutions obtained using finite element analysis.

Santare, Michael H.↗

Flutter, vibration, and buckling of truncated orthotropic conical shells with generalized elastic edge restraint

A theoretical investigation has been made of the flutter, vibration, and buckling of truncated conical shells with generalized elastic edge restraint. The shell analysis is of the classical Donnell type, in-plane inertias and structural damping are neglected, and the aerodynamic loading is represented by the inviscid two-dimensional quasi-steady approximation. An approximate solution is obtained by the generalized Galerkin method. The accuracy and limitations of the analysis are illustrated by comparing numerical results for buckling an vibration with results of other investigations for various boundary conditions, applied loads, and shell geometries and stiffness. Sufficient numerical results are presented to permit the determination of the flutter condition for simply supported isotropic conical shells for a wide range of cone angle, length-radius ration, and radius-thickness ratio. Results are also presented to indicate some effects of variations in edge restraint, applied loads, and ring or stringer stiffening.

Flutter↗

Scaling laws for homogeneous turbulent shear flows in a rotating frame

The scaling properties of plane homogeneous turbulent shear flows in a rotating frame are examined mathematically by a direct analysis of the Navier-Stokes equations. It is proved that two such shear flows are dynamically similar if and only if their initial dimensionless energy spectrum E star (k star, 0), initial dimensionless shear rate SK sub 0/epsilon sub 0, initial Reynolds number K squared sub 0/nu epsilon sub 0, and the ration of the rotation rate to the shear rate omega/S are identical. Consequently, if universal equilibrium states exist, at high Reynolds numbers, they will only depend on the single parameter omega/S. The commonly assumed dependence of such equilibrium states on omega/S through the Richardson number Ri=-2(omega/S)(1-2 omega/S) is proven to be inconsistent with the full Navier-Stokes equations and to constitute no more than a weak approximation. To be more specific, Richardson number similarity is shown to only rigorously apply to certain low-order truncations of the Navier-Stokes equations (i.e., to certain second-order closure models) wherein closure is achieved at the second-moment level by assuming that the higher-order moments are a small perturbation of their isotropic states. The physical dependence of rotating turbulent shear flows on omega/S is discussed in detail along with the implications for turbulence modeling.

Speziale, Charles G.↗

Scaling laws for homogeneous turbulent shear flows in a rotating frame

The scaling properties of plane homogeneous turbulent shear flows in a rotating frame are examined mathematically by a direct analysis of the Navier-Stokes equations. It is proved that two such shear flows are dynamically similar if and only if their initial dimensionless energy spectrum E star (k star, 0), initial dimensionless shear rate SK sub 0/epsilon sub 0, initial Reynolds number K squared sub 0/nu epsilon sub 0, and the ration of the rotation rate to the shear rate omega/S are identical. Consequently, if universal equilibrium states exist, at high Reynolds numbers, they will only depend on the single parameter omega/S. The commonly assumed dependence of such equilibrium states on omega/S through the Richardson number Ri = -2(omega/S)(1-2 omega/S) is proven to be inconsistent with the full Navier-Stokes equations and to constitute no more than a weak approximation. To be more specific, Richardson number similarity is shown to only rigorously apply to certain low-order truncations of the Navier-Stokes equations (i.e., to certain second-order closure models) wherein closure is achieved at the second-moment level by assuming that the higher-order moments are a small perturbation of their isotropic states. The physical dependence of rotating turbulent shear flows on omega/S is discussed in detail along with the implications for turbulence modeling.

Speziale, Charles G.↗

On implementation of the extended interior penalty function

The extended interior penalty function formulation is implemented. A rational method for determining the transition between the interior and extended parts is set forth. The formulation includes a straightforward method for avoiding design points with some negative components, which are physically meaningless in structural analysis. The technique, when extended to problems involving parametric constraints, can facilitate closed form integration of the penalty terms over the most important parts of the parameter interval. The method lends itself well to the use of approximation concepts, such as design variable linking, constraint deletion and Taylor series expansions of response quantities in terms of design variables. Examples demonstrating the algorithm, in the context of planar orthogonal frames subjected to ground motion, are included.

Cassis, J. H.↗

An investigation of the SNS Josephson junction as a three-terminal device

A particular phenomenon of the SNS Josephson junction was investigated; i.e., control by a current entering the normal region and leaving through one of the superconducting regions. The effect of the control current on the junction was found to be dependent upon the ration of the resistances of the two halves of the N layer. A low frequency, lumped, nonlinear model was proposed to describe the electrical characteristics of the device, and a method was developed to plot the dynamic junction resistance as a function of junction current. The effective thermal noise temperature of the sample was determined. Small signal linearized analysis of the device suggests its use as an impedance transformer, although geometric limitations must be overcome. Linear approximation indicates that it is reciprocal and no power gain is possible. It is felt that, with suitable metallurgical and geometrical improvements, the device has promise to become a superconducting transistor.

Meissner, H.↗

A laser flash photolysis-resonance fluorescence kinetics study of the reaction Cl/2P/ + CH4 yields CH3 + HCl

The technique of laser flash photolysis-resonance fluorescence is employed to study the kinetics of the reaction Cl(2P) + CH4 yields CH3 + HCl over the temperature range 221-375 K. At temperatures less than or equal to 241 K the apparent bimolecular rate constant is found to be dependent upon the identity of the chemically inert gases in the reaction mixture. For Cl2/CH4/He reaction mixtures (total pressure = 50 torr) different bimolecular rate constants are measured at low and high methane concentrations. For Cl2/CH4/CCl/He and Cl2/CH4/Ar reaction mixtures, the bimolecular rate constant is independent of methane concentration, being approximately equal to the rate constant measured at low methane concentrations for Cl2/CH4/He mixtures. These rate constants are in good agreement with previous results obtained using the discharge flow-resonance fluorescence and competitive chlorination techniques. At 298 K the measured bimolecular rate constant is independent of the identity of the chemically inert gases in the reaction mixture and in good agreement with all previous investigations. The low-temperature results obtained in this investigation and all previous investigations can be rationalized in terms of a model which assumes that the Cl(2P 1/2) state reacts with CH4 much faster than the Cl(2P 3/2) state. Extrapolation of this model to higher temperatures, however, is not straightforward.

Ravishankara, A. R.↗

A degree of controllability definition - Fundamental concepts and application to modal systems

Starting from basic physical considerations, this paper develops a concept of the degree of controllability of a control system, and then develops numerical methods to generate approximate values of the degree of controllability for any linear time-invariant system. In many problems, such as the control of future, very large, flexible spacecraft and certain chemical process control problems, the question of how to choose the number and locations of the control system actuators is an important one. The results obtained here offer the control system designer a tool which allows him to rank the effectiveness of alternative actuator distributions, and hence to choose the actuator locations on a rational basis. The degree of controllability is shown to take a particularly simple form when the dynamic equations of a satellite are in second-order modal form. The degree of controllability concept has still other fundamental uses - it allows one to study the system structural relations between the various inputs and outputs of a linear system, which has applications to decoupling and model reduction.

Viswanathan, C. N.↗

Data reduction using cubic rational B-splines

A geometric method is proposed for fitting rational cubic B-spline curves to data that represent smooth curves including intersection or silhouette lines. The algorithm is based on the convex hull and the variation diminishing properties of Bezier/B-spline curves. The algorithm has the following structure: it tries to fit one Bezier segment to the entire data set and if it is impossible it subdivides the data set and reconsiders the subset. After accepting the subset the algorithm tries to find the longest run of points within a tolerance and then approximates this set with a Bezier cubic segment. The algorithm uses this procedure repeatedly to the rest of the data points until all points are fitted. It is concluded that the algorithm delivers fitting curves which approximate the data with high accuracy even in cases with large tolerances.

Chou, Jin J.↗

The Development Of Enabling Technologies For Submillimeter-Wave Remote Sensing of Ice Clouds From Space

This paper presents the chronological development of technologies and techniques that have led to a satellite mission concept aimed at quantifying the temporal and spatial distributions of upper tropospheric ice clouds. The Submillimeter-wave and Infrared Ice Cloud Experiment (SIRICE) is an Earth System Science Pathfinder mission concept designed to improve our understanding of the upper tropospheric water cycle and its coupling to the Earth s radiation budget. Ice outflow from convective storm systems is known to play an important role in regional energy budgets; however, ice generation and subsequent precipitation and sublimation are poorly quantified. SIRICE will provide measurements of ice cloud distributions and microphysical properties which are needed for understanding the crucial link between the hydrologic and energy cycles. The SIRICE measurement platform is comprised of two integrated instruments, the Submillimeter/millimeter-wave radiometer (SM4) and the Infrared Cloud Ice Radiometer (IRCIR). The primary instrument is the SM4, a conical scanner that provides a 1600 km swath of the Earth's surface at 53 degree incidence. The SM4 has 6 linearly polarized receivers measuring 12 spectral bands centered at 183 GHz, 325 GHz, 448 GHz, 643 GHz and 874 GHz; two receivers at 643 GHz measure horizontal and vertical polarizations. Submillimeter-wavelengths are well suited to the remote sensing of ice clouds due to the relative size of the wavelengths to particle sizes. Upwelling emission from lower tropospheric water vapor is scattered by the ice clouds thus causing a brightness temperature depression at submillimeter wavelengths. The IRCIR is a push broom imager with approximately 1500 km swath and spectral channels at 11 and 12 micrometers. This combination of coincident infrared and submillimeter-wavelength measurements were chosen because of its ability to provide retrieval of ice water path and median particle size for a wide range of ice clouds from thin cirrus to thick anvil structures. Over the past decade there has been a parallel development of submillimeter-wave technologies, demonstration instruments, and remote sensing techniques that have led to the present SIRICE mission concept. Mapping of these developmental paths reveals the origins, rational and maturity of features of the SIRICE payload such as its channel selection, compact design, and multipoint calibration. This presentation traces the evolution of the SIRICE mission concept from the early 1990's to its present status.

Racette, Paul↗

Estimation of the vortex length scale and intensity from two-dimensional samples

A method is proposed for estimating flow features that influence flame wrinkling in reciprocating internal combustion engines, where traditional statistical measures of turbulence are suspect. Candidate methods were tested in a computed channel flow where traditional turbulence measures are valid and performance can be rationally evaluated. Two concepts are tested. First, spatial filtering is applied to the two-dimensional velocity distribution and found to reveal structures corresponding to the vorticity field. Decreasing the spatial-frequency cutoff of the filter locally changes the character and size of the flow structures that are revealed by the filter. Second, vortex length scale and intensity is estimated by computing the ensemble-average velocity distribution conditionally sampled on the vorticity peaks. The resulting conditionally sampled 'average vortex' has a peak velocity less than half the rms velocity and a size approximately equal to the two-point-correlation integral-length scale.

Reuss, D. L.↗

Effects of Convective Transport of Solute and Impurities on Defect-Causing Kinetics Instabilities

For in-situ studies of the formation and evolution of step patterns during the growth of protein crystals, we have designed and assembled an experimental setup based on Michelson interferometry with the surface of the growing protein crystal as one of the reflective surfaces. The crystallization part of the device allows optical monitoring of a face of a crystal growing at temperature stable within 0.05 C in a developed solution flow of controlled direction and speed. The reference arm of the interferometer contains a liquid-crystal element that allows controlled shifts of the phase of the interferograms. We employ an image processing algorithm which combines five images with a pi/2 phase difference between each pair of images. The images are transferred to a computer by a camera capable of capturing 6-8 frames per second. The device allows data collection data regarding growth over a relatively large area (approximately .3 sq. mm) in-situ and in real time during growth. The estimated dept resolution of the phase shifting interferometry is about 100 A. The lateral resolution, depending on the zoom ratio, varies between 0.3 and 0.6 micrometers. We have now collected quantitative results on the onset, initial stages and development of instabilities in moving step trains on vicinal crystal surfaces at varying supersaturation, position on the facet, crystal size and temperature with the proteins ferritin, apoferritin and thaumatin. Comparisons with theory, especially with the AFM results on the molecular level processes, see below, allow tests of the rational for the effects of convective flows and, as a particular case, the lack thereof, on step bunching.

Vekilov, Peter G.↗

Chemical Thermodynamics and the Mathematical Integration of Reaction Kinetics

Key advances in the development of numerical methods for non-reacting compressible flows have been enabled by translating physical requirements into concrete numerical guidelines, such as the satisfaction of entropy inequalities for shock-capturing techniques [Lax, Contributions to Nonlinear Functional Analysis (1971) 603-634]. In the present work, we present nonlinear numerical analysis tools that draw from Chemical Thermodynamics , the branch of Nonequilibrium Thermodynamics that deals with chemical reactions. Through Gibbs formalism, chemical thermodynamics provides a well-known theoretical expression for the chemical equilibrium constant of a reaction in terms of reduced chemical potentials. A less-known, yet extremely valuable result, due to [Krambeck, Arch. Ration. Mech. Anal. , 38 (1970) 317], states that when this expression is implemented, mass-action kinetic models are consistent with the dynamical prescriptions of the 2nd law of thermodynamics. For fixed-temperature ordinary differential equations modeling constant-volume reacting gas mixtures, this leads to a decreasing Helmholtz free energy. If the temperature is allowed to vary in accordance with conservation of energy (1st law), this leads to the statement of increasing entropy. These nonlinear prescriptions can, and should be, used to further develop temporal integration techniques for reaction kinetics. We demonstrate that Krambeck's result holds even when the equilibrium constants are approximated from data. We prove this result by constructing the implicit free energy and the implicit entropy inherent to a given approximation. This is first done for a 5-species, 17-reaction model problem for air. With this structure established, elements of discrete entropy-stability theory [Tadmor, Acta Numer. , 12 (2003) 451] are leveraged to examine the consistency of time-integration schemes with these prescriptions. Using chemical potentials, one can compute the respective contributions of the kinetics model and of the temporal scheme to free energy/entropy variations. We introduce a nonlinear-stable version of the Discontinuous-Galerkin (DG) scheme in time which shows robustness improvements over the standard linearly-stable version. Most notably, the maximum timestep that can be resolved with the nonlinearly-stable variant tends to grow with polynomial order, in contrast to the linearly-stable variant. We generalize our constructions to arbitrary systems of reversible chemical reactions, ultimately showing that the compressible reacting Euler system admits the opposite of the implicitly constructed thermodynamic entropy as a mathematical entropy . This lays important theoretical foundations towards robust scheme development [Harten, J. Comput. Phys. 49 (1983) 151-164].

STMD↗

Continuum models for static and dynamic analysis of repetitive lattices

A simple, rational approach is presented for developing continuum models for large repetitive lattice structures subjected to static and dynamic loadings. The procedure involves introducing kinematic assumptions to reduce the dimensionality of the lattice, and equating the strain and kinetic energies of the continuum model to those of the original lattice structure. The proposed procedure is applied to obtain effective elastic and dynamic characteristics of continuum plate models for double-layered grids. Numerical results are presented of stress and free vibration problems for two double-layered grids. These problems demonstrate the high accuracy of the continuum plate models developed. Also, an assessment is made of different approximations to these models.

Noor, A. K.↗

MER Field Geologic Traverse in Gusev Crater, Mars: Initial Results From the Perspective of Spirit

This report casts the initial results of the traverse and science investigations by the Mars Exploration Rover (MER) Spirit at Gusev crater [1] in terms of data sets commonly used in field geologic investigations: Local mapping of geologic features, analyses of selected samples, and their location within the local map, and the regional context of the field traverse in terms of the larger geologic and physiographic region. These elements of the field method are represented in the MER characterization of the Gusev traverse by perspective-based geologic/morphologic maps, the placement of the results from Mossbauer, APXS, Microscopic Imager, Mini-TES and Pancam multispectral studies in context within this geologic/ morphologic map, and the placement of the overall traverse in the context of narrow-angle MOC (Mars Orbiter Camera) and descent images. A major campaign over a significance fraction of the mission will be the first robotic traverse of the ejecta from a Martian impact crater along an approximate radial from the crater center. The Mars Exploration Rovers have been conceptually described as 'robotic field geologists', that is, a suite of instruments with mobility that enables far-field traverses to multiple sites located within a regional map/image base at which in situ analyses may be done. Initial results from MER, where the field geologic method has been used throughout the initial course of the investigation, confirm that this field geologic model is applicable for remote planetary surface exploration. The field geologic method makes use of near-field geologic characteristics ('outcrops') to develop an understanding of the larger geologic context through continuous loop of rational steps focused on real-time hypothesis identification and testing. This poster equates 'outcrops' with the locations of in situ investigations and 'regional context' with the geology over distance of several kilometers. Using this fundamental field geologic method, we have identified the basic local geologic materials on the floor of Gusev at this site, their compositions and likely lithologies, origins, processes that have modified these materials, and their potential significance in the interpretation of the regional geology both spatially and temporally.

Crumpler, L.↗

Aubrite basalt vitrophyres: High sulfur silicate melts and a snapshot of aubrite formation

Two aubrite basalt vitrophyre clasts have been found within AMNH thin sections from the Parsa EH3 chondrite and the Khor Temiki aubrite. Polished sections of the Parsa Aubrite Inclusion (PAI) and the Khor Temiki Inclusion (KTI) were studied by optical, electron probe microanalysis (EPMA), and scanning electron microscopy (SEM) techniques with broad-beam and low absorbed EPMA currents used to minimize glass volatile loss. Some data have previously been reported for PAI and KTI may possibly correlate to a previously reported inclusion in Khor Tiimiki. In polished sections, PAI and KTI are approximately equal 4 mm in diameter and contain a large volume of glass. The clasts have similar textural characteristics and are akin to lunar vitrophyre textures. The glasses have high alkali rhyodacitic compositions Al-though PAI is peraluminous, KTI is significantly peralkaline. Additionally, the glasses have elevated sulfur concentrations that are extremely high by geochemical standards. SEM examination for beam overlap of microscopic CaS, FeS, and (Mg, Mn, Fe) S inclusions showed no such contamination. Furthermore, homogeneity of glass S content and low FeO contents help rule out contamination. Materials research data show that under reducing conditions alumino-silicate melts can dissolve up to several weight percent sulfur in the absence of Fe. The high S and alkali contents, the lack of associated high shock features, and the rationalized phase equilibria suggest that PAI and KTI are igneous melting products of an E-chondrite-like source material. Although large-scale impact melting cannot totally be ruled out, the above observations eliminate the possibility of in-situ shock melting.

Fogel, R. A.↗

Symmetry as Bias: Rediscovering Special Relativity

This paper describes a rational reconstruction of Einstein's discovery of special relativity, validated through an implementation: the Erlanger program. Einstein's discovery of special relativity revolutionized both the content of physics and the research strategy used by theoretical physicists. This research strategy entails a mutual bootstrapping process between a hypothesis space for biases, defined through different postulated symmetries of the universe, and a hypothesis space for physical theories. The invariance principle mutually constrains these two spaces. The invariance principle enables detecting when an evolving physical theory becomes inconsistent with its bias, and also when the biases for theories describing different phenomena are inconsistent. Structural properties of the invariance principle facilitate generating a new bias when an inconsistency is detected. After a new bias is generated. this principle facilitates reformulating the old, inconsistent theory by treating the latter as a limiting approximation. The structural properties of the invariance principle can be suitably generalized to other types of biases to enable primal-dual learning.

Lowry, Michael R.↗

A Satellite Borne Cadmium Sulfide Total Corpuscular Energy Detector

The properties of single crystals of cadmium sulfide as radiation detectors are described. It has been found possible to select crystals such that: (a) The ratio of increase of conductivity under irradiation to the rate of absorption of energy in the crystal is substantially independent of particle energy (over the examined ranges of 500 ev to 80 kev for electrons and 5 kev to 180 kev for protons) and of the magnitude of energy flux (over the range from.005 to 10 ergs/cm(sup 2 -sec); and (b) The above ration is substantially the same for protons, electrons, alpha particles, x-rays, and gamma-rays. For a driving voltage of 100 volts, typical crystal yield currents of 10(sup -7) to 10(sup- 6) amperes for each erg/cm(sup 2-sec) of energy absorbed by the crystal. The threshold of such crystal detectors (resulting from dark currents of the order of 10(sup 10 amp) is typically 10(sup -3) ergs/cm(sup 2- sec). For the selected crystals a response-temperature coefficient of -0.25% per degree centigrade is found for the temperature range -50 deg C to + 50 deg C. A description is given of a complete CdS total corpuscular energy detector for the study of geomagnetically trapped radiation by means of a satellite. The detector described has a dynamic range great than 10(sup 4), a solid angle of 10(exp -3) steradian, and a detection threshold of approximately 1 erg/cm(sup 2-sec-sterad). A similar detector employing a small magnet for the selective exclusion of electrons is also described. Noteworthy practical features of these detectors for satellite and space probe experiments are: (a) Use of bare crystals, without covering foils, in order to detect charged particles having energies as low as hundreds of electron volts. (b) Simplicity of electronic auxiliaries. (c) Compactness, lightweight and nechanical ruggedness. (d) Low electrical power requirements; and (e) Conversion of conduction current to the rate of a twostate relaxation oscillator in order to facilitate telemetric transmission of data. A pair of such detectors was flown as part of the s-46 satellite payload on March 23, 1960, but due to vehicular failure an orbit was not achieved and the operation of the CdS detectors was observed for only, six minutes of flight.

Freeman, John W.↗