Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shape analytics”

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 1,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

A Process Model of Interpersonal Relationship Formation in Isolated and Confined Environments

Future space crews will face several challenges such as living and working in a confined environment, isolated from others. These circumstances increase the importance of interpersonal compatibility, teamwork, and team performance. The interpersonal compatibility of space crews has been and continues to be of interest to both NASA and the Institute of Biomedical Problems (IBMP), whose research informs operations for Roscosmos. In NASA-sponsored research, interpersonal compatibility has been examined in terms of how the combination of team members’ personality traits, values, and demographics shape team member relationships and team performance overtime. IBMP-sponsored research mostly has moved away from trait-based approaches toward an idiographic (in-depth, heavily descriptive) approach to researching crew interpersonal relations. This research uses software such as Personal Self-Perception and Attitudes (PSPA), network approaches to team member relations, and content analysis of interactions to assess interpersonal compatibility, and infer states and team dynamics. Our research program integrated these ideas. Our primary research aim was to develop and empirically test a process model of interpersonal relationship formation in isolated and confinement environments. We created a model, collected data from teams in an isolated and confined environment, and applied a novel analytical strategy that combines trait, state, and interaction data (i.e., relational events). We focused specifically on the formation of strained relationships in isolation and confinement—or with whom fellow crewmembers find it difficult to work.

S T Bell↗
NASA NTRS1990

Analytical and experimental study of the vibration of bonded beams with a lap joint

A theoretical model to study the flexural vibration of a bonded lap joint system is described in this paper. First, equations of motion at the joint region are derived using a differential element approach. The transverse displacements of the upper and lower beam are considered to be different. The adhesive is assumed to be linearly viscoelastic and the widely used Kelvin-Voight model is used to represent the viscoelastic behavior of the adhesive. The shear force at the interface between the adhesive and the beam is obtained from the simple bending motion equations of the two beams. The resulting equations of motion are combined with the equations of transverse vibration of the beams in the unjointed regions. These are later solved as a boundary value problem to obtain the eigenvalues and eigenvectors of the system. The model can be used to predict the natural frequencies, modal damping ratios, and mode shapes of the system for free vibration. Good agreement between numerical and experimental results was obtained for a system of graphite epoxy beams lap-jointed by an epoxy adhesive.

Rao, M. D.↗
NASA NTRS1975

A computerized symbolic integration technique for development of triangular and quadrilateral composite shallow-shell finite elements

Computerized symbolic integration was used in conjunction with group-theoretic techniques to obtain analytic expressions for the stiffness, geometric stiffness, consistent mass, and consistent load matrices of composite shallow shell structural elements. The elements are shear flexible and have variable curvature. A stiffness (displacement) formulation was used with the fundamental unknowns consisting of both the displacement and rotation components of the reference surface of the shell. The triangular elements have six and ten nodes; the quadrilateral elements have four and eight nodes and can have internal degrees of freedom associated with displacement modes which vanish along the edges of the element (bubble modes). The stiffness, geometric stiffness, consistent mass, and consistent load coefficients are expressed as linear combinations of integrals (over the element domain) whose integrands are products of shape functions and their derivatives. The evaluation of the elemental matrices is divided into two separate problems - determination of the coefficients in the linear combination and evaluation of the integrals. The integrals are performed symbolically by using the symbolic-and-algebraic-manipulation language MACSYMA. The efficiency of using symbolic integration in the element development is demonstrated by comparing the number of floating-point arithmetic operations required in this approach with those required by a commonly used numerical quadrature technique.

Anderson, C. M.↗
NASA NTRS1993

Experimental unsteady pressures at flutter on the Supercritical Wing Benchmark Model

This paper describes selected results from the flutter testing of the Supercritical Wing (SW) model. This model is a rigid semispan wing having a rectangular planform and a supercritical airfoil shape. The model was flutter tested in the Langley Transonic Dynamics Tunnel (TDT) as part of the Benchmark Models Program, a multi-year wind tunnel activity currently being conducted by the Structural Dynamics Division of NASA Langley Research Center. The primary objective of this program is to assist in the development and evaluation of aeroelastic computational fluid dynamics codes. The SW is the second of a series of three similar models which are designed to be flutter tested in the TDT on a flexible mount known as the Pitch and Plunge Apparatus. Data sets acquired with these models, including simultaneous unsteady surface pressures and model response data, are meant to be used for correlation with analytical codes. Presented in this report are experimental flutter boundaries and corresponding steady and unsteady pressure distribution data acquired over two model chords located at the 60 and 95 percent span stations.

Dansberry, Bryan E.↗
NASA NTRS2010

High-Precision Shape Control of In-Space Deployable Large Membrane/Thin-Shell Reflectors

This innovation has been developed to improve the resolutions of future spacebased active and passive microwave antennas for earth-science remote sensing missions by maintaining surface figure precisions of large membrane/thin-shell reflectors during orbiting. The intention is for these sensing instruments to be deployable at orbit altitudes one or two orders of magnitude higher than Low Earth Orbit (LEO), but still being able to acquire measurements at spatial resolution and sensitivity similar to those of LEO. Because active and passive microwave remote sensors are able to penetrate through clouds to acquire vertical profile measurements of geophysical parameters, it is desirable to elevate them to the higher orbits to obtain orbital geometries that offer large spatial coverage and more frequent observations. This capability is essential for monitoring and for detailed understanding of the life cycles of natural hazards, such as hurricanes, tropical storms, flash floods, and tsunamis. Major components of this high-precision antenna-surface-control system include a membrane/thin shell reflector, a metrology sensor, a controller, actuators, and corresponding power amplifier and signal conditioning electronics (see figure). Actuators are attached to the back of the reflector to produce contraction/ expansion forces to adjust the shape of the thin-material reflector. The wavefront-sensing metrology system continuously measures the surface figure of the reflector, converts the surface figure to digital data and feeds the data to the controller. The controller determines the control parameters and generates commands to the actuator system. The flexible, piezoelectric polymer actuators are thus activated, providing the control forces needed to correct any distortions that exist in the reflector surface. Piezoelectric polymer actuators are very thin and flexible. They can be implemented on the back of the membrane/thin-shell reflector without introducing significant amounts of mass or stiffness to the reflector. They can be rolled up or folded to accommodate the packaging needed for launch. An analytical model of the system, which includes the membrane reflector, actuator, and controller has been developed to investigate the functionality of this control system on a 35-meter-diameter membrane reflector. The performance of this system under external disturbances such as in space thermal loads and W-error due to inflation has been investigated. A subscale breadboard has been developed, and the functionality of this control concept has been demonstrated by this breadboard.

Watkins, Ronald↗
NASA NTRS2000

Applications of Fractal Analytical Techniques in the Estimation of Operational Scale

The observational scale and the resolution of remotely sensed imagery are essential considerations in the interpretation process. Many atmospheric, hydrologic, and other natural and human-influenced spatial phenomena are inherently scale dependent and are governed by different physical processes at different spatial domains. This spatial and operational heterogeneity constrains the ability to compare interpretations of phenomena and processes observed in higher spatial resolution imagery to similar interpretations obtained from lower resolution imagery. This is a particularly acute problem, since longterm global change investigations will require high spatial resolution Earth Observing System (EOS), Landsat 7, or commercial satellite data to be combined with lower resolution imagery from older sensors such as Landsat TM and MSS. Fractal analysis is a useful technique for identifying the effects of scale changes on remotely sensed imagery. The fractal dimension of an image is a non-integer value between two and three which indicates the degree of complexity in the texture and shapes depicted in the image. A true fractal surface exhibits self-similarity, a property of curves or surfaces where each part is indistinguishable from the whole, or where the form of the curve or surface is invariant with respect to scale. Theoretically, if the digital numbers of a remotely sensed image resemble an ideal fractal surface, then due to the self-similarity property, the fractal dimension of the image will not vary with scale and resolution, and the slope of the fractal dimension-resolution relationship would be zero. Most geographical phenomena, however, are not self-similar at all scales, but they can be modeled by a stochastic fractal in which the scaling properties of the image exhibit patterns that can be described by statistics such as area-perimeter ratios and autocovariances. Stochastic fractal sets relax the self-similarity assumption and measure many scales and resolutions to represent the varying form of a phenomenon as the pixel size is increased in a convolution process. We have observed that for images of homogeneous land covers, the fractal dimension varies linearly with changes in resolution or pixel size over the range of past, current, and planned space-borne sensors. This relationship differs significantly in images of agricultural, urban, and forest land covers, with urban areas retaining the same level of complexity, forested areas growing smoother, and agricultural areas growing more complex as small pixels are aggregated into larger, mixed pixels. Images of scenes having a mixture of land covers have fractal dimensions that exhibit a non-linear, complex relationship to pixel size. Measuring the fractal dimension of a difference image derived from two images of the same area obtained on different dates showed that the fractal dimension increased steadily, then exhibited a sharp decrease at increasing levels of pixel aggregation. This breakpoint of the fractal dimension/resolution plot is related to the spatial domain or operational scale of the phenomenon exhibiting the predominant visible difference between the two images (in this case, mountain snow cover). The degree to which an image departs from a theoretical ideal fractal surface provides clues as to how much information is altered or lost in the processes of rescaling and rectification. The measured fractal dimension of complex, composite land covers such as urban areas also provides a useful textural index that can assist image classification of complex scenes.

Emerson, Charles W.↗
NASA NTRS2016

Cost-Benefit Analysis for the Advanced Near Net Shape Technology (ANNST) Method for Fabricating Stiffened Cylinders

Low Technology Readiness Levels (TRLs) and high levels of uncertainty make it challenging to develop cost estimates of new technologies in the R&D phase. It is however essential for NASA to understand the costs and benefits associated with novel concepts, in order to prioritize research investments and evaluate the potential for technology transfer and commercialization. This paper proposes a framework to perform a cost-benefit analysis of a technology in the R&D phase. This framework was developed and used to assess the Advanced Near Net Shape Technology (ANNST) manufacturing process for fabricating integrally stiffened cylinders. The ANNST method was compared with the conventional multi-piece metallic construction and composite processes for fabricating integrally stiffened cylinders. Following the definition of a case study for a cryogenic tank cylinder of specified geometry, data was gathered through interviews with Subject Matter Experts (SMEs), with particular focus placed on production costs and process complexity. This data served as the basis to produce process flowcharts and timelines, mass estimates, and rough order-of-magnitude cost and schedule estimates. The scalability of the results was subsequently investigated to understand the variability of the results based on tank size. Lastly, once costs and benefits were identified, the Analytic Hierarchy Process (AHP) was used to assess the relative value of these achieved benefits for potential stakeholders. These preliminary, rough order-of-magnitude results predict a 46 to 58 percent reduction in production costs and a 7-percent reduction in weight over the conventional metallic manufacturing technique used in this study for comparison. Compared to the composite manufacturing technique, these results predict cost savings of 35 to 58 percent; however, the ANNST concept was heavier. In this study, the predicted return on investment of equipment required for the ANNST method was ten cryogenic tank barrels when compared with conventional metallic manufacturing. The AHP study results revealed that decreased final cylinder mass and improved quality assurance were the most valued benefits of cylinder manufacturing methods, therefore emphasizing the relevance of the benefits achieved with the ANNST process for future projects.

Ivanco, Marie L.↗
NASA NTRS1992

Crack-mouth displacements for semielliptical surface cracks subjected to remote tension and bending loads

The exact analytical solution for an embedded elliptical crack in an infinite body subjected to arbitrary loading was used in conjunction with the finite element alternating method to obtain crack-mouth-opening displacements (CMOD) for surface cracks in finite plates subjected to remote tension. Identical surface-crack configurations were also analyzed with the finite element method using 20-noded element for plates subjected to both remote tension and bending. The CMODs from these two methods generally agreed within a few percent of each other. Comparisons made with experimental results obtained from surface cracks in welded aluminum alloy specimens subjected to tension also showed good agreement. Empirical equations were developed for CMOD for a wide range of surface-crack shapes and sizes subjected to tension and bending loads. These equations were obtained by modifying the Green-Sneddon exact solution for an elliptical crack in an infinite body to account for finite boundary effects. These equations should be useful in monitoring surface-crack growth in tests and in developing complete crack-face-displacement equations for use in three-dimensional weight-function methods.

Raju, Ivatury S.↗
NASA NTRS1998

Dynamics of Liquids in Edges and Corners (DYLCO): IML-2 Experiment for the BDPU

Knowledge of the behavior of fluids possessing free surfaces is important to many fluid systems, particularly in space, where the normally subtle effects of surface wettability play a more dramatic and often surprising role. DYLCO for the IML-2 mission was proposed as a simple experiment to probe the particular behavior of capillary surfaces in containers of irregular cross section. Temperature control was utilized to vary the fluid-solid contact angle, a questionable thermodynamic parameter of the system, small changes in which can dramatically influence the configuration, stability, and flow of a capillary surface. Container shapes, test fluid, and temperature ranges were selected for observing both local changes in interface curvature as well as a global change in fluid orientation due to a critical wetting phenomenon. The experiment hardware performed beyond what was expected and fluid interfaces could be readily digitized post flight to show the dependence of the interface curvature on temperature. For each of the containers tested surfaces were observed which did not satisfy the classic equations for the prediction of interface shape with constant contact angle boundary condition. This is explained by the presence of contact angle hysteresis arising from expansion and contraction of the liquid during the heating and cooling steps of the test procedure. More importantly, surfaces exceeding the critical surface curvature required for critical wetting were measured, yet no wetting was observed. These findings are indeed curious and pose key questions concerning the role of hysteresis for this critical wetting phenomena. The stability of such surfaces was determined numerically and it is shown that stability is enhance (reduced) when a surface is in its 'advancing' ('receding') state, The analysis shows complete instability as the critical wetting condition is reached. The case of ideal dynamic wetting is addressed analytically in detail with results of significant flow characteristics presented in closed form. The solutions indicate a square root of T dependence of the capillary 'rise' rate which is corroborated by drop tower tests. The analysis clearly shows that infinite time is necessary for surfaces to reorient at the critical wetting transition.

Langbein, D.↗
NASA NTRS1985

Calculation of particle dynamics effects on laser velocimeter data

In order to analytically assess the effects of particle dynamics on laser velocimeter measurements, the particle response to the fluid in which it is entrained must be modeled. The equation of particle motion derived at Arnold Engineering and Development Center (AEDC), is used to model the motion of the particle. The equation is derived from the Navier-Stokes equation under the assumption of unsteady Stokes flow. It should be noted that the Stokes flow eliminates the need for a far field boundary condition, and hence eliminates any normal (lift) forces. To extend the equation to higher particle slip velocities, the Stokes drag is replaced with the particle drag law, modified to include heat transfer effects. The drag coefficient accounts for inertial, rarefaction, compressibility, shape, and heat transfer effects.

Nichols, R. H.↗
NASA NTRS1986

Evaluating space station applications of automation and robotics technologies from a human productivity point of view

The role that automation, robotics, and artificial intelligence will play in Space Station operations is now beginning to take shape. Although there is only limited data on the precise nature of the payoffs that these technologies are likely to afford there is a general consensus that, at a minimum, the following benefits will be realized: increased responsiveness to innovation, lower operating costs, and reduction of exposure to hazards. Nevertheless, the question arises as to how much automation can be justified with the technical and economic constraints of the program? The purpose of this paper is to present a methodology which can be used to evaluate and rank different approaches to automating the functions and tasks planned for the Space Station. Special attention is given to the impact of advanced automation on human productivity. The methodology employed is based on the Analytic Hierarchy Process. This permits the introduction of individual judgements to resolve the confict that normally arises when incomparable criteria underly the selection process. Because of the large number of factors involved in the model, the overall problem is decomposed into four subproblems individually focusing on human productivity, economics, design, and operations, respectively. The results from each are then combined to yield the final rankings. To demonstrate the methodology, an example is developed based on the selection of an on-orbit assembly system. Five alternatives for performing this task are identified, ranging from an astronaut working in space, to a dexterous manipulator with sensory feedback. Computational results are presented along with their implications. A final parametric analysis shows that the outcome is locally insensitive to all but complete reversals in preference.

Bard, J. F.↗
NASA NTRS1990

The method of lines in analyzing solids containing cracks

A semi-numerical method is reviewed for solving a set of coupled partial differential equations subject to mixed and possibly coupled boundary conditions. The line method of analysis is applied to the Navier-Cauchy equations of elastic and elastoplastic equilibrium to calculate the displacement distributions in various, simple geometry bodies containing cracks. The application of this method to the appropriate field equations leads to coupled sets of simultaneous ordinary differential equations whose solutions are obtained along sets of lines in a discretized region. When decoupling of the equations and their boundary conditions is not possible, the use of a successive approximation procedure permits the analytical solution of the resulting ordinary differential equations. The use of this method is illustrated by reviewing and presenting selected solutions of mixed boundary value problems in three dimensional fracture mechanics. These solutions are of great importance in fracture toughness testing, where accurate stress and displacement distributions are required for the calculation of certain fracture parameters. Computations obtained for typical flawed specimens include that for elastic as well as elastoplastic response. Problems in both Cartesian and cylindrical coordinate systems are included. Results are summarized for a finite geometry rectangular bar with a central through-the-thickness or rectangular surface crack under remote uniaxial tension. In addition, stress and displacement distributions are reviewed for finite circular bars with embedded penny-shaped cracks, and rods with external annular or ring cracks under opening mode tension. The results obtained show that the method of lines presents a systematic approach to the solution of some three-dimensional mechanics problems with arbitrary boundary conditions. The advantage of this method over other numerical solutions is that good results are obtained even from the use of a relatively coarse grid.

Gyekenyesi, John P.↗
NASA NTRS1979

A review and assessment of the stress-intensity factors for surface cracks

The stress-intensity factor solutions proposed for a surface crack in a finite plate subjected to uniform tension are reviewed. Fourteen different solutions obtained over the past 16 years using approximate analytical methods, experimental methods, and engineering estimates are compared. The accuracy of the various solutions is assessed by correlating fracture data on surface-cracked tension specimens made of a brittle epoxy material. Fracture of the epoxy material was characterized by a constant value of stress-intensity factor as failure. Thus, the correctness of various solutions is judged by the variations in the stress-intensity factors at failure. The solutions were ranked in order of minimum standard deviation. The highest ranking solutions correlated 95% of data analyzed within + or - 10% whereas the lowest ranking solutions correlated 95% of data analyzed within + or - 20%. Some solutions could be applied to all data considered, whereas others were limited with respect to crack shapes and crack sizes that could be analyzed.

Newman, J. C., Jr.↗
NASA NTRS1995

FTIR instrumentation to monitor vapors from Shuttle tile waterproofing materials

The Space Shuttle Thermal Protection System (TPS) tiles and blankets are waterproofed using DimethylEthoxySilane (DMEX) in the Orbiter Processing Facilities (OPF). DMES has a Threshold Limit Value (TLV) for exposure of personnel to vapor concentration in air of 0.5 ppm. The OPF high bay cannot be opened for normal work after a waterproofing operation until the DMES concentration is verified by measurement to be below the TLV. On several occasions the high bay has been kept closed for up to 8 hours following waterproofing operations due to high DMES measurements. In addition, the Miran 203 and Miran 1 BX infrared analyzers calibrated at different wavelengths gave different readings under the same conditions. There was reason to believe that some of the high DMES concentration readings were caused by interference form water and ethanol vapors. The Toxic Vapor Detection Laboratory (TVDL) was asked to test the existing DMES instruments and identify the best qualified instrument. In addition the TVDL was requested to develop instrumentation to ensure the OPF high bay could be opened safely as soon as possible after a waterproofing operation. A Fourier Transform Infrared (FTIR) spectrophotometer instrument developed for an earlier project was reprogrammed to measure DMES vapor along with ethanol, water, and several common solvent vapors. The FTIR was then used to perform a series of laboratory and field tests to evaluate the performance of the single wavelength IR instruments in use. The results demonstrated that the single wavelength IR instruments did respond to ethanol and water vapors, more or less depending on the analytical IR wavelength selected. The FTIR was able to separate the responses to DMES, water and ethanol, and give consistent readings for the DMES vapor concentration. The FTIR was then deployed to the OPF to monitor real waterproofing operations. The FTIR was also used to measure the time for DMES to evaporate from TPS tile under a range of humidity conditions in controlled laboratory tests. The combination of laboratory and field tests with the FTIR instrument demonstrated superior sensitivity, ability to reject interference from water and ethanol vapors, ruggedness to be transported from the lab to the OPF and set up without special procedures or degradation of performance. The multiple component vapor analysis algorithm was developed at KSC and incorporates automatic baseline correction and shape fitting of the spectra. The analysis for DMES, TetraMethylDiSiloxane (TMDS), ethanol, methanol, isopropanol, and baseline parameters uses 161 points per sample at 4 cm(exp -1) resolution, and processes an eight scan sample every ten seconds. The standard deviation of the measurements is 0.013 ppm and the upper linear limit is 125 ppm DMES. Based on successful demonstration of capabilities we produced three mobile instrument carts to be used in each OPF to support future waterproofing operations. The design and building of the 'DMES Carts' were accomplished in Fiscal year 1995.

Mattson, C. B.↗
NASA NTRS1996

Turbulence Scales, Rise Times, Caustics, and the Simulation of Sonic Boom Propagation

The general topic of atmospheric turbulence effects on sonic boom propagation is addressed with especial emphasis on taking proper and efficient account of the contributions of the portion oi the turbulence that is associated with extremely high wavenumber components. The recent work reported by Bart Lipkens in his doctoral thesis is reexamined to determine whether the good agreement between his measured rise times with the 1971 theory of the author is fortuitous. It is argued that Lipken's estimate of the distance to the first caustic was a gross overestimate because of the use of a sound speed correlation function shaped like a gaussian curve. In particular, it is argued that the expected distance to the first caustic varies with the kinematic viscosity nu and the energy epsilon dissipated per unit mass per unit time, and the sound speed c as : d(sub first caustic) = nu(exp 7/12) c(exp 2/3)/ epsilon(exp 5/12)(nu x epsilon/c(exp 4))(exp a), where the exponent a is greater than -7/12 and can be argued to be either O or 1/24. In any event, the surprising aspect of the relationship is that it actually goes to zero as the viscosity goes to zero with s held constant. It is argued that the apparent overabundance of caustics can be grossly reduced by a general computational and analytical perspective that partitions the turbulence into two parts, divided by a wavenumber k(sub c). Wavenumbers higher than kc correspond to small-scale turbulence, and the associated turbulence can be taken into account by a renormalization of the ambient sound speed so that the result has a small frequency dependence that results from a spatial averaging over of the smaller-scale turbulent fluctuations. Selection of k(sub c). can be made so large that only a very small number of caustics are encountered if one adopts the premise that the frequency dispersion of pulses is caused by that part of the turbulence spectrum which lies in the inertial range originally predicted by Kolmogoroff. The acoustic propagating wave's dispersion relation has the acoustic wavenumber being of the form k = (omega/c) + F(omega), where c is a spatially averaged sound speed and where, for mechanical turbulence, the extra term F(omega) must depend on only the angular frequency omega, the sound speed c, and the turbulent energy dissipation epsilon per unit fluid mass and per unit time. If the turbulence is weak, then the quantity F(omega) has to be of second order in the portions of the turbulent fluid velocity in the inertial range, so, following Kolmogoroff's reasoning, it must vary with epsilon as epsilon(exp 2/3). Simple dimensional analysis then reveals that F(omega) is K epsilon(exp 2/3) c(exp -7/3) omega(exp l/3), K being a universal dimensionless complex constant.

Pierce, Allan D.↗
NASA NTRS1979

Minimizing pattern effects in semiconductor lasers at high rate pulse modulation

This paper deals with analytical and experimental work related to modulation of a semiconductor laser used in high bit rate communication. The approach is based upon minimizing the charge storage effect by a proper choice of the area of the modulating pulses and the bias current. The concept of using additional current pulses to probe for variations in electron density between pulses is investigated. The primary limitation on bit rate is found to be the ability to generate laser drive pulses free of ringing or similar transients. This and the 300 ps pulsewidth, an experimental constraint, limit the bit rate to about 1 Gbit/s. However, by using this approach it appears that bit rates considerably higher than 2 Gbits/s could be reached with sufficiently accurate control of drive pulse shape. It is found that the laser bias and the current pulse area had to be controlled within 1 percent and 10 percent, respectively.

Torphammar, P.↗
NASA NTRS1993

Laplace's equation and the Dirichlet-Neumann map in multiply connected domains

A variety of problems in material science and fluid dynamics require the solution of Laplace's equation in multiply connected domains. Integral equation methods are natural candidates for such problems, since they discretize the boundary alone, require no special effort for free boundaries, and achieve superalgebraic convergence rates on sufficiently smooth domains in two space dimensions, regardless of shape. Current integral equation methods for the Dirichlet problem, however, require the solution of M independent problems of dimension N, where M is the number of boundary components and N is the total number of points in the discretization. In this paper, we present a new boundary integral equation approach, valid for both interior and exterior problems, which requires the solution of a single linear system of dimension N + M. We solve this system by making use of an iterative method (GMRES) combined with the last multipole method for the rapid calculation of the necessary matrix vector products. For a two-dimensional system with 200 components and 100 points on each boundary, we gain a speedup of a factor of 100 from the new analytic formulation and a factor of 50 from the fast multipole method. The resulting scheme brings large scale calculations in extremely complex domains within practical reach.

Greenbaum, A.↗
NASA NTRS2007

White Light Demonstration of One Hundred Parts per Billion Irradiance Suppression in Air by New Starshade Occulters

A new mission concept for the direct imaging of exo-solar planets called the New Worlds Observer (NWO) has been proposed. The concept involves flying a meter-class space telescope in formation with a newly-conceived, specially-shaped, deployable star-occulting shade several meters across at a separation of some tens of thousands of kilometers. The telescope would make its observations from behind the starshade in a volume of high suppression of incident irradiance from the star around which planets orbit. The required level of irradiance suppression created by the starshade for an efficacious mission is of order 0.1 to 10 parts per billion in broadband light. This paper discusses the experimental setup developed to accurately measure the suppression ratio of irradiance produced at the null position behind candidate starshade forms to these levels. It also presents results of broadband measurements which demonstrated suppression levels of just under 100 parts per billion in air using the Sun as a light source. Analytical modeling of spatial irradiance distributions surrounding the null are presented and compared with photographs of irradiance captured in situ behind candidate starshades.

Levinton, Douglas B.↗