Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “science & technology”

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 325 records · Page 18

Many-Body Benchmark of Electronic Charge and Spin Densities for Li 1–x NiO 2

Accurate benchmarks are particularly important for highly correlated oxides as mean-field approximations often fail to describe the subtle balance of charge transfer and magnetism in these materials with an accuracy comparable to experimental needs. Here we present accurate diffusion Monte Carlo (DMC) results of the electronic charge and spin densities for the tunable highly correlated oxide Li 1–x NiO 2 for x = 0, 1/2, and 1. To enable quantitative comparisons, we introduce a robust density-partitioning scheme, extending Voronoi analysis to assign atomic charges from spatially noisy DMC densities. We then benchmark common approximations used in density functional theory (DFT). Comparison against DMC shows that r 2 SCAN delivers the most balanced performance across charge, spin, and radial density descriptors, nearly reproducing DMC results for LiNiO 2 and apical Ni sites in Li 0.5 NiO 2 . Hybrid functionals (PBE0, SCAN0) perform unexpectedly poorly, and PBE + U + V yields inconsistent trends between charge and spin densities. Therefore, the r 2 SCAN functional minimizes errors relative to DMC while capturing the variable valence of the Ni ion and also retaining the computational efficiency of DFT for large-scale simulations of the tunable structural and electronic phases of Li1−xNiO2. Our study highlights the importance of accurate benchmarking of the fundamental quantities involved in DFT to select appropriate DFT approximations in order to advance the predictive modeling of charge-transfer-driven phenomena in correlated electron systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Spatial Mapping of Chiral-Induced Spin Selectivity in Chiral Perovskite via Spin-Schottky Junction

Chiral halide perovskite (c-HP) semiconductors exhibit on average a large chiral-induced spin selectivity (CISS) effect. Nevertheless, the microscopic details of CISS and its integration in opto-spintronic constructs remain nascent. Reliable reporting of CISS performance characteristics represents a significant challenge in providing the necessary design rules. We show a Kelvin probe force microscopy (KPFM) method that can quantitatively evaluate and spatially map the chirality-dependent surface contact potential difference resulting from the formation of a spin-Schottky junction. We revealed inhomogeneity in the CISS response, where low-CISS regions in the c-HP films reduce the overall macroscopic average, likely serving as a key factor in optimizing macroscopic performance. We also observed that although c-HP films made from higher precursor concentrations lead to thicker films and higher carrier concentrations with subsequent larger barrier heights in the Schottky junction, stronger spin relaxation due to non-ideal film quality reduces spin polarization.

14 SOLAR ENERGY↗

Formation of carbon homonuclear bonds in β-SiC under neutron irradiation at various temperatures and neutron doses

To elucidate radiation defect processes in SiC, Raman spectroscopy was systematically applied to high-purity, polycrystalline β-SiC that was neutron irradiated at a range of temperature and dose conditions. The analysis specifically focused on formation of carbon homonuclear bonds by irradiation; these bonds were indicated by D and G bands and amorphous carbon peaks. Intensity of the carbon peaks relative to SiC peaks significantly decreased in the case of high temperature and/or high neutron dose of 500 °C to 29 displacements per atom (dpa) and about 800 °C to 1.38 and 29 dpa. The absence of carbon bond peaks under those conditions was explained by growth of stoichiometric defect clusters, consistent with previous atomistic simulations on SiC defect stability. The lack of Raman bands associated with carbon clusters under high-temperature and high-dose radiation conditions accounts for the resistance of SiC to phase separation under irradiation. The findings further suggest that material compositions and chemical properties that are inherently resistant to chemical disordering under high-dose radiation conditions are indicative of the long-term durability of ceramic compounds in radiation environments.

36 MATERIALS SCIENCE↗

Applications of visualization technology in the structural sciences

The structural sciences are undergoing a transformation driven by advancements in visualization technologies that aid researchers in understanding and communicating experimental data from complex molecular systems. New applications of integrative structural biological and biophysical approaches add a wide variety of complementary information from a broad range of scientific disciplines. These approaches extend structural biophysical methodologies to enable research by the incorporation of a variety of data streams and utilization of tools like molecular graphics, virtual reality, and machine learning. To redefine how structural data—particularly from cryo-electron microscopy and x-ray crystallography—are fed forward for scientific exploration and communication, the advances in tools for data visualization and interpretation have been critical. By bringing molecular systems into an interactive three-dimensional space, these novel technologies enhance research workflows, facilitate structure-based drug design, and create engaging educational experiences. Taken together, these visualization innovations are essential tools for advancing the field by making concepts more accessible and compelling.

Eng, Edward T. [New York Structural Biology Center↗

Directed energy deposition of functionally graded V-4Cr-4Ti to Fe-9Cr transition for fusion power systems

This study proposes a graded structure via additive manufacturing for divertor and first wall blanket applications in fusion reactors. Materials were selected based on thermodynamic calculations to operate from 1100 °C at the plasma-facing level to 550 °C at the structural steel level. Conventional joining methods often lead to failures due to discrete reaction layers with significant mechanical property differences. Using laser beam-directed energy deposition (LB-DED), this study demonstrates the fabrication of a VCrTi-Gr91 steel functionally graded component through a novel process parameter optimization framework. A systematic approach included powder characterization, single-track depositions, and construction of printability maps. Near full-density specimens of each interlayer were additively manufactured, and a transition from V-based alloys to reduced activation ferritic martensitic steels was achieved. Computational material selection of interlayer alloys and thermodynamic/diffusion kinetics simulations prevented most interface incompatibilities. A brittle intermetallic formed at one interface, causing cracking, which was not predicted by current thermodynamic models. Transition alloy design approach was updated with a more recent database and a mitigation strategy has been proposed to eliminate the formation of deleterious intermetallic phases. Ultimately, LB-DED has proven effective for producing multi-material graded systems for fusion applications, with the demonstrated process parameter optimization framework applicable to various materials.

Additive manufacturing↗

A Geometric Volume of Fluid-Based Multiphase Flow Solver Extension to the Reacting Flow Solver, PeleLM

A new algorithm is presented to simulate multiphase flows with surface tension in a pathway for spray combustion simulation. The algorithm combines capabilities from two open-source packages, including the interface reconstruction library (IRL), a library of computational geometry routines to enable the volume of fluid (VOF) method, and PeleLM, a solver for the reacting Navier-Stokes equations. Additionally, surface tension is implemented using the continuum surface force (CSF) model with an improved height function technique in the volume fraction field. Spurious errors in volume fraction arising from our combined strategy are corrected through a topology-based method that improves both numerical stability and accuracy. Multiple validation simulations are conducted, including (i) translations and rotations of Zalesak's disk, (ii) a stationary circular droplet with surface tension, (iii) an oscillating elliptical droplet, and (iv) three-dimensional deformation of a spherical droplet. Results indicate that the combined scheme retains the favorable properties of each of the component algorithms.

42 ENGINEERING↗

Equilipy

There is high demand for high-throughput calculations of phase equilibria based on the CALPHAD method (CALculation of PHAse Diagram). High-throughput calculations are possible through HPC, currently by a commercial program called Thermo-Calc. However, the number of nodes/processors are limited to the number of purchased license (16 processors per license). This program provides a toolkit for high-throughput calculations of phase equilibria by the CALPHAD method. The program is prepared for Python environment, so that it can be easily installed and used together with other open-source programs. The program can be run in supercomputer.

Kwon, Sunyong [Oak Ridge National Laboratory (ORNL↗

Expanding the homogeneity range of UFe 2 : Formation of Laves phase UFe 2 +Ni

Uranium-based intermetallic Laves phases (UM 2 ) provide valuable insight into the interplay of 5f orbital hybridization and magnetic properties. However, ferromagnetism in these compounds is rare. UFe 2 is a cubic C15 Laves phase with greatest Curie temperature (T C ) of all such compounds, at 162 K. In this study, we report the synthesis and characterization of a novel U-based Laves phase, UFe 2 +Ni, effectively forcing excess Ni in to the UFe 2 structure. UFe 2 +Ni crystallizes in the cubic C15 Laves phase, appearing to deviate from the typical UM 2 formula and exhibiting a significantly extended homogeneity range. Annealing this material leads to an increased ferromagnetic transition temperature, between 35 and 155 K, and decreased lattice parameter with higher annealing temperatures. This variation is attributed to changes in Fe and Ni solubility with annealing temperature. In conclusion, our findings highlight the potential for tuning magnetic properties in U-based Laves phases through compositional and thermal modifications, expanding the understanding of magnetically ordered intermetallics.

Annealing effects↗

High-performance permanent-magnet materials based on CeFe 12

We perform first-principles calculations on the Mo/Zr-alloyed ThMn 12 -structure prototype CeFe 12 to explore its electronic and magnetic properties and potential for permanent magnet applications. A recent experiment has shown that the Mo/Zr-alloyed CeFe 12 phase with the symmetry 𝐼⁢4/𝑚⁢𝑚⁢𝑚 can be stabilized in bulk. It encourages an investigation into the intrinsic magnetic characteristics of the base structure CeFe 12 and the effect of transition metal alloying in it. Our calculations reveal that the base structure prototype CeFe 12 in its tetragonal phase exhibits uniaxial magnetic anisotropy with a magnetic anisotropy energy (𝐾 1 ) ∼ 1 MJ/m −3 and a substantial magnetic moment of ∼ 1.8T, showing its promise as a permanent magnet. The introduction of the stabilizing transition metals Zr, Mo, and W reduces these values but still keeps the materials very promising. Notably, when 50% of Ce is substituted by Sm in Mo-alloyed CeFe 12 , the 𝐾 1 is significantly enhanced, accompanied by a fairly decent magnetic moment of ∼ 1.2T. These findings establish this alloy as a strong candidate for high-performance permanent magnet applications. In conclusion, this study emphasizes the importance of strategic element substitution and site selection in transition metal-alloyed CeFe 12 to achieve both structural stability and remarkable magnetic properties.

Density of states↗

Screening of Coulombic Interactions To Achieve a Higher Power Factor in Conjugated Polymers

Thermoelectric properties of conducting polymers typically suffer from molecular chain disordering, as charge transport is predominantly controlled by morphology. This is especially more problematic when counterions are introduced to tune the carrier concentration for optimal thermoelectric performance, which disturbs the morphology further. Here, in this work, we introduce a new avenue for enhancing thermoelectric properties without needing to regulate the morphology, namely, by controlling the coulombic interaction between polarons and counterions. We perform in situ de-doping thermoelectric experiments over 3 orders of magnitude change in electrical conductivity of three distinct thermoelectric polymers, namely, poly(3-hexylthiophene-2,5-diyl) (P3HT), poly[2,5-bis(3-dodecylthiophen-2-yl)thieno[3,2-b]thiophene] (PBTTT-C 12 ), and poly[2,5-(2-octyldodecyl)-3,6-diketopyrrolopyrrole-alt-5,5-(2,5-di(thien-2-yl)thieno[3,2-b]thiophene)] (OD-PDPP2T-TT) conjugated polymers, followed by grazing-incidence wide-angle X-ray scattering (GIWAXS) to study their respective morphologies. We demonstrate a 9-fold enhancement in the thermoelectric power factor in OD-PDPP2T-TT compared to PBTTT-C 12 and link it to the coulombic screening of charge carriers, including in the optimally doped regime. We support this hypothesis using Boltzmann transport equations and show that, in both P3HT and PBTTT-C 12 , as the polymer is doped, impurity scattering remains the dominant scattering mechanism, while in OD-PDPP2T-TT, the scattering mechanism changes from impurity to acoustic phonon limited, resulting in more effective screening of ionized counterions. Our results provide an additional knob to enhance the fundamental understanding of thermoelectric physics of conducting polymers and provide a pathway to achieve higher performance in the field of organic thermoelectrics.

36 MATERIALS SCIENCE↗

Conceptual study on using Doppler backscattering to measure magnetic pitch angle in tokamak plasmas

We introduce a new approach to measure the magnetic pitch angle profile in tokamak plasmas with Doppler backscattering (DBS), a technique traditionally used for measuring flows and density fluctuations. The DBS signal is maximised when its probe beam's wavevector is perpendicular to the magnetic field at the cutoff location, independent of the density fluctuations [Hillesheim \emph{et al} 2022 \emph{Nucl. Fusion} \textbf{55} 073024]. Hence, if one could isolate this effect, DBS would then yield information about the magnetic pitch angle. By varying the toroidal launch angle, the DBS beam reaches cutoff with different angles with respect to the magnetic field, but with other properties remaining similar. Hence, the toroidal launch angle which gives maximum backscattered power is thus that which is matched to the pitch angle at the cutoff location, enabling inference of the magnetic pitch angle. We performed systematic scans of the DBS toroidal launch angle for repeated DIII-D tokamak discharges. Experimental DBS data from this scan were analysed and combined with Gaussian beam-tracing simulations using the Scotty code [Hall-Chen \emph{et al} 2022 \emph{Plasma Phys. Control. Fusion} \textbf{64} 095002]. The pitch-angle inferred from DBS is consistent with that from magnetics-only and motional-Stark-effect-constrained (MSE) equilibrium reconstruction in the edge. In the core, the pitch angles from DBS and magnetics-only reconstructions differ by one to two degrees, while simultaneous MSE measurements were not available. The uncertainty in these measurements was under a degree; we show that this uncertainty is primarily due to the error in toroidal steering, the number of toroidally separated measurements, and shot-to-shot repeatability. We find that the error of pitch-angle measurements can be reduced by optimising the poloidal launch angle and initial beam properties. Since DBS has high spatial and temporal resolutions, is non-perturbative, does not require neutral beams, and is likely robust to neutron damage of and debris on the first mirrors, using DBS to measure the pitch angle in future fusion energy systems is especially appealing.

beam tracing↗

Geometrical optics without singularities: using the ray time as the coordinate space

Geometrical optics (GO) is widely used for reduced modelling of waves in plasmas, but it fails near reflection points, where it predicts a spurious singularity of the wave amplitude. We show how to avoid this singularity by adopting a different representation of the wave equation. Instead of the physical coordinate 𝑥 and the wavevector 𝑘, we use the ray time 𝜏 as the new canonical coordinate and the ray energy ℎ as the associated canonical momentum. To derive the envelope equation in the 𝜏-representation, we construct the Weyl symbol calculus on the (𝜏,ℎ) space and show that the corresponding Weyl symbols are related to their (𝑥,𝑘) counterparts by the Airy transform. This allows us to express the coefficients in the envelope equation through the known properties of the original dispersion operator. When necessary, solutions of this equation can be mapped to the 𝑥-space using a generalised metaplectic transform. However, the field per se might not even be needed in practice. Instead, knowing the corresponding Wigner function usually suffices for linear and quasilinear calculations. As a Weyl symbol itself, the Wigner function can be mapped analytically, using the aforementioned Airy transform. We show that the standard Airy patterns that form in regions where conventional GO fails are successfully reproduced within metaplectic GO (MGO) simply by remapping the field from the 𝜏-space to the 𝑥-space. An extension to mode-converting waves is also presented. This formulation, which we call generalised MGO, can be particularly useful, for example, for reduced modelling of the O–X conversion in inhomogeneous plasma near the critical density, an effect that is important for fusion applications and also occurs in the ionosphere. Overall, MGO can replace GO for any practical purposes, because it better handles cutoffs and is similar otherwise.

plasma waves↗

Aerosol effective radiative forcing increases Earth’s energy imbalance in recent decades

Earth’s energy imbalance (EEI), a key driver of climate change, has risen markedly over the last two decades. Greenhouse gas forcing, aerosol forcing, and cloud feedback all contribute to this increase. However, the role of aerosol forcing, particularly effective radiative forcing through aerosol-cloud interactions (ERF ACI ), remains highly uncertain and closely intertwined with cloud feedback. Here, we estimate shortwave ERF ACI using satellite observations and show its importance as a contributor to the EEI increase between 2003 and 2023. ERF ACI exhibits a significant increasing trend of 0.33 Wm −2 /decade averaged over oceans between 60°S and 60°N. The increasing trend of ERF ACI stems from a global decline in cloud droplet number concentration driven by decreasing anthropogenic aerosol emissions. It is of similar magnitude as the combined instantaneous forcing trend from greenhouse gases and aerosol-radiation interactions. Our results may reduce the gap between the simulated and observed shortwave contribution to EEI trends while implying a weak shortwave cloud feedback. They also suggest a stronger aerosol forcing than the multi-model mean.

Yuan, Tianle [NASA Goddard Space Flight Center (GS↗

Application of Accelerator Technology to Quantum Information Science

The intersection of accelerator and quantum information science (QIS) offers a unique platform to advance both fields through shared technology and infrastructure. This talk will discuss the synergies which exist between these two vastly different but complementary domains. We demonstrate how we leverage pre-existing infrastructure and knowledge to perform research and development which helps to realize dramatic improvements in both 10 km long accelerators and 10 cm large quantum processors. We will explore niobium superconducting radio-frequency (SRF) cavities, a highly advanced technology that excels in efficiently storing electromagnetic energy, enabling ultra-long photon lifetimes critical for quantum processors and facilitating the characterization of quantum materials with parts-per-billion precision. We will also discuss how advancements in superconducting materials, cryogenic systems, and control techniques help to reduce cost and improve performance for both quantum systems and particle accelerators. Moreover, we will discuss cross-disciplinary applications such as dark-matter searches and demonstrate the convergence of these fields in addressing fundamental scientific questions.

Bafia, Daniel [Fermilab]↗

LANL: Missions and Student Opportunities [Slides]

At Los Alamos, we deliver science and technology to protect our nation and promote world stability. Our mission began by applying science and technology to address an international crisis. Today, we are responsible for a wide range of programs focused on national security that rely on our strong science and engineering capabilities. We offer unparalleled career opportunities in science, engineering, manufacturing, business, and more.

99 GENERAL AND MISCELLANEOUS↗

Elimination of NO Pollutant in Semi-enclosed Spaces over Sodium-promoted Cobalt Oxyhydroxide (CoOOH) by Oxidation and Adsorption Mechanism

Low concentrations of NOx accumulated in air in semi-enclosed spaces is a severe health hazard, but its efficient removal at low cost is a grand challenge. A simple method in preparing CoOOH for low-concentration NO removal at room temperature is described here. The samples are synthesized by oxidizing Co(OH)2 to CoOOH with different oxidants (Na2S2O8, O2 and H2O2). Among them, CoOOH prepared by Na2S2O8 (Co-NS) shows remarkable performance for NO abatement, while the samples synthesized by and H2O2 (Co-HO) and O2 (Co-OG) have no apparent NO removal capacity. It is found that the Co-NS material contains a higher concentration of residual Na, which contributes to the decreased crystallinity and increased concentrations of surface active species (Co2+, active oxygen) for NO storage. Following NO storage, the material can be readily regenerated by water washing. The authors from East China University of Science and Technology acknowledge support from National Basic Research Program of China (2013CB933200), National Natural Science Foundation of China (21577035, 21577034), Commission of Science and Technology of Shanghai Municipality (15DZ1205305) and 111 Project (B08021). Aiyong Wang acknowledges the China Scholarship Council for the Joint-Training Scholarship Program with the Pacific Northwest National Laboratory (PNNL). PNNL is operated by Battelle for the US DOE under contract DE-AC05-76RL01830. FG is supported by the U.S. DOE/Office of Energy Efficiency and Renewable Energy, Vehicle Technologies Office.

NOx abatement, NOx adsorber, ambient temperature, ↗

CCN data from TAMU TRACER campaign in the Houston TX region from July to September 2022

During TRACER, the Texas A&M Rapid Onsite Atmospheric Measurements Van (ROAM-V) was deployed to capture airmasses behind (maritime) and ahead (continental) of the passage of the sea-breeze front through Houston. On select sampling days, ROAM-V sampled in the morning/mid-day on the coast and then transited to a second inland site for the afternoon/evening. The suite of instruments deployed on ROAM-V included a Condensation Particle Counter (CPC; GRIMM Model 5.403 CPC), Scanning Mobility Particle Sizer (SMPS; TSI 3750 detector, TSI 3082 classifier, TSI 3088 neutralizer, TSI 3081A Differential Mobility Analyzer), Cloud Condensation Nuclei counter (Droplet Measurement Technologies CCN Counter), micro pulse lidar (Droplet Measurement Technologies Micro Pulse LiDAR (miniMPL)), and a Davis Rotating Uniform size-cut Monitor (DRUM; DRUMAir 4-DRUM). Before sampling at each location, the latitude and longitude were recorded using the GPS on the phone application “My Altitude”. Onboard the ROAM-V, aerosol samples are drawn through a shared isokinetic inlet at a flow rate ranging from 3.5 to 7.0 LPM. A portion of this flow is directed through a cyclone impactor (Brechtel, Inc. Model SCC 0.732) and 0.5 LPM is directed to the CCN. To calculate particle losses, we used a two-step method. First, the measured SMPS size distributions were used to calculate particle loss through the inlet during sampling. Second, the corrected SMPS data was used to calculate the average of the total losses per scan down the CCN line. Then, the correction was applied to the CCN data. This calculation was done separately for each deployment location due to changes in the measured size distributions between locations. Particle loss from diffusion (based on Kesten, 1991 and Gormley, 1949), inertial impaction in 90-degree bends (based on Aerosol Measurement, 2011 and Crane, 1977), and cyclone impactor efficiency (based on Dirgo, 1985) were included in the loss calculation. When the SMPS was not sampling at a location (in the case of an instrument malfunction or operator error), the reported CPC data was corrected with an average of the total losses for the entire campaign at the specified deployment location (e.g., if we needed to correct Galveston data, then the average of all calculated losses at Galveston was taken). These flatline corrections were used for all data on 22/07/13, 22/07/20, 22/07/22, and the data from Galveston on 22/08/09. The supersaturation uncertainty is estimated conservatively at +/- 0.03%, where variation in the inlet temperature, pressure, and calibration technique prevents a more accurate measurement. Confidence in the reported supersaturation measurements is based on a pre-campaign calibration (following the methods from our previous work and Deng, 2014 based on Rose, 2008) in addition to inter-comparisons with the DOE for two days (22/08/18 and 22/09/01) where TAMU was co-located with AMF1. The inter-comparisons show good agreement between our instrument and the DOEs instrument on both days at all supersaturations. After the last inter-comparison on 22/09/01, there was no indication of a malfunction by our instrument through the rest of the campaign. Unfortunately, the instrument was dropped during demobilization. A post-campaign calibration was conducted, which showed a substantial departure from the pre-campaign calibration. The drop may have damaged the instrument’s ability to produce the desired supersaturations. Therefore, we do consider the data after 22/09/01 to be correct, but it should be used with caution. The CCN counter sampled for 3 minutes at each supersaturation setpoint (0.2, 0.4, 0.6, 0.8, 1.0, and 1.2%). At the end of a cycle, the instrument was set to 0.01% supersaturation for 5 minutes. The data is comprised of the last 60 seconds of each supersaturation set point (0.2, 0.4, 0.6, 0.8, 1.0, and 1.2%) to ensure the instrument stabilized and was able to reach thermal equilibrium. We removed the data during the periods where there were operational difficulties, setup, or maintenance. This data was collected for ARM Field Campaign AFC07055 and supported by DOE ASR grant DE-SC0021047. For any further questions, please feel free to contact the instrument PI, Sarah D. Brooks, sbrooks@tamu.edu. Rose et. al. Calibration and Measurement Uncertainties of a Continuous-Flow Cloud Condensation Nuclei Counter (DMT-CCNC): CCN Activation of Ammonium Sulfate and Sodium Chloride Aerosol Particles in Theory and Experiment. Atmos. Chem. Phys., 8, 1153-1179, 2008. Deng et. al. Using Raman Microspectroscopy to Determine Chemical Composition and Mixing State of Airborne Marine Aerosols over the Pacific Ocean. Aerosol Science and Technology, Vol 48, Issue 2, 2014. Kesten et. al. Calibration of a TSI Model 3025 Ultrafine Condensation Particle Counter. Aerosol Science and Technology, 15:2, 107-111, 1991. Gormley et. al. Diffusion from a Stream Flowing through a Cylindrical Tube. Proceedings of the Royal Irish Academy, Vol 52, 163-169, 1948. Aerosol Measurement: Principles, Techniques, and Applications, Third Edition. John Wiley & Sons, Inc, 2011. Crane et. al. Inertial Deposition of Particles in a Bent Pipe. Journal of Aerosol Science, Vol 8, 161-170, 1977. Dirgo et. al. Cyclone Collection Efficiency: Comparison of Experimental Results with Theoretical Predictions. Aerosol Science and Technology, 4:4, 401-415, 1985.

54 ENVIRONMENTAL SCIENCES↗