Search NASA⌕ Search

SEARCH · Search NASA

Results for “POINT MATCHING METHOD”

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 343 records · Page 19

Planning Satellite Swarm Measurements for Climate Models: Comparing Dynamic Constraint Processing and MILP Methods

We present D-SHIELD, a challenging climate science application to plan coordinated measurements (observations) for a constellation of satellites, each containing two different sensors, each with 61 pointing angle options. The L-band and P-band radar sensors collect data fed into a soil moisture model which tracks and predicts soil moisture across 1.67 million Ground Positions (GP). Soil moisture is an important predictor of wildfires, and then a predictor of floods, landslides and debris flow after a fire. Each measurement covers multiple GP due to the sensor footprint. Each GP has a "model error" which represents the uncertainty of the the soil moisture state prediction. Model error changes at different rates for each GP as the time since last observation increases and after significant events like rain. The planner's goal is to select measurements which maximize soil moisture model improvement (reduce model uncertainty). This problem is combinatorically explosive, involving many degrees of freedom for planner choices. Good domain heuristics can find solutions within a reasonable time for our application needs but cannot be proven optimal. In this paper we compare two different planning approaches to this problem: Dynamic Constraint Processing (DCP) and Mixed Integer Linear Programming (MILP). We match inputs and metrics for both DCP and MILP algorithms to enable a direct apples-to-apples comparison. We demonstrate and discuss the trades between DCP flexibility and performance vs. MILP's promise of provable optimality.

Rich Levinson↗

MANN: A program to transfer designs for diffractive optical elements to a MANN photolithographic mask generator

There are two basic areas of interest for diffractive optics. In the first, the property of wavefront division is exploited for achieving optical fanout, analogous to the more familiar electrical fanout of electronic circuitry. The basic problem here is that when using a simple uniform diffraction grating the energy input is divided unevenly among the output beams. The other area of interest is the use of diffractive elements to replace or supplement standard refractive elements such as lenses. Again, local grating variations can be used to control the amount of bending imparted to optical rays, and the efficiency of the diffractive element will depend on how closely the element can be matched to the design requirements. In general, production restrictions limit how closely the element approaches the design, and for the common case of photolithographic production, a series of binary masks is required to achieve high efficiency. The actual design process is much more involved than in the case of elements for optical fanout, as the desired phase of the optical wavefront over some reference plane must be specified and the phase alteration to be introduced at each point by the diffraction element must be known. This generally requires the utilization of a standard optical design program. Two approaches are possible. In the first approach, the diffractive element is treated as a special type of lens and the ordinary optical design equations are used. Optical design programs tend to follow a second approach, namely, using the equations of optical interference derived from holographic theory and then allowing the introduction of phase front corrections in the form of polynomial equations. By using either of these two methods, diffractive elements can be used not only to compensate for distortions such as chromatic or spherical aberration, but also to perform the work of a variety of other optical elements such as null correctors, beam shapers, etc. The main focus of the project described in this report is how the design information from the lens design program is incorporated into the photolithographic process. It is shown that the MANN program, a photolithographic mask generator, fills the need for a link between lens design programs and mask generation controllers.The generated masks can be used to expose a resist-coated substrate which is etched and then must be re-coated, re-exposed, and re-etched for making copies, just as in the electronics industry.

Matthys, Donald R.↗

A Match Made in Space

Just before the space shuttle reaches orbit, its three main engines shut down so that it can achieve separation from the massive external tank that provided the fuel required for liftoff and ascent. In jettisoning the external tank, which is completely devoid of fuel at this point in the flight, the space shuttle fires a series of thrusters, separate from its main engines, that gives the orbiter the maneuvering ability necessary to safely steer clear of the descending tank and maintain its intended flight path. These thrusters make up the space shuttle s Reaction Control System. While the space shuttle s main engines only provide thrust in one direction (albeit a very powerful thrust), the Reaction Control System engines allow the vehicle to maneuver in any desired direction (via small amounts of thrust). The resulting rotational maneuvers are known as pitch, roll, and yaw, and are very important in ensuring that the shuttle docks properly when it arrives at the International Space Station and safely reenters the Earth s atmosphere upon leaving. To prevent the highly complex Reaction Control System from malfunctioning during space shuttle flights, and to provide a diagnosis if such a mishap were to occur, NASA turned to a method of artificial intelligence that truly defied the traditional laws of computer science.

Source record↗

Measurement of Residual Flexibility for Substructures Having Prominent Flexible Interfaces

Verification of a dynamic model of a constrained structure requires a modal survey test of the physical structure and subsequent modification of the model to obtain the best agreement possible with test data. Constrained-boundary or fixed-base testing has historically been the most common approach for verifying constrained mathematical models, since the boundary conditions of the test article are designed to match the actual constraints in service. However, there are difficulties involved with fixed-base testing, in some cases making the approach impractical. It is not possible to conduct a truly fixed-base test due to coupling between the test article and the fixture. In addition, it is often difficult to accurately simulate the actual boundary constraints, and the cost of designing and constructing the fixture may be prohibitive. For use when fixed-base testing proves impractical or undesirable, alternate free-boundary test methods have been investigated, including the residual flexibility technique. The residual flexibility approach has been treated analytically in considerable detail and has had limited frequency response measurements for the method. This concern is well-justified for a number of reasons. First, residual flexibilities are very small numbers, typically on the order of 1.0E-6 in/lb for translational diagonal terms, and orders of magnitude smaller for off-diagonal values. This poses difficulty in obtaining accurate and noise-free measurements, especially for points removed from the excitation source. A second difficulty encountered in residual measurements lies in obtaining a clean residual function in the process of subtracting synthesized modal data from a measured response function. Inaccuracies occur since modes are not subtracted exactly, but only to the accuracy of the curve fits for each mode; these errors are compounded with increasing distance from the excitation point. In this paper, the residual flexibility method is applied to a simple structure in both test and analysis. Measured and predicted residual functions are compared, and regions of poor data in the measured curves are described. It is found that for accurate residual measurements, frequency response functions having prominent stiffness lines in the acceleration/force format are needed. The lack of such stiffness lines increases measurement errors. Interface drive point frequency respose functions for shuttle orbiter payloads exhibit dominant stiffness lines, making the residual test approach a good candidate for payload modal tests when constrained tests are inappropriate. Difficulties in extracting a residual flexibility value from noisy test data are discussed. It is shown that use of a weighted second order least-squares curve fit of the measured residual function allows identification of residual flexibility that compares very well with predictions for the simple structure. This approach also provides an estimate of second order residual mass effects.

Tinker, Michael L.↗

Galaxy assembly bias and large-scale distribution: a comparison between IllustrisTNG and a semi-analytic model

In this work, we compare large scale structure observables for stellar mass selected samples at z = 0, as predicted by two galaxy models, the hydrodynamical simulation IllustrisTNG and the Santa-Cruz semi-analytic model (SC-SAM). Although both models have been independently calibrated to match observations, rather than each other, we find good agreement between the two models for two-point clustering and galaxy assembly bias signatures. The models also show a qualitatively similar response of occupancy and clustering to secondary halo parameters other than mass, such as formation history and concentration, although with some quantitative differences. Thus, our results demonstrate that the galaxy–halo relationships in SC-SAM and TNG are quite similar to first order. However, we also find areas in which the models differ. For example, we note a strong correlation between halo gas content and environment in TNG, which is lacking in the SC-SAM, as well as differences in the occupancy predictions for low-mass haloes. Moreover, we show that higher order statistics, such as cumulants of the density field, help us to accurately describe the galaxy distribution and discriminate between models that show degenerate behaviour for two-point statistics. Our results suggest that SAMs are a promising cost-effective and intuitive method for generating mock catalogues for next generation cosmological surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

Nanoflare Heating Frequency of an X-Ray Bright Point Observed By MaGIXS

Nanoflares are thought to be one of the prime candidates that can heat the solar corona to its multimillion kelvin temperature. Individual nanoflares are difficult to detect with the present generation instruments, however their presence can be inferred by comparing simulated nanoflare-heated plasma emissions with the observed emission. Using HYDRAD coronal loop simulations, we model the emission from an X-ray bright point observed by the Marshall Grazing Incidence X-ray Spectrometer (MaGIXS), along with concurrent observations from the Atmospheric Imaging Assembly (AIA) onboard Solar Dynamics Observatory (SDO) and X-Ray Telescope (XRT) onboard Hinode observatory. The length and magnetic field strength of the coronal loops are derived from the potential field extrapolation of the observed photospheric magnetogram by Helioseismic and Magnetic Imager (HMI) onboard SDO. Each loop is assumed to be heated by random nanoflares, whose magnitude and frequency are determined by the loop length and magnetic field strength. The simulation results are then compared and matched against the measured intensity from AIA, XRT, and MaGIXS. Our model results indicate the loop morphology and emissions from the XBP under study could be well matched by a distribution of nanoflares with average delay times 400 s to 800 s, which strongly suggest that the heating is dominated by high-frequency events. Further, we demonstrate the high sensitivity of MaGIXS and XRT to diagnose the heating frequency using this method, while AIA passbands are found to be the least sensitive.

coronal heating↗

Determining The Nanoflare Heating Frequency of an X-Ray Bright Point Observed by MaGIXS

Nanoflares are thought to be one of the prime candidates that can heat the solar corona to its multi-million kelvin temperature. Individual nanoflares are difficult to detect with the present generation instruments, however their presence can be inferred by comparing simulated nanoflare-heated plasma emissions with the observed emission. Using HYDRAD coronal loop simulations, we model the emission from an X-ray bright point (XBP) observed by the Marshall Grazing Incidence X-ray Spectrometer (MaGIXS), along with nearest-available observations from the Atmospheric Imaging Assembly (AIA) onboard Solar Dynamics Observatory (SDO) and X-Ray Telescope (XRT) onboard Hinode observatory. The length and magnetic field strength of the coronal loops are derived from the linear-force-free extrapolation of the observed photospheric magnetogram by Helioseismic and Magnetic Imager (HMI) onboard SDO. Each loop is assumed to be heated by random nanoflares, whose magnitude and frequency are determined by the loop length and magnetic field strength. The simulation results are then compared and matched against the measured intensity from AIA, XRT, and MaGIXS. Our model results indicate the observed emissions from the XBP under study could be well matched by a distribution of nanoflares with average delay times 1500 s to 3000 s, which suggest that the heating is dominated by high-frequency events. Further, we demonstrate the high sensitivity of MaGIXS and XRT to diagnose the heating frequency using this method, while AIA passbands are found to be the least sensitive.

coronal heating↗

Models of bars − II. Exponential profiles

ABSTRACT We present a new model for galactic bars with exponentially falling major axis luminosity profiles and Gaussian cross-sections. This is based on the linear superposition of Gaussian potential–density pairs with an exponential weight function, using an extension of the method originally introduced by Long & Murali. We compute the density, potential, and forces, using Gaussian quadrature. These quantities are given as explicit functions of position. There are three independent scaled bar parameters that can be varied continuously to produce bespoke bars of a given mass and shape. We categorize the effective potential by splitting a reduced parameter space into six regions. Unusually, we find bars with three stable Lagrange points on the major axis are possible. Our model reveals a variety of unexpected orbital structure, including a bifurcating x1 orbit coexisting with a stable x4 orbit. Propeller orbits are found to play a dominant role in the orbital structure, and we find striking similarities between our bar configuration and the model of Kaufmann & Contopoulos. We find a candidate orbital family, sired from the propeller orbits, that may be responsible for the observed high-velocity peaks in the Milky Way’s bar. As a cross-check, we inspect, for the first time, the proper motions of stars in the high-velocity peaks, which also match our suggested orbital family well. This work adds to the increasing body of evidence that real galactic bars may be supported at least partly by propeller orbits rather than solely by elliptical-like orbits of the x1 family.

McGough, D. P.↗

The Spatial Vision Tree: A Generic Pattern Recognition Engine- Scientific Foundations, Design Principles, and Preliminary Tree Design

New foundational ideas are used to define a novel approach to generic visual pattern recognition. These ideas proceed from the starting point of the intrinsic equivalence of noise reduction and pattern recognition when noise reduction is taken to its theoretical limit of explicit matched filtering. This led us to think of the logical extension of sparse coding using basis function transforms for both de-noising and pattern recognition to the full pattern specificity of a lexicon of matched filter pattern templates. A key hypothesis is that such a lexicon can be constructed and is, in fact, a generic visual alphabet of spatial vision. Hence it provides a tractable solution for the design of a generic pattern recognition engine. Here we present the key scientific ideas, the basic design principles which emerge from these ideas, and a preliminary design of the Spatial Vision Tree (SVT). The latter is based upon a cryptographic approach whereby we measure a large aggregate estimate of the frequency of occurrence (FOO) for each pattern. These distributions are employed together with Hamming distance criteria to design a two-tier tree. Then using information theory, these same FOO distributions are used to define a precise method for pattern representation. Finally the experimental performance of the preliminary SVT on computer generated test images and complex natural images is assessed.

Rahman, Zia-ur↗

Development and Applications of a Stage Stacking Procedure

The preliminary design of multistage axial compressors in gas turbine engines is typically accomplished with mean-line methods. These methods, which rely on empirical correlations, estimate compressor performance well near the design point, but may become less reliable off-design. For land-based applications of gas turbine engines, off-design performance estimates are becoming increasingly important, as turbine plant operators desire peaking or load-following capabilities and hot-day operability. The current work develops a one-dimensional stage stacking procedure, including a newly defined blockage term, which is used to estimate the off-design performance and operability range of a 13-stage axial compressor used in a power generating gas turbine engine. The new blockage term is defined to give mathematical closure on static pressure, and values of blockage are shown to collapse to curves as a function of stage inlet flow coefficient and corrected shaft speed. In addition to these blockage curves, the stage stacking procedure utilizes stage characteristics of ideal work coefficient and adiabatic efficiency. These curves are constructed using flow information extracted from computational fluid dynamics (CFD) simulations of groups of stages within the compressor. Performance estimates resulting from the stage stacking procedure are shown to match the results of CFD simulations of the entire compressor to within 1.6% in overall total pressure ratio and within 0.3 points in overall adiabatic efficiency. Utility of the stage stacking procedure is demonstrated by estimation of the minimum corrected speed which allows stable operation of the compressor. Further utility of the stage stacking procedure is demonstrated with a bleed sensitivity study, which estimates a bleed schedule to expand the compressors operating range.

Kulkarni, Sameer↗

Inner Front Pole Cover Erosion in the 12.5 kW HERMeS Hall Thruster Over a Range of Operating Conditions

Erosion characteristics on the cover of the inner front pole in a 12.5 kW Hall thruster were measured over a wide range of operating conditions in tests of 6 to 14 hours duration using an accelerated test method and a very sensitive radioactive tracer-based erosion diagnostic. The operating points included the nominal 300 - 600 V conditions on a constant 20.8 A throttle curve, but included additional conditions at other currents spanning the throttling envelope and measurements at varying magnetic eld strength, facility pressure, and discharge voltage oscillation amplitude. The results show that the 300 V condition produces the highest wear rates on the 20.8 A throttle curve, but that rates actually increase with decreasing current. The wear rate was insensitive to discharge voltage ripple, but increased monotonically with magnetic eld strength, particularly near the inner radius of the pole cover. The inner region was also sensitive to facility pressure, showing lower rates at a higher pressure level. Separate experiments in which the energy distributions of ions generated by the hollow cathode were measured suggest that the cathode plume may be a source of energetic ions responsible for some of the erosion trends, in addition to ions originating in the thruster plume. The Hall thruster simulation code Hall2De is able to reproduce the erosion characteristics observed at 600 V, 20.8 A, but cannot currently match the rates measured at lower voltages and currents.

Lopez Ortega, Alejandro↗

Handheld White Light Interferometer for Measuring Defect Depth in Windows

Accurate quantification of defects (scratches and impacts) is vital to the certification of flight hardware and other critical components. The amount of damage to a particular component contributes to the performance, reliability, and safety of a system, which ultimately affects the success or failure of a mission or test. The launch-commit criteria on a Space Shuttle Orbiter window are governed by the depth of the defects that are identified by a visual inspection. This measurement of a defect is not easy to obtain given the environment, size of the defect, and location of the window(s). The determination of depth has typically been performed by taking a mold impression and measuring the impression with an optical profiling instrument. Another method of obtaining an estimate of the depth is by using a refocus microscope. To use a refocus microscope, the surface of the glass and bottom of the defect are, in turn, brought into focus by the operator. The amount of movement between the two points corresponds to the depth of the defect. The refocus microscope requires a skilled operator and has been proven to be unreliable when used on Orbiter windows. White light interferometry was chosen as a candidate to replace the refocus microscope. The White Light Interferometer (WLI) was developed to replace the refocus microscope as the instrument used for measuring the depth of defects in Orbiter windows. The WLI consists of a broadband illumination source, interferometer, detector, motion control, displacement sensor, mechanical housing, and support electronics. The illumination source for the WLI is typically a visible light emitting diode (LED) or a near-infrared superluminescent diode (SLD) with power levels of less than a milliwatt. The interferometer is a Michelson configuration consisting of a 1-in. (2.5-cm) cube beam splitter, a 0.5-in. (1.3-cm) optical window as a movable leg (used to closely match the return intensity of the fixed leg from the window), and a mirrored prism to fold the optics into the mechanical housing. The detector may be one of many C-mount CCD (charge-coupled device) cameras. Motion is provided by a commercial nanostepping motor with a serial interface. The displacement sensor is a custom device specifically designed for this application. The mechanical housing and support electronics were designed to integrate the various components into an instrument that could be physically handled by a technician and easily transported.

Youngquist, Robert↗

TPSAS-NF1676L-16988-DND

There has been renewed interest to uniformly recalibrate historical geostationary (GEO) data records to aid in climate monitoring. GEO sensors have annual repeatable angular sampling over a given location. The view angle is fixed and the imaging schedule is usually constant through its lifetime. Given the fact that colocated GEOs are always share the same sub-satellite point and maintain their imaging schedules provides repeatable angular sampling over decades. One of the biggest challenges in transferring a reference sensor calibration using invariant desert targets to another sensor is the accuracy of the bidirectional reflectance distribution function (BDRF). A well-calibrated GEO can be used to predict the daily exoatmospheric radiance model (DERM) over a desert target for a given GMT that is valid for any GEO sensor at the same location. The advantage of this method is that a BRDF is not needed. Another challenge of invariant target calibration is the unique spectra signature of the desert. However, since most GEOs are built in batches, the spectral response functions (SRF) are very similar for most historical GEOs, the spectral band adjustment factor (SBAF) between GEO sensors is much smaller than for MODIS and GEO sensors. Since the water vapor burden over the desert is seasonal, both the TOA and desert surface can be considered invariant for a given day of the year. The reference GEO can be inter-calibrated with MODIS or VIIRS, which have onboard visible calibration using solar diffusers, using other methods, such as ray-matching or deep convective clouds. Also the next generation GEOs will have onboard visible calibration, which will increase the accuracy of this method. Three Meteosats over the Libyan desert will be used to illustrate the DERM method. The reference Meteosat will be inter-calibrated against Aqua-MODIS. The reference GEO DERM will be constructed and used to calibrate the remaining Meteosats. The DERM calibration will be validated by comparing the calibration using Aqua-MODIS ray-matching. Similarly, two GOES sensors using the Sonoran desert will also be highlighted. An uncertainty analysis will also be performed with emphasis on the SBAF, derived over the desert targets using both SCIAMACHY and Hyperion hyper-spectral radiances.

David Doelling↗

Non-Contact Thrust Stand Calibration Method for Repetitively-Pulsed Electric Thrusters

A thrust stand calibration technique for use in testing repetitively-pulsed electric thrusters for in-space propulsion has been developed and tested using a modified hanging pendulum thrust stand. In the implementation of this technique, current pulses are applied to a solenoidal coil to produce a pulsed magnetic field that acts against the magnetic field produced by a permanent magnet mounted to the thrust stand pendulum arm. The force on the magnet is applied in this non-contact manner, with the entire pulsed force transferred to the pendulum arm through a piezoelectric force transducer to provide a time-accurate force measurement. Modeling of the pendulum arm dynamics reveals that after an initial transient in thrust stand motion the quasisteady average deflection of the thrust stand arm away from the unforced or zero position can be related to the average applied force through a simple linear Hooke s law relationship. Modeling demonstrates that this technique is universally applicable except when the pulsing period is increased to the point where it approaches the period of natural thrust stand motion. Calibration data were obtained using a modified hanging pendulum thrust stand previously used for steady-state thrust measurements. Data were obtained for varying impulse bit at constant pulse frequency and for varying pulse frequency. The two data sets exhibit excellent quantitative agreement with each other as the constant relating average deflection and average thrust match within the errors on the linear regression curve fit of the data. Quantitatively, the error on the calibration coefficient is roughly 1% of the coefficient value.

Wong, Andrea R.↗

Status of the NASA GMAO Observing System Simulation Experiment

An Observing System Simulation Experiment (OSSE) is a pure modeling study used when actual observations are too expensive or difficult to obtain. OSSEs are valuable tools for determining the potential impact of new observing systems on numerical weather forecasts and for evaluation of data assimilation systems (DAS). An OSSE has been developed at the NASA Global Modeling and Assimilation Office (GMAO, Errico et al 2013). The GMAO OSSE uses a 13-month integration of the European Centre for Medium- Range Weather Forecasts 2005 operational model at T511/L91 resolution for the Nature Run (NR). Synthetic observations have been updated so that they are based on real observations during the summer of 2013. The emulated observation types include AMSU-A, MHS, IASI, AIRS, and HIRS4 radiance data, GPS-RO, and conventional types including aircraft, rawinsonde, profiler, surface, and satellite winds. The synthetic satellite wind observations are colocated with the NR cloud fields, and the rawinsondes are advected during ascent using the NR wind fields. Data counts for the synthetic observations are matched as closely as possible to real data counts, as shown in Figure 2. Errors are added to the synthetic observations to emulate representativeness and instrument errors. The synthetic errors are calibrated so that the statistics of observation innovation and analysis increments in the OSSE are similar to the same statistics for assimilation of real observations, in an iterative method described by Errico et al (2013). The standard deviations of observation minus forecast (xo-H(xb)) are compared for the OSSE and real data in Figure 3. The synthetic errors include both random, uncorrelated errors, and an additional correlated error component for some observational types. Vertically correlated errors are included for conventional sounding data and GPS-RO, and channel correlated errors are introduced to AIRS and IASI (Figure 4). HIRS, AMSU-A, and MHS have a component of horizontally correlated error. The forecast model used by the GMAO OSSE is the Goddard Earth Observing System Model, Version 5 (GEOS-5) with Gridpoint Statistical Interpolation (GSI) DAS. The model version has been updated to v. 5.13.3, corresponding to the current operational model. Forecasts are run on a cube-sphere grid with 180 points along each edge of the cube (approximately 0.5 degree horizontal resolution) with 72 vertical levels. The DAS is cycled at 6-hour intervals, with 240 hour forecasts launched daily at 0000 UTC. Evaluation of the forecasting skill for July and August is currently underway. Prior versions of the GMAO OSSE have been found to have greater forecasting skill than real world forecasts. It is anticipated that similar forecast skill will be found in the updated OSSE.

OSSE↗

Development And Test of A Digitally Steered Antenna Array for The Navigator GPS Receiver

Global Positioning System (GPS)-based navigation has become common for low-Earth orbit spacecraft as the signal environment is similar to that on the Earth s surface. The situation changes abruptly, however, for spacecraft whose orbital altitudes exceed that of the GPS constellation. Visibility is dramatically reduced and signals that are present may be very weak and more susceptible to interference. GPS receivers effective at these altitudes require increased sensitivity, which often requires a high-gain antenna. Pointing such an antenna can pose a challenge. One efficient approach to mitigate these problems is the use of a digitally steered antenna array. Such an antenna can optimally allocate gain toward desired signal sources and away from interferers. This paper presents preliminary results in the development and test of a digitally steered antenna array for the Navigator GPS research program at NASA s Goddard Space Flight Center. In particular, this paper highlights the development of an array and front-end electronics, the development and test of a real-time software GPS receiver, and implementation of three beamforming methods for combining the signals from the array. Additionally, this paper discusses the development of a GPS signal simulator which produces digital samples of the GPS L1C/A signals as they would be received by an arbitrary antenna array configuration. The simulator models transmitter and receiver dynamics, near-far and multipath interference, and has been a critical component in both the development and test of the GPS receiver. The GPS receiver system was tested with real and simulated GPS signals. Preliminary results show that performance improvement was achieved in both the weak signal and interference environments, matching analytical predictions. This paper summarizes our initial findings and discusses the advantages and limitations of the antenna array and the various beamforming methods.

Pinto, Heitor David↗

Identification of major moisture sources across the Mediterranean Basin

We employ a Lagrangian based moisture back trajectory method on an ensemble of four reanalysis datasets to provide a comprehensive understanding of moisture sources over the Mediterranean land region (30° N–49.5° N and 9.75° W–61.5° E) at seasonal timescales for 1980–2013 period. Using a source region between 10° S–71.35° N along the latitude and 80° W–84.88° E along the longitude that is subdivided into ten complimentary sub-regions, our analyses is able to backtrack up to > 90% of seasonal precipitation at each grid point within the target region. Our results indicate a significant role of moisture advected from the North Atlantic and Mediterranean Sea, and locally recycled moisture over the target region in shaping the spatial organization of seasonal precipitation. However, a clear east–west contrast is witnessed in determining the relative importance of each of these major moisture sources where the North Atlantic dictates the moisture supply over the western Mediterranean while moisture from Mediterranean Sea and local recycling play a key role over the eastern Mediterranean. Our analyses also demonstrate a major footprint of the North Atlantic Oscillation (NAO) on precipitation variability over the Mediterranean land as dynamic and thermodynamic anomalies during the negative phase of NAO match with those during wet years and vice versa. The findings reported here are generally consistent across the four reanalysis datasets. Overall, this study establishes the relative roles of adjacent and far-off oceanic and terrestrial evaporative sources over the Mediterranean land and should help in understanding the drivers of precipitation variability and change at varying timescales.

54 ENVIRONMENTAL SCIENCES↗

Scalability of Parallel Spatial Direct Numerical Simulations on Intel Hypercube and IBM SP1 and SP2

The implementation and performance of a parallel spatial direct numerical simulation (PSDNS) approach on the Intel iPSC/860 hypercube and IBM SP1 and SP2 parallel computers is documented. Spatially evolving disturbances associated with the laminar-to-turbulent transition in boundary-layer flows are computed with the PSDNS code. The feasibility of using the PSDNS to perform transition studies on these computers is examined. The results indicate that PSDNS approach can effectively be parallelized on a distributed-memory parallel machine by remapping the distributed data structure during the course of the calculation. Scalability information is provided to estimate computational costs to match the actual costs relative to changes in the number of grid points. By increasing the number of processors, slower than linear speedups are achieved with optimized (machine-dependent library) routines. This slower than linear speedup results because the computational cost is dominated by FFT routine, which yields less than ideal speedups. By using appropriate compile options and optimized library routines on the SP1, the serial code achieves 52-56 M ops on a single node of the SP1 (45 percent of theoretical peak performance). The actual performance of the PSDNS code on the SP1 is evaluated with a "real world" simulation that consists of 1.7 million grid points. One time step of this simulation is calculated on eight nodes of the SP1 in the same time as required by a Cray Y/MP supercomputer. For the same simulation, 32-nodes of the SP1 and SP2 are required to reach the performance of a Cray C-90. A 32 node SP1 (SP2) configuration is 2.9 (4.6) times faster than a Cray Y/MP for this simulation, while the hypercube is roughly 2 times slower than the Y/MP for this application. KEY WORDS: Spatial direct numerical simulations; incompressible viscous flows; spectral methods; finite differences; parallel computing.

Joslin, Ronald D.↗