SEARCH · Engineering Papers

Results for “Variational approach”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 667 records · Page 37

All sources NASA NTRS DOE OSTI

NASA NTRS1993

Using experimental modal data to detect damage in a space truss

The effects of damage on a model space truss are investigated using experimentally measured modes along with a procedure that uses the vibratory residual forces and a weighted sensitivity analysis to estimate the extent of mass and/or stiffness variations, where damage is characterized as a reduction in structural stiffness. It is shown that this approach yields excellent results for a four-bay space truss under consideration and was able to accurately predict not only the location of the damage region but also estimate the relative severity for three different configurations. In addition, the method was able to accurately detect and differentiate a shift in structural mass.

Tavares, R.↗

NASA NTRS1995

Unified approach for incompressible flows

A unified approach for solving incompressible flows has been investigated in this study. The numerical CTVD (Centered Total Variation Diminishing) scheme used in this study was successfully developed by Sanders and Li for compressible flows, especially for the high speed. The CTVD scheme possesses better mathematical properties to damp out the spurious oscillations while providing high-order accuracy for high speed flows. It leads us to believe that the CTVD scheme can equally well apply to solve incompressible flows. Because of the mathematical difference between the governing equations for incompressible and compressible flows, the scheme can not directly apply to the incompressible flows. However, if one can modify the continuity equation for incompressible flows by introducing pseudo-compressibility, the governing equations for incompressible flows would have the same mathematical characters as compressible flows. The application of the algorithm to incompressible flows thus becomes feasible. In this study, the governing equations for incompressible flows comprise continuity equation and momentum equations. The continuity equation is modified by adding a time-derivative of the pressure term containing the artificial compressibility. The modified continuity equation together with the unsteady momentum equations forms a hyperbolic-parabolic type of time-dependent system of equations. Thus, the CTVD schemes can be implemented. In addition, the physical and numerical boundary conditions are properly implemented by the characteristic boundary conditions. Accordingly, a CFD code has been developed for this research and is currently under testing. Flow past a circular cylinder was chosen for numerical experiments to determine the accuracy and efficiency of the code. The code has shown some promising results.

Chang, Tyne-Hsien↗

NASA NTRS1994

On the Interaction of CO and NH3 with BH3 and BF3

BH3CO, BF3CO, BH3NH3, and BF3NH3 have been studied using density functional theory (DFT). The computed geometries are in reasonable agreement with experiment. BF3CO is found to be weakly bound. The binding energies of BH3CO, BH3NH3 , and BF3NH3 have also been computed using the coupled cluster singles and doubles level of theory, including a perturbational estimate of the connected triples [CCSD(T)]. These CCSD(T) calculations are performed using correlation consistent polarized-valence triple-zeta basis sets. These are expected to be the most accurate binding energies to date. The DFT binding energies are in good agreement with those obtained at the CCSD(T) level. The bonding is analyzed using the constrained space orbital variation (CSOV) technique. In addition to analyzing the bonding in BH3CO, BH3NH3, and BF3NH3, the CSOV approach is used to study why the bonding in BF3CO is so weak.

Bauschlicher, Charles W., Jr.↗

NASA NTRS2000

NASA/ARMY/BELL XV-15 Tiltrotor Low Noise Terminal Area Operations Flight Research Program

A series of three XV-15 acoustic flight tests have been conducted over a five year period by a NASA/Army/Bell Helicopter team to evaluate the noise reduction potential for tiltrotor aircraft during terminal area operations. Lower hemispherical noise characteristics for a wide range of steady-state terminal area type operating conditions were measured during the phase 1 test and indicated that the takeoff and level flight conditions were not significant contributors to the total noise of tiltrotor operations. Phase 1 results were used to design low noise approach profiles that were tested during the phase 2 and phase 3 tests, which used large area microphone arrays to directly measure the ground noise footprints. Approach profile designs emphasized noise reduction while maintaining handling qualities sufficient for tiltrotor commercial passenger ride comfort and flight safety under Instrument Flight Rules (IFR) conditions. This paper will discuss the weather, aircraft, tracking, guidance, and acoustic instrumentation systems, as well as the approach profile design philosophy, and the overall test program philosophy. Acoustic results are presented documenting the variation in tiltrotor noise due to changes in operating condition, indicating the potential for significant noise reduction using the unique tiltrotor capability of nacelle tilt.

Conner, David A.↗

NASA NTRS2002

Nonlinear Propagation of Planet-Generated Tidal Waves

The propagation and evolution of planet-generated density waves in protoplanetary disks is considered. The evolution of waves, leading to shock formation and wake dissipation, is followed in the weakly nonlinear regime. The 2001 local approach of Goodman and Rafikov is extended to include the effects of surface density and temperature variations in the disk as well as the disk cylindrical geometry and nonuniform shear. Wave damping due to shocks is demonstrated to be a nonlocal process spanning a significant fraction of the disk. Torques induced by the planet could be significant drivers of disk evolution on timescales of approx. 10(exp 6)-10(exp 7) yr, even in the absence of strong background viscosity. A global prescription for angular momentum deposition is developed that could be incorporated into the study of gap formation in a gaseous disk around the planet.

Rafikov, R. R.↗

NASA NTRS2004

Implicit and Explicit Spacecraft Gyro Calibration

This paper presents a comparison between two approaches to sensor calibration. According to one approach, called explicit, an estimator compares the sensor readings to reference readings, and uses the difference between the two to estimate the calibration parameters. According to the other approach, called implicit, the sensor error is integrated to form a different entity, which is then compared with a reference quantity of this entity, and the calibration parameters are inferred from the difference. In particular this paper presents the comparison between these approaches when applied to in-flight spacecraft gyro calibration. Reference spacecraft rate is needed for gyro calibration when using the explicit approach; however, such reference rates are not readily available for in-flight calibration. Therefore the calibration parameter-estimator is expanded to include the estimation of that reference rate, which is based on attitude measurements in the form of attitude-quaternion. A comparison between the two approaches is made using simulated data. It is concluded that the performances of the two approaches are basically comparable. Sensitivity tests indicate that the explicit filter results are essentially insensitive to variations in given spacecraft dynamics model parameters.

Bar-Itzhack, Itzhack Y.↗

NASA NTRS2009

Estimating Accurate Relative Spacecraft Angular Position from DSN VLBI Phases Using X-Band Telemetry or DOR Tones

At present spacecraft angular position with Deep Space Network (DSN) is determined using group delay estimates from very long baseline interferometer (VLBI) phase measurements employing differential one way ranging (DOR) tones. As an alternative to this approach, we propose estimating position of a spacecraft to half a fringe cycle accuracy using time variations between measured and calculated phases as the Earth rotates using DSN VLBI baseline(s). Combining fringe location of the target with the phase allows high accuracy for spacecraft angular position estimate. This can be achieved using telemetry signals of at least 4-8 MSamples/sec data rate or DOR tones.

x-band telemetry↗

NASA NTRS2018

Implementation of J-A Methodology Elastic-Plastic Crack Instability Analysis Capability into the WARP-3D Code

Characterization of the near crack-tip stress/strain fields is the foundation of fracture mechanics. The description of the near tip stress field and the prediction of when fracture occurs is well established for brittle materials that exhibit linear elastic behavior. However, in ductile materials or conditions that violate linear elastic assumptions (Aluminum alloys, Al 2024-T3, Al 2024- T351 etc.), the elastic-plastic crack-tip stress fields are characterized by the Hutchison-Rice-Rosengren (HRR) field. The J-integral is commonly used to characterize amplitude of the HRR field under elastic-plastic conditions. The J-integral has been demonstrated for crack-tip fields that are under high constraint conditions (i.e., small-scale plasticity where the J-dominance is maintained). However, as the external load increases, yielding changes from small- to largescale plasticity and usually a loss of constraint (i.e., reduction in the triaxial stress field along the crack front). The loss of constraint leads to the deviation of the crack-tip stress fields from that given by the HRR field. Hence, the J-dominance will be gradually lost and additional parameter(s) are required to quantify the crack-tip stress fields and predict fracture behavior. The assessment objectives were to: 1) implement a two-parameter (i.e., J-A) fracture criterion into an elastic-plastic three-dimensional (3D) finite element analysis (FEA), 2) validate the implementation by comparison with the A parameter from literature data, 3) conduct material characterization tests to quantify the material behavior and provide fracture data for validation of the J-A fracture criteria, and (4) perform evaluations to establish if the J-A criteria can be used to predict fracture in a ductile metallic material (e.g., aluminum alloys). The A parameter in these criteria is the second parameter in a three-term elastic-plastic asymptotic expansion of the neartip stress behavior. A series of extensive FEAs were performed using WARP3D software package to obtain solutions for the A parameter for different specimen configurations. The methodology needed for the estimation of the A parameter in the asymptotic expansion was developed and implemented using Matlab®. A user material (UMAT) routine was used to model the material stress-strain response using a Ramberg-Osgood power law with a hardening exponent (n) and a material coefficient (alpha). This UMAT routine was successfully implemented in WARP3D software and validated through comparison with the experimental data. Three configurations were extracted from published results: 1) center cracked plate (CCP), 2) single edge-cracked plate (SECP), and 3) double edge-cracked plate (DECP). These configurations and four other configurations (three-hole tension (THT)), three-point bend (3PTB), three-hole compact tension (3PCT), and compact tension (CT)) were analyzed to verify the methodology that was developed and implemented into WARP3D. Solutions of the A parameter were obtained for remote tension loading conditions that started with small-scale yielding and continued into the large-scale plasticity regime. The results indicate that the methodology developed can be used to calculate the elastic-plastic J-A parameters for test specimens with a range of crack geometries, material strain hardening behaviors, and loading conditions. The J-A parameters were implemented as fracture criteria and used to predict the test results. For comparison, other fracture criteria were used to predict the same test results. Major findings include: The A constraint parameter A varies with specimen type and applied load thus accurate determination is crucial in predicting the failure load, and the A parameter is asymptotic as the failure load is approached, making an accurate determination difficult (i.e., small differences in the A parameter can cause large variations in failure load) for materials exhibiting elastic-plastic behavior. The failure predictions from J-A methodology were more accurate than the traditionally used KC and J methods, and have comparable scatter to that observed when using the crack-tip opening angle (CTOA) method. However, the J-A methodology requires considerable effort (expertise level and labor) to implement and to evaluate the A parameter for different specimen types and materials, or to apply this methodology to part-through crack (e.g., 3D problems) structural applications.

Hamm, Kenneth R., Jr.↗

NASA NTRSDate not supplied

Simulating a Low-Weight Encoding of the Fermi-Hubbard Model Feedback and Support with QAOA

Recent work [1] has proposed a compact fermion-to-qubit mapping for the Fermi-Hubbard model (FHM), potentially lowering the overhead needed to simulate such a system on quantum hardware. In prior work [2], we numerically investigated how well a quantum annealer could prepare ground states of the FHM encoded using the resource-efficient mapping. In order to compare these results to gate-model systems, we turn to a QAOA-style approach and replace the annealing schedule with both Trotterizations thereof and also more general Hamiltonian-variational ansatze based on the terms of the mapped FH Hamiltonian. We investigate the achievable ground state fidelities of these gate sequences, consider various optimizations of the ansatze structure, and offer resource estimations of the optimized QAOA protocols for near-term hardware.

quantum computing↗

NASA NTRS1999

Ultrasound Thermal Field Imaging of Opaque Fluids

We have initiated an experimental program to develop an ultrasound system for non-intrusively imaging the thermal field in opaque fluids under an externally imposed temperature gradient. Many industrial processes involve opaque fluids, such as molten metals, semiconductors, and polymers, often in situations in which thermal gradients are important. For example, one may wish to understand semiconductor crystal growth dynamics in a Bridgman apparatus. Destructive testing of the crystal after the process is completed gives only indirect information about the fluid dynamics of the formation process. Knowledge of the coupled thermal and velocity fields during the growth process is then essential. Most techniques for non-intrusive velocity and temperature measurement in fluids are optical in nature, and hence the fluids studied must be transparent. In some cases (for example, LDV (laser Doppler velocimetry) and PIV (particle imaging velocimetry)) the velocities of small neutrally buoyant seed particles suspended in the fluid, are measured. Without particle seeding one can use the variation of the index of refraction of the fluid with temperature to visualize, through interferometric, Schlieren or shadowgraph techniques, the thermal field. The thermal field in turn gives a picture of the pattern existing in the fluid. If the object of study is opaque, non-optical techniques must be used. In this project we focus on the use of ultrasound, which propagates easily through opaque liquids and solids. To date ultrasound measurements have almost exclusively relied on the detection of sound scattered from density discontinuities inside the opaque material of interest. In most cases it has been used to visualize structural properties, but more recently the ultrasound Doppler velocimeter has become available. As in the optical case, it relies on seed particles that scatter Doppler shifted sound back to the detector. Doppler ultrasound techniques are, however, not useful for studying convective fluid flow in crystal growth, because particle seeding is unacceptable and flow velocities are typically too low to be resolved, and may be even lower in microgravity conditions where buoyancy forces are negligible. We will investigate a different use of ultrasound to probe the flows of opaque fluids non-intrusively and without the use of seed particles: our goal is to ultrasonically visualize the thermal field of opaque fluids with relatively high spatial resolution. The proposed technique relies upon the variation of sound speed with temperature of the fluid. A high frequency ultra-sound pulse passing through a fluid-filled chamber will traverse the chamber in a time determined by the relevant chamber dimension and the temperature of the fluid through which the pulse passes. With high time-resolution instrumentation that compares the excitation signal with the received pulse we can detect the influence of the fluid temperature on the pulse travel time. This is effectively an interferometric system, which in its optical form is an extremely sensitive approach to measuring thermal fields in fluids. Moreover, the temperature dependence of sound velocity in liquid metals is comparable to the temperature dependence of the speed of light required for accurate interferometric thermal images in transparent fluids. With an array of transducers scanned electronically a map of the thermal field over the chamber could be produced. An alternative approach would be to use the ultrasound analog of the shadowgraph technique. In the optical version, collimated light passes through the fluid, where it is focused or defocused locally by temperature field induced variations of the index of refraction. The resulting image reveals the thermal field through the spatial pattern of light intensity variations. By analogy, an ultrasound plane wave traversing an opaque fluid sample would be also locally focused or defocused depending on the speed of sound variations, giving rise to spatial variations in sound intensity that will reveal the thermal field pattern. These approaches could be applied to any situation in which temperature differences are expected to occur, and will rapidly provide information about the flow that simply cannot be obtained by any current intrusive or non-intrusive diagnostic technique. As materials processing in microgravity matures it will become increasingly important to have available simple and versatile diagnostic tools, such as we will develop, for studying the flows of opaque fluids under thermal forcing.

Andereck, C. David↗

NASA NTRS1994

Northern middle-latitude ozone profile features and trends observed by SBUV and Umkehr, 1979-1990

A comparison of Umkehr ozone profile data with the reprocessed solar backscatter ultraviolet (SBUV) ozone data in the northern middle-latitude region, 30 deg to 50 deg N, is reported. Although significant biases exist between the two types of observations, the long-term variations and least squares linear regression trends agree remarkably well over the comparison period of 1979 to 1990. The ozone trend in the upper stratosphere is of the order of -0.9%/yr. Near 25 km, little if any trend appears, but a larger negative trend is seen in the lower stratosphere near 15 km. Comparisons show that the average annual ozone cycles in the profiles also agree well. The upper stratospheric ozone results are consistent with photochemical model predictions of ozone depletion near 40 km that are due to the release of anthropogenically produced chlorofluorocarbons. The lower stratospheric ozone trend results are in reasonable agreement with published ozonesonde data trends. It is shown that the ozone trends in the lower stratospheric layers impact significantly on the total ozone trend of the order of -0.47%/yr. The good agreement now seen between the two types of observations suggest that the combined ground-based and satellite approach could provide a valuable database for long-term monitoring of stratospheric ozone for trends and extraordinary variations.

Deluisi, J. J.↗

NASA NTRS2004

Laser Sounder for Global Measurement of CO2 Concentrations in the Troposphere from Space: Progress

We describe progress toward developing a laser-based technique for the remote measurement of the tropospheric CO2 concentrations from orbit. Our goal is to demonstrate a lidar technique and instrument technology that will permit measurements of the CO2 column abundance in the lower troposphere from aircraft at the few ppm level, with a capability of scaling to permit global CO2 measurements from orbit. Accurate measurements of the tropospheric CO2 mixing ratio from space are challenging due to the many potential error sources. These include possible interference from other trace gas species, the effects of temperature, clouds, aerosols & turbulence in the path, changes in surface reflectivity, and variability in dry air density caused by changes in atmospheric pressure, water vapor and topographic height. Some potential instrumental errors include frequency drifts in the transmitter, small transmission and sensitivity drifts in the instrument. High signal-to-noise ratios and measurement stability are needed for mixing ratio estimates at the few ppm level. We have been developing a laser sounder approach as a candidate for a future space mission. It utilizes multiple different laser transmitters to permit simultaneous measurement of CO2 and O2 extinction, and aerosol backscatter in the same measurement path. It directs the narrow co-aligned laser beams from the instrument's fiber lasers toward nadir, and measures the energy of the strong laser echoes reflected from the Earth's land and water surfaces. During the measurement its narrow linewidth lasers are rapidly tuned on- and off- selected CO2 line near 1572 nm and an O2 absorption line near 770 nm. The receiver measures the energies of the laser echoes from the surface and any clouds and aerosols in the path with photon counting detectors. Ratioing the on- to off-line echo pulse energies for each gas permits the column extinction and column densities of CO2 and O2 to be estimated simultaneously via the differential absorption lidar technique. For the on-line wavelengths, the side of the selected absorption lines are used, which due to pressure broadening, weights the measurements to the lower troposphere, where CO2 variations caused by surface sources and sinks are largest. Simultaneous measurements of O2 column abundance are made using an identical approach using an O2 line. The laser backscatter profiles from clouds and aerosols are measured with other lidar channels, which permits identifying measurements influenced by clouds and/or aerosol scattering in the path. For space use, our lidar would continuously measure at nadir in near polar circular orbit. Using dawn and dusk measurements made over the same region will make it possible to sample the diurnal variations in CO2 mixing ratios. A 1-m diameter telescope is used for the receiver for all wavelengths. When averaging over 50 seconds, our calculations show a SNR of approximately 1500 is achievable for each gas at each on- and off-line measurement. Measurements from such a mission can be used to generate monthly global maps of the lower tropospheric CO2 column abundance. Our calculations show global coverage with an accuracy of a few ppm with a spatial resolution of approximately 50,000 sq. km are achievable each month. We have demonstrated some key elements of the laser, detector and receiver approaches in the laboratory and with measurements over a 206 m horizontal path. These include stable measurements of CO2 line shapes in an absorption cell using a fiber laser amplifier seeded by a tunable diode laser, measurement of small amplitude changes at low optical signal levels with the PMT receiver, and comparison of the horizontal path measurements of CO2 against those from an in-situ instrument.

Abshire, J. B.↗

NASA NTRS2010

Architecture-Based Unit Testing of the Flight Software Product Line

This paper presents an analysis of the unit testing approach developed and used by the Core Flight Software (CFS) product line team at the NASA GSFC. The goal of the analysis is to understand, review, and reconunend strategies for improving the existing unit testing infrastructure as well as to capture lessons learned and best practices that can be used by other product line teams for their unit testing. The CFS unit testing framework is designed and implemented as a set of variation points, and thus testing support is built into the product line architecture. The analysis found that the CFS unit testing approach has many practical and good solutions that are worth considering when deciding how to design the testing architecture for a product line, which are documented in this paper along with some suggested innprovennents.

Ganesan, Dharmalingam↗

NASA NTRS1975

Time-dependent hydrodynamic phenomena in the envelopes of the Be stars

Temporal and angular variations in the motion and distribution of circumstellar material about Be stars are studied by means of a simple hydrodynamical approach. The flow variables are each written as the sum of a known, steady-state axisymmetric term and a smaller component which contains both temporal and angular dependence. Substitution of these dichotomous expressions for the flow variables into Euler's equation and the continuity equation leads to a set of linearized equations for the terms which depend on time and angle. Solutions can be found to a somewhat approximated version of the linearized equations, and these solutions explain the qualitative features of the observed short-time variations.

Morgan, T. H.↗

NASA NTRS2004

Human-Induced Climate Variations Linked to Urbanization: From Observations to Modeling

The goal of this session is to bring together scientists from interdisciplinary backgrounds to discuss the data, scientific approaches and recent results focusing on the impact of urbanization on the climate. The discussion will highlight current observational and modeling capabilities being employed for investigating the urban environment and its linkage to the change in the Earth's climate system. The goal of the session is to identify our current stand and the future direction on the topic. Urbanization is one of the extreme cases of land use change. Most of population of the world has moved to urban areas. By 1995, more than 70% of population of North America and Europe were living in cities. By 2025, the United Nations estimates that 60% of the worlds population will live in cities. Although currently only 1.2% of the land is urban, better understanding of how the atmosphere-ocean-land-biosphere components interact as a coupled system and the influence of human activities on this system is critical. Our understanding of urbanization effect is incomplete, partly because human activities induce new changes on climate in addition to the original natural variations, and partly because previously few data available for study urban effect globally. Urban construction changes surface roughness, albedo, heat capacity and vegetation coverage. Traffic and industry increase atmospheric aerosol. It is suggested that urbanization may modify rainfall processes through aerosol-cloud interactions or dynamic feedbacks. Because urbanization effect on climate is determined by many factors including land cover, the city's microscale features, population density, and human lifestyle patterns, it is necessary to study urban areas over globe.

Shepherd, J. Marshall↗

NASA NTRS2007

Laser Sounder for Global Measurement of CO2 Concentrations in the Troposphere from Space

Measurements of tropospheric CO2 abundance with global-coverage, a few hundred km spatial and monthly temporal resolution are needed to quantify processes that regulate CO2 storage by the land and oceans. The Orbiting Carbon Observatory (OCO) is the first space mission focused on atmospheric CO2 for measuring total column CO, and O2 by detecting the spectral absorption in reflected sunlight. The OCO mission is an essential step, and will yield important new information about atmospheric CO2 distributions. However there are unavoidable limitations imposed by its measurement approach. These include best accuracy only during daytime at moderate to high sun angles, interference by cloud and aerosol scattering, and limited signal from CO2 variability in the lower tropospheric CO2 column. We have been developing a new laser-based technique for the remote measurement of the tropospheric CO2 concentrations from orbit. Our initial goal is to demonstrate a lidar technique and instrument technology that will permit measurements of the CO2 column abundance in the lower troposphere from aircraft. Our final goal is to develop a space instrument and mission approach for active measurements of the CO2 mixing ratio at the 1-2 ppmv level. Our technique is much less sensitive to cloud and atmospheric scattering conditions and would allow continuous measurements of CO2 mixing ratio in the lower troposphere from orbit over land and ocean surfaces during day and night. Our approach is to use the 1570nm CO2 band and a 3-channel laser absorption spectrometer (i.e. lidar used an altimeter mode), which continuously measures at nadir from a near polar circular orbit. The approach directs the narrow co-aligned laser beams from the instrument's lasers toward nadir, and measures the energy of the laser echoes reflected from land and water surfaces. It uses several tunable fiber laser transmitters which allowing measurement of the extinction from a single selected CO2 absorption line in the 1570 nm band. This band is free from interference from other gases and has temperature insensitive absorption lines. During the measurement the lasers are tuned on- and off- a selected CO2 line near 1572 nm and a selected O2 line near 768 nm in the Oxygen A band at kHz rates. The lasers use tunable diode seed lasers followed by fiber amplifiers, and have spectral widths much narrower than the gas absorption lines. The receiver uses a 1-m diameter telescope and photon counting detectors and measures the background light and energies of the laser echoes from the surface. The extinction and column densities for the CO2 and O2 gases are estimated from the ratio of the on and offline surface echo via the differential optical absorption technique. Our technique rapidly alternates between several on-line wavelengths set to the sides of the selected gas absorption lines. It exploits the atmospheric pressure broadening of the lines to weight the measurement sensitivity to the atmospheric column below 5 km. This maximizes sensitivity to CO2 in the boundary layer, where variations caused by surface sources and sinks are largest. Simultaneous measurements of O2 column will use an identical approach with an O2 line. Thee laser frequencies are tunable and have narrow (MHz) line widths. In combination with sensitive photon counting detectors these enables much higher spectral resolution and precision than is possible with passive spectrometer. 1aser backscatter profiles are also measured, which permits identifying measurements made to cloud tops and through aerosol layers. The measurement approach using lasers in common-nadir-zenith path allows retrieving CO2 column mixing ratios in the lower troposphere irrespective of sun angle. Pulsed laser signals, time gated receiver and a narrow receiver field-of-view are used to isolate the surface laser echo signals and to exclude photons scattered from clouds and aerosols. Nonetheless, the optical absorption change due to a change of a few ppO2 is small, <1 % which makes achieving the needed measurement sensitivities and stabilities quite challenging. Measurement SNRs and stabilities of >600:1 are needed to estimate CO2 mixing ratio at the 1-2 ppm level. We have calculated characteristics of the technique and have demonstrated aspects of the laser, detector and receiver approaches in th e laboratory We have also measured O2 in an absorption cell, and made C02 measurements over a 400 m long (one way) horizontal path using a sensor breadboard. We will describe these and more details of our approach in the paper.

Abshire, James B.↗

NASA NTRS2005

Issues and Methods Concerning the Evaluation of Hypersingular and Near-Hypersingular Integrals in BEM Formulations

It is known that higher order modeling of the sources and the geometry in Boundary Element Modeling (BEM) formulations is essential to highly efficient computational electromagnetics. However, in order to achieve the benefits of hIgher order basis and geometry modeling, the singular and near-singular terms arising in BEM formulations must be integrated accurately. In particular, the accurate integration of near-singular terms, which occur when observation points are near but not on source regions of the scattering object, has been considered one of the remaining limitations on the computational efficiency of integral equation methods. The method of singularity subtraction has been used extensively for the evaluation of singular and near-singular terms. Piecewise integration of the source terms in this manner, while manageable for bases of constant and linear orders, becomes unwieldy and prone to error for bases of higher order. Furthermore, we find that the singularity subtraction method is not conducive to object-oriented programming practices, particularly in the context of multiple operators. To extend the capabilities, accuracy, and maintainability of general-purpose codes, the subtraction method is being replaced in favor of the purely numerical quadrature schemes. These schemes employ singularity cancellation methods in which a change of variables is chosen such that the Jacobian of the transformation cancels the singularity. An example of the sin,oularity cancellation approach is the Duffy method, which has two major drawbacks: 1) In the resulting integrand, it produces an angular variation about the singular point that becomes nearly-singular for observation points close to an edge of the parent element, and 2) it appears not to work well when applied to nearly-singular integrals. Recently, the authors have introduced the transformation u(x(prime))= sinh (exp -1) x(prime)/Square root of ((y prime (exp 2))+ z(exp 2) for integrating functions of the form I = Integral of (lambda(r(prime))((e(exp -jkR))/(4 pi R) d D where A (r (prime)) is a vector or scalar basis function and R = Square root of( (x(prime)(exp2) + (y(prime)(exp2) + z(exp 2)) is the distance between source and observation points. This scheme has all of the advantages of the Duffy method while avoiding the disadvantages listed above. In this presentation we will survey similar approaches for handling singular and near-singular terms for kernels with 1/R(exp 2) type behavior, addressing potential pitfalls and offering techniques to efficiently handle special cases.

Fink, P. W.↗

NASA NTRS2004

Mars Express: The Exploration of Phobos

The origin, structure and composition of Phobos are essentially still unresolved issue, with major outcomes for understanding the origin and evolution of the Solar system, and of the Mars system. Several missions have been designed with this unique body the only target. Is Phobos a captured small body (and in this case, can we identify the degree of differentiation such a small object has undergone ?), or is there any co-generic process with Mars involved ? Are there still volatile species trapped, and if so where, and in what form ? What are the reciprocal contributions of Mars and Phobos material accreted in the other body? The number of key clues for deciphering the early solar system processes, the Mars-Phobos binary evolution, and the role Phobos could play in the future of Mars exploration (including human expeditions), are numerous, and exciting for a large community. Part of the answers are in the high resolution coverage of this object (optical and IR, for composition variation), and in its full spectral analysis, from UV to thermal IR : Mars Express is the unique mission to provide this opportunity, and orbit 756, August 22, a totally unique opportunity to approach Phobos so closely. OMEGA should provide all across Stickney, the deepest accessible material, and over more than 20 km along, a spatial resolution less than 200 m, and assess the composition (and possibly their variations) wrt silicates (mafic materials), oxides, and eventually hydrated minerals and organics. By no means we should miss this opportunity: no one would understand the rationale not to turn the instruments ON while so close, and pointing them towards the Phobos surface.

Phobos↗