Search NASA⌕ Search

SEARCH · Search NASA

Results for “First principles calculations”

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

The scattering of electromagnetic pulses by a slit in a conducting screen

A direct method for calculating the impulse response of a slit in a conducting screen is presented which is derived specifically for the analysis of transient scattering by two-dimensional objects illuminated by a plane incident wave. The impulse response is obtained by assuming that the total response is composed of two sequences of diffracted waves. The solution is determined for the first two waves in one sequence by using Green's functions and the equivalence principle, for additional waves in the sequence by iteration, and for the other sequence by a transformation of coordinates. The cases of E-polarization and H-polarization are considered.

Ackerknecht, W. E., III↗

Applications of a General Finite-Difference Method for Calculating Bending Deformations of Solid Plates

This paper reports the findings of an investigation of a finite - difference method directly applicable to calculating static or simple harmonic flexures of solid plates and potentially useful in other problems of structural analysis. The method, which was proposed in doctoral thesis by John C. Houbolt, is based on linear theory and incorporates the principle of minimum potential energy. Full realization of its advantages requires use of high-speed computing equipment. After a review of Houbolt's method, results of some applications are presented and discussed. The applications consisted of calculations of the natural modes and frequencies of several uniform-thickness cantilever plates and, as a special case of interest, calculations of the modes and frequencies of the uniform free-free beam. Computed frequencies and nodal patterns for the first five or six modes of each plate are compared with existing experiments, and those for one plate are compared with another approximate theory. Beam computations are compared with exact theory. On the basis of the comparisons it is concluded that the method is accurate and general in predicting plate flexures, and additional applications are suggested. An appendix is devoted t o computing procedures which evolved in the progress of the applications and which facilitate use of the method in conjunction with high-speed computing equipment.

Walton, William C., Jr.↗

Thermochemical production of hydrogen via multistage water splitting processes

This paper presents and reviews the fundamental thermodynamic principles underlying thermochemical water splitting processes. The overall system is considered first and the temperature limitation in process thermal efficiency is developed. The relationship to an ideal water electrolysis cell is described and the nature of efficient multistage reaction processes is discussed. The importance of the reaction entropy change and the relation of the reaction free energy change to the work of separation is described. A procedure for analyzing thermochemical water splitting processes is presented and its use to calculate individual stage efficiency is demonstrated. A number of processes are used to illustrate the concepts and procedures.

Funk, J. E.↗

FD-TD modeling of 2-D dielectric waveguides for propagation and scattering of femtosecond optical solitons

Experimentalists have produced all-optical switches capable of 100-fs responses. To adequately model such switches, nonlinear effects in optical materials (both instantaneous and dispersive) must be included. In principle, the behavior of electromagnetic fields in nonlinear dielectrics can be determined by solving Maxwell's equations subject to the assumption that the electric polarization has a nonlinear relation to the electric field. However, until our previous work, the resulting nonlinear Maxwell's equations have not been solved directly. Rather, approximations have been made that result in a class of generalized nonlinear Schrodinger equations (GNLSE) that solve only for the envelope of the optical pulses. In this paper, we present first-time calculations from the vector nonlinear Maxwell's equations of femtosecond soliton propagation and scattering, including carrier waves, in two-dimensional systems of dielectric waveguides exhibiting the Kerr and Raman quantum effects. We use the finite-difference time-domain (FD-TD) method in an extension of our 1-D work. There, in a fundamental innovation, we treated the linear and nonlinear convolutions for the electric polarization as new dependent variables. By differentiating these convolutions in the time domain, we derived an equivalent system of coupled, nonlinear second-order ODE's. These equations together with Maxwell's equations form the system that is solved to determine the electromagnetic fields in inhomogeneous nonlinear dispersive media. Backstorage in time is limited to only that needed by the time-integration algorithm for the ODE's, rather than that needed to store the time-history of the kernel functions of the convolutions (1000-10,000 time steps). Thus, a 2-D nonlinear optics model from Maxwell's equations is now feasible.

Joseph, Rose↗

Synthesis of line profiles from models of structured winds

On the basis of a careful analysis of resonance line formation (both for singlets and doublets) in structured winds, present time dependent models of the line driven winds of hot stars are shown to be able to explain a number of observational features with respect to variability and structure: they are (in principle) able to reproduce the black and broad troughs (without any artificial 'turbulence velocity') and the 'blue edge variability' observed in saturate resonance lines: they might explain the 'long lived narrow absorption components' often observed in unsaturated lines at high velocities; they predict a relation between the 'edge velocity' of UV-lines and the radiation temperature of the observed X-ray emission. As a first example of the extent to which theoretical models can be constrained by comparisons between observations and profiles calculated by spectrum synthesis from structured winds, we show here that models with deep-seated onset of structure formation (approximately greater than 1.1 R(sub *)) produce resonance lines which agree qualitatively with observational findings; in contrast, the here presented models with structure formation only well out in the wind (approximately greater than 1.6 R(sub *) fail in this respect.

Puls, J.↗

Symbolic derivation of material property matrices in finite element analysis

The principles and operation of MMAX, a symbolic-computation program which automates the process of generating property matrices for structural materials, are briefly described and illustrated with sample analyses of a rubberlike material and an elastoplastic material. MMAX is written in LISP under the symbolic finite-element generator FINGER and the general symbolic manipulator MACSYMA; it first derives the formulas required by mathematical manipulation, and then translates the formulas into FORTRAN code, adapted to the particular type of machine to be used for the numerical calculations. This approach is shown to combine efficiently the advantages of symbolic and numerical computation for engineering applications.

Tan, H. Q.↗

Free Vibration of Uncertain Unsymmetrically Laminated Beams

Monte Carlo Simulation and Stochastic FEA are used to predict randomness in the free vibration response of thin unsymmetrically laminated beams. For the present study, it is assumed that randomness in the response is only caused by uncertainties in the ply orientations. The ply orientations may become random or uncertain during the manufacturing process. A new 16-dof beam element, based on the first-order shear deformation beam theory, is used to study the stochastic nature of the natural frequencies. Using variational principles, the element stiffness matrix and mass matrix are obtained through analytical integration. Using a random sequence a large data set is generated, containing possible random ply-orientations. This data is assumed to be symmetric. The stochastic-based finite element model for free vibrations predicts the relation between the randomness in fundamental natural frequencies and the randomness in ply-orientation. The sensitivity derivatives are calculated numerically through an exact formulation. The squared fundamental natural frequencies are expressed in terms of deterministic and probabilistic quantities, allowing to determine how sensitive they are to variations in ply angles. The predicted mean-valued fundamental natural frequency squared and the variance of the present model are in good agreement with Monte Carlo Simulation. Results, also, show that variations between plus or minus 5 degrees in ply-angles can affect free vibration response of unsymmetrically and symmetrically laminated beams.

Kapania, Rakesh K.↗

Computational aerodynamics and artificial intelligence

The general principles of artificial intelligence are reviewed and speculations are made concerning how knowledge based systems can accelerate the process of acquiring new knowledge in aerodynamics, how computational fluid dynamics may use expert systems, and how expert systems may speed the design and development process. In addition, the anatomy of an idealized expert system called AERODYNAMICIST is discussed. Resource requirements for using artificial intelligence in computational fluid dynamics and aerodynamics are examined. Three main conclusions are presented. First, there are two related aspects of computational aerodynamics: reasoning and calculating. Second, a substantial portion of reasoning can be achieved with artificial intelligence. It offers the opportunity of using computers as reasoning machines to set the stage for efficient calculating. Third, expert systems are likely to be new assets of institutions involved in aeronautics for various tasks of computational aerodynamics.

Mehta, U. B.↗

Jet shielding of jet noise

An experimental and theoretical study was conducted to develop a validated first principle analysis for predicting the jet noise reduction achieved by shielding one jet exhaust flow with a second, closely spaced, identical jet flow. A generalized fuel jet noise analytical model was formulated in which the acoustic radiation from a source jet propagates through the velocity and temperature discontinuity of the adjacent shielding jet. Input variables to the prediction procedure include jet Mach number, spacing, temperature, diameter, and source frequency. Refraction, diffraction, and reflection effects, which control the dual jet directivity pattern, are incorporated in the theory. The analysis calculates the difference in sound pressure level between the dual jet configuration and the radiation field based on superimposing two independent jet noise directivity patterns. Jet shielding was found experimentally to reduce noise levels in the common plane of the dual jet system relative to the noise generated by two independent jets.

Simonich, J. C.↗

Sensitivity of Cirrus Bidirectional Reflectance at MODIS Bands to Vertical Inhomogeneity of Ice Crystal Habits and Size Distribution

A common assumption in satellite imager-based cirrus retrieval algorithms is that the radiative properties of a cirrus cloud may be represented by those associated with a specific ice crystal shape (or habit) and a single particle size distribution. However, observations of cirrus clouds have shown that the shapes and sizes of ice crystals may vary substantially with height within the clouds. In this study we investigate the sensitivity of the top-of-atmosphere bidirectional reflectances at two MODIS bands centered at 0.65 micron and 2.11 micron to the cirrus models assumed to be either a single homogeneous layer or three distinct but contiguous, layers. First, we define the single- and three-layer cirrus cloud models with respect to ice crystal habit and size distribution on the basis of in situ replicator data acquired during the First ISCCP Regional Experiment (FIRE-II), held in Kansas during the fall of 1991. Subsequently, fundamental light scattering and radiative transfer theory is employed to determine the single scattering and the bulk radiative properties of the cirrus cloud. Regarding the radiative transfer computations, we present a discrete form of the adding/doubling principle by introducing a direct transmission function, which is computationally straightforward and efficient an improvement over previous methods. For the 0.65 micron band, at which absorption by ice is negligible, there is little difference between the bidirectional reflectances calculated for the one- and three-layer cirrus models, suggesting that the vertical inhomogeneity effect is relatively unimportant. At the 2.11 micron band, the bidirectional reflectances computed for both optically thin (tau = 1) and thick (tau = 10) cirrus clouds show significant differences between the results for the one- and three-layer models. The reflectances computed for the three-layer cirrus model are substantially larger than those computed for the single-layer cirrus. Finally, we find that cloud reflectance is very sensitive to the optical properties of the small crystals that predominate in the top layer of the three-layer cirrus model. It is critical to define the most realistic geometric shape for the small "quasi-spherical" ice crystals in the top layer for obtaining reliable single-scattering parameters and bulk radiative properties of cirrus.

Yang, P.↗

ATD-2 Surface Scheduling and Metering Concept

This presentation describes the concept of ATD-2 tactical surface scheduling and metering. The concept is composed of several elements, including data exchange and integration; surface modeling; surface scheduling; and surface metering. The presentation explains each of the elements. Surface metering is implemented to balance demand and capacity• When surface metering is on, target times from surface scheduler areconverted to advisories for throttling demand• Through the scheduling process, flights with CTOTs will not get addedmetering delay (avoids potential for ‘double delay’)• Carriers can designate certain flights as exempt from metering holds• Demand throttle in Phase 1 at CLT is through advisories sent to rampcontrollers for pushback instructions to the flight deck– Push now– Hold for an advised period of time (in minutes)• Principles of surface metering can be more generally applied to otherairports in the NAS to throttle demand via spot-release times (TMATs Strong focus on optimal use of airport resources• Flexibility enables stakeholders to vary the amount of delay theywould like transferred to gate• Addresses practical aspects of executing surface metering in aturbulent real world environment• Algorithms designed for both short term demand/capacityimbalances (banks) or sustained metering situations• Leverage automation to enable surface metering capability withoutrequiring additional positions• Represents first step in Tactical/Strategic fusion• Provides longer look-ahead calculations to enable analysis ofstrategic surface metering potential usage

Departure↗

Structural, dynamical and electronic properties of CaCuO2

The scalar relativistic version of an accurate first principles full potential self-consistent linearized muffin tin orbital (LMTO) method has been employed for describing the physical properties of the parent system of the high-T(sub c) oxide superconductors, i.e., CaCuO2. The presently employed modified version of the LMTO method is quite fast and goes beyond the usual LMTO ASA method in the sense that it permits a completely general shape of the potential and the charge density. Also, in contrast to LMTO ASA, the present method is also capable of treating distorted lattice structures accurately. The calculated values of the lattice parameters of pure CaCuO2 lie within 3% of the experimentally measured values for the Sr-doped system Ca(0.86)Sr(0.14)CuO(2). The computed electronic structures and the density of states is quite similar to those of the other oxide superconductors, except of their three- dimensional character because of the presence of strong coupling between the closely spaced CuO2 layers. The van Hove singularity peak appears slightly below the Fermi level and a small concentration of oxygenation /or/ substitutional doping may pin it at the Fermi level. The calculated frequencies for some symmetric frozen phonons for undoped CaCuO2 are quite near to the measured data for the Sr-doped CaCuO2.

Agrawal, Bal K.↗

The Impacts of Bowtie Effect and View Angle Discontinuity on MODIS Swath Data Gridding

We have analyzed two effects of the MODIS viewing geometry on the quality of gridded imagery. First, the fact that the MODIS scans a swath of the Earth 10 km wide at nadir, causes abrupt change of the view azimuth angle at the boundary of adjacent scans. This discontinuity appears as striping of the image clearly visible in certain cases with viewing geometry close to principle plane over the snow of the glint area of water. The striping is a true surface Bi-directional Reflectance Factor (BRF) effect and should be preserved during gridding. Second, due to bowtie effect, the observations in adjacent scans overlap each other. Commonly used method of calculating grid cell value by averaging all overlapping observations may result in smearing of the image. This paper describes a refined gridding algorithm that takes the above two effects into account. By calculating the grid cell value by averaging the overlapping observations from a single scan, the new algorithm preserves the measured BRF signal and enhances sharpness of the image.

Wang, Yujie↗

The rainout parameterization in a photochemical model

Photochemical models of tropospheric trace species typically simulate the effects of wet removal of soluble species from the atmosphere through a first-order rainout parameterization, which is independent of the species solubility and constant in time. In this work a modified first-order rainout parameterization is proposed that is a function of the species' solubility and takes into account the intermittent nature of atmospheric wet removal. In the case of a highly soluble species, such as HNO3, the inclusion of the effects of intermittent rainout in the parameterization leads to averaged rainout lifetimes considerably higher than ones used in most photochemical models. For illustrative purposes this parameterization is applied to a one-dimensional steady state model with parameterized eddy diffusion, although in principle it can be used in more sophisticated multidimensional models. The model results are compared to analagous calculations obtained with a 3-D tracer transport model.

Giorgi, F.↗

High Performance Parallel Methods for Space Weather Simulations

This is the final report of our NASA AISRP grant entitled 'High Performance Parallel Methods for Space Weather Simulations'. The main thrust of the proposal was to achieve significant progress towards new high-performance methods which would greatly accelerate global MHD simulations and eventually make it possible to develop first-principles based space weather simulations which run much faster than real time. We are pleased to report that with the help of this award we made major progress in this direction and developed the first parallel implicit global MHD code with adaptive mesh refinement. The main limitation of all earlier global space physics MHD codes was the explicit time stepping algorithm. Explicit time steps are limited by the Courant-Friedrichs-Lewy (CFL) condition, which essentially ensures that no information travels more than a cell size during a time step. This condition represents a non-linear penalty for highly resolved calculations, since finer grid resolution (and consequently smaller computational cells) not only results in more computational cells, but also in smaller time steps.

Hunter, Paul↗

First Use of Synoptic Vector Magnetograms for Global Nonlinear, Force-Free Coronal Magnetic Field Models

Context. The magnetic field permeating the solar atmosphere is generally thought to provide the energy for much of the activity seen in the solar corona, such as flares, coronal mass ejections (CMEs), etc. To overcome the unavailability of coronal magnetic field measurements, photospheric magnetic field vector data can be used to reconstruct the coronal field. Currently, there are several modelling techniques being used to calculate three-dimensional field lines into the solar atmosphere. Aims. For the first time, synoptic maps of a photospheric-vector magnetic field synthesized from the vector spectromagnetograph (VSM) on Synoptic Optical Long-term Investigations of the Sun (SOLIS) are used to model the coronal magnetic field and estimate free magnetic energy in the global scale. The free energy (i.e., the energy in excess of the potential field energy) is one of the main indicators used in space weather forecasts to predict the eruptivity of active regions. Methods. We solve the nonlinear force-free field equations using an optimization principle in spherical geometry. The resulting threedimensional magnetic fields are used to estimate the magnetic free energy content E(sub free) = E(sub nlfff) − E(sub pot), which is the difference of the magnetic energies between the nonpotential field and the potential field in the global solar corona. For comparison, we overlay the extrapolated magnetic field lines with the extreme ultraviolet (EUV) observations by the atmospheric imaging assembly (AIA) on board the Solar Dynamics Observatory (SDO). Results. For a single Carrington rotation 2121, we find that the global nonlinear force-free field (NLFFF) magnetic energy density is 10.3% higher than the potential one. Most of this free energy is located in active regions.

Sun: corona↗

The determination of exhaust cloud dimensions from films of space shuttle launches

Principles of photogrammetry are used to calculate the dimension of ground clouds produced from the Space Shuttle launch system. For each of three launches (Missions 41C, 41D, and 51A), a 16 mm camera recorded the ground cloud from three different locations. Measurements were made from outlines of the cloud and other features of interest which were traced onto paper at one-minute intervals using a 16 mm stop-action projector. Cloud characteristics such as top, max width at top, average width, and base are presented as a function of time. A temperature inversion was present each launch day and the cloud responded by first rising into the inversion and then descending to about the inversion height. Max tops were achieved in about three to five minutes and ranged from 2200 m for Mission 41C to 3500 m for Mission 41D. Cloud bases rose steadily to between 700 and 1000 m after 10 minutes. Average widths ranged from 500 to 1700 m depending on mission and camera. Photographs of digitizations of selected film frames are included in the report to show the irregular cloud shapes. Error sources for this analysis are also discussed.

Zak, R. A.↗

Localized basis functions and other computational improvements in variational nonorthogonal basis function methods for quantum mechanical scattering problems involving chemical reactions

The Generalized Newton Variational Principle for 3D quantum mechanical reactive scattering is briefly reviewed. Then three techniques are described which improve the efficiency of the computations. First, the fact that the Hamiltonian is Hermitian is used to reduce the number of integrals computed, and then the properties of localized basis functions are exploited in order to eliminate redundant work in the integral evaluation. A new type of localized basis function with desirable properties is suggested. It is shown how partitioned matrices can be used with localized basis functions to reduce the amount of work required to handle the complex boundary conditions. The new techniques do not introduce any approximations into the calculations, so they may be used to obtain converged solutions of the Schroedinger equation.

Schwenke, David W.↗