Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “large deviations”

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 973 records · Page 54

Exploring the parameter space of MagLIF implosions using similarity scaling. II. Current scaling

Magnetized liner inertial fusion (MagLIF) is a magneto-inertial-fusion (MIF) concept, which is presently being studied on the Z pulsed power facility. The MagLIF platform has achieved interesting plasma conditions at stagnation and produced significant fusion yields in the laboratory. Given the relative success of MagLIF, there is a strong interest to scale the platform to higher peak currents. However, scaling MagLIF is not entirely straightforward due to the large dimensionality of the experimental input parameter space and the numerous physical processes involved in MIF implosions. Here, in this work, we propose a novel method to scale MagLIF loads to higher currents. Our method is based on similarity (or similitude) scaling and attempts to preserve much of the physics regimes already known or being studied on today's Z pulsed-power driver. By avoiding significant deviations into unexplored and/or less well-understood regimes, the risk of unexpected outcomes on future scaled-up experiments is reduced. Using arguments based on similarity scaling, we derive the scaling rules for the experimental input parameters characterizing a MagLIF load (as functions of the characteristic current driving the implosion). We then test the estimated scaling laws for various metrics measuring performance against results of 2D radiation–magneto-hydrodynamic hydra simulations. Agreement is found between the scaling theory and the simulation results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Much of zero emissions commitment occurs before reaching net zero emissions

Abstract We explore the response of the Earth’s coupled climate and carbon system to an idealized sequential addition and removal of CO 2 to the atmosphere, following a symmetric and continuous emissions pathway, in contrast to the discontinuous emissions pathways that have largely informed our understanding of the climate response to net zero and net negative emissions to date. We find, using both an Earth system model and an ensemble of simple climate model realizations, that warming during the emissions reduction and negative emissions phases is defined by a combination of a proportionality of warming to cumulative emissions characterized by the transient climate response to emissions (TCRE), and a deviation from that proportionality that is governed by the zero emissions commitment (ZEC). About half of the ZEC is realized before reaching zero emissions, and the ZEC thus also controls the timing between peak cumulative CO 2 emissions and peak temperature, such that peak temperature may occur before peak cumulative emissions if ZEC is negative, underscoring the importance of ZEC in climate policies aimed to limit peak warming. Thus we argue that ZEC is better defined as the committed warming relative to the expected TCRE proportionality, rather than as the additional committed warming that will occur after reaching net zero CO 2 emissions. Once established, the combined TCRE and ZEC relationship holds almost to complete removal of prior cumulative CO 2 emissions. As cumulative CO 2 emissions approach zero through negative CO 2 emissions, CO 2 concentrations drop below preindustrial values, while residual long-term climate change continues, governed by multicentennial dynamical processes.

54 ENVIRONMENTAL SCIENCES↗

The nonaxisymmetric baroclinic instability in thin accretion disks

The instability of linear, nonaxisymmetric baroclinic waves in the radiative regions of a thin accretion disk is investigated with simplifying assumptions used in the representation of the structure of the basic state. As in stars with rotation laws deviating from rotation on cylinders, the departure from barotropy ('baroclinicity') resulting from such rotation feeds potential energy from the basic state into the kinetic energy of the baroclinic motions. The baroclinic instability is considered as a possible mechanism for driving turbulence in disks. Estimations of the Reynolds stress from baroclinic waves indicate that the degree of baroclinicity in stationary disk models is insufficient to account for the large viscosity inferred from observations of, e.g., cataclysmic variables.

Cabot, W.↗

Comets, impacts, and atmospheres

Studies of element abundances and values of D/H in the atmospheres of the giant planets and Titan have emphasized the important role of icy planetesimals in the formation of these bodies. In these atmospheres, C/H and D/H increase as the relative masses of the 'cores' of the planets increase. N/H appears to deviate from this trend in an interesting way. In the inner solar system, the traditional approach of using carbonaceous chondrites as the source of planetary volatiles is in serious trouble because of the depletion of xenon and the unusual pattern of xenon isotopes found in the atmospheres of Earth and Mars, and because of the solar-type abundance ratios of argon, krypton and xenon and the large amounts of neon and argon on Venus. Recent studies of elemental abundances in comets, especially P/Halley, coupled with laboratory studies of the trapping of gas in ice formed at low temperatures by A. Bar-Nun et al. provide a consistent interpretation of all of these results. This interpretation emphasizes the fundamental importance of icy planetesimals (comets) and the randomness of early impacts in the formation of planetary systems. Cometary delivery by itself will not explain the noble gas abundances on the inner planets. There is good evidence for at least one additional source, which presumably consists of the rocky material making up the bulk of the planets. The existence of this rocky reservoir is manifested in the nucleogenic isotopes and in the neon which is found in all these atmospheres and is also present in the Earth's mantle. This neon may well be a relic of the planets' earliest, accretional atmospheres.

Owen, Tobias↗

M2-SCREAM: A Stratospheric Composition Reanalysis of Aura MLS Data With MERRA-2 Transport

MERRA-2 Stratospheric Composition Reanalysis of Aura Microwave Limb Sounder (M2-SCREAM) is a new reanalysis of stratospheric ozone, water vapor, hydrogen chloride (HCl), nitric acid (HNO 3 ) and nitrous oxide (N 2 O) between 2004 and the present (with a latency of several months). The assimilated fields are provided at a 50-km horizontal resolution and at a three-hourly frequency. M2-SCREAM assimilates version 4.2 Microwave Limb Sounder (MLS) profiles of the five constituents alongside total ozone column from the Ozone Monitoring Instrument. Dynamics and tropospheric water vapor are constrained by the MERRA-2 reanalysis. The assimilated species are in excellent agreement with the MLS observations, except for HNO 3 in polar night, where data are not assimilated. Comparisons against independent observations show that the reanalysis realistically captures the spatial and temporal variability of all the assimilated constituents. In particular, the standard deviations of the differences between M2-SCREAM and constituent mixing ratio data from The Atmospheric Chemistry Experiment Fourier Transform Spectrometer are much smaller than the standard deviations of the measured constituents. Evaluation of the reanalysis against aircraft data and balloon-borne frost point hygrometers indicates a faithful representation of small-scale structures in the assimilated water vapor, HNO 3 and ozone fields near the tropopause. Comparisons with independent observations and a process-based analysis of the consistency of the assimilated constituent fields with the MERRA-2 dynamics and with large-scale stratospheric processes demonstrate the utility of M2-SCREAM for scientific studies of chemical and transport variability on time scales ranging from hours to decades. Analysis uncertainties and guidelines for data usage are provided.

MERRA-2↗

BOSS DR12 full-shape cosmology: Λ CDM constraints from the large-scale galaxy power spectrum and bispectrum monopole

Here, we present a full Λ CDM analysis of the BOSS DR12 dataset, including information from the power spectrum multipoles, the real-space power spectrum, the reconstructed power spectrum and the bispectrum monopole. This is the first analysis to feature a complete treatment of the galaxy bispectrum, including a consistent theoretical model and without large-scale cuts. Unlike previous works, the statistics are measured using window-free estimators: this greatly reduces computational costs by removing the need to window-convolve the theory model. Our pipeline is tested using a suite of high-resolution mocks and shown to be robust and precise, with systematic errors far below the statistical thresholds. Inclusion of the bispectrum yields consistent parameter constraints and shrinks the σ 8 posterior by 13% to reach < 5 % precision; less conservative analysis choices would reduce the error bars further. Our constraints are broadly consistent with Planck: in particular, we find H 0 = 69.6$^{+1.1}_{–1.3}$ km s –1 Mpc –1 , σ 8 = 0.692$^{+0.035}_{–0.041}$ and n s = 0.870$^{+0.067}_{–0.064}$, including a BBN prior on the baryon density. When n s is set by Planck, we find H 0 = 68.31$^{+0.83}_{–0.86}$ km s –1 Mpc –1 and σ 8 = 0.722$^{+0.032}_{–0.036}$. Our S 8 posterior, 0.751 ± 0.039, is consistent with weak lensing studies, but lower than Planck. Constraints on the higher-order bias parameters are significantly strengthened from the inclusion of the bispectrum, and we find no evidence for deviation from the dark matter halo bias relations. These results represent the most complete full-shape analysis of BOSS DR12 to-date, and the corresponding spectra will enable a variety of beyond-Λ CDM analyses, probing phenomena such as the neutrino mass and primordial non-Gaussianity.

79 ASTRONOMY AND ASTROPHYSICS↗

Assessing the thermal maturity of black shales using vitrinite reflectance: Insights from Devonian black shales in the eastern United States

Thermal maturity of source rocks indicated by vitrinite reflectance (R o ) is an important and reliable parameter to determine the petroleum potential of sedimentary basins. Originally, R o was used to characterize the degree of coalification of coals and thus works well for coal seams, but it has limitations when used to determine the thermal maturity of black shale successions. A suite of six coal-shale pairs from the Upper Devonian black shale formations in the eastern United States (New Albany Shale, Ohio Shale, Dunkirk Shale, and Rhinestreet Shale) was selected to study the applicability of R o in assessing the thermal maturity of black shale successions. The results show that vitrinite in the studied coal samples is dominated by collotelinite, whereas vitrinite in black shales occurs as small dispersed particles (~5 μm) in the mineral matrix. When comparing the size and morphology of dispersed vitrinite particles in shales and collotelinite fragments in coals, vitrinite in shales and adjacent coals should have the same origin. The measured mean random R o of vitrinite in coals ranges from 0.51–0.68%, and is 0.04–0.11% (average 0.07%) lower than that of dispersed vitrinite particles in enclosing shales. This observation contrasts with previously reported R o suppression in liptinite-rich black shales. A high standard deviation of R o measurements indicates a highly heterogeneous nature of vitrinite in black shales. Three mechanisms can be envisioned to contribute to the higher R o of dispersed vitrinite particles in black shales relative to enclosed coals. First, small vitrinite particles may become more oxidized because small particles are more frequently suspended during transport to the site of deposition and reworking of bottom sediments after deposition than large pieces of driftwood that will turn into enclosed coal lenses during burial diagenesis. Second, misidentification of zooclast (e.g., chitinozoan) fragments as vitrinite due to loss of diagnostic morphology. Third, various macerals in the vitrinite group (e.g., corpogelinite and collotelinite) may have had different original reflectance values. Based on the petrographic characteristics of vitrinite and vitrinite-like particles, caution should be applied when using the R o values of dispersed vitrinite and vitrinite-like particles in black shales as indicators of thermal maturity.

58 GEOSCIENCES↗

Revealing the short-range structure of the mirror nuclei 3 H and 3 He

When protons and neutrons (nucleons) are bound into atomic nuclei, they are close enough to feel significant attraction, or repulsion, from the strong, short-distance part of the nucleon-nucleon interaction. These strong interactions lead to hard collisions between nucleons, generating pairs of highly energetic nucleons referred to as short-range correlations (SRCs). SRCs are an important but relatively poorly understood part of nuclear structure, and mapping out the strength and the isospin structure (neutron-proton (np) versus proton-proton (pp) pairs) of these virtual excitations is thus critical input for modelling a range of nuclear, particle and astrophysics measurements. Two-nucleon knockout or 'triple coincidence' reactions have been used to measure the relative contribution of np-SRCs and pp-SRCs by knocking out a proton from the SRC and detecting its partner nucleon (proton or neutron). These measurements have shown that SRCs are almost exclusively np pairs, but they had limited statistics and required large model-dependent final-state interaction corrections. Here we report on measurements using inclusive scattering from the mirror nuclei hydrogen-3 and helium-3 to extract the np/pp ratio of SRCs in systems with a mass number of three. In this work we obtain a measure of the np/pp SRC ratio that is an order of magnitude more precise than previous experiments, and find a marked deviation from the near-total np dominance observed in heavy nuclei. This result implies an unexpected structure in the high-momentum wavefunction for hydrogen-3 and helium-3. Understanding these results will improve our understanding of the short-range part of the nucleon-nucleon interaction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Computational investigations of Dienes defect- and vacancy-induced changes in the electronic and vibrational properties of carbon fiber structural units

Carbon fiber (CF) is a promising lightweight alternative to steel and is of significant interest for energy applications. As CF continues to find new uses and is exposed to new external conditions, a noninvasive method of monitoring its structural integrity is critical. Raman spectroscopy is a commonly used method for this monitoring; however, it is highly inferential, and the interpretation of the data is not always straightforward. In this work, we perform density functional theory (DFT) calculations to investigate changes in the vibrational properties of CF structural units (i.e., graphene and graphite) caused by monovacancy and Dienes defects as a foundation for modeling more complex defects that move our model toward that of realistic CF. Using large computational supercells, we can understand how these defects change the electronic structure and vibrational properties of graphene and graphite for interdefect distances near those of the lower experimental limit. The monovacancy opens an electronic bandgap at the K point. Although no such electronic gap is opened by the Dienes defect, both defects introduce flat defect bands near the Fermi energy. The Dienes defect creates long-range deviations of the phonons, leading to substantial broadening of the highest frequency optical modes in the band structure compared to that of the pristine material. In contrast, the phonon changes caused by the monovacancy are short range, and only minor changes in the band structure or phonon density of states were observed. These findings can assist in the interpretation of experimental results by providing atomic-scale insight into key electronic and vibrational features.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Local and global approaches to treat the torsional barriers of 4-methylacetophenone using microwave spectroscopy

The Fourier transform microwave spectrum of 4-methylacetophenone recorded from 8 GHz to 18 GHz under jet-cooled conditions has revealed large tunneling splittings arising from a low barrier to internal rotation of the ring methyl group and small splittings from a high torsional barrier of the acetyl methyl group. The large splittings are especially challenging to model, while the small splittings are difficult to analyze due to the resolution limit of 120 kHz. The combination of two methyl groups undergoing internal rotations caused each rotational transition to split into five torsional species, which were resolved and fitted using a modified version of the XIAM code and the newly devel- oped ntop code to a root-mean-square deviation close to measurement accuracy, providing an estimate of the V 3 potential barriers of about 22 cm -1 and 584–588 cm -1 for the ring and the acetyl methyl groups, respectively. The assignment was aided by separately fitting the five torsional species using odd-power order operators. Only one conformer in which all heavy atoms are located on a symmetry plane could be identified in the spectrum, in agreement with results from conformation analysis using quantum chemical calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Rain attenuation measurements: Variability and data quality assessment

Year to year variations in the cumulative distributions of rain rate or rain attenuation are evident in any of the published measurements for a single propagation path that span a period of several years of observation. These variations must be described by models for the prediction of rain attenuation statistics. Now that a large measurement data base has been assembled by the International Radio Consultative Committee, the information needed to assess variability is available. On the basis of 252 sample cumulative distribution functions for the occurrence of attenuation by rain, the expected year to year variation in attenuation at a fixed probability level in the 0.1 to 0.001 percent of a year range is estimated to be 27 percent. The expected deviation from an attenuation model prediction for a single year of observations is estimated to exceed 33 percent when any of the available global rain climate model are employed to estimate the rain rate statistics. The probability distribution for the variation in attenuation or rain rate at a fixed fraction of a year is lognormal. The lognormal behavior of the variate was used to compile the statistics for variability.

Crane, Robert K.↗

Kinetic Alfven waves as a source of plasma transport at the dayside magnetopause

As the shocked solar wind with variable plasma density and magnetic field impinges on the dayside magnetopause, it is likely to generate large-scale Alfven waves at the solar wind magnetosphere interface. However, large gradients in the density and magnetic field at the magnetopause boundary effectively couple large-scale Alfven waves with kinetic Alfven waves may play an important role in plasma transport at the dayside magnetopause and in electron acceleration along field lines. The transport can occur because, unlike the magnetohydrodynamic (MHD) shear Alfven wave, the kinetic Alfven wave has an associated parallel electric field which breaks down the 'frozen-in' condition and decouples the plasma from field lines. We calculate the average deviation of the plasma from the field line from which we estimate the diffusion coefficient associated with these 'bundles' of lead to acceleration of electrons along field lines in the magnetopause boundary and may possibly provide an explanation for observed counterstreaming electron beams characterized by energies of 50-200 eV.

Lee, L. C.↗

Influence of Subpixel Scale Cloud Top Structure on Reflectances from Overcast Stratiform Cloud Layers

Recent observational studies have shown that satellite retrievals of cloud optical depth based on plane-parallel model theory suffer from systematic biases that depend on viewing geometry, even when observations are restricted to overcast marine stratus layers, arguably the closest to plane parallel in nature. At moderate to low sun elevations, the plane-parallel model significantly overestimates the reflectance dependence on view angle in the forward-scattering direction but shows a similar dependence in the backscattering direction. Theoretical simulations are performed that show that the likely cause for this discrepancy is because the plane-parallel model assumption does not account for subpixel, scale variations in cloud-top height (i.e., "cloud bumps"). Monte Carlo simulation, comparing ID model radiances to radiances from overcast cloud field with 1) cloud-top height variation, but constant cloud volume extinction; 2) flat tops but horizontal variations in cloud volume extinction; and 3) variations in both cloud top height and cloud extinction are performed over a approximately equal to 4 km x 4 km domain (roughly the size of an individual GAC AVHRR pixel). The comparisons show that when cloud-top height variations are included, departures from 1D theory are remarkably similar (qualitatively) to those obtained observationally. In contrast, when clouds are assumed flat and only cloud extinction is variable, reflectance differences are much smaller and do not show any view-angle dependence. When both cloud-top height and cloud extinction variations are included, however, large increases in cloud extinction variability can enhance reflectance difference. The reason 3D-1D reflectance differences are more sensitive to cloud-top height variations in the forward-scattering direction (at moderate to low, sun elevations) is because photons leaving the cloud field in that direction experience fewer scattering events (low-order scattering) and are restricted to the topmost portions of the cloud. While reflectance deviations from 1D theory are much larger for bumpy clouds than for flat clouds with variable cloud extinction, differences in cloud albedo are comparable for these two cases.

Loeb, N. G.↗

Electric Field Double Probe Measurements on Satellites with Very Low Perigees

Measurements of DC and AC electric Fields at low altitudes (less than 150 km) in the Earth's ionosphere address a wide array of important scientific questions including Joule Heating, the closure of field aligned currents, and the deviation of E x B electron velocities from ion velocity vectors influenced by collisions. Double probes represent a well-proven technique for gathering high quality DC and AC electric field measurements for which the design of the boom system is of critical importance for the success of the experiment. In general, ionospheric DC electric field instruments that achieve accuracies of 0.1 mV/m or better, place sensors at large distances from the spacecraft body in order to extend well beyond the spacecraft wake and sheath and to achieve the large signal-to-noise ratios for DC and long wavelength measurements. Additional sets of sensors inboard of the primary, outermost sensors provide useful additional information, both for diagnostics of the plasma contact potentials, which particularly enhance the DC electric field measurements on non-spinning spacecraft, and for wavelength and phase velocity measurements that use the spaced receiver or "Interferometer" technique. Accurate attitude knowledge enables V x B contributions to be subtracted from the measured potentials, and permits the measured components to be rotated into meaningful geophysical reference frames. This presentation discusses the expected performance of electric field double probe experiments and their boom mechanisms on both spinning and non-spinning satellite platforms with very low perigees. Careful selection of probe surface materials, such as titanium nitride, for the low perigee environment, as well as thermal considerations are also discussed.

Pfaff, R. F., Jr.↗

Observational and Numerical Diagnostics of Galaxy Cluster Outer Regions

Aims. We present the analysis of a local (z = 0.04 - 0.2) sample of 31 galaxy clusters with the aim of measuring the density of the X-ray emitting gas in cluster outskirts. We compare our results with numerical simulations to set constraints on the azimuthal symmetry and gas clumping in the outer regions of galaxy clusters. Methods. We exploit the large field-of-view and low instrumental background of ROSAT/PSPC to trace the density of the intracluster gas out to the virial radius. We perform a stacking of the density profiles to detect a signal beyond r(sub 200) and measure the typical density and scatter in cluster outskirts. We also compute the azimuthal scatter of the profiles with respect to the mean value to look for deviations from spherical symmetry. Finally, we compare our average density and scatter profiles with the results of numerical simulations. Results. As opposed to several recent results, we observe a steepening of the density profiles beyond approximately 0.3r(sub 500). Comparing our density profiles with simulations, we find that non-radiative runs predict too steep density profiles, whereas runs including additional physics and/or gas clumping are in better agreement with the observed gas distribution. We note a systematic difference between cool-core and non-cool core clusters beyond approximately 0.3r(sub 200), which we explain by a different distribution of the gas in the two classes. Beyond approximately r(sub 500), galaxy clusters deviate significantly from spherical symmetry, with only little differences between relaxed and disturbed systems. We find good agreement between the observed and predicted scatter profiles, but only when the 1% densest clumps are filtered out in the simulations. Conclusions. The general trend of steepening density around the virial radius indicates that the shallow density profiles found in several recent works were probably obtained along particular directions (e.g., filaments) and are not representative of the typical behavior of clusters. Comparing our results with numerical simulations, we find that non-radiative simulations fail to reproduce the gas distribution, even well outside cluster cores. Therefore, a detailed treatment of gas cooling, star formation, clumping, and AGN feedback is required to construct realistic models of cluster outer regions.

Eckert, D.↗

Dust Optical Properties Over North Africa and Arabian Peninsula Derived from the AERONET Dataset

Dust optical properties over North Africa and the Arabian Peninsula are extracted from the quality assured multi-year datasets obtained at 14 sites of the Aerosol Robotic Network (AERONET). We select the data with (a) large aerosol optical depth (AOD >= 0.4 at 440 nm) and (b) small Angstrom exponent (A(sub ext)<= 0.2) for retaining high accuracy and reducing interference of non-dust aerosols. The result indicates that the major fraction of high aerosol optical depth days are dominated by dust over these sites even though it varies depending on location and time. We have found that the annual mean and standard deviation of single scattering albedo, asymmetry parameter, real refractive index, and imaginary refractive index for Saharan and Arabian desert dust is 0.944 +/- 0.005, 0.752 +/- 0.014, 1.498 +/- 0.032, and 0.0024 +/- 0.0034 at 550 nm wavelength, respectively. Dust aerosol selected by this method is less absorbing than the previously reported values over these sites. The weaker absorption of dust from this study is consistent with the studies using remote sensing techniques from satellite. These results can help to constrain uncertainties in estimating global dust shortwave radiative forcing.

Kim, D.↗

Stagnation Point Convective Heating Correlations for Entry Into H2/He Atmospheres

New convective heating correlations have been developed over a broad range of entry conditions into H2/He atmospheres. The correlations are based on calculations of stagnation point convective heat fluxes over a large range of velocities, densities, and effective nose radii. The heat fluxes and corresponding flowfields are computed using NASA’s hypersonic computational fluid dynamics code, DPLR. Cases in this study have been computed using both a temperaturelimited, radiative equilibrium wall and a cold wall temperature boundary condition. The results for the scaled heat flux, as a function of the difference in freestream and wall enthalpies, show a clear deviation from the linear trend predicted by the method of Sutton and Graves. Therefore, updated non-linear correlations are proposed for the prediction of stagnation point convective heat fluxes. Generally, these updated correlations predict lower heat loads and therefore smaller required TPS thicknesses when compared to previous correlations for a range of test cases in the Neptune, Saturn, and Uranus atmospheres.

Alex T. Carroll↗

Agriculture, forestry, range resources

In the area of crop specie identification, it has been found that temporal data analysis, preliminary stratification, and unequal probability analysis were several of the factors that contributed to high identification accuracies. Single data set accuracies on fields of greater than 80,000 sq m (20 acres) are in the 70- to 90-percent range; however, with the use of temporal data, accuracies of 95 percent have been reported. Identification accuracy drops off significantly on areas of less than 80,000 sq m (20 acres) as does measurement accuracy. Forest stratification into coniferous and deciduous areas has been accomplished to a 90- to 95-percent accuracy level. Using multistage sampling techniques, the timber volume of a national forest district has been estimated to a confidence level and standard deviation acceptable to the Forest Service at a very favorable cost-benefit time ratio. Range specie/plant community vegetation mapping has been accomplished at various levels of success (69- to 90-percent accuracy). However, several investigators have obtained encouraging initial results in range biomass (forage production) estimation and range readiness predictions. Soil association map correction and soil association mapping in new area appear to have been proven feasible on large areas; however, testing in a complex soil area should be undertaken.

Crea, W. J.↗