Search NASA⌕ Search

SEARCH · Search NASA

Results for “Closed form solution”

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 775 records · Page 43

Benchmark solutions for the galactic heavy-ion transport equations with energy and spatial coupling

Nontrivial benchmark solutions are developed for the galactic heavy ion transport equations in the straightahead approximation with energy and spatial coupling. Analytical representations of the ion fluxes are obtained for a variety of sources with the assumption that the nuclear interaction parameters are energy independent. The method utilizes an analytical LaPlace transform inversion to yield a closed form representation that is computationally efficient. The flux profiles are then used to predict ion dose profiles, which are important for shield design studies.

Ganapol, Barry D.↗

The effect of dissipative inhomogeneous medium on the statistics of the wave intensity

One of the main theoretical points in the theory of wave propagation in random medium is the derivation of closed form equations to describe the statistics of the propagating waves. In particular, in one dimensional problems, the closed form representation of the multiple scattering effects is important since it contributes in understanding such problems like wave localization, backscattering enhancement, and intensity fluctuations. In this the propagation of plane waves in a layer of one-dimensional dissipative random medium is considered. The medium is modeled by a complex permittivity whose real part is a constant representing the absorption. The one dimensional problem is mathematically equivalent to the analysis of a transmission line with randomly perturbed distributed parameters and a single mode lossy waveguide and the results can be used to study the propagation of radio waves through atmosphere and the remote sensing of geophysical media. It is assumed the scattering medium consists of an ensemble of one-dimensional point scatterers randomly positioned in a layer of thickness L with diffuse boundaries. A Poisson impulse process with density lambda is used to model the position of scatterers in the medium. By employing the Markov properties of this process an exact closed form equation of Kolmogorov-Feller type was obtained for the probability density of the reflection coefficient. This equation was solved by combining two limiting cases: (1) when the density of scatterers is small; and (2) when the medium is weakly dissipative. A two variable perturbation method for small lambda was used to obtain solutions valid for thick layers. These solutions are then asymptotically evaluated for small dissipation. To show the effect of dissipation, the mean and fluctuations of the reflected power are obtained. The results were compared with a lossy homogeneous medium and with a lossless inhomogeneous medium and the regions where the effect of absorption is not essential were discussed.

Saatchi, Sasan S.↗

Filtered Mass Density Function for Design Simulation of High Speed Airbreathing Propulsion Systems

The objective of this research is to develop and implement new methodology for large eddy simulation of (LES) of high-speed reacting turbulent flows. We have just completed two (2) years of Phase I of this research. This annual report provides a brief and up-to-date summary of our activities during the period: September 1, 2000 through August 31, 2001. In the work within the past year, a methodology termed "velocity-scalar filtered density function" (VSFDF) is developed and implemented for large eddy simulation (LES) of turbulent flows. In this methodology the effects of the unresolved subgrid scales (SGS) are taken into account by considering the joint probability density function (PDF) of all of the components of the velocity and scalar vectors. An exact transport equation is derived for the VSFDF in which the effects of the unresolved SGS convection, SGS velocity-scalar source, and SGS scalar-scalar source terms appear in closed form. The remaining unclosed terms in this equation are modeled. A system of stochastic differential equations (SDEs) which yields statistically equivalent results to the modeled VSFDF transport equation is constructed. These SDEs are solved numerically by a Lagrangian Monte Carlo procedure. The consistency of the proposed SDEs and the convergence of the Monte Carlo solution are assessed by comparison with results obtained by an Eulerian LES procedure in which the corresponding transport equations for the first two SGS moments are solved. The unclosed SGS convection, SGS velocity-scalar source, and SGS scalar-scalar source in the Eulerian LES are replaced by corresponding terms from VSFDF equation. The consistency of the results is then analyzed for a case of two dimensional mixing layer.

Drozda, T. G.↗

Time-domain analysis of planar microstrip devices using a generalized Yee-algorithm based on unstructured grids

The generalized Yee-algorithm is presented for the temporal full-wave analysis of planar microstrip devices. This algorithm has the significant advantage over the traditional Yee-algorithm in that it is based on unstructured and irregular grids. The robustness of the generalized Yee-algorithm is that structures that contain curved conductors or complex three-dimensional geometries can be more accurately, and much more conveniently modeled using standard automatic grid generation techniques. This generalized Yee-algorithm is based on the the time-marching solution of the discrete form of Maxwell's equations in their integral form. To this end, the electric and magnetic fields are discretized over a dual, irregular, and unstructured grid. The primary grid is assumed to be composed of general fitted polyhedra distributed throughout the volume. The secondary grid (or dual grid) is built up of the closed polyhedra whose edges connect the centroid's of adjacent primary cells, penetrating shared faces. Faraday's law and Ampere's law are used to update the fields normal to the primary and secondary grid faces, respectively. Subsequently, a correction scheme is introduced to project the normal fields onto the grid edges. It is shown that this scheme is stable, maintains second-order accuracy, and preserves the divergenceless nature of the flux densities. Finally, for computational efficiency the algorithm is structured as a series of sparse matrix-vector multiplications. Based on this scheme, the generalized Yee-algorithm has been implemented on vector and parallel high performance computers in a highly efficient manner.

Gedney, Stephen D.↗

Enabling attractive-repulsive potentials in binary-collision-approximation monte-carlo codes for ion-surface interactions

Abstract Binary Collision Approximation (BCA) codes for ion-material interactions, such as SRIM, Tridyn, F-TRIDYN, and SDtrimSP, have historically been limited to screened Coulomb potentials even at low energies due to the difficulty in numerically solving the Distance of Closest Approach (DOCA) problem for attractive-repulsive potentials. Techniques such as direct n-body simulation or modifications to Newton’s method are either prohibitively costly or not guaranteed to work for all potentials. Advanced rootfinding techniques, such as companion matrix solvers, offer a solution. For many attractive-repulsive potentials, however, a companion matrix cannot be used directly, because there is no way to put the associated functions into a monomial basis form. A complementary technique is proxy rootfinding—by finding the best-fit polynomial approximant of a function, the zeros of the approximant can be guaranteed to be close to the zeros of the function. Using the Chebyshev basis and grid offers additional guarantees with regards to the quality of the approximation, the speed of convergence, and the avoidance of Runge’s phenomenon. By finding Chebyshev interpolants and using the Chebyshev-Frobenius companion matrix, the zeros of any real function on a bounded domain can be found. Here we show that using an Adaptive Chebyshev Proxy Rootfinder with Automatic Subdivision (ACPRAS) with appropriate scaling functions, numerical issues presented by attractive-repulsive potentials, including those of scale, can be handled. Using these techniques, we show that it is possible to include any physically reasonable interatomic potential in a BCA code, and to guarantee correctness of the resulting scattering angle calculations.

Materials Science↗

Transient nucleation in condensed systems

Using classical nucleation theory we consider transient nucleation occurring in a one-component, condensed system under isothermal conditions. We obtain an exact closed-form expression for the time dependent cluster populations. In addition, a more versatile approach is developed: a numerical simulation technique which models directly the reactions by which clusters are produced. This simulation demonstrates the evolution of cluster populations and nucleation rate in the transient regime. Results from the simulation are verified by comparison with exact analytical solutions for the steady state. Experimental methods for measuring transient nucleation are assessed, and it is demonstrated that the observed behavior depends on the method used. The effect of preexisting cluster distributions is studied. Previous analytical and numerical treatments of transient nucleation are compared to the solutions obtained from the simulation. The simple expressions of Kashchiev are shown to give good descriptions of the nucleation behavior.

Kelton, K. F.↗

The effects of viscosity on the stability of a trailing-line vortex in compressible flow

We consider the effects of viscosity on the inviscid stability of the Batchelor vortex in a compressible flow. The problem is tackled asymptotically, in the limit of large (streamwise and azimuthal) wavenumbers, together with large Mach numbers. Previous studies, with viscous effects neglected, found that the nature of the solution passes through different regimes as the Mach number increases, relative to the wavenumber. This structure persists when viscous effects are included in the analysis. In the present study the mode present in the incompressible case ceases to be unstable at high Mach numbers and a center mode forms, whose stability characteristics are determined primarily by conditions close to the vortex axis. We find generally that viscosity has a stabilizing influence on the flow, while in the case of center modes, viscous effects become important at much larger Reynolds numbers than for the first class of disturbance.

Stott, Jillian A. K.↗

Structural optimization with approximate sensitivities

Computational efficiency in structural optimization can be enhanced if the intensive computations associated with the calculation of the sensitivities, that is, gradients of the behavior constraints, are reduced. Approximation to gradients of the behavior constraints that can be generated with small amount of numerical calculations is proposed. Structural optimization with these approximate sensitivities produced correct optimum solution. Approximate gradients performed well for different nonlinear programming methods, such as the sequence of unconstrained minimization technique, method of feasible directions, sequence of quadratic programming, and sequence of linear programming. Structural optimization with approximate gradients can reduce by one third the CPU time that would otherwise be required to solve the problem with explicit closed-form gradients. The proposed gradient approximation shows potential to reduce intensive computation that has been associated with traditional structural optimization.

Patnaik, S. N.↗

Helicopter Human Factors

Even under optimal conditions, helicopter flight is a most demanding form of human-machine interaction, imposing continuous manual, visual, communications, and mental demands on pilots. It is made even more challenging by small margins for error created by the close proximity of terrain in NOE flight and missions flown at night and in low visibility. Although technology advances have satisfied some current and proposed requirements, hardware solutions alone are not sufficient to ensure acceptable system performance and pilot workload. However, human factors data needed to improve the design and use of helicopters lag behind advances in sensor, display, and control technology. Thus, it is difficult for designers to consider human capabilities and limitations when making design decisions. This results in costly accidents, design mistakes, unrealistic mission requirements, excessive training costs, and challenge human adaptability. NASA, in collaboration with DOD, industry, and academia, has initiated a program of research to develop scientific data bases and design principles to improve the pilot/vehicle interface, optimize training time and cost, and maintain pilot workload and system performance at an acceptable level. Work performed at Ames, and by other research laboratories, will be reviewed to summarize the most critical helicopter human factors problems and the results of research that has been performed to: (1) Quantify/model pilots use of visual cues for vehicle control; (2) Improve pilots' performance with helmet displays of thermal imagery and night vision goggles for situation awareness and vehicle control; (3) Model the processes by which pilots encode maps and compare them to the visual scene to develop perceptually and cognitively compatible electronic map formats; (4) Evaluate the use of spatially localized auditory displays for geographical orientation, target localization, radio frequency separation; (5) Develop and flight test control/display concepts; (6) Quantify, model, predict, and improve pilots, workload-management strategies; and (7) Design computer-game trainers to reduce training time and cost.

Hart, Sandra G.↗

Modal element method for potential flow in non-uniform ducts: Combining closed form analysis with CFD

An analytical procedure is presented, called the modal element method, that combines numerical grid based algorithms with eigenfunction expansions developed by separation of variables. A modal element method is presented for solving potential flow in a channel with two-dimensional cylindrical like obstacles. The infinite computational region is divided into three subdomains; the bounded finite element domain, which is characterized by the cylindrical obstacle and the surrounding unbounded uniform channel entrance and exit domains. The velocity potential is represented approximately in the grid based domain by a finite element solution and is represented analytically by an eigenfunction expansion in the uniform semi-infinite entrance and exit domains. The calculated flow fields are in excellent agreement with exact analytical solutions. By eliminating the grid surrounding the obstacle, the modal element method reduces the numerical grid size, employs a more precise far field boundary condition, as well as giving theoretical insight to the interaction of the obstacle with the mean flow. Although the analysis focuses on a specific geometry, the formulation is general and can be applied to a variety of problems as seen by a comparison to companion theories in aeroacoustics and electromagnetics.

Baumeister, Kenneth J.↗

The determination of transport parameters of minority carriers in n-p junctions by means of an electron microscope - Critique of recent developments

It has recently been shown that amplitude modulated electron beams provide a novel means for the determination of minority carrier lifetimes, diffusion lengths, etc., in n-p junctions. In this paper it is shown that: (1) a recently published analysis based on a cylindrically symmetric configuration is incorrect, (2) the correct approach leads to a system of dual integral equations for which the formal solution is given, (3) in general, the short circuit current can only be determined by means of extensive computer calculations except in the case of large front surface recombination velocities, and (4) the difficulties encountered with cylindrically symmetric configurations (circular ohmic contacts and the like) are completely avoided with a choice of a planar geometry since simple closed form expressions for the short circuit current are readily available in this case.

Von Roos, O.↗

Inspection with Robotic Microscopic Imaging

Future Mars rover missions will require more advanced onboard autonomy for increased scientific productivity and reduced mission operations cost. One such form of autonomy can be achieved by targeting precise science measurements to be made in a single command uplink cycle. In this paper we present an overview of our solution to the subproblems of navigating a rover into place for microscopic imaging, mapping an instrument target point selected by an operator using far away science camera images to close up hazard camera images, verifying the safety of placing a contact instrument on a sample or finding nearby safe points, and analyzing the data that comes back from the rover. The system developed includes portions used in the Multiple Target Single Cycle Instrument Placement demonstration at NASA Ames in October 2004, and portions of the MI Toolkit delivered to the Athena Microscopic Imager Instrument Team for the MER mission still operating on Mars today. Some of the component technologies are also under consideration for MSL mission infusion.

Pedersen, Liam↗

Aqueous Alteration on Mars

Aqueous alteration is the change in composition of a rock, produced in response to interactions with H2O-bearing ices, liquids, and vapors by chemical weathering. A variety of mineralogical and geochemical indicators for aqueous alteration on Mars have been identified by a combination of surface and orbital robotic missions, telescopic observations, characterization of Martian meteorites, and laboratory and terrestrial analog studies. Mineralogical indicators for aqueous alteration include goethite (lander), jarosite (lander), kieserite (orbiter), gypsum (orbiter) and other Fe-, Mg-, and Ca-sulfates (landers), halides (meteorites, lander), phyllosilicates (orbiter, meteorites), hematite and nanophase iron oxides (telescopic, orbiter, lander), and Fe-, Mg-, and Ca-carbonates (meteorites). Geochemical indicators (landers only) for aqueous alteration include Mg-, Ca-, and Fe-sulfates, halides, and secondary aluminosilicates such as smectite. Based upon these indicators, several styles of aqueous alteration have been suggested on Mars. Acid-sulfate weathering (e.g., formation of jarosite, gypsum, hematite, and goethite), may occur during (1) the oxidative weathering of ultramafic igneous rocks containing sulfides, (2) sulfuric acid weathering of basaltic materials, and (3) acid fog (i.e., vapors rich in H2SO4) weathering of basaltic or basaltic-derived materials. Near-neutral or alkaline alteration occurs when solutions with pH near or above 7 move through basaltic materials and form phases such as phyllosilicates and carbonates. Very low water:rock ratios appear to have been prominent at most of the sites visited by landed missions because there is very little alteration (leaching) of the original basaltic composition (i.e., the alteration is isochemical or in a closed hydrologic system). Most of the aqueous alteration appears to have occurred early in the history of the planet (3 to 4.5 billion years ago); however, minor aqueous alteration may be occurring at the surface even today (e.g., in thin films of water or by acid fog).

Ming, Douglas W.↗

A Taylor-Galerkin finite element algorithm for transient nonlinear thermal-structural analysis

A Taylor-Galerkin finite element solution algorithm for transient nonlinear thermal-structural analysis of large, complex structural problems subjected to rapidly applied thermal-structural loads is described. The two-step Taylor-Galerkin algorithm is an application of an algorithm recently developed for problems in compressible fluid dynamics. The element integrals that appear in the algorithm can be evaluated in closed form for two and three dimensional elements.

Thornton, Earl A.↗

Csi-star: a Low-cost CSI Orbital Testbed

The topics are presented in viewgraph form and include the following: rationale for an on-orbit control-structurs interaction (CSI) test facility; CSI flight experiment objectives; feasibility study objectives; CSI free-flyer solution; feasibility study technical status summary; CSI-Star - a low-cost CSI free flyer; conceptual experiment design - option 2 configuration; Delta 2 - Quickstar Interface clambband capability; open and closed loop response of baselined truss with active struts; experiment weight baseline (option 1) configuration; experiment weight option 2 configuration; experiment power baseline (option 1) configuration; experiment power option 2 configuration; CSI Quickstar capabilities/requirements; and remaining work.

Edberg, D.↗

Non-Fermi Liquid Behavior in the Three-Dimensional Hubbard Model

We present a numerical study on the non-Fermi liquid behavior of a three-dimensional strongly correlated system. The Hubbard model in a simple cubic lattice is simulated by the dynamical cluster approximation; in particular, the quasi-particle weight is calculated at finite dopings for a range of temperatures. By fitting the quasi-particle weight to the marginal Fermi liquid form at finite doping near the putative quantum critical point, we find evidence of a separatrix between Fermi liquid and non-Fermi liquid regions. Our results suggest that a marginal Fermi liquid and possibly a quantum critical point exist in the non-symmetry broken solution of the three-dimensional interacting electron systems. We also calculate the spectral function, close to the half-filling, and we obtain evidence of pseudogap.

Kellar, Samuel↗

Synthetically Reversible, Proton-Mediated Nitrite N–O Bond Cleavage at a Dicopper Site

A monocationic dicopper(I,I) nitrite complex [Cu 2 (μ-κ 1 :κ 1 -O 2 N)DPFN][NTf 2 ] (2) (DPFN = 2,7-bis(fluoro-di(2-pyridyl)methyl)-1,8-naphthyridine, NTf 2 – = N(SO 2 CF 3 ) 2 – ), was synthesized by treatment of a dicopper acetonitrile complex, [Cu 2 (μ-MeCN)DPFN][NTf 2 ] 2 (1), with tetrabutylammonium nitrite ([nBu 4 N][NO 2 ]). DFT calculations indicate that 2 is one of three linkage isomers that are close in energy and presumably accessible in solution. Reaction of the μ-κ 1 :κ 1 -O 2 N complex with p-TolSH produces nitrous acid (HONO) and the corresponding dicopper thiolate species via an acid–base exchange reaction. Notably, treatment of 2 with HNTf 2 results in N–O bond cleavage in the putative, HONO-ligated complex to form the more thermodynamically favorable nitrosyl-bridged dicopper complex [Cu 2 (μ-NO)(μ-OH)DPFN][NTf 2 ] 2 (4). This scission can be reversed via deprotonation of the hydroxy ligand with KO t Bu. X-ray diffraction studies confirmed the solid-state molecular structures of 2 and 4. DFT calculations were used to construct a reaction coordinate diagram detailing formation of the μ-NO complex and to describe its electronic structure. The nitrosyl ligand in 4 is chemically labile, as demonstrated by its ready displacement in reactions with CO or NO 2 – .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Periodic motions (close to stationary) of an axisymmetric satellite with magnetic damping

Close to stationary periodic motions of an axisymmetric satellite in a circular orbit are considered. The satellite was equipped with a spherical magnetic damper. The investigation was conducted on the assumption that a strong magnet was installed on the damper float. Stationary rotations of the satellite around the axis of symmetry are selected as the generating solutions. The solutions are constructed in the form of power series of the smaller parameter, and they are extended numerically to the region of random values of the damping coefficient. The stability of the resulting solutions was investigated.

Ovchinnikov, M. Y.↗