SEARCH · Engineering Papers
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.
Thermal radiation of solids - calculation of optical constants, radiation characteristics of rough surfaces of inhomogeneous materials, and complex systems
The purpose of the paper is to present the steps that have been taken towards an algorithm theory of stability and robustness.
A method of using ellipsometry or polarization analysis of light in total internal reflection of a surface to determine the number density of gold nanoparticles on a smooth substrate has been developed. The method can be modified to enable determination of densities of sparse distributions of nanoparticles in general, and is expected to be especially useful for measuring gold-nanoparticle-labeled biomolecules on microarrays. The method is based on theoretical calculations of the ellipsometric responses of gold nanoparticles. Elements of the calculations include the following: For simplicity, the gold nanoparticles are assumed to be spherical and to have the same radius. The distribution of gold nanoparticles is assumed to be a sub-monolayer (that is, sparser than a monolayer). The optical response of the sub-monolayer is modeled by use of a thin-island-film theory, according to which the polarizabilities parallel and perpendicular to the substrate are functions of the wavelength of light, the dielectric functions (permittivities expressed as complex functions of frequency or wavelength) of the gold and the suspending medium (in this case, the suspending medium is air), the fraction of the substrate area covered by the nanoparticles, and the radius of the nanoparticles. For the purpose of the thin-island-film theory, the dielectric function of the gold nanoparticles is modeled as the known dielectric function of bulk gold plus a correction term that is necessitated by the fact that the mean free path length for electrons in gold decreases with decreasing radius, in such a manner as to cause the imaginary part of the dielectric function to increase with decreasing radius (see figure). The correction term is a function of the nanoparticle radius, the wavelength of light, the mean free path and the Fermi speed of electrons in bulk gold, the plasma frequency of gold, and the speed of light in a vacuum. These models are used to calculate ellipsometric responses for various concentrations of gold nanoparticles having an assumed radius. The modeled data indicates distinct spectral features for both the real and the imaginary part of the dielectric function. An ellipsometric measurement would determine this distinct feature and thus can be used to measure nanoparticle concentration. By "ellipsometric responses" is meant the intensities of light measured in various polarization states as functions of the angle of incidence and the polarization states of the incident light. These calculated ellipsometric responses are used as calibration curves: Data from subsequent ellipsometric measurements on real specimens are compared with the calibration curves. The concentration of the nanoparticles on a specimen is assumed to be that of the calibration curve that most closely matches the data pertaining to that specimen.
From the viewpoint of a project manager responsible for the V&V (verification and validation) of a software system, mathematical verification techniques provide a possibly useful orthogonal dimension to otherwise standard empirical testing. However, the value they add to an empirical testing regime both in terms of coverage and in fault detection has been difficult to quantify. Furthermore, potential cost savings from replacing testing with mathematical verification techniques cannot be realized until the tradeoffs and synergies can be formulated. Integration of formal verification with empirical testing is also difficult because the idealized view of mathematical verification providing a correctness proof with total coverage is unrealistic and does not reflect the limitations imposed by computational complexity of mathematical techniques. This paper first describes a framework based on software reliability and formalized fault models for a theory of software design fault detection - and hence the utility of various tools for debugging. It then describes a utility model for integrating mathematical and empirical techniques with respect to fault detection and coverage analysis. It then considers the optimal combination of black-box testing, white-box (structural) testing, and formal methods in V&V of a software system. Using case studies from NASA software systems, it then demonstrates how this utility model can be used in practice.
The theory of galactic heavy ion fragmentation is furthered by incorporating a T matrix approach into the description of the three step process of abrasion, ablation, and final state interations. The connection between this T matrix and the interaction potential is derived. For resonant states, the substitution of complex energies for real energies in the transition rate is formerly justified for up to third order processes. The previously developed abrasion-ablation fragmentation theory is rederived from first principles and is shown to result from time ordering, classical probability, and zero width resonance approximations. Improvements in the accuracy of the total fragmentation cross sections require an alternative to the latter two approximations. A Lorentz invariant differential abrasion-ablation cross section is derived which explicitly includes the previously derived abrasion total cross sections. It is demonstrated that spectral and angular distributions can be obtained from the general Lorentz invariant form.
The design of large-scale complex engineered systems (LaCES) such as an aircraft is inherently interdisciplinary. Multiple engineering disciplines, drawing from a team of hundreds to thousands of engineers and scientists, are woven together throughout the research, development, and systems engineering processes to realize one system. Though research and development (R&D) is typically focused in single disciplines, the interdependencies involved in LaCES require interdisciplinary R&D efforts. This study investigates the interdisciplinary interactions that take place during the R&D and early conceptual design phases in the design of LaCES. Our theoretical framework is informed by both engineering practices and social science research on complex organizations. This paper provides preliminary perspective on some of the organizational influences on interdisciplinary interactions based on organization theory (specifically sensemaking), data from a survey of LaCES experts, and the authors experience in the research and design. The analysis reveals couplings between the engineered system and the organization that creates it. Survey respondents noted the importance of interdisciplinary interactions and their significant benefit to the engineered system, such as innovation and problem mitigation. Substantial obstacles to interdisciplinarity are uncovered beyond engineering that include communication and organizational challenges. Addressing these challenges may ultimately foster greater efficiencies in the design and development of LaCES and improved system performance by assisting with the collective integration of interdependent knowledge bases early in the R&D effort. This research suggests that organizational and human dynamics heavily influence and even constrain the engineering effort for large-scale complex systems.
The theoretical foundation and formulation of a numerical method for predicting the viscous flowfield in and about isolated three dimensional nozzles of geometrically complex configuration are presented. High Reynolds number turbulent flows are of primary interest for any combination of subsonic, transonic, and supersonic flow conditions inside or outside the nozzle. An alternating-direction implicit (ADI) numerical technique is employed to integrate the unsteady Navier-Stokes equations until an asymptotic steady-state solution is reached. Boundary conditions are computed with an implicit technique compatible with the ADI technique employed at interior points of the flow region. The equations are formulated and solved in a boundary-conforming curvilinear coordinate system. The curvilinear coordinate system and computational grid is generated numerically as the solution to an elliptic boundary value problem. A method is developed that automatically adjusts the elliptic system so that the interior grid spacing is controlled directly by the a priori selection of the grid spacing on the boundaries of the flow region.
The normal dispersion analysis for linear adiabatic wave propagation in stratified atmospheres adopts a real frequency and solves for the complex vertical wavenumber. We show that an exponentially stratified atmosphere does not have any spatially bounded normal modes for real frequencies. The usual treatment involves a representation where the imaginary part of the vertical wavenumber yields a rho(sup -1/2) dependence of the velocity amplitude which diverges as the absolute value of z approaches infinity. This solution includes a cutoff frequency below which acoustic modes cannot propagate. The standard dispersion analysis is a local representation of the wave behavior in both space and time but which is assumed to represent the motion throughout - infinity is less than t is less than infinity and 0 is less than infinity. However, any solution which has a purely sinusoidal time dependence extends through this full domain and is divergent due to the rho(sup -1/2) dependence. We show that a proper description is in terms of a near field of a boundary piston which is driven arbitrarily as a function of space and time. The atmosphere which responds to this piston is a semi-infinite layer which has an initially constant sound speed but which has the usual gravitational stratification. In a restricted domain of space and time above this boundary, the wavelike behavior of the medium may be described by frequencies and vertical wavenumbers which are both complex. When both parameters are allowed to have imaginary components, a new range of solutions is found for which there is virtually no cutoff frequency. We show that vertical energy propagation can take place through the solar atmosphere as a result of oscillations below the nominal cutoff frequency. Previously, the largest amplitude oscillations which generally have low frequencies were dropped from the calculation of energy flux becuase their frequencies are below the cutoff frequency. This new family of near-field waves permits these modes to carry energy vertically outward and raises the possibility that the largest amplitude 5 minute oscillations play a substantial role in the transport of acoustic energy to the chromosphere.
In this paper we describe the numerical analysis underlying our efforts to develop an accurate and reliable code for simulating flame propagation using complex physical and chemical models. We discuss our spatial and temporal discretization schemes, which in our current implementations range in order from two to six. In space we use staggered meshes to define discrete divergence and gradient operators, allowing us to approximate complex diffusion operators while maintaining ellipticity. Our temporal discretization is based on the use of preconditioning to produce a highly efficient linearly implicit method with good stability properties. High order for time accurate simulations is obtained through the use of extrapolation or deferred correction procedures. We also discuss our techniques for computing stationary flames. The primary issue here is the automatic generation of initial approximations for the application of Newton's method. We use a novel time-stepping procedure, which allows the dynamic updating of the flame speed and forces the flame front towards a specified location. Numerical experiments are presented, primarily for the stationary flame problem. These illustrate the reliability of our techniques, and the dependence of the results on various code parameters.
A systematic method is proposed for the development, optimization, and comparison of controller-models for the human operator. This is suitable for any designed model, even multiparameter systems. A random search technique is chosen for the parameter optimization. As valuation criteria for the quality of the model development the criterion function - the comparison between the input and output functions of the human operator and those of the model - and the most important characteristic values and functions of the statistical signal theory are used. A nonlinear multiparameter model for the human operator is being designed which considers the complex input information rate per time in a single display. The nonlinear features of the model are effected by a modified threshold element and a decision algorithm. Different display-configurations as well as various transfer functions of the controlled element are explained by different optimized parameter-combinations.
A study was designed to investigate the effects of behavioral and cognitive control on post-stress performance. Half of the subjects exposed to loud unpredictable noise bursts were given 'behavioral' control (a button which would terminate the noise), while the other half had no behavioral control. In addition, subjects were provided with one of three levels of feedback (success, failure, or no feedback) regarding their performance during the noise. It was expected that information about performance would provide subjects with an increased sense of 'cognitive' control which would affect their appraisal of stressful events and their later performance. The results indicated that subjects given feedback performed better on subsequent tasks than those given no feedback. Perceived behavioral control had little effect on performance. The causal attributions made by subjects were used to interpret these effects. These findings were viewed as supportive of Averill's (1973) notion that various types of control are related to stress in a complex fashion. The data may also support the reformulation by Abramson et al. (1978) of learned helplessness theory.
Due to geometric nonlinearities and complex dynamics, a decentralized technique for adaptive control for multilink robot arms is attractive. Lyapunov-function theory for stability analysis provides an approach to robust stabilization. Each joint of the arm is treated as a component subsystem. The adaptive controller is made locally stable with servo signals including proportional and integral gains. This results in the bound on the dynamical interactions with other subsystems. A nonlinear controller which stabilizes the system with uniform boundedness is used to improve the robustness properties of the overall system. As a result, the robot tracks the reference trajectories with convergence. This strategy makes computation simple and therefore facilitates real-time implementation.
In this work, we utilize fuzzy sets theory to evaluate and make predictions of flexural strength and density of NASA 6Y silicon nitride ceramic. Processing variables of milling time, sintering time, and sintering nitrogen pressure are used as an input to the fuzzy system. Flexural strength and density are the output parameters of the system. Data from 273 Si3N4 modulus of rupture bars tested at room temperature and 135 bars tested at 1370 C are used in this study. Generalized mean operator and Hamming distance are utilized to build the fuzzy predictive model. The maximum test error for density does not exceed 3.3 percent, and for flexural strength 7.1 percent, as compared with the errors of 1.72 percent and 11.34 percent obtained by using neural networks, respectively. These results demonstrate that fuzzy sets theory can be incorporated into the process of designing materials, such as ceramics, especially for assessing more complex relationships between the processing variables and parameters, like strength, which are governed by randomness of manufacturing processes.
The journal issue consists of abstracts presented at the International Society for the Study of the Origins of Life (ISSOL) conference. Topics include research on biological and chemical evolution including prebiotic evolution: cosmic and terrestrial; mechanisms of abiogenesis including synthesis and reactions of biomonomers; and analysis of cometary matter and its possible relationship to organic compounds on Earth. Theories and research on origins of ribonucleic acids (RNA), deoxyribonucleic acid (DNA), and other amino acids and complex proteins including their autocatalysis, replication, and translation are presented. Abiotic synthesis of biopolymers, mechanisms of the Genetic Code, precellular membrane systems and energetics are considered. Earth planetary evolution including early microfossils and geochemical conditions and simulations to study these conditions are discussed. The role of chirality in precellular evolution and the taxonomy and phylogeny of very simple organisms are reported. Past and future explorations in exobiology and space research directed toward study of the origins of life and solar system evolution are described.
A perturbative correction to the method of configuration interaction with single substitutions (CIS) is presented. This CIS(D) correction approximately introduces the effect of double substitutions which are absent in CIS excited states. CIS(D) is a second-order perturbation expansion of the coupled-cluster excited state method, restricted to single and double substitutions, in a series in which CIS is zeroth order, and the first-order correction vanishes. CIS (D) excitation energies are size consistent and the calculational complexity scales with the fifth power of molecular size, akin to second-order Moller-Plesset theory for the ground state. Calculations on singlet excited states of ethylene, formaldehyde, acetaldehyde, butadiene and benzene show that CIS (D) is a uniform improvement over CIS. CIS(D) appears to be a promising method for examining excited states of large molecules, where more accurate methods are not feasible.
Ionizing radiation produces several radiolytic species such as.OH, e-aq, and H. when interacting with biological matter. Following their creation, radiolytic species diffuse and chemically react with biological molecules such as DNA. Despite years of research, many questions on the DNA damage by ionizing radiation remains, notably on the indirect effect, i.e. the damage resulting from the reactions of the radiolytic species with DNA. To simulate DNA damage by ionizing radiation, we are developing a step-by-step radiation chemistry code that is based on the Green's functions of the diffusion equation (GFDE), which is able to follow the trajectories of all particles and their reactions with time. In the recent years, simulations based on the GFDE have been used extensively in biochemistry, notably to simulate biochemical networks in time and space and are often used as the "gold standard" to validate diffusion-reaction theories. The exact GFDE for partially diffusion-controlled reactions is difficult to use because of its complex form. Therefore, the radial Green's function, which is much simpler, is often used. Hence, much effort has been devoted to the sampling of the radial Green's functions, for which we have developed a sampling algorithm This algorithm only yields the inter-particle distance vector length after a time step; the sampling of the deviation angle of the inter-particle vector is not taken into consideration. In this work, we show that the radial distribution is predicted by the exact radial Green's function. We also use a technique developed by Clifford et al. to generate the inter-particle vector deviation angles, knowing the inter-particle vector length before and after a time step. The results are compared with those predicted by the exact GFDE and by the analytical angular functions for free diffusion. This first step in the creation of the radiation chemistry code should help the understanding of the contribution of the indirect effect in the formation of DNA damage and double-strand breaks.
Application of three-dimensional (3D) woven composites is growing as an alternative to the use of ply-based composite materials. However, the design, analysis, modeling, and optimization of these materials is more challenging due to their complex and inherently multiscale geometries. Herein, a multiscale modeling procedure, based on efficient, semi-analytical micromechanical theories rather than the traditional finite element approach, is presented and applied to a 3D woven carbon–epoxy composite. A crack-band progressive damage model was employed for the matrix constituent to capture the globally observed nonlinear response. Realistic microstructural dimensions and tow-fiber volume fractions were determined from detailed X-ray computed tomography (CT) and scanning electron microscopy data. Pre-existing binder-tow disbonds and weft-tow waviness, observed in X-ray CT scans of the composite, were also included in the model. The results were compared with experimental data for the in-plane tensile and shear behavior of the composite. The tensile predictions exhibited good correlations with the test data. While the model was able to capture the less brittle nature of the in-plane shear response, quantitative measures were underpredicted to some degree.
Tribology, the study of friction and wear of materials, has achieved a new interest because of the need for energy conservation. Fundamental understanding of this field is very complex and requires a knowledge of solid-state physics, material science, chemistry, and mechanical engineering. This paper is meant to be didactic in nature and outlines some of the considerations needed for a tribology research program. The approach is first to present a simple model, a field emission tip in contact with a flat surface, in order to elucidate important considerations, such as contact area, mechanical deformations, and interfacial bonding. Then examples from illustrative experiments are presented. Finally, the current status of physical theories concerning interfacial bonding are presented.