Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Smooth approximation”

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 595 records · Page 33

Longitudinal-control design approach for high-angle-of-attack aircraft

This paper describes a control synthesis methodology that emphasizes a variable-gain output feedback technique that is applied to the longitudinal channel of a high-angle-of-attack aircraft. The aircraft is a modified F/A-18 aircraft with thrust-vectored controls. The flight regime covers a range up to a Mach number of 0.7; an altitude range from 15,000 to 35,000 ft; and an angle-of-attack (alpha) range up to 70 deg, which is deep into the poststall region. A brief overview is given of the variable-gain mathematical formulation as well as a description of the discrete control structure used for the feedback controller. This paper also presents an approximate design procedure with relationships for the optimal weights for the selected feedback control structure. These weights are selected to meet control design guidelines for high-alpha flight controls. Those guidelines that apply to the longitudinal-control design are also summarized. A unique approach is presented for the feed-forward command generator to obtain smooth transitions between load factor and alpha commands. Finally, representative linear analysis results and nonlinear batch simulation results are provided.

Ostroff, Aaron J.↗

Multicolor eclipse studies of UU Aquarii. 1: Observations and system parameters

A study of the eclipses in UU Aqr from multicolor high-speed photometry is presented. A revised ephemeris for the times of minimum and an upper limit for orbital period variations are obtained. We use measurements of contact phases in the eclipse light curve to derive the binary geometry and to estimate masses and relevant dimensions. We find a mass ratio of q = 0.30 +/- 0.07 and an inclination of i = 78 deg +/- 2 deg. The masses of the component stars are M(sub 1) = 0.67 +/- 0.14 solar mass and M(sub 2) = 0.20 +/- 0.07 solar mass. Our photometric model predicts K(sub 1) = 84 +/- 26 km/s, which is approximately 30% smaller than the velocity amplitude obtained from the emission lines. From the white dwarf fluxes we estimate T(sub wd) approximately = 34,000 K and a distance of d = 270 +/- 50 pc if the inner disk is opaque. UU Aqr has long term brightness variations of approximately = 0.3 m on timescales of approximately 4 yr. The system was in a 'high' state in 1989 and 1990 and in a 'low' state in 1988 and 1992. The high state results from an increase in the brightness of the outer and cooler parts of the disk, mainly due to the appearance of a bright spot at disk rim. Based on the smooth and gradual eclipse shape and on the absence of a prominent hump in the light curve we suggest that UU Aqr is a high mass-transfer nova-like system with a relatively bright and optically thick accretion disk. We find no perceptible eclipse in the H-alpha emission line. The fluxes at mid-eclipse can be fitted by a compostion of a late-type spectrum plus an optically thin hydrogen emission-line spectrum. These evidences suggest that the emission lines are formed in an extended region only partially occulted during eclipse.

Baptista, R.↗

Combined Modeling of Acceleration, Transport, and Hydrodynamic Response in Solar Flares: The Numerical Model - 1

Acceleration and transport of high-energy particles and fluid dynamics of atmospheric plasma are interrelated aspects of solar flares, but for convenience and simplicity they were artificially separated in the past. We present here self consistently combined Fokker-Planck modeling of particles and hydrodynamic simulation of flare plasma. Energetic electrons are modeled with the Stanford unified code of acceleration, transport, and radiation, while plasma is modeled with the Naval Research Laboratory flux tube code. We calculated the collisional heating rate directly from the particle transport code, which is more accurate than those in previous studies based on approximate analytical solutions. We repeated the simulation of Mariska et al. with an injection of power law, downward-beamed electrons using the new heating rate. For this case, a -10% difference was found from their old result. We also used a more realistic spectrum of injected electrons provided by the stochastic acceleration model, which has a smooth transition from a quasi-thermal background at low energies to a non thermal tail at high energies. The inclusion of low-energy electrons results in relatively more heating in the corona (versus chromosphere) and thus a larger downward heat conduction flux. The interplay of electron heating, conduction, and radiative loss leads to stronger chromospheric evaporation than obtained in previous studies, which had a deficit in low-energy electrons due to an arbitrarily assumed low-energy cutoff. The energy and spatial distributions of energetic electrons and bremsstrahlung photons bear signatures of the changing density distribution caused by chromospheric evaporation. In particular, the density jump at the evaporation front gives rise to enhanced emission, which, in principle, can be imaged by X-ray telescopes. This model can be applied to investigate a variety of high-energy processes in solar, space, and astrophysical plasmas.

STOCHASTIC ACCELERATION MODEL↗

High-Performance Semiempirical Excited-State Molecular Dynamics Powered by Graphics Processing Units

Here, this Letter introduces excited-state molecular dynamics in PYSEQM, a GPU-accelerated semiempirical quantum chemistry engine implemented in PyTorch. The new module enables Born–Oppenheimer molecular dynamics (BOMD) using configuration-interaction singles and random phase approximation for excited states, allowing long trajectories and large statistical ensembles to be simulated efficiently on a single GPU. We also implement an extended Lagrangian excited-state BOMD (XL-ESMD) scheme that propagates auxiliary electronic variables, enabling relaxed ground and excited-state convergence thresholds without compromising energy conservation. The excited-state BOMD implementation scales smoothly from small chromophores to a nearly 900-atom dendrimer (taking 6.5 s per MD step). PYSEQM also supports batched execution, allowing many geometries or trajectories to be evaluated in a single GPU launch, substantially increasing throughput and making ensemble-based protocols routine. As a demonstration, we compute absorption, emission, and infrared spectra from trajectories propagated on the ground and first excited states. The XL-ESMD scheme yields identical spectra at significantly lower computational cost, establishing the role of extended Lagrangian based dynamics for efficient excited-state BOMD simulations. Beyond raw performance, PYSEQM’s PyTorch foundation provides automatic differentiation for forces, efficient GPU batching, and seamless interfacing with machine learning models. These capabilities position PYSEQM as a practical platform for machine learning-augmented excited-state dynamics and lay the foundation for future data-driven nonadiabatic excited-state dynamics modeling of ultrafast spectroscopic probes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Carbon and nitrogen abundance determinations from transition layer lines

For red giants a smooth increase in the nitrogen to carbon abundance ratio for increasing B-V as is expected for the first dredge up phase when the outer convection zone deepens is found. An average increase in the nitrogen to silicon ratio for B-V = 0.6 which goes back to almost solar values for cool giants with B - V approximately 1.0 is reported. It looks as if Si would be enriched for deeper mixing contrary to expectations from standard evolution theory.

Boehm-Vitense, Erika↗

Effect of nose shape on three-dimensional stagnation region streamlines and heating rates

A new method for calculating the three-dimensional inviscid surface streamlines and streamline metrics using Cartesian coordinates and time as the independent variable of integration has been developed. The technique calculates the streamline from a specified point on the body to a point near the stagnation point by using a prescribed pressure distribution in the Euler equations. The differential equations, which are singular at the stagnation point, are of the two point boundary value problem type. Laminar heating rates are calculated using the axisymmetric analog concept for three-dimensional boundary layers and approximate solutions to the axisymmetric boundary layer equations. Results for elliptic conic forebody geometries show that location of the point of maximum heating depends on the type of conic in the plane of symmetry and the angle of attack, and that this location is in general different from the stagnation point. The new method was found to give smooth predictions of heat transfer in the nose region where previous methods gave oscillatory results.

Hassan, Basil↗

Rapid Structured Volume Grid Smoothing and Adaption Technique

A rapid, structured volume grid smoothing and adaption technique, based on signal processing methods, was developed and applied to the Shuttle Orbiter at hypervelocity flight conditions in support of the Columbia Accident Investigation. Because of the fast pace of the investigation, computational aerothermodynamicists, applying hypersonic viscous flow solving computational fluid dynamic (CFD) codes, refined and enhanced a grid for an undamaged baseline vehicle to assess a variety of damage scenarios. Of the many methods available to modify a structured grid, most are time-consuming and require significant user interaction. By casting the grid data into different coordinate systems, specifically two computational coordinates with arclength as the third coordinate, signal processing methods are used for filtering the data [Taubin, CG v/29 1995]. Using a reverse transformation, the processed data are used to smooth the Cartesian coordinates of the structured grids. By coupling the signal processing method with existing grid operations within the Volume Grid Manipulator tool, problems related to grid smoothing are solved efficiently and with minimal user interaction. Examples of these smoothing operations are illustrated for reduction in grid stretching and volume grid adaptation. In each of these examples, other techniques existed at the time of the Columbia accident, but the incorporation of signal processing techniques reduced the time to perform the corrections by nearly 60%. This reduction in time to perform the corrections therefore enabled the assessment of approximately twice the number of damage scenarios than previously possible during the allocated investigation time.

Alter, Stephen J.↗

Rapid Structured Volume Grid Smoothing and Adaption Technique

A rapid, structured volume grid smoothing and adaption technique, based on signal processing methods, was developed and applied to the Shuttle Orbiter at hypervelocity flight conditions in support of the Columbia Accident Investigation. Because of the fast pace of the investigation, computational aerothermodynamicists, applying hypersonic viscous flow solving computational fluid dynamic (CFD) codes, refined and enhanced a grid for an undamaged baseline vehicle to assess a variety of damage scenarios. Of the many methods available to modify a structured grid, most are time-consuming and require significant user interaction. By casting the grid data into different coordinate systems, specifically two computational coordinates with arclength as the third coordinate, signal processing methods are used for filtering the data [Taubin, CG v/29 1995]. Using a reverse transformation, the processed data are used to smooth the Cartesian coordinates of the structured grids. By coupling the signal processing method with existing grid operations within the Volume Grid Manipulator tool, problems related to grid smoothing are solved efficiently and with minimal user interaction. Examples of these smoothing operations are illustrated for reductions in grid stretching and volume grid adaptation. In each of these examples, other techniques existed at the time of the Columbia accident, but the incorporation of signal processing techniques reduced the time to perform the corrections by nearly 60%. This reduction in time to perform the corrections therefore enabled the assessment of approximately twice the number of damage scenarios than previously possible during the allocated investigation time.

Alter, Stephen J.↗

Initial On-Orbit Spatial Resolution Characterization of OrbView-3 Panchromatic Images

Characterization was conducted under the Memorandum of Understanding among Orbital Sciences Corp., ORBIMAGE, Inc., and NASA Applied Sciences Directorate. Acquired five OrbView-3 panchromatic images of the permanent Stennis Space Center edge targets painted on a concrete surface. Each image is available at two processing levels: Georaw and Basic. Georaw is an intermediate image in which individual pixels are aligned by a nominal shift in the along-scan direction to adjust for the staggered layout of the panchromatic detectors along the focal plane array. Georaw images are engineering data and are not delivered to customers. The Basic product includes a cubic interpolation to align the pixels better along the focal plane and to correct for sensor artifacts, such as smile and attitude smoothing. This product retains satellite geometry - no rectification is performed. Processing of the characterized images did not include image sharpening, which is applied by default to OrbView-3 image products delivered by ORBIMAGE to customers. Edge responses were extracted from images of tilted edges in two directions: along-scan and cross-scan. Each edge response was approximated with a superposition of three sigmoidal functions through a nonlinear least-squares curve-fitting. Line Spread Functions (LSF) were derived by differentiation of the analytical approximation. Modulation Transfer Functions (MTF) were obtained after applying the discrete Fourier transform to the LSF.

Blonski, Slawomir↗

First NuSTAR Observations of the BL Lac-Type Blazar PKS 2155-304: Constraints on the Jet Content and Distribution of Radiating Particles

We report the first hard X-ray observations with NuSTAR of the BL Lac-type blazar PKS2155-304, augmented with soft X-ray data from XMM-Newton and gamma-ray data from the Fermi Large Area Telescope, obtained in 2013April when the source was in a very low flux state. A joint NuSTAR and XMM spectrum, covering the energy range 0.5-60 keV, is best described by a model consisting of a log-parabola component with curvature Beta = -0.3(+0.2 -0.1) and a (local) photon index 3.04 +/- 0.15 at photon energy of 2 keV, and a hard power-law tail with photon index 2.2 +/- 0.4. The hard X-ray tail can be smoothly joined to the quasi-simultaneous gamma-ray spectrum by a synchrotron self-Compton component produced by an electron distribution with index p 2.2. Assuming that the power-law electron distribution extends down to gamma (sub min) = 1 and that there is one proton per electron, an unrealistically high total jet power of Lp approximately 10 (exp 47) erg s(sub -1) is inferred. This can be reduced by two orders of magnitude either by considering a significant presence of electron-positron pairs with lepton-to-proton ratio n(sub e+e-/n(sub p) approx. 30, or by introducing an additional, low-energy break in the electron energy distribution at the electron Lorentz factor gamma br1 approx. 100. In either case, the jet composition is expected to be strongly matter-dominated

gamma ray↗

Bayesian sparse learning with preconditioned stochastic gradient MCMC and its applications

Deep neural networks have been successfully employed in an extensive variety of research areas, including solving partial differential equations. Despite its significant success, there are some challenges in effectively training DNN, such as avoiding overfitting in over-parameterized DNNs and accelerating the optimization in DNNs with pathological curvature. Here, we propose a Bayesian type sparse deep learning algorithm. The algorithm utilizes a set of spike-and-slab priors for the parameters in the deep neural network. The hierarchical Bayesian mixture will be trained using an adaptive empirical method. That is, one will alternatively sample from the posterior using preconditioned stochastic gradient Langevin Dynamics (PSGLD), and optimize the latent variables via stochastic approximation. The sparsity of the network is achieved while optimizing the hyperparameters with adaptive searching and penalizing. A popular SG-MCMC approach is Stochastic gradient Langevin dynamics (SGLD). However, considering the complex geometry in the model parameter space in nonconvex learning, updating parameters using a universal step size in each component as in SGLD may cause slow mixing. To address this issue, we apply a computationally manageable preconditioner in the updating rule, which provides a step-size parameter to adapt to local geometric properties. Moreover, by smoothly optimizing the hyperparameter in the preconditioning matrix, our proposed algorithm ensures a decreasing bias, which is introduced by ignoring the correction term in the preconditioned SGLD. According to the existing theoretical framework, we show that the proposed algorithm can asymptotically converge to the correct distribution with a controllable bias under mild conditions. Numerical tests are performed on both synthetic regression problems and learning solutions of elliptic PDE, which demonstrate the accuracy and efficiency of the present work.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Aerodynamic Characteristics of a 0.04956-Scale Model of the Convair F-102A Airplane at Transonic Speeds

Tests have been conducted in the Langley 8-foot transonic tunnel on a 0.04956-scale model of the Convair F-102A airplane which employed an indented and extended fuselage, cambered wing leading edges, and deflected wing tips. Force and moment characteristics were obtained for Mach numbers from 0.60 to 1.135 at angles of attack up to 20 . In addition, tests were made over a limited angle-of-attack range to determine the effects of the cambered leading edges, deflected tips, and a nose section with a smooth area distribution. Fuselage modifications employed on the F-102A were responsible for a 25.percent reduction in the minimum drag-coefficient rise between the Mach numbers of 0.85 and 1.075 when compared with that for the earlier versions of the F-102. Although the wing modifications increased the F-102A subsonic minimum drag-coefficient level approximately 0.0020, they produced large decreases in drag at lifting conditions over that for the original (plane-wing) F-102. The F-102A had 15 to 25 percent higher maximum lift-drag ratios than did the original F-102. The F-102A had about 15 percent lower maximum lift-drag ratios at Mach numbers below 0.95 and slightly higher maximum lift-drag ratios at supersonic speeds when compared with those ratios for sn earlier modified-wing version of the F-102. Chordwise wing fences which provided suitable longitudinal stability for the original F-102 were not adequate for the cambered-wing F-102A The pitching-moment curves indicated a region of near neutral stability with possible pitch-up tendencies for the F-102A at high subsonic Mach numbers for lift coefficients between about 0.4 and 0.5.

Tempelmeyer, Kenneth E.↗

Microstructure and Thermal History of Metal Particles in CH Chondrites

Fe-Ni metal particles with smooth Ni, Co, and Cr zoning patterns, 8-13 wt.% Ni in the center of the particle to 3-5 wt% Ni at the rim, have been identified in several CR-clan (CH, Bencubbinlike, and CR) chondrites. These zoning patterns are broadly consistent with an origin by gas-solid condensation in the solar nebula at temperatures between approximately 1500 to 1300 K and fast cooling rates, 2 to 25 K/day. Apparently, this condensate metal was not melted during chondrule formation or affected significantly in the solid-state by alteration during parent body processing. Consideration of diffusional redistribution of Ni, Co, Cr and siderophile elements have further constrained the calculated condensation temperatures and cooling rates of the zoned condensates. These condensate metals have irregular shapes and vary in size from 50 to 350 m as revealed in some detail by optical and SEM techniques. In addition to zoned condensate particles, other types of metal particles have been observed. These include zoned condensates with exsolution-precipitates, unzoned homogeneous metal with no exsolution precipitates, unzoned metal exhibiting exsolution precipitates and high Ni metal grains.

Goldstein, J. I.↗

A multi-center quadrature scheme for the molecular continuum

We report a common way to evaluate electronic integrals for polyatomic molecules is to use Becke’s partitioning scheme (Becke and Chem, 1988) in conjunction with overlapping grids centered at each atomic site. The Becke scheme was designed for integrands that fall off rapidly at large distances, such as those approximating bound electronic states. When applied to states in the electronic continuum, however, Becke scheme exhibits slow convergence and it is highly redundant. Here, we present a modified version of Becke scheme that is applicable to functions of the electronic continuum, such as those involved in molecular photoionization and electron–molecule scattering, and which ensures convergence and efficiency comparable to those realized in the calculation of bound states. In this modified scheme, the atomic weights already present in Becke’s partition are smoothly switched off within a range of few bond lengths from their respective nuclei, and complemented by an asymptotically unitary weight. The atomic integrals are evaluated on small spherical grids, centered on each atom, with size commensurate to the support of the corresponding atomic weight. The residual integral of the interstitial and long-range region is evaluated with a central master grid. The accuracy of the method is demonstrated by evaluating integrals involving integrands containing Gaussian Type Orbitals and Yukawa potentials, on the atomic sites, as well as spherical Bessel functions centered on the master grid. These functions are representative of those encountered in realistic electron-scattering and photoionization calculations in polyatomic molecules.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Anomalies in the topology of the temperature fluctuations in the cosmic microwave background: An analysis of the NPIPE and FFP10 data releases

We present a topological analysis of the temperature fluctuation maps from the Planck 2020 Data Release 4 NPIPE dataset and the Planck 2018 Data Release 3 FFP10 dataset. We performed a multiscale analysis in terms of the homology characteristics of the maps, invoking relative homology to account for the analysis in the presence of masks. We performed our analysis for a range of smoothing scales spanning sub- and super-horizon scales corresponding to a full width at half maximum (FWHM) of 5',10',20',40',80',160',320', and 640', and employed simulations based on the standard model for comparison, which assumes the initial fluctuation field to be an isotropic and homogeneous Gaussian random field. Examining the behavior of topological components, represented by the 0D homology group, we find the observations to be approximately 2σ or less deviant from the simulations for all resolutions and scales for the NPIPE dataset. For the FFP10 dataset, we detect a 2.96σ deviation between the observations and simulations at N = 128, FWHM = 80'. For the topological loops, represented by the first homology group, the simulations and observations are consistent within 2σ for most resolutions and scales for both the datasets. However, for the NPIPE dataset, we observe a high deviation between the observation and simulations in the number of loops at FWHM = 320', but at a low dimensionless threshold ν = –2.5. Under a Gaussian assumption, this would amount to a deviation of ~4σ. However, the distribution in this bin is manifestly non-Gaussian and does not obey Poisson statistics either. In the absence of a true theoretical understanding, we simply note that the significance is higher than what may be resolved by 600 simulations, yielding an empirical p-value of at most 0.0016. Specifically in this case, our tests indicate that the numbers arise from a statistically stable regime, despite being based on small numbers. For the FFP10 dataset, the differences are not as strong as for the NPIPE dataset, indicating a 2.77σ deviation at this resolution and threshold. The Euler characteristic, which is the alternating sum of the ranks of relative homology groups, reflects the deviations in the components and loops. To assess the significance of combined levels for a given scale, we employed the empirical and theoretical versions of the χ 2 test as well as the nonparametric Tukey depth test. Although all statistics exhibit a stable distribution, we favor the empirical version of the χ 2 test in the final interpretation, as it indicates the most conservative differences. For the NPIPE dataset, we find that the components and loops differ at more than 95%, but agree within the 99% confidence level with respect to the base model at N = 32, FWHM = 320'. The Euler characteristic at this resolution displays a per mil deviation. In contrast, the FFP10 dataset shows that the observations are consistent with the base model within the 95% confidence level, at this and smaller scales. This is consistent with the observations of the Planck analysis pipeline via Minkowski functionals. For the largest smoothing scale, N = 16, FWHM = 640', both datasets exhibit an anomalous behavior of the loops, where FFP10 data exhibit a deviation that is larger by an order of magnitude than that of the NPIPE dataset. In contrast, the values for the topological components and the Euler characteristic agree between observations and model to within a confidence level of 99%. However, for the largest scales, the statistics are based on low numbers and may have to be regarded with caution. Even though both datasets exhibit mild to significant discrepancies, they also exhibit contrasting behaviors at various instances. Therefore, we do not find it feasible to convincingly accept or reject the null hypothesis. Disregarding the large-scale anomalies that persist at similar scales in WMAP and Planck, observations of the cosmic microwave background are largely consistent with the standard cosmological model within 2σ.

79 ASTRONOMY AND ASTROPHYSICS↗

Chandra Observations of the X-ray Environs of SN 1998 bw/GRB 980425

We report X-ray studies of the environs of SN 1998bw and GRB 980425 using the Chandra X-Ray Observatory 1281 days after the GRB. Eight X-ray point sources were localized, three and five each in the original error boxes--S1 and S2--assigned for variable X-ray counterparts to the GRB by BeppoSAX. The sum of the discrete X-ray sources plus continuous emission in S2 observed by CXO on day 1281 is within a factor of 1.5 of the maximum and the upper limits seen by BeppoSAX. We conclude that S2 is the sum of several variable sources that have not disappeared, and therefore is not associated with the GRB. Within S1, clear evidence is seen for a decline of approximately a factor of 12 between day 200 and day 1281. One of the sources in S1, S1a, is coincident with the well-determined radio location of SN 1998bw, and is certainly the remnant of that explosion. The nature of the other sources is also discussed. Combining our observation of the supernova with others of the GRB afterglow, a smooth X-ray light curve, spanning (approx) 1300 days, is obtained by assuming the burst and supernova were coincident at 35. 6 Mpc. When this X-ray light curve is compared with those of the X-ray 'afterglows' of ordinary GRBs, X-ray Flashes, and ordinary supernovae, evidence emerges for at least two classes of lightcurves, perhaps bounding a continuum.

X-RAY SOURCES↗

Stability and Control Flight Tests of a Vertically Rising Airplane Model Similar to the Lockheed XFV-1 Airplane

This paper presents the results of an investigation of the dynamic stability and controllability of a model which approximately represents the Lockheed XFV-1 airplane to a 1/8 scale. The investigation consisted of hovering flights in still air at a considerable height above the ground, hovering flights very close to the ground, vertical take-offs and landings, flights through the transition range from hovering to normal forward flight, and sideways translational flights. The model could be flown smoothly and easily in hovering flight despite the fact that the uncontrolled pitching and yawing motions were unstable oscillations. There was a noticeable reduction in the controllability of the model when hovered very close to the ground but take-offs could be made easily and landings on a g,ven spot could be made accurately in spite of this adverse ground effect. Flights through the transition range from hovering to normal forward flight could be performed fairly easily. The model seemed to have stability of angle of attack and angle of roll over most of the transition range. The yawing motion was divergent in the very high angle-of-attack range but could be controlled easily. At the lower angles of attack, the model seemed to become stable in yaw. In sideways flight there was an increasingly strong tendency to diverge in roll as the speed was increased and finally, at a speed of about 25 knots (full scale), the model rolled off despite efforts of the pilot to control it.

Kirby, Robert H.↗

The Herschel ATLAS: Evolution of the 250 Micrometer Luminosity Function Out to z = 0.5

We have determined the luminosity function of 250 micrometer-selected galaxies detected in the approximately equal to 14 deg(sup 2) science demonstration region of the Herschel-ATLAS project out to a redshift of z = 0.5. Our findings very clearly show that the luminosity function evolves steadily out to this redshift. By selecting a sub-group of sources within a fixed luminosity interval where incompleteness effects are minimal, we have measured a smooth increase in the comoving 250 micrometer luminosity density out to z = 0.2 where it is 3.6(sup +1.4) (sub -0.9) times higher than the local value.

Dye, S.↗