Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Well Log”

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 217 records · Page 12

Growth and Morphology of Supercritical Fluids, a Fluid Physics Experiment Conducted on Mir, Complete

The Growth and Morphology of Supercritical Fluids (GMSF) is an international experiment facilitated by the NASA Glenn Research Center and under the guidance of U.S. principal investor Professor Hegseth of the University of New Orleans and three French coinvestigators: Daniel Beysens, Yves Garrabos, and Carole Chabot. The GMSF experiments were concluded in early 1999 on the Russian space station Mir. The experiments spanned the three science themes of near-critical phase separation rates, interface dynamics in near-critical boiling, and measurement of the spectrum of density fluctuation length scales very close to the critical point. The fluids used were pure CO2 or SF6. Three of the five thermostats used could adjust the sample volume with the scheduled crew time. Such a volume adjustment enabled variable sample densities around the critical density as well as pressure steps (as distinct from the usual temperature steps) applied to the sample. The French-built ALICE II facility was used for these experiments. It allows tightly thermostated (left photograph) samples (right photograph) to be controlled and viewed/measured. Its diagnostics include interferometry, shadowgraph, high-speed pressure measurements, and microscopy. Data were logged on DAT tapes, and PCMCIA cards and were returned to Earth only after the mission was over. The ground-breaking near critical boiling experiment has yielded the most results with a paper published in Physical Review Letters (ref. 1). The boiling work also received press in Science Magazine (ref. 2). This work showed that, in very compressible near-critical two-phase pure fluids, a vapor bubble was induced to temporarily overheat during a rapid heating of the sample wall. The temperature rise in the vapor was 23-percent higher than the rise in the driving container wall. The effect is due to adiabatic compression of the vapor bubble by the rapid expansion of fluid near the boundary during heatup. Thermal diffusivity is low near the critical point, so getting heat out of the compressed bubble is observably slow. This gives the appearance of a backward heat flow, or heat flow from a cold surface to a warm fluid.

Wilkinson, R. Allen↗

The Kinematics of Turbulent Boundary Layer Structure

The long history of research into the internal structure of turbulent boundary layers has not provided a unified picture of the physics responsible for turbulence production and dissipation. The goals of the present research are to: (1) define the current state of boundary layer structure knowledge; and (2) utilize direct numerical simulation results to help close the unresolved issues identified in part A and to unify the fragmented knowledge of various coherent motions into a consistent kinematic model of boundary layer structure. The results of the current study show that all classes of coherent motion in the low Reynolds number turbulent boundary layer may be related to vortical structures, but that no single form of vortex is representative of the wide variety of vortical structures observed. In particular, ejection and sweep motions, as well as entrainment from the free-streem are shown to have strong spatial and temporal relationships with vortical structures. Disturbances of vortex size, location, and intensity show that quasi-streamwise vortices dominate the buffer region, while transverse vortices and vortical arches dominate the wake region. Both types of vortical structure are common in the log region. The interrelationships between the various structures and the population distributions of vortices are combined into a conceptual kinematic model for the boundary layer. Aspects of vortical structure dynamics are also postulated, based on time-sequence animations of the numerically simulated flow.

Robinson, Stephen Kern↗

SCDU Testbed Automated In-Situ Alignment, Data Acquisition and Analysis

In the course of fulfilling its mandate, the Spectral Calibration Development Unit (SCDU) testbed for SIM-Lite produces copious amounts of raw data. To effectively spend time attempting to understand the science driving the data, the team devised computerized automations to limit the time spent bringing the testbed to a healthy state and commanding it, and instead focus on analyzing the processed results. We developed a multi-layered scripting language that emphasized the scientific experiments we conducted, which drastically shortened our experiment scripts, improved their readability, and all-but-eliminated testbed operator errors. In addition to scientific experiment functions, we also developed a set of automated alignments that bring the testbed up to a well-aligned state with little more than the push of a button. These scripts were written in the scripting language, and in Matlab via an interface library, allowing all members of the team to augment the existing scripting language with complex analysis scripts. To keep track of these results, we created an easily-parseable state log in which we logged both the state of the testbed and relevant metadata. Finally, we designed a distributed processing system that allowed us to farm lengthy analyses to a collection of client computers which reported their results in a central log. Since these logs were parseable, we wrote query scripts that gave us an effortless way to compare results collected under different conditions. This paper serves as a case-study, detailing the motivating requirements for the decisions we made and explaining the implementation process.

Automation↗

Analysis of CO2 Composite Spray Cleaning System Results

In a previous study, the cleaning efficiency of a CO2 composite cleaning system for particulate removal was tested. The study covered particles from spores to fluorescent particles of different sizes as well as a variety of substrate surfaces, including aluminum, titanium, stainless steel, and nitinol. Particles were deposited using aerosol (dry) and droplet (wet) deposition. Results from the previous study show that the CO2 composite spray system is capable of cleaning to sterility for aerosol deposited spores and is capable of cleaning a minimum of a 4-log reduction for droplet deposited spores. This minimum 4-log reduction matches current Planetary Protection dry heat microbial reduction requirements. In this paper we will present new data to further correlate the cleaning efficiency with how contamination was introduced to the surface, the surface roughness, and particle size. Possible causes for such correlations will be discussed.

Chung, Shirley↗

JPL control/structure interaction test bed real-time control computer architecture

The Control/Structure Interaction Program is a technology development program for spacecraft that exhibit interactions between the control system and structural dynamics. The program objectives include development and verification of new design concepts - such as active structure - and new tools - such as combined structure and control optimization algorithm - and their verification in ground and possibly flight test. A focus mission spacecraft was designed based upon a space interferometer and is the basis for design of the ground test article. The ground test bed objectives include verification of the spacecraft design concepts, the active structure elements and certain design tools such as the new combined structures and controls optimization tool. In anticipation of CSI technology flight experiments, the test bed control electronics must emulate the computation capacity and control architectures of space qualifiable systems as well as the command and control networks that will be used to connect investigators with the flight experiment hardware. The Test Bed facility electronics were functionally partitioned into three units: a laboratory data acquisition system for structural parameter identification and performance verification; an experiment supervisory computer to oversee the experiment, monitor the environmental parameters and perform data logging; and a multilevel real-time control computing system. The design of the Test Bed electronics is presented along with hardware and software component descriptions. The system should break new ground in experimental control electronics and is of interest to anyone working in the verification of control concepts for large structures.

Briggs, Hugh C.↗

Retrieval of Normalized Gamma Size Distribution Parameters Using Precipitation Imaging Package (Pip) Snowfall Observations During Ice-Pop 2018

Parameters of the normalized gamma particle size distribution (PSD) have been retrieved from the Precipitation Image Package (PIP) snowfall observations collected during the International Collaborative Experiment - PyeongChang Olympics and Paralympic (ICE-POP 2018). Two of the gamma PSD parameters, the mass weighted particle diameter (Dmass) and the normalized intercept parameter NW, have median values of 1.15-1.31 mm and 2.84-3.04 log(mm-1 m-3), respectively. This range arises from the choice of the relationship between the maximum versus equivalent diameter, Dmx-Deq, and the relationship between the Reynolds and Best numbers, Re-X. Normalization of snow water equivalent rate (SWER) and ice water content (W) by NW reduces the range in NW resulting in well fitted power law relationship, between SWER/NW and Dmass and between W/NW and Dmass. The bulk descriptors of snowfall are calculated from PIP observations and from the gamma PSD with values of the shape parameter (μ) ranging from -2 to 10. NASA’s Global Precipitation Measurement (GPM) mission, which adopted the normalized gamma PSD, assumes μ = 2 and μ = 3 in its two separate algorithms. The mean fractional bias (MFB) of the snowfall parameters changes with μ, where the functional dependence on μ depends on the specific snowfall parameter of interest. The MFB of the total concentration was underestimated by 0.23-0.34 when μ = 2 and by 0.29-0.40 when μ = 3, while the MFB of SWER had a much narrower range (-0.03 to 0.04) for the same μ values.

Snowfall↗

Efficient parallel algorithms for string editing and related problems

The string editing problem for input strings x and y consists of transforming x into y by performing a series of weighted edit operations on x of overall minimum cost. An edit operation on x can be the deletion of a symbol from x, the insertion of a symbol in x or the substitution of a symbol x with another symbol. This problem has a well known O((absolute value of x)(absolute value of y)) time sequential solution (25). The efficient Program Requirements Analysis Methods (PRAM) parallel algorithms for the string editing problem are given. If m = ((absolute value of x),(absolute value of y)) and n = max((absolute value of x),(absolute value of y)), then the CREW bound is O (log m log n) time with O (mn/log m) processors. In all algorithms, space is O (mn).

Apostolico, Alberto↗

HR 1362 - A test case for stellar dynamo theories

An analysis of 11 yr of photometry of HR 1362 = EK Eri has shown a well-established photometric period of 335 days. This confirms an earlier period determination of about 310 days from 3 yr of APT data and is consistent with a new v sin i determination of 2 + or - 2 km/s. Recent Ca II H and K observations show moderately strong H and K emission lines at a flux level of log F prime (K) of about 6.5 ergs sq m/s. If the photometric period is interpreted as the rotation period, this surface flux is more than one order of magnitude larger than what can be expected from empirical activity versus rotation relationships. The spectrum and the U-B and B-V colors are consistent with a single G8 III-IV spectral classification, while the V-R and V-I values indicate a small color excess, typical of chromospherically active stars.

Strassmeier, Klaus G.↗

Nitrogen Oxides and Ozone from B-747 Measurements (NOXAR) During POLINAT-2 and SONEX: Overview and Case Studies on Continental and Marine Convection

In the framework of the project POLINAT 2 (Pollution in the North Atlantic Flight Corridor) we measured NO(x) (NO and NO2) and ozone on 85 flights through the North Atlantic Flight Corridor (NAFC) with a fully automated system permanently installed aboard an in-service Swissair B-747 airliner in the period of August to November 1997. The averaged NO(x) concentrations both in the NAFC and at the U.S. east coast were similar to that measured in autumn 1995 with the same system. The patchy occurrence of NO(x) enhancements up to 3000 pptv over several hundred kilometers (plumes), predominately found over the U.S. east coast lead to a log-normal NO(x) probability density function. In three case studies we examine the origins of such plumes by combining back-trajectories with brightness temperature enhanced (IR) satellite imagery, lightning observations from the U.S. National Lightning Detection Network (NLDN) and the Optical Transient Detector (OTD) satellite. We demonstrate that the location of NO(x) plumes can be well explained with maps of convective influence. We show that the number of lightning flashes in cluster of marine thunderstorms is proportional to the NO(x) concentrations observed several hundred kilometers downwind of the anvil outflows. From the fact that in autumn the NO(x) maximum was found several hundred kilometers off the U.S. east coast, it can be inferred that thunderstorms triggered over the warm Gulf Stream current are major sources for the regional upper tropospheric NO(x) budget in autumn.

Jeker, Dominique P.↗

A Possible Explanation for the Z-R Parameter Inconsistency when Comparing Stratiform and Convective Rainfall

The well-known Z-R power law Z A R(sup b) uses two parameters, A and b, in order to relate rainfall rate R to measured weather radar reflectivity Z. A common method used by researchers is to compute Z and R from disdrometer data and then extract the A-b parameter pair from a log-linear line fit to a scatter plot of Z-R pairs. Even though it may seem far more truthful to extract the parameter pair from a fit of radar Z(sub R) versus gauge rainfall rate R(sub G), the extreme difference in spatial and temporal sampling volumes between radar and rain gauge creates a slew of problems that can generally only be solved by using rain gauge arrays and long sampling averages. Disdrometer derived A-b parameters are easily obtained and can provide information for the study of stratiform versus convective rainfall. However, an inconsistency appears when comparing averaged A-b pairs from various researchers. Values of b range from 1.26 to 1.51 for both stratiform and convective events. Paradoxically the values of A fall into three groups: 150 to 200 for convective; 200 to 400 for stratiform; and 400 to 500 again for convective. This apparent inconsistency can be explained by computing the A-b pair using the gamma DSD coupled with a modified drop terminal velocity model, v(D) alpha D(sup beta) - w, where w is a somewhat artificial constant vertical velocity of the air above the disdrometer. This model predicts three regions of A, corresponding to w less than 0, w = 0, and w greater than 0, which approximately matches observed data.

gamma DSD↗

Interpolation and Extrapolation of Creep Rupture Data by the Minimum Commitment Method. Part 3: Analysis of Multiheats

The Minimum Commitment Method was applied to two sets of data for which multiple heat information was available. For one alloy, a 304 stainless steel studied in Japan, data on nine well characterized heats were used, while for a proprietary low alloy carbon steel studied in the United Kingdom data were available on seven heats - in many cases to very long rupture times. For this preliminary study no instability factors were used. It was discovered that heat-to-heat variations would be accounted for by introducing heat identifiers in the form A + B log sigma where sigma is the stress and the constants A and B depend only on the heat. With these identifiers all the data could be collapsed onto a single master curve, even though there was considerable scatter among heats. Using these identifiers together with the average behavior of all heats made possible the determination of an accurate constitutive equation for each individual heat. Two basic approaches are discussed for applying the results of the analysis.

Manson, S. S.↗

Modification of Eccentric Gaze-Holding

Clear vision and accurate localization of objects in the environment are prerequisites for reliable performance of motor tasks. Space flight confronts the crewmember with a stimulus rearrangement that requires adaptation to function effectively with the new requirements of altered spatial orientation and motor coordination. Adaptation and motor learning driven by the effects of cerebellar disorders may share some of the same demands that face our astronauts. One measure of spatial localization shared by the astronauts and those suffering from cerebellar disorders that is easily quantified, and for which a neurobiological substrate has been identified, is the control of the angle of gaze (the "line of sight"). The disturbances of gaze control that have been documented to occur in astronauts and cosmonauts, both in-flight and postflight, can be directly related to changes in the extrinsic gravitational environment and intrinsic proprioceptive mechanisms thus, lending themselves to description by simple non-linear statistical models. Because of the necessity of developing robust normal response populations and normative populations against which abnormal responses can be evaluated, the basic models can be formulated using normal, non-astronaut test subjects and subsequently extended using centrifugation techniques to alter the gravitational and proprioceptive environment of these subjects. Further tests and extensions of the models can be made by studying abnormalities of gaze control in patients with cerebellar disease. A series of investigations were conducted in which a total of 62 subjects were tested to: (1) Define eccentric gaze-holding parameters in a normative population, and (2) explore the effects of linear acceleration on gaze-holding parameters. For these studies gaze-holding was evaluated with the subjects seated upright (the normative values), rolled 45 degrees to both the left and right, or pitched back 30 and 90 degrees. In a separate study the further effects of acceleration on gaze stability was examined during centrifugation (+2 G (sub x) and +2 G (sub z) using a total of 23 subjects. In all of our investigations eccentric gaze-holding was established by having the subjects acquire an eccentric target (+/-30 degrees horizontal, +/- 15 degrees vertical) that was flashed for 750 msec in an otherwise dark room. Subjects were instructed to hold gaze on the remembered position of the flashed target for 20 sec. Immediately following the 20 sec period, subjects were cued to return to the remembered center position and to hold gaze there for an additional 20 sec. Following this 20 sec period the center target was briefly flashed and the subject made any corrective eye movement back to the true center position. Conventionally, the ability to hold eccentric gaze is estimated by fitting the natural log of centripetal eye drifts by linear regression and calculating the time constant (G) of these slow phases of "gaze-evoked nystagmus". However, because our normative subjects sometimes showed essentially no drift (tau (sub c) = m), statistical estimation and inference on the effect of target direction was performed on values of the decay constant theta = 1/(tau (sub c)) which we found was well modeled by a gamma distribution. Subjects showed substantial variance of their eye drifts, which were centrifugal in approximately 20 % of cases, and > 40% for down gaze. Using the ensuing estimated gamma distributions, we were able to conclude that rightward and leftward gaze holding were not significantly different, but that upward gaze holding was significantly worse than downward (p<0.05). We also concluded that vertical gaze holding was significantly worse than horizontal (p<0.05). In the case of left and right roll, we found that both had a similar improvement to horizontal gaze holding (p<0.05), but didn't have a significant effect on vertical gaze holding. For pitch tilts, both tilt angles significantly decreased gaze-holding ility in all directions (p<0.05). Finally, we found that hyper-g centrifugation significantly decreased gaze holding ability in the vertical plane. The main findings of this study are as follows: (1) vertical gaze-holding is less stable than horizontal, (2) gaze-holding to upward targets is less stable than to downward targets, (3) tilt affects gaze holding, and (4) hyper-g affects gaze holding. This difference between horizontal and vertical gaze-holding may be ascribed to separate components of the velocity-to-position neural integrator for eye movements, and to differences in orbital mechanics. The differences between upward and downward gaze-holding may be ascribed to an inherent vertical imbalance in the vestibular system. Because whole body tilt and hyper-g affects gaze-holding, it is implied that the otolith organs have direct connections to the neural integrator and further studies of astronaut gaze-holding are warranted. Our statistical method for representing the range of normal eccentric gaze stability can be readily applied to normals who maybe exposed to environments which may modify the central integrator and require monitoring, and to evaluate patients with gaze-evoked nystagmus by comparing to the above established normative criteria.

Reschke, M. F.↗

SIAM-SERVIR: An Environmental Monitoring and Decision Support System for Mesoamerica

In 2002/2003 NASA, the World Bank and the United States Agency for International Development (USAID) joined with the Central American Commission for Environment and Development (CCAD) to develop an advanced decision support system for Mesoamerica (named SERVIR) as part of the Mesoamerican Environmental Information System (SIAM). Mesoamerica, composed of the seven Central American countries and the five southernmost states of Mexico, make up only a small fraction of the world's land surface. However, the region is home to seven to eight percent of the planet's biodiversity (14 biosphere reserves, 31 Ramsar sites, 8 world heritage sites, 589 protected areas) and 45 million people including more than 50 different ethnic groups. Today Mesoamerica's biological and cultural diversity is severely threatened by extensive deforestation, illegal logging, water pollution, and uncontrolled slash and burn agriculture. Additionally, Mesoamerica's distinct geology and geography result in disproportionate vulnerability to natural disasters such as earthquakes, hurricanes, drought, and volcanic eruptions. NASA Marshall Space Flight Center, together with the University of Alabama in Huntsville (UAH) and the SIAM-SERVIR partners are developing state-of-the-art decision support tools for environmental monitoring as well as disaster prevention and mitigation in Mesoamerica. These partners are contributing expertise in space-based observation with information management technologies and intimate knowledge of local ecosystems to create a system that is being used by scientists, educators, and policy makers to monitor and forecast ecological changes, respond to natural disasters and better understand both natural and human induced effects. In its first year of development and operation, the SIAM-SERVIR project has already yielded valuable information on Central American fires, weather conditions, and the first ever real-time data on red tides. This paper presents the progress thus far in the development of SIAM-SERVIR and the plans for the future.

Irwin, D. E.↗

SIAM-SERVIR: An Environmental Monitoring and Decision Support System for Mesoamerica

In 2002/2003 NASA, the World Bank and the United States Agency for International Development (USAID) joined with the Central American Commission for Environment and Development (CCAD) to develop an advanced decision support system for Mesoamerica (named SERVIR) as part of the Mesoamerican Environmental Information System (SIAM). Mesoamerica - composed of the seven Central American countries and the five southernmost states of Mexico - make up only a small fraction of the world s land surface. However, the region is home to seven to eight percent of the planet s biodiversity (14 biosphere reserves, 31 Ramsar sites, 8 world heritage sites, 589 protected areas) and 45 million people including more than 50 different ethnic groups. Today Mesoamerica s biological and cultural diversity is severely threatened by extensive deforestation, illegal logging, water pollution, and uncontrolled slash and burn agriculture. Additionally, Mesoamerica's distinct geology and geography result in disproportionate vulnerability to natural disasters such as earthquakes, hurricanes, drought, and volcanic eruptions. NASA Marshall Space Flight Center, together with the University of Alabama in Huntsville (UAH) and the SIAM-SERVIR partners are developing state-of-the-art decision support tools for environmental monitoring as well as disaster prevention and mitigation in Mesoamerica. These partners are contributing expertise in space-based observation with information management technologies and intimate knowledge of local ecosystems to create a system that is being used by scientists, educators, and policy makers to monitor and forecast ecological changes, respond to natural disasters and better understand both natural and human induced effects. In its first year of development and operation, the SIAM-SERVIR project has already yielded valuable information on Central American fires, weather conditions, and the first ever real-time data on red tides. This paper presents the progress thus far in the development of SIAM-SERVIR and the plans for the future.

Irwin, Daniel E.↗

An improved logarithmic amplifier circuit for PDS microdensitometers

A high speed logarithmic amplifier circuit for the Perkin-Elmer PDS microdensitometer is discussed. The circuit is designed around an Analog Devices 757p log-amp module. The amplifier produces undistorted profiles of rapidly changing, dense ( 4) images over the available range of scanning speeds. The circuit board was designed to replace directly the manufacturer supplied unit; neither electrical nor mechanical modifications of the basic PDS are required. The performance of the circuit is illustrated through its effects on the overall modulation transfer function of the instrument and by scans of a well exposed stellar image. Circuit diagrams and parts lists a mechanism for obtaining either circuit board layout art work or finished circuit boards are presented.

Anderson, C. M.↗

Nitrogen Oxides and Ozones from B-747 Measurements (NOXAR) during POLINAT 2 and SONEX: Overview and Case-Studies on Continental and Marine Convection

In the framework of the project POLINAT 2 (Pollution in the North Atlantic Flight Corridor) we measured NO(x) (NO and NO2) and ozone on 98 flights through the North Atlantic Flight Corridor (NAFC) with a fully automated system permanently installed aboard an in-service Swissair B-747 airliner in the period of August to November 1997. The averaged NO, concentrations both in the NAFC and at the U.S. east coast were similar to that measured in autumn 1995 with the same system. The patchy occurrence of NO(x), enhancements up to 3000 pptv over several hundred kilometers (plumes), predominately found over the U.S. east coast lead to a log-normal NO(x) probability density function. In three case-studies we examine the origins of such plumes by combining back-trajectories with brightness temperature enhanced (IR) satellite imagery, with lightning observations from the U.S. National Lightning Detection Network (NLDN) or with the Optical Transient Detector (OTD) satellite. For frontal activity above the continental U.S., we demonstrate that the location of NO(x) plumes can be well explained with maps of convective influence. For another case we show that the number of lightning flashes in a cluster of marine thunderstorms is proportional to the NO(x) concentrations observed several hundred kilometers downwind of the anvil outflows and suggest that lightning was the dominant source. From the fact that in autumn the NO, maximum was found several hundred kilometers off the U.S. east coast, it can be inferred that thunderstorms triggered over the warm Gulf Stream current are an important source for the regional upper tropospheric NO(x) budget in autumn.

Jeker, Dominique↗

Gamma-Ray Burst Class Properties

Guided by the Supervised pattern recognition algorithm C4.5, we examine the three gamma-ray burst classes identified by Mukherjee et al. C4.5 provides strong statistical support for this classification. However, with C4.5 and our knowledge of the BATSE instrument, we demonstrate that Class 3 (intermediate fluence, intermediate duration, soft) does not have to be a distinct source population: statistical/systematic errors in measuring burst attributes combined with the well-known hardness/intensity correlation can cause low peak flux Class I (high fluence, long, intermediate hardness) bursts to take on Class 3 characteristics naturally. Based on our hypothesis that the third class is not a distinct one, we provide rules so that future events can be placed in either Class I or Class 2 (low fluence, short, hard). Using classified bursts from the BATSE 4B Catalog, we plot log(N>P) vs. log(P) curves and study spectral features of each class. We find that the two classes are relatively distinct on the basis of spectral parameters, alpha, Beta, and E(sub peak) alone. Although this does not indicate a better basis for classification, it does suggest that different physical conditions exist for Class I and Class 2 bursts. In the process of studying burst class characteristics, we identify a new bias that affects measurement of burst fluences and durations. Using a simple model of how burst duration can be underestimated, we generally characterize how this fluence duration bias affects BATSE measurements, and demonstrate the type of effect it can have on the BATSE fluence vs. peak flux diagram.

Hakkila, Jon↗

The Herschel* PEP-HERMES Luminosity Function- I. Probing the Evolution of PACS Selected Galaxies to z approx. equal to 4

We exploit the deep and extended far-IR data sets (at 70, 100 and 160 μm) of the Herschel Guaranteed Time Observation (GTO) PACS Evolutionary Probe (PEP) Survey, in combination with the Herschel Multi-tiered Extragalactic Survey data at 250, 350 and 500 μm, to derive the evolution of the rest-frame 35-, 60-, 90- and total infrared (IR) luminosity functions (LFs) up to z ∼ 4.We detect very strong luminosity evolution for the total IR LF (LIR ∝ (1 + z)(sup 3.55 +/- 0.10) up to z ∼ 2, and ∝ (1 + z)(sup 1.62 +/- 0.51) at 2 less than z less than approximately 4) combined with a density evolution (∝ (1 + z)(sup −0.57 +/- 0.22) up to z ∼ 1 and ∝ (1 + z)(sup −3.92 +/- 0.34) at 1 less than z less than approximately 4). In agreement with previous findings, the IR luminosity density (ρIR) increases steeply to z ∼ 1, then flattens between z ∼ 1 and z ∼ 3 to decrease at z greater than approximately 3. Galaxies with different spectral energy distributions, masses and specific star formation rates (SFRs) evolve in very different ways and this large and deep statistical sample is the first one allowing us to separately study the different evolutionary behaviours of the individual IR populations contributing to ρIR. Galaxies occupying the well-established SFR-stellar mass main sequence (MS) are found to dominate both the total IR LF and ρIR at all redshifts, with the contribution from off-MS sources (≥0.6 dex above MS) being nearly constant (∼20 per cent of the total ρIR) and showing no significant signs of increase with increasing z over the whole 0.8 < z <2.2 range. Sources with mass in the range 10 ≤ log(M/solar mass) ≤ 11 are found to dominate the total IR LF, with more massive galaxies prevailing at the bright end of the high-z (greater than approximately 2) LF. A two-fold evolutionary scheme for IR galaxies is envisaged: on the one hand, a starburst-dominated phase in which the Super Massive Black Holes (SMBH) grows and is obscured by dust (possibly triggered by a major merging event), is followed by an AGN-dominated phase, then evolving towards a local elliptical. On the other hand, moderately star-forming galaxies containing a low-luminosity AGN have various properties suggesting they are good candidates for systems in a transition phase preceding the formation of steady spiral galaxies.

Herschel↗