Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Radiative capture”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
NASA NTRS2001

Large Eddy Simulation of a Near Sonic Turbulent Jet and Its Radiated Noise

In this paper numerical simulations are used to calculate the turbulence dynamics simultaneously with the sound field for a high-speed near-sonic (Ma=0.9) compressible jet at two Reynolds numbers of 3,600 and 72,000. LES (Large Eddy Simulation) in conjunction with accurate numerical schemes is used to calculate the unsteady flow and sound in the near field of the jet. It is shown that the jet mean parameters, mean velocity fields and turbulence statistics are in good agreement with experimental data and results from other simulations. The sound in the near-field is calculated directly from the simulations. The calculations are shown to capture the peak in the dilatation and pressure spectra around a Strouhal number St=0.25-0.3, in agreement with typical jet-noise spectra measured in experiments. Dilatation contours in the near-field show the formation of acoustic waves with a dominant wavelength of 3.2-4 jet diameters, corresponding to the peak in the dilatation spectra. As expected, the non-compact noise sources are found to be most dominant in the region corresponding to the end of the potential core. The contribution of the LES model to the radiated noise appears to be weak and does not contaminate the sound field with spurious high-frequency noise. However, the frequency spectra of the sound show a rapid falloff away from the peak frequency. This is attributed to the quasi-laminar state of the shear-layers in the region prior to potential core closure, and a possible effect of insufficient azimuthal resolution at the observed location. Further analysis of the effect of the LES model, especially at high frequencies, is needed.

Constantinescu, G. S.↗
NASA NTRS2006

Standardized Radiation Shield Design Methods: 2005 HZETRN

Research committed by the Langley Research Center through 1995 resulting in the HZETRN code provides the current basis for shield design methods according to NASA STD-3000 (2005). With this new prominence, the database, basic numerical procedures, and algorithms are being re-examined with new methods of verification and validation being implemented to capture a well defined algorithm for engineering design processes to be used in this early development phase of the Bush initiative. This process provides the methodology to transform the 1995 HZETRN research code into the 2005 HZETRN engineering code to be available for these early design processes. In this paper, we will review the basic derivations including new corrections to the codes to insure improved numerical stability and provide benchmarks for code verification.

Wilson, John W.↗
NASA NTRSDate not supplied

Benchmarking GOCART-2G in the Goddard Earth Observing System (GEOS)

The Goddard Chemistry Aerosol Radiation and Transport (GOCART) model, which controls the sources sinks and chemistry within the Goddard Earth Observing System, recently underwent a major refactoring and update to the representation of physical processes. The code refactoring increases flexibility such multiple instances of an aerosol species can be run and interact with radiation and cloud microphysics, in addition to the output of multiple wavelength aerosol optical properties in support of data assimilation. From a science perspective, a new radiatively active tracer, brown carbon, was added to distinguish smoke from other sources of organic aerosol thereby improving optical properties entering the radiative calculations. A four-year benchmark simulation was evaluated using in situ and space borne measurements to develop a baseline and prioritize future development. A comparison of simulated aerosol optical depth between GOCART-2G and MODIS retrievals indicates the model captures the overall spatial pattern and seasonal cycle of aerosol optical depth but overestimates aerosol extinction over dusty regions and underestimates aerosol extinction over northern hemisphere boreal forests, requiring further tuning of emissions. This MODIS-based analysis is corroborated by comparisons to MISR and selected AERONET stations. Despite the underestimate of aerosol optical depth in biomass burning regions in GEOS, there is an overestimate in the surface mass of organic carbon in the United States, especially during the summer months.

Allison Collow↗
NASA NTRS2024

Benchmarking GOCART-2G in the Goddard Earth Observing System (GEOS)

The Goddard Chemistry Aerosol Radiation and Transport (GOCART) model, which controls the sources sinks and chemistry within the Goddard Earth Observing System, recently underwent a major refactoring and update to the representation of physical processes. This paper serves to document code changes that were included in GOCART 2nd Generation (GOCART-2G) and establishes a benchmark simulation that is to be used for future development of the system. The code refactoring increases flexibility such multiple instances of an aerosol species can be run and interact with radiation and cloud microphysics, in addition to the output of multiple wavelength aerosol optical properties in support of data assimilation. From a science perspective, a new radiatively active tracer, brown carbon, was added to distinguish smoke from other sources of organic aerosol thereby improving optical properties entering the radiative calculations. A four-year benchmark simulation was evaluated using in situ and space borne measurements to develop a baseline and prioritize future development. A comparison of simulated aerosol optical depth between GOCART-2G and MODIS retrievals indicates the model captures the overall spatial pattern and seasonal cycle of aerosol optical depth but overestimates aerosol extinction over dusty regions and underestimates aerosol extinction over northern hemisphere boreal forests, requiring further tuning of emissions. This MODIS-based analysis is corroborated by comparisons to MISR and selected AERONET stations. Despite the underestimate of aerosol optical depth in biomass burning regions in GEOS, there is an overestimate in the surface mass of organic carbon in the United States, especially during the summer months.

Allison B Collow↗
NASA NTRS1980

Climatic effects of atmospheric aerosols

It is shown quantitatively how radiative climatic forcing by aerosols depends on the physical properties of the aerosols. The special case of atmospheric aerosols produced by volcanic explosions is considered, and evidence is presented which indicates that even the simple climate models available today may be able to capture some of the basic effects of aerosols on global climate. Possible effects of anthropogenic aerosols are briefly examined. It is noted that additional evidence for the significant influence of aerosols on climate may be found in data obtained from recent spacecraft missions to Venus and Mars.

Hansen, J. E.↗
NASA NTRS2019

El Niño Related Tropical Land Surface Water and Energy Response in MERRA-2

The El Niño Southern Oscillation (ENSO) is a coupled Earth system circulation phenomena that reaches all around the globe. The heat added to the atmosphere by increased precipitation produces circulation changes that have global reach and, over time, warms the entire tropical band, and much of the Earth. Many studies have noted that El Niño causes warm and dry (and sometimes drought) conditions over tropical land masses. We develop a composite analysis of El Niño to identify the predominate features of tropical land response. This analysis shows that the land’s lagged response is related to a reduction clouds that leads to increase surface shortwave radiation that increases the surface temperature. The precipitation lag is somewhat longer, and then leads to a reduction in soil water and, in concert with increased SW induced surface warming, leads to increased sensible heating of the atmosphere above. The M2AMIP simulation generally captures these features, but the response is strongest with increased temporal and spatial proximity to the El Niño peak warming. The regionality of these features is also discussed, and it is noted that even the strongest individual El Niño events can vary from this composite mechanistic paradigm.

Bosilovich, Michael G.↗
NASA NTRS2022

Vulnerability of Passive Microwave Snowfall Retrievals to Physical Properties of Snowpack: A Perspective from Dense Media Radiative Transfer Theory

The uncertainty of passive microwave retrievals of snowfall is notoriously high where the high-frequency surface emissivity is significantly reduced and varies markedly in response to changes of snowpack physical properties. Using the dense media radiative transfer theory, this article studies the potential effects of terrestrial snow-cover depth, density, and6grain size on high-frequency channels 89 and 166 GHz of the radiometer onboard the Global Precipitation Measurement (GPM) core satellite, which are commonly used to capture the snowfall scattering signals. Integrating the inference across all feasible grain sizes, ranges of snowpack density and depth are identified over which the snowfall scattering signatures can be time varying and potentially obscured. Using 10 years of reanalysis data, the seasonal vulnerability of snowfall retrievals to changes of snowpack emissivity in the Northern Hemisphere is mapped in a probabilistic sense and connections are made with uncertainties of the GPM passive microwave snowfall retrievals. It is found that among different snow classes, relatively light Arctic tundra snow in fall, with a density below 260kg m−3, and shallow prairie snow during the winter, with a depth of less than 40cm, can reduce the surface emissivity and obscure the snowfall passive microwave signatures. It is demonstrated that, during the winter, the highly vulnerable areas are over Kazakhstan, and Mongolia with taiga and prairie snow. In the fall, these areas are largely over tundra and taiga snow in north of Russia and the Arctic Archipelagos as well as prairies in Canada and the Great Plains in the United States.

Reyhaneh Rahimi↗
NASA NTRS1975

Dynamic analysis of the long-term zonal earth energy budget experiment (LZEEBE) spacecraft

The motions of the flexible LZEEBE spacecraft after injection into a circular orbit were analyzed. The spacecraft consists of three spherical balloons connected to a central hub by three long flexible booms which lie in a plane such that the angle between any pair of booms is 120 degrees. The major torques acting on the spacecraft are gravity-gradient torques and torques due to solar radiation pressure acting on the balloons which have different reflective properties. A development of the equations of motion is presented. Computer simulations indicate that the spacecraft will have random motion, provided the injection into orbit does not create conditions favorable to gravity-gradient capture.

Meirovitch, L.↗
NASA NTRS1994

Some considerations on velocity vector accuracy in dust trajectory analysis

The relative contributions of comets and asteroids to the reservoir of dust in the interplanetary medium is not known. There are direct observations of dust released from comets and there is evidence to associate the IRAS dust bands with possible collisions of asteroids in the main belt. A means towards sorting out the parent sources has been proposed in the establishment of a dust collector in orbit about the Earth. The purpose of such a facility would be to collect not only cosmic dust particles intact but also the state vectors, as they arrive at the detector, the idea being that one may combine analytical laboratory analysis of the physics and chemistry of the captured particles with orbital data in order to help distinguish between bodies and identify parent bodies. The theoretical study of dust particle orbits in the solar system takes on greatly more importance if we use collected trajectory data. The orbital motion of dust when radiation and forces alone are acting is well understood. When gravitational forces due to the planets are included, the motion can become quite complex. In order to characterize the orbits of particles as they crossed the Earth's orbits, a study of the long-time dust orbital evolution was undertaken. We have considered various parameters associated with these dust orbits to see if one may in a general way discriminate between particles evolved from comets and asteroids. We proceed in this study as we have done previously. That is, we considered the dust particles as ideal black bodies, of density 1 gm/cc, spherical, with radii 10-100 microns. Particles of this size are affected by radiation forces, photon pressure, and Poynting-Robertson drag. Account was also taken of solar wind drag, which amounts to about 30 percent of the Poynting-Robertson drag negligible. The gravitational forces due to the planets are included, unlike in our previous study; the planetary orbits are those of true n-body interaction so that the possibility of secular resonance is included. Our method was to calculate explicitly by a numerical procedure the orbits of dust particles after they left their parent bodies. The motion is determined numerically with the implicit Runge-Kutta integrator using Gauss-Radau spacings.

Jackson, A. A.↗
NASA NTRS2018

Multi-Sensor Integration of Vegetation Index Products for Long-Term Monitoring of Vegetation Dynamics: A Case Study from MODIS to VIIRS

Spectral vegetation index (VI) time series data from moderate resolution sensors, such as Earth Observing System (EOS) Moderate Resolution Imaging Spectroradiometer (MODIS), have widely been used to identify "hot spot" areas of vegetation changes and to characterize long-term trends of vegetation changes. The Visible Infrared Imaging Radiometer Suite (VIIRS) sensor series of the Joint Polar Satellite System (JPSS) program is slated to continue the highly calibrated data stream initiated with MODIS. The first VIIRS sensor has operationally been acquiring Earth-reflected radiation since March 2012. The second VIIRS sensor (NOAA-20) was launched in November 2017, which is to continue the measurements made by the first VIIRS sensor. In this study, we examined an integrated use of MODIS and VIIRS VI time series data on capturing vegetation dynamics in the Asia-Pacific region using their overlapped period of observations (2013-2017). Three VIs, the "top-of-canopy (TOC)" normalized difference vegetation index (NDVI), TOC enhanced vegetation index (EVI), and TOC two-band enhanced vegetation index (EVI2), were investigated. For all the three VIs, VIIRS VIs were systematically higher than the MODIS counterparts due mainly to their spectral bandpass differences. However, both VIIRS and MODIS VIs showed the comparable spatial patterns in their temporal variations. Empirical spectral corrections allowed to merge the two data streams, which slightly improved the temporal resolution of the VI temporal profiles. These results suggest the suitability of VIIRS data to extend and merge into the MODIS VI record for long-term vegetation dynamics studies.

MODIS↗
NASA NTRS2019

Flash Drought as Captured by Reanalysis Data: Disentangling the Contributions of Precipitation Deficit and Excess Evapotranspiration

Flash droughts – uncharacteristically rapid dryings of the land system – are naturally associated with extreme precipitation deficits. Such precipitation deficits, however, do not tell the whole story, for land surface drying can be exacerbated by anomalously high evapotranspiration (ET) rates driven by anomalously high temperatures (e.g., during heatwaves), anomalously high incoming radiation (e.g., from reduced cloudiness), and other meteorological anomalies. In this study, the relative contributions of precipitation and ET anomalies to flash drought generation in the Northern Hemisphere are quantified through the analysis of diagnostic fields contained within the MERRA-2 reanalysis product. Unique to the approach is the explicit treatment of soil moisture impacts on ET through relationships diagnosed from the reanalysis data; under this treatment, an ET anomaly that is negative relative to the local long-term climatological mean is still considered positive in terms of its contribution to a flash drought if it is high for the concurrent value of soil moisture. Maps produced in the analysis show the fraction of flash drought production stemming specifically from ET anomalies and illustrate how ET anomalies for some droughts are related to temperature and radiation anomalies. While ET is found to have an important impact on flash drought production in the central US and in parts of Russia known from past studies to be prone to heatwave-related drought, and while this impact does appear stronger during the onset (first several days) of flash droughts, overall the contribution of ET to these droughts is small relative to the contribution of precipitation deficit.

R D Koster↗
NASA NTRS1992

Induced activation study of LDEF

Analysis of the induced radioactivity of LDEF is continuing with extraction of specific activities for various spacecraft materials. Data and results of activation measurements from eight national facilities are being collected for interpretation. The major activation mechanisms in LDEF components is the proton flux in the South Atlantic Anomaly (SAA) inner radiation belt. This flux is highly anisotropic, and exposes the west side of the spacecraft to higher radiation doses. The directionally dependent activation due to these protons has clearly been observed in the data from Al experiment tray clamps, steel trunnions, and is also indicated by the presence of a variety of radioisotopes in other materials. A secondary production mechanism, thermal neutron capture, was observed in two materials having large capture cross sections, Co and Ta. The neutrons could be thermalized in nearby low Z material, although this has yet to be verified. Specific activities are presented for a number of materials which show SAA effects and thermal neutron capture. Trends are examined in the measured results that show the effect of shielding and non-SAA related activation.

Harmon, B. A.↗
NASA NTRS1968

The Use of Large Lithium-Drifted, Germanium Diodes for Gamma-Ray Spectral Measurements at Balloon Altitudes

A series of balloon-borne experiments has been initiated at the Oak Ridge National Laboratory (ORNL) in cooperation with the Space Sciences Laboratory of the NASA G. C. Marshall Space Flight Center to determine the feasibility of using large (greater than 20 cc) lithium-drifted, germanium diodes to investigate the spectrum of atmospheric gamma rays (60 keV < E < 1.5 MeV) at altitudes corresponding to a few grams per cm (sup 2) residual atmosphere. The results of these measurements will also provide a basis for designing an appropriate shield for a multi-diode, highly directional gamma-ray spectrometer for use in astrophysical measurements. Two flights were accomplished during 1967 using the same two diodes and basic instrument package. Both flights were launched at the NCAR Scientific Balloon Flight Base, Palestine, Texas, and attained a float altitude of 117, 000 ft. The measured atmospheric gamma-ray spectrum shows at least one distinct line at 511 keV (annihilation radiation) superimposed on a continuous distribution of gamma rays attributable to both energy-degraded gamma rays and brems-strahlung photons. A second diode, encased in a passive shield of plastic and lead, shows the expected lines resulting from both neutron inelastic scattering and capture in the germanium. Data acquisition was accomplished on an event-by-event basis through the use of 512-channel, on-board ADC and word-generator circuits and a ground-based, on-line telemetry decoder. The decoder makes it possible to store the data in the memory of a modified pulse-height analyzer simultaneously with storage on magnetic tape. This provides a real-time visual observation of the data as it is accumulated and greatly facilitates preflight calibrations.

Chapman, G. T.↗
NASA NTRS2003

Ultrasail

Ultrasail is a complete sail system for the launch, deployment, stabilization and control of very large solar sails enabling reduced mission times for interplanetary and deep space spacecraft. Ultrasail is an innovative, non-traditional approach to propulsion technology achieved by combining propulsion and control systems developed for formation-flying microsatellites with an innovative solar sail architecture to achieve sq km-class controllable sail areas, sail subsystem area densities of 1 gm per sq m, and thrust levels equivalent to 400 kW ion thruster systems used for comparable deep space missions. Ultrasail can conceivably even achieve outer planetary rendezvous, a deep space capability now reserved for high-mass nuclear and chemical systems. Ultrasail is a Delta IV-launched multi-blade spin-stabilized system with blade lengths as long as 50 km, reminiscent of the MacNeal Heliogyro. The primary innovation is the near-elimination of sail supporting structures by attaching the sail tip to a rigid formation-flying microsatellite truss which deploys the sail blade, and which then articulates the blade to provide attitude control, including spin stabilization and precession of the spin axis. These tip microsatellites are controlled by a solar-powered 3-axis microthruster system (electric or cold gas) to maintain proper sail film tension during deployment and spin-up. The satellite mass also provides a stabilizing centrifugal force on the blade while in rotation. Understanding the dynamics of individual blades is key to the overall dynamics of Ultrasail. Forces and torques that must be modeled include those due to solar pressure, those generated by the microsatellite at the blade tip and by torques applied at the blade root. Centrifugal forces also play a significant role in the deployment and maintenance of the sail configuration. To capture the dynamics of the overall system, the equations of motion for the blades have been derived. Using these differential equations, a control law will be derived to maneuver Ultrasail. This law involves the pitching of the individual blades thereby moving the distribution of the radiation pressure on each individual blade and inducing a resultant torque on the system. The direction of the angular momentum vector and its rate of precession can be controlled through the pitch angle of the blades. The Ultrasail trajectory is also being studied. Optimal or near-optimal trajectories are being generated to showcase Ultrasail performance. Various missions, e.g. outer planet and solar polar missions for observation of the Sun, are currently being investigated to demonstrate the performance enhancements generated by Ultrasail technology. Calculus-of-variations-based optimization software is used to produce optimal Ultrasail trajectories. The performance of these trajectories is being compared to optimal results generated with other propulsion models, including chemical propulsion, ion propulsion, and competing solar sail concepts. Results of these studies will quantify the performance of Ultrasail compared to existing solar sail concepts for high energy missions.

Burton, R.↗
NASA NTRS2018

Flash Drought as Captured by MERRA-2: Disentangling the Contributions of Precipitation Deficit and Excess Energy

Flash droughts uncharacteristically rapid dryings of the land system are naturally associated with extreme precipitation deficits. Such precipitation deficits, however, do not tell the whole story. Land surface drying can be exacerbated by anomalously high evapotranspiration (ET) rates driven, for example, by anomalously high temperatures (e.g., during heatwaves) and/or anomalously high incoming radiation (e.g., from reduced cloudiness). In this study the relative contributions of precipitation and ET anomalies to flash drought generation are quantified through the analysis of diagnostic fields contained within the MERRA-2 reanalysis product. Unique to the overall approach is the explicit treatment of soil moisture impacts on ET an ET anomaly that is negative relative to the local long-term climatological mean is still considered positive in terms of its contribution to a flash drought, if it is high for the concurrent value of soil moisture. Maps produced in the analysis show the fraction of flash drought production stemming specifically from ET anomalies and the degree to which these ET anomalies are related to temperature and radiation anomalies. ET is found to have a large impact, for example, on flash drought production in the east-central US and in parts of Russia known from past studies to be prone to heatwave-related drought.

Koster, Randal↗
NASA NTRS2016

A Quasi-Global Approach to Improve Day-Time Satellite Surface Soil Moisture Anomalies through the Land Surface Temperature Input

Passive microwave observations from various spaceborne sensors have been linked to the soil moisture of the Earth's surface layer. A new generation of passive microwave sensors are dedicated to retrieving this variable and make observations in the single theoretically optimal L-band frequency (1-2 GHz). Previous generations of passive microwave sensors made observations in a range of higher frequencies, allowing for simultaneous estimation of additional variables required for solving the radiative transfer equation. One of these additional variables is land surface temperature, which plays a unique role in the radiative transfer equation and has an influence on the final quality of retrieved soil moisture anomalies. This study presents an optimization procedure for soil moisture retrievals through a quasi-global precipitation-based verification technique, the so-called Rvalue metric. Various land surface temperature scenarios were evaluated in which biases were added to an existing linear regression, specifically focusing on improving the skills to capture the temporal variability of soil moisture. We focus on the relative quality of the day-time (01:30 pm) observations from the Advanced Microwave Scanning Radiometer for Earth Observing System (AMSR-E), as these are theoretically most challenging due to the thermal equilibrium theory, and existing studies indicate that larger improvements are possible for these observations compared to their night-time (01:30 am) equivalent. Soil moisture data used in this study were retrieved through the Land Parameter Retrieval Model (LPRM), and in line with theory, both satellite paths show a unique and distinct degradation as a function of vegetation density. Both the ascending (01:30 pm) and descending (01:30 am) paths of the publicly available and widely used AMSR-E LPRM soil moisture products were used for benchmarking purposes. Several scenarios were employed in which the land surface temperature input for the radiative transfer was varied by imposing a bias on an existing regression. These scenarios were evaluated through the Rvalue technique, resulting in optimal bias values on top of this regression. In a next step, these optimal bias values were incorporated in order to re-calibrate the existing linear regression, resulting in a quasi-global uniform LST relation for day-time observations. In a final step, day-time soil moisture retrievals using the re-calibrated land surface temperature relation were again validated through the Rvalue technique. Results indicate an average increasing Rvalue of 16.5%, which indicates a better performance obtained through the re-calibration. This number was confirmed through an independent Triple Collocation verification over the same domain, demonstrating an average root mean square error reduction of 15.3%. Furthermore, a comparison against an extensive in situ database (679 stations) also indicates a generally higher quality for the re-calibrated dataset. Besides the improved day-time dataset, this study furthermore provides insights on the relative quality of soil moisture retrieved from AMSR-E's day- and night-time observations.

soil moisture↗
NASA NTRS1993

Orbital resonances and Poynting-Robertson drag

The phenomenon of resonance trapping with Poynting-Robertson drag in the simplest case - the circular restricted three-body problem - is elucidated. Attention is given to what determines whether a grain of a given size passes through a given resonance or is trapped there, to how and why a trapped particle's orbit evolves with time, and to why Poynting-Robertson drag resonances are only temporary, while gas-drag resonances appear to be stable. The possibility of trapping a grain into resonance with a planet depends on the combination of the following parameters: the ratio of radiation pressure force to solar gravity, the mass of the perturbing planet normalized to the solar mass, an integer, and eccentricity. In general, the peak eccentricity and sometimes the threshold value are large enough so that crossing orbits and close approaches to the planet can inhibit capture and aid escape from resonance.

Weidenschilling, S. J.↗
NASA NTRS2001

Toward a Unified Parameterization of the Boundary Layer and Moist Convection: Simulations of Clear and Cloudy Convection - Part 3

A model that employs a new form of mass-flux closure (described in Part I of this paper) is applied to a variety of clear and cloudy planetary boundary layers (PBLs) including dry convection from the Wangara Experiment, trade wind cumulus from the Barbados Oceanographic and Meteorological Experiment (BOMEX), and marine stratocumulus from the Atlantic Stratocumulus Experiment (ASTEX). For Wangara, the simulated variances and fluxes match that expected from similarity arguments, while the mean state is a little less mixed than the observations. In the BOMEX simulation. the shape and magnitude of the fluxes and the turbulence kinetic energy budget agree with LES results and observations. However, the liquid water mixing ratio is too large. This is attributed to an underprediction of the skewness. In agreement with observations from the ASTEX experiment, many of the model-simulated fields distinctly reflect a regime in transition between the trade wind cumulus and the classic stratocumulus-topped boundary layers. In general, the simulated entrainment rate tends to be a little underpredicted in regimes where there is little cloud-top radiative cooling (Wangara and BOMEX), while it is overpredicted in regimes where this process is more critical (e.g., ASTEX). Prior work suggests that this may be related to the manner in which the pressure terms are parameterized in the model. Overall, the model is able to capture some key physical features of these PBL regimes, and appears to have the potential to represent both cloud and boundary layer processes. Thus, this approach is a first step toward unifying these processes in large-scale models.

Lappen, Cara-Lyn↗