SEARCH · Search NASA
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.
Huygen principle for electromagnetic field in moving, isotropic, homogeneous and linear medium, using Maxwell-Minkowski equations and Green function for field equation
If a set of N points or nodes with a nonnegative cost associated with each ordered pair is known, it is desired to find a path from one given node to another given node which minimizes the cost sum. An algorithm is presented which yields a global minimum solution after at most N - 1 iterations or on a typical large third-generation computer, after 1 hour of computation time for a 10,000-node problem. The rapid-access data storage capacity demanded by the algorithm is approximately 3N words for costs read in from slow-access storage or 2N words for calculable costs. The time-storage requirements of the algorithm known to the authors. When the problem is viewed as a discretized optimal control problem, after N-1 iterations, an optimal control or node transition is established for each of the N nodes or states; thus, the algorithm can be applied to situations were there may be errors in the control that necessitate a closed loop control that necessitate a closed loop control philosophy.
Huygen principle for analysis of electromagnetic field in moving isotropic homogeneous and linear medium
The shape of a wave carries all information about the spatial and temporal structure of its source, given that the medium and its properties are known. Most modern imaging methods seek to utilize this nature of waves originating from Huygens' principle. We discuss the retrieval of the complete kinetic energy distribution from the acoustic trace that is recorded when a short ion bunch deposits its energy in water. This novel method, which we refer to as Ion-Bunch Energy Acoustic Tracing (I-BEAT), is a refinement of the ionoacoustic approach. With its capability of completely monitoring a single, focused proton bunch with prompt readout and high repetition rate, I-BEAT is a promising approach to meet future requirements of experiments and applications in the field of laser-based ion acceleration. We demonstrate its functionality at two laser-driven ion sources for quantitative online determination of the kinetic energy distribution in the focus of single proton bunches.
The acoustic fields and levitation forces produced along the axis of a single-axis resonance system were measured. The system consisted of a St. Clair generator and a planar reflector. The levitation force was measured for bodies of various sizes and geometries (i.e., spheres, cylinders, and discs). The force was found to be roughly proportional to the volume of the body until the characteristic body radius reaches approximately 2/k (k = wave number). The acoustic pressures along the axis were modeled using Huygens principle and a method of imaging to approximate multiple reflections. The modeled pressures were found to be in reasonable agreement with those measured with a calibrated microphone.
Flow fields involving turbulent flames in premixed gases under a variety of conditions are modeled by the use of a numerical technique based on the random vortex method to solve the Navier-Stokes equations and a flame propagation algorithm to trace the motion of the front and implement the Huygens principle, both due to Chorin. A successive over-relaxation hybrid method is applied to solve the Euler equation for flows in an arbitrarily shaped domain. The method of images, conformal transformation, and the integral-equation technique are also used to treat flows in special cases, according to their particular requirements. Salient features of turbulent flame propagation in premixed gases are interpreted by relating them to the aerodynamic properties of the flow field. Included among them is the well-known cellular structure of flames stabilized by bluff bodies, as well as the formation of the characteristic tulip shape of flames propagating in ducts. In its rudimentary form, the mechanism of propagation of a turbulent flame is shown to consist of: (1) rotary motion of eddies at the flame front, (2) self-advancement of the front at an appropriate normal burning speed, and (3) dynamic effects of expansion due to exothermicity of the combustion reaction. An idealized model is used to illustrate these fundamental mechanisms and to investigate basic aerodynamic features of flames in premixed gases. The case of a confined flame stabilized behind a rearward-facing step is given particular care and attention. Solutions are shown to be in satisfactory agreement with experimental results, especially with respect to global properties such as the average velocity profiles and reattachment length.
Space exploration and research require large amounts of power. Solar pumped laser systems have been shown to have the potential for meeting the performance requirements necessary for power transmission in space. The successful design of laser systems for these applications will depend, in part, on having a clear understanding of the development of the dynamical processes in the laser cavity and on the effects that changes in physical and design parameters have on laser performance. In particular, it is necessary to know the amplitude and phase distributions of the laser beam at the output aperture when steady state operation is achieved in order to determine the far-field power distribution. The output from the laser will depend on the active medium, the optical environment of the active material and on the gain distribution in the active region as laser action builds up and reaches steady state. An important component in the design process will be a realistic model of the active laser cavity. A computer model of the laser cavity, based on Huygens' principle was developed. The code calculates the amplitude and phase of an optical wave reflected back and forth between the mirrors of a laser cavity. The original code assumes a gain distribution which does not change with the buildup of oscillations in the cavity. A step in the direction of realism is the inclusion of a saturable gain medium in the cavity. The objective is to incorporate saturation effects into the existing computer model.
In this paper, a unified computational technique is presented which allows the incorporation of the curved frequency-selective surface (FSS) geometry in the computation of the reflector antenna radiation pattern. The scattered fields from an illuminated FSS reflector are formalized using Huygens' principle in such a way that the 'reflecting' and the 'transparent' FSS subreflector cases are treated identically and the thickness of the FSS subreflector remains arbitrary. The analysis utilizes local surface coordinates to describe the reflection/transmission matrices of the FSS subreflector where it is assumed that these matrices are available. In most cases one may use the local tangent plane for approximating the plane of the FSS in the local coordinate surface of the reflector. The paper demonstrates how the local curved coordinate system can be introduced in the diffraction modeling of FSS reflectors and its importance in accurately predicting the side-lobe and crosspolarization levels. Results of numerical simulations are presented for several FSS subreflector configurations.
Extreme ultraviolet (EUV) waves are large-scale propagating disturbances observed in the solar corona, frequently associated with coronal mass ejections and flares. They appear as faint, extended structures propagating from a source region across the structured solar corona. To measure these waves, we have constructed the Automated Wave Analysis and REduction (AWARE) algorithm. AWARE is implemented in two stages. In the first stage, we use simple image processing techniques to isolate the propagating, brightening wave fronts as they move across the corona. In the second stage, AWARE measures the distance, velocity and acceleration of that wave front across the Sun. We explore the use of the Huygens principle, dynamic time warping and a simple parametric representation of the wave front as potential methods for detecting and tracking non-radial EUV wave propagation.
Recent studies have used the superposition principle (extended Huygens-Fresnel principle) to characterize completely the statistics of a field that has propagated through a thick slab of turbulent air in terms of the statistics for spherical-wave sources. In this paper, we consider the normal-mode decomposition associated with this linear system propagation model. In particular, we use the statistics of the atmospheric impulse response (Green's function) to show that the atmospheric mode decomposition exhibits far-field and near-field regimes very similar to those of free-space propagation. The significance of these results for optical communication through the atmosphere is briefly discussed.
The current effort addresses two issues important to the research conducted by the Thermal Radiation Group at Virginia Tech. The first research topic involves the development of a method which can properly model the diffraction of radiation as it enters an instrument aperture. The second topic involves the study of a potential next-generation space-borne radiometric instrument concept. Presented are multiple modeling efforts to describe the diffraction of monochromatic radiant energy passing through an aperture for use in the Monte-Carlo ray-trace environment. Described in detail is a deterministic model based upon Heisenberg's uncertainty principle and the particle theory of light. This method is applicable to either Fraunhofer or Fresnel diffraction situations, but is incapable of predicting the secondary fringes in a diffraction pattern. Also presented is a second diffraction model, based on the Huygens-Fresnel principle with a correcting obliquity factor. This model is useful for predicting Fraunhofer diffraction, and can predict the secondary fringes because it keeps track of phase. NASA is planning for the next-generation of instruments to follow CERES (Clouds and the Earth's Radiant Energy System), an instrument which measures components of the Earth's radiant energy budget in three spectral bands. A potential next-generation concept involves modification of the current CERES instrument to measure in a larger number of wavelength bands. This increased spectral partitioning would be achieved by the addition of filters and detectors to the current CERES geometry. The capacity of the CERES telescope to serve for this purpose is addressed in this thesis.
The paper examines the operation of an optical system forming an image of an incoherent source object under conditions of atmospheric noise. The effect of atmospheric turbulence is considered. The calculations are performed in the Fresnel – Huygens approximation. The characteristics of the image for ‘long’ and ‘short’ exposures are investigated. The issues related to the adaptive correction of the general light wavefront slopes are discussed. The limits of applicability of phase adaptive correction using a reference source are shown. (paper)
Improved correlation function in Kirchhoff- Huygen method of solution of scattering of waves from statistically rough surfaces
The structures of the Fresnel zones of an occulting body with a substantial atmosphere are analyzed. The first Fresnel zone is defined as the set of all contiguous points connected to the transmitter and receiver via paths with a combined length within one half wavelength of the stationary value, and calculated by the application of the Huygens-Fresnel principle to a wave front frozen in time as it traverses the region of interest, taking into account wave front curvature of waves undergoing refraction. Calculations of the Fresnel zones of a planet determined by the meeting of wave fronts represented as tori reveal that the zone along the propagation path traversing the nearest portion of the limb to the transmitter has a shape similar to an ellipse, while that of the far zone resembles the letter X. The phase delay along the farside path is found to be a stationary point rather than a minimum, since delay is distributed in the form of a hyperbolic paraboloid near the ray and in the form of a torus intersection for greater separations. Examples are presented for S band propagation deep in the atmospheres of Venus and Jupiter.
Two boundary element formulations of acoustical behavior exist; the Direct Boundary Element Methods (DBEM) based on the Helmholtz Integral Equations and the Indirect Boundary Element Methods (IBEM) based on Huygen's Principle. In this investigation, both methods are implemented utilizing a simple linear superparametric element. In addition the IBEM is studied using a quadratic isoparametric element. The accuracy and relative efficiency of the various techniques are examined. In order to properly model aircraft interior cavities the additional capability to model wall treatments and internal point sources is added to the methods. The procedures are verified for several well-understood cavity problems. The relative merits of each boundary element method and the finite element method are examined.
Based on three formulations of the Huygen's principle, explicit expressions is given for the far field contribution from a small ray tube. This expression is useful in shooting and bouncing rays for solving complex scattering problems.
Numerical techniques such as Finite Difference Time Domain (FDTD) computer programs, which were first developed to analyze the external electromagnetic environment of an aircraft during a wave illumination, a lightning event, or any kind of current injection, are now very powerful investigative tools. The program called GORFF-VE, was extended to compute the inner electromagnetic fields that are generated by the penetration of the outer fields through large apertures made in the all metallic body. Then, the internal fields can drive the electrical response of a cable network. The coupling between the inside and the outside of the helicopter is implemented using Huygen's principle. Moreover, the spectacular increase of computer resources, as calculations speed and memory capacity, allows the modellization structures as complex as these of helicopters with accuracy. This numerical model was exploited, first, to analyze the electromagnetic environment of an in-flight helicopter for several injection configurations, and second, to design a coaxial return path to simulate the lightning aircraft interaction with a strong current injection. The E field and current mappings are the result of these calculations.
Interferometric synthetic aperture radar (IFSAR) measurements at X-, C-, L-, and P-band are used to derive ground topography at meter level resolution. Interpretation of the derived topography requires attention due to the complex interaction of the radar signal with ground cover. The presence of penetrable surfaces such as vegetation, and tree canopies poses a challenge since the depth of penetration depends on a number of parameters such as the operating radar frequency, polarization, incident angle, as well as terrain structure. The dependence of the reconstructed topography on polarization may lead to the characterization of the ground cover. Simulation of interferometric measurements is useful for interpretation of the derived topography (B. Houshmand, Proceedings of URSI, 314, 1997). In this talk , time domain simulations for interferometric measurement for TE- and TM- polarization are presented. Time domain simulation includes the effects of the surface material property as well geometry comparable the radar signal wavelength (B. Houshmand, Proceedings of the URSI, 25, 1998). The IFSAR simulation is carried out in two steps. First, the forward scattering data is generated based on full wave analysis. Next, the electromagnetic information is inverted to generate surface topography. This inversion is based on the well known IFSAR processing technique which is composed of signal compression, and formation of an interferogram. The full wave forward scattering data is generated by the scattered-field formulation of the FDTD algorithm. The simulation is carried out by exciting the computational domain by a radar signal. The scattered field is then computed and translated to the receiving interferometric antennas using the time-domain Huygen's principle. The inversion process starts by compressing the time-domain data. The range compressed data from both receivers are then coregistered to form an interferogram. The resulting interferogram is then related to the ground topography using the radar imaging geometry. In this talk, the simulation results are compared with the C-band TM IFSAR derived topography, and the TE/TM SAR images at L-Band.