Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Conditional diffusion models”

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 1,495 records · Page 83
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Q2 Report for FY25 Theory and Simulation Performance Target: Development of an integrated modeling framework for fusion reactor design and assessment

This report describes the work and activities carried out towards the completion of each of the following milestones in FY25 Q2: 1. Demonstrate workflow for generating self-consistent CESOL plasma profiles + first wall and divertor loading prediction and generate the CAT plasma and neutron loading needed for further engineering analysis: $\circ$ Run CESOL with BOUT++/Hermes-3 and immersed boundary condition to directly map to wall: • Run BOUT++/Hermes-3 through the IPS workflow to find radial particle and energy diffusivities to match either the Eich or the physics-based scaling of the SOL heat flux width, and • Expand source of first wall heat flux to include charged particles, neutrals, and radiation from the core+edge. 2. Generate medium fidelity parametrized CAD: $\circ$ Develop the TRACER tool to read an existing CAD, regenerate the geometry based on vertex location and connectivity information, define vertex translation and parameters needed for scaling the CAD, and $\circ$ Utilize the FreeGS code to determine CAT PF coil placement, including minimizing the number of coils, coil current, and electromechanical stresses. 3. Utilize plasma loading for engineering analysis: $\circ$ Couple the plasma loading to input for OpenFOAM and demonstrate initial test of thermal analysis of CAT first wall loading with typical DCLL blanket component cooling boundary conditions. 4. Demonstrate nuclear analysis: $\circ$ Apply initial analysis of tritium transport in DCLL blanket by evaluating spatially resolved tritium generation rates, tritium diffusion and convection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2007

Multi-Scale Modeling of Magnetospheric Reconnection

One of the major challenges in modeling the magnetospheric magnetic reconnection is to quantify the interaction between large-scale global magnetospheric dynamics and microphysical processes in diffusion regions near reconnection sites. There is still considerable debate as to what degree microphysical processes on kinetic scales affect the global evolution and how important it is to substitute numerical dissipation and/or ad hoc anomalous resistivity by a physically motivated model of dissipation. Comparative studies of magnetic reconnection in small scale geometries demonstrated that MHD simulations that included non-ideal processes in terms of a resistive term $\eta J$ did not produce the fast reconnection rates observed in kinetic simulations. For a broad range of physical parameters in collisionless magnetospheric plasma, the primary mechanism controlling the dissipation in the vicinity of the reconnection site is non-gyrotropic effects with spatial scales comparable with the particle Larmor radius. We utilize the global MHD code BATSRUS and incorporate nongyrotropic effects in diffusion regions in terms of corrections to the induction equation. We developed an algorithm to search for magnetotail reconnection sites, specifically where the magnetic field components perpendicular to the local current direction approaches zero and form an X-type configuration. Spatial scales of the diffusion region and magnitude of the reconnection electric field are calculated selfconsistently using MHD plasma and field parameters in the vicinity of the reconnection site. The location of the reconnection sites is updated during the simulations. To clarify the role of nongyrotropic effects in diffusion region on the global magnetospheric dynamic we perform simulations with steady southward IMF driving of the magnetosphere. Ideal MHD simulations with magnetic reconnection supported by numerical resistivity produce steady configuration with almost stationary near-earth neutral line (NENL). Simulations with non-gyrotropic corrections demonstrate dynamic quasi-periodic response to the steady driving condition. The loading/unloading cycle in non-gyrotropic MHD results has a non-stationary reconnection site in the magnetotail, with the retreating during the stretching phase and then a new NENL forming in the resulting thin plasma sheet. We expect that this model will lead to improved representations of space weather event in the magnetosphere.

Kuznetsova, M. M.↗
DOE OSTI2025

Corrosion of 316 L stainless steel under the natural circulation of molten NaCl-MgCl 2 salt

The corrosion behavior of 316 L stainless steel (SS) was studied via the natural circulation of molten eutectic NaCl-MgCl 2 salt through a microloop. The post-corrosion tested 316 L SS microloop sections were characterized with microscopy techniques to determine the microstructural and microchemical changes that occurred at the alloy/salt interface. It was found that 316 L SS showed heterogeneous dissolution at the hot-leg, whereas deposition of corrosion products occurred at the cold-leg. For the first time, experimentally obtained molten salt flow-induced corrosion of 316 L SS results were combined with computational thermodynamic-kinetic models to validate the dissolution and deposition in terms of elemental compositional changes at the alloy/salt interface. The thermodynamic-kinetic modeling predicted that the heterogeneous dissolution of Cr and Fe from the hot-leg section of 316 L SS persisted throughout the salt circulation. The model also estimated that, despite Cr deposition starting earlier than Fe, the total redeposition of Fe is expected to be significantly greater than that of Cr over the circulation of salt. Furthermore, the modeling accurately predicted the subsurface enrichment of Mo which is attributed to the reduced Cr activity and the relatively higher diffusion rate of Mo within the alloy matrix. Here, the agreement between modeling and experimental results confirms that Fe chlorides dissolve at the hot-leg and subsequently deposit at the cold-leg due to activity changes driven by the thermal gradient. By contrast, Cr was not detected in the cold-leg deposits, which is attributed to its weaker temperature dependence on activity, limiting its redeposition under these conditions.

36 - MATERIALS SCIENCE↗
DOE OSTI2019

On Thermodynamic and Kinetic Mechanisms for Stabilizing Surface Solid Solutions

Many processes for energy storage rely on transformations between phases with strong separation tendencies. In these systems, performance limitations can arise from undesirable chemical and mechanical factors associated with the phase separation behavior. Solid solutions represent a desirable alternative, provided the conditions for their formation are known. In the present work, we invoke linear stability theory and diffuse-interface mesoscopic simulations to demonstrate that solid solutions can be stabilized near surface layers of phase-separating systems. Two factors are found to drive surface solid-solution formation: surface relaxation of solution self-strain energy and anisotropy of diffusion mobility. Using a strongly phase-separating LiXFePO4 particle as a model system, we show that the relaxation of the solution self-strain energy competes against the relaxation of the coherency strain energy to stabilize surface solid solutions. Our theoretical understanding also suggests that highly anisotropic diffusion mobility can provide an alternative kinetic route to achieve the same aim, with stabilizing behavior strongly dependent on the specific alignment of the surface orientation. Our findings provide fundamental guidance for manipulating solid-solution behavior in nanoscale structures, in which surface effects become especially significant. Beyond energy storage materials, our findings have important implications for understanding solid-solution formation in other phase-separating systems from metal alloys to ceramics.

25 ENERGY STORAGE↗
DOE OSTI2018

Improving accuracy of electrochemical capacitance and solvation energetics in first-principles calculations

Here, reliable first-principles calculations of electrochemical processes require accurate prediction of the interfacial capacitance, a challenge for current computationally efficient continuum solvation methodologies. We develop a model for the double layer of a metallic electrode that reproduces the features of the experimental capacitance of Ag(100) in a non-adsorbing, aqueous electrolyte, including a broad hump in the capacitance near the potential of zero charge and a dip in the capacitance under conditions of low ionic strength. Using this model, we identify the necessary characteristics of a solvation model suitable for first-principles electrochemistry of metal surfaces in non-adsorbing, aqueous electrolytes: dielectric and ionic nonlinearity, and a dielectric-only region at the interface. The dielectric nonlinearity, caused by the saturation of dipole rotational response in water, creates the capacitance hump, while ionic nonlinearity, caused by the compactness of the diffuse layer, generates the capacitance dip seen at low ionic strength. We show that none of the previously developed solvation models simultaneously meet all these criteria. We design the nonlinear electrochemical soft-sphere solvation model which both captures the capacitance features observed experimentally and serves as a general-purpose continuum solvation model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2012

AB INITIO SIMULATIONS FOR MATERIAL PROPERTIES ALONG THE JUPITER ADIABAT

We determine basic thermodynamic and transport properties of hydrogen-helium-water mixtures for the extreme conditions along Jupiter's adiabat via ab initio simulations, which are compiled in an accurate and consistent data set. In particular, we calculate the electrical and thermal conductivity, the shear and longitudinal viscosity, and diffusion coefficients of the nuclei. We present results for associated quantities like the magnetic and thermal diffusivity and the kinematic shear viscosity along an adiabat that is taken from a state-of-the-art interior structure model. Furthermore, the heat capacities, the thermal expansion coefficient, the isothermal compressibility, the Grueneisen parameter, and the speed of sound are calculated. We find that the onset of dissociation and ionization of hydrogen at about 0.9 Jupiter radii marks a region where the material properties change drastically. In the deep interior, where the electrons are degenerate, many of the material properties remain relatively constant. Our ab initio data will serve as a robust foundation for applications that require accurate knowledge of the material properties in Jupiter's interior, e.g., models for the dynamo generation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Electrode Mesoscale as a Collection of Particles: Coupled Electrochemical and Mechanical Analysis of NMC Cathodes

Battery electrodes are composed of polydisperse particles and a porous, composite binder domain. These materials are arranged into a complex mesostructure whose morphology impacts both electrochemical performance and mechanical response. We present image-based, particle-resolved, mesoscale finite element model simulations of coupled electrochemical-mechanical performance on a representative NMC electrode domain. Beyond predicting macroscale quantities such as half-cell voltage and evolving electrical conductivity, studying behaviors on a per-particle and per-surface basis enables performance and material design insights previously unachievable. Voltage losses are primarily attributable to a complex interplay between interfacial charge transfer kinetics, lithium diffusion, and, locally, electrical conductivity. Mesoscale heterogeneities arise from particle polydispersity and lead to material underutilization at high current densities. Particle-particle contacts, however, reduce heterogeneities by enabling lithium diffusion between connected particle groups. While the porous composite binder domain (CBD) may have slower ionic transport and less available area for electrochemical reactions, its high electrical conductivity makes it the preferred reaction site late in electrode discharge. Mesoscale results are favorably compared to both experimental data and macrohomogeneous models. This work enables improvements in materials design by providing a tool for optimization of particle sizes, CBD morphology, and manufacturing conditions.

42 ENGINEERING↗
NASA NTRS1979

Mathematical model investigation of long-term transport of ocean-dumped sewage sludge related to remote sensing

An existing, three-dimensional, Eulerian-Lagrangian finite-difference model was modified and used to examine the transport processes of dumped sewage sludge in the New York Bight. Both in situ and laboratory data were utilized in an attempt to approximate model inputs such as mean current speed, horizontal diffusion coefficients, particle size distributions, and specific gravities. The results presented are a quantitative description of the fate of a negatively buoyant sewage sludge plume resulting from continuous and instantaneous barge releases. Concentrations of the sludge near the surface were compared qualitatively with those remotely sensed. Laboratory study was performed to investigate the behavior of sewage sludge dumping in various ambient density conditions.

Kuo, C. Y.↗
NASA NTRS1993

Comment on 'Ionospheric evidence of hot oxygen in the upper atmosphere of Venus'

The conclusion of Mahajan et al. (1992) that 'the existence of O(+) as dominant at (Venusian) ionopause altitudes in excess of 500-1000 km can only be explained if atomic oxygen is the major neutral constituent' is argued to be incorrect. It is suggested that at a transition region of about 200 km, thermal atomic oxygen is the dominant neutral gas, and hot oxygen is a minor species; thus the O(+) to H(+) ratio at high altitudes is not an indicator of the presence of hot oxygen at these altitudes. A 1D model for H(+) and O(+) appropriate for the dayside ionosphere of Venus shows that within hot atomic oxygen density values from 1000 to 10 exp 6/ cu cm at 150 km, the calculated H(+) and O(+) densities did not change in any meaningful way, because the hot oxygen population remained a minor neutral constituent below 200 km, which is the approximate height of the transition between chemical and diffusive equilibrium conditions for the ions.

Nagy, Andrew F.↗
NASA NTRS1994

Gamma-ray constraints on supernova nucleosynthesis

Gamma-ray spectroscopy holds great promise for probing nucleosynthesis in individual supernova explosions via short-lived radioactivity, and for measuring current global Galactic supernova nucleosynthesis with longer-lived radioactivity. It was somewhat surprising that the former case was realized first for a Type II supernova, when both Co-56 and Co-57 were detected in SN 1987A. These provide unprecedented constraints on models of Type II explosions and nucleosynthesis. Live Al-26 in the Galaxy might come from Type II supernovae, and if it is eventually shown to be so, can constrain massive star evolution, supernova nucleosynthesis, and the Galactic Type II supernova rate. Type Ia supernovae, thought to be thermonuclear explosions, have not yet been detected in gamma-rays. This is somewhat surprising given current models and recent Co-56 detection attempts. Ultimately, gamma-ray measurements can confirm their thermonuclear nature, probe the nuclear burning conditions, and help evaluate their contributions to Galactic nucleosynthesis. Type Ib/c supernovae are poorly understood. Whether they are core collapse or thermonuclear events might be ultimately settled by gamma-ray observations. Depending on details of the nuclear processing, any of these supernova types might contribute to a detectable diffuse glow of Fe-60 gamma-ray lines. Previous attempts at detection have come very close to expected emission levels. Remnants of any type of age less that a few centuries might be detectable as individual spots of Ti-44 gamma-ray line emission. It is in fact quite surprising that previous surveys have not discovered such spots, and the constraints on the combination of nucleosynthesis yields and supernova rates are very interesting. All of these interesting limits and possibilities mean that the next mission, International Gamma-Ray Astrophysics Laboratory (INTEGRAL), if it has sufficient sensitivity, is very likely to lead to the realization of much of the great potential of gamma-ray spectroscopy for understanding supernovae.

Leising, Mark D.↗
NASA NTRS1999

Determination of Creep Behavior of Thermal Barrier Coatings Under Laser Imposed High Thermal and Stress Gradient Conditions

A laser sintering/creep technique has been established to determine the creep behavior of thermal barrier coatings under steady-state high heat flux conditions. For a plasma sprayed zirconia-8 wt. % yttria coating, a significant primary creep strain and a low apparent creep activation energy were observed. Possible creep mechanisms involved include stress induced mechanical sliding and temperature and stress enhanced cation diffusion through the splat and grain boundaries. The elastic modulus evolution, stress response, and total accumulated creep strain variation across the ceramic coating are simulated using a finite difference approach. The modeled creep response is consistent with experimental observations.

Zhu, Dongming↗
DOE OSTI2023

Kernel Buffer Volume Fraction Margin of the AGR Designed Fuel Particle

Modeling results used to assess the fuel performance of the TRISO-coated fuel particles as a function of kernel/buffer volume fraction include SiC tangential stress, formation of the buffer/IPyC gap, particle temperature profile, internal particle pressure, fission gas released from the kernel, probability of fuel particle failure, and fission product diffusion. These results were evaluated at two burnup levels and irradiation temperatures to bound expected steady-state irradiation conditions. In general, increasing the kernel/buffer volume fraction increases the SiC stress and subsequently the failure probability of a fuel particle when compared to the AGR designed particle. There was little impact on the fission product diffusion through the particle as the kernel/buffer volume fraction increased.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Neural Network Enhanced RKPM for Electrochemical-Mechanical Coupled Damage Modeling of Energy Storage Materials

Energy storage materials undergo significant charge cycling, which makes understanding their reliability and durability fundamental in predicting performance and service life. Strong electrochemical-mechanical coupling and highly anisotropic material properties contribute to the formation and propagation of micro-cracking, largely along material interfaces and grain boundaries. For Li-ion batteries, for example, lithium moving between electrodes during charging and discharging process causes expansion and contraction of grains, and the strongly anisotropic and nonlinearly [Li]-dependent grain material properties can cause grains to expand into and contract away from each other, leading to chemo-mechanical cracking. In the first part of this work, a RKPM based computational framework for solving the coupled solid-phase lithium conservation with Fickian diffusion and the lithium concentration dependent anisotropic mechanical problem subjected to a highly nonlinear Butler-Volmer boundary condition is introduced. The choice of RKPM completeness conditions for lithium concentration and mechanical deformation fields, and the variational consistency condition for the domain integration of the coupled problem is first determined. In the second part of this work, a neural network-enhanced reproducing kernel particle method (NN-RKPM) [1] is leveraged to accurately capture damage and crack propagation throughout the material, by learning the location, orientation, and sharpness of discontinuity while allowing for a coarser nodal distribution than that is necessary for capturing sharp solution transitions using traditional mesh-based methods. NN-RKPM is used to inform how crack opening and closure in turn affect the coupled chemical equations and material microstructure.

damage modeling↗
NASA NTRS1995

The neutral oxygen spectrum. 2: Pumping by hydrogen Lyman-beta under the optically thin condition: A first application to the classical novae

A calculation, employing a detailed model of neutral oxygen, is carried out to give fluorescent line intensities expected in a long-proposed photoexcitation by accidental resonance (PAR) process in which hydrogen Lyman-beta photoexcites the oxygen spectrum. The results pertain to the optically thin case but provide an upper limit to the fluorescent intensities which can be attained. They are applied to analyze line ratios involving the strong 8446 A line observed in classical novae during the diffusion-enhanced and Orion phases. Operation of the PAR process in the novae is verified. It is found that photoexcitation rates in the ejecta reach values greater than 0.1/sec, corresponding to hydrogen Lyman-beta radiation field intensities greater than 1250 ergs/cm/sec/sr.

Kastner, S. O.↗
DOE OSTI2024

Simulations of flaming combustion and flaming-to-smoldering transition in wildland fire spread at flame scale

Our objective in the present study is to provide basic insights into the coupling between external-gas and solid biomass vegetation processes that control the dynamics of flame spread in wildland fire problems. Here we focus on a modeling approach that resolves processes occurring at vegetation and flame scales, i.e., the formation of flammable vapors due to the thermal degradation of the solid biomass, the subsequent combustion in ambient air, the thermal feedback to the biomass through radiative and convective heat transfer, and the possible transition from flaming combustion (taking place outside of the solid biomass) to smoldering combustion (taking place inside the solid biomass). The capability uses a multiphase combustion framework and treats external-gas processes through a Large Eddy Simulation solver and solid biomass processes through a discrete particle model. The discrete particle model adopts a one-dimensional porous medium formulation, includes descriptions of drying, thermal pyrolysis, oxidative pyrolysis, and char oxidation, as well as a description of the external-gas-to-solid-biomass diffusion of oxygen mass; the discrete particle model thereby provides a treatment of in-depth oxidative processes and allows the simulation of smoldering combustion. The modeling capability is applied to the simulation of fire spread across a surrogate biomass vegetation bed corresponding to a discrete array of cylindrical-shaped, vertically-oriented, pine wood sticks, characterized by a monomodal size distribution, in horizontal flat terrain and under wind-aided conditions. The numerical results demonstrate that the model can simulate successful flaming-to-smoldering transition followed by complete biomass consumption.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2014

Oxygen distribution in tumors: A qualitative analysis and modeling study providing a novel Monte Carlo approach

Purpose: To construct a Monte Carlo (MC)-based simulation model for analyzing the dependence of tumor oxygen distribution on different variables related to tumor vasculature [blood velocity, vessel-to-vessel proximity (vessel proximity), and inflowing oxygen partial pressure (pO{sub 2})]. Methods: A voxel-based tissue model containing parallel capillaries with square cross-sections (sides of 10 μm) was constructed. Green's function was used for diffusion calculations and Michaelis-Menten's kinetics to manage oxygen consumption. The model was tuned to approximately reproduce the oxygenational status of a renal carcinoma; the depth oxygenation curves (DOC) were fitted with an analytical expression to facilitate rapid MC simulations of tumor oxygen distribution. DOCs were simulated with three variables at three settings each (blood velocity, vessel proximity, and inflowing pO{sub 2}), which resulted in 27 combinations of conditions. To create a model that simulated variable oxygen distributions, the oxygen tension at a specific point was randomly sampled with trilinear interpolation in the dataset from the first simulation. Six correlations between blood velocity, vessel proximity, and inflowing pO{sub 2} were hypothesized. Variable models with correlated parameters were compared to each other and to a nonvariable, DOC-based model to evaluate the differences in simulated oxygen distributions and tumor radiosensitivities for different tumor sizes. Results: For tumors with radii ranging from 5 to 30 mm, the nonvariable DOC model tended to generate normal or log-normal oxygen distributions, with a cut-off at zero. The pO{sub 2} distributions simulated with the six-variable DOC models were quite different from the distributions generated with the nonvariable DOC model; in the former case the variable models simulated oxygen distributions that were more similar to in vivo results found in the literature. For larger tumors, the oxygen distributions became truncated in the lower end, due to anoxia, but smaller tumors showed undisturbed oxygen distributions. The six different models with correlated parameters generated three classes of oxygen distributions. The first was a hypothetical, negative covariance between vessel proximity and pO{sub 2} (VPO-C scenario); the second was a hypothetical positive covariance between vessel proximity and pO{sub 2} (VPO+C scenario); and the third was the hypothesis of no correlation between vessel proximity and pO{sub 2} (UP scenario). The VPO-C scenario produced a distinctly different oxygen distribution than the two other scenarios. The shape of the VPO-C scenario was similar to that of the nonvariable DOC model, and the larger the tumor, the greater the similarity between the two models. For all simulations, the mean oxygen tension decreased and the hypoxic fraction increased with tumor size. The absorbed dose required for definitive tumor control was highest for the VPO+C scenario, followed by the UP and VPO-C scenarios. Conclusions: A novel MC algorithm was presented which simulated oxygen distributions and radiation response for various biological parameter values. The analysis showed that the VPO-C scenario generated a clearly different oxygen distribution from the VPO+C scenario; the former exhibited a lower hypoxic fraction and higher radiosensitivity. In future studies, this modeling approach might be valuable for qualitative analyses of factors that affect oxygen distribution as well as analyses of specific experimental and clinical situations.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Model-informed, Adaptive Physical Vapor Deposition to Fabricate Hierarchical Binary-alloy Thin-films

Designing next generation thin films, tailor-made for specific applications, relies on the availability of robust processing-structure-property relationships. Traditional structure zone diagrams are limited to low-dimensional mappings, with machine-learning methods only recently attempting to relate multiple processing parameters to the final microstructure. Despite this progress, structure-processing relationships are unknown for processing conditions that vary during thin-film deposition, limiting the range of microstructures and properties achievable. In this project, we employed a phase-field computational model combined with a genetic algorithm (GA) to identify and design time-dependent processing protocols that achieve tailor-made microstructures. We simulate the physical vapor deposition of a binary-alloy thin film by employing a phase-field model, where deposition rates and diffusivities are controlled via the genetic algorithm. Our GA-guided protocols achieve targeted microstructures with lateral and vertical concentration modulations, as well as more complex, hierarchical microstructures previously not described in simple structure zone diagrams. Our algorithm provides insight to experimentalists looking for additional avenues to design novel thin-film microstructures.

36 MATERIALS SCIENCE↗
NASA NTRS1987

Monitoring biological impacts of space shuttle launches from Vandenberg Air Force Base: Establishment of baseline conditions

Space shuttle launches produce environmental impacts resulting from the formation of an exhaust cloud containing hydrogen chloride aerosols and aluminum oxide particulates. Studies have shown that most impacts occur near-field (within 1.5 km) of the launch site while deposition from launches occurs far-field (as distant as 22 km). In order to establish baseline conditions of vegetation and soils in the areas likely to be impacted by shuttle launches from Vandenberg Air Force Base (VAFB), vegetation and soils in the vicinity of Space Launch Complex-6 (SLC-6) were sampled and a vegetation map prepared. The areas likely to be impacted by launches were determined considering the structure of the launch complex, the prevailing winds, the terrain, and predictions of the Rocket Exhaust Effluent Diffusion Model (REEDM). Fifty vegetation transects were established and sampled in March 1986 and resampled in September 1986. A vegetation map was prepared for six Master Planning maps surrounding SLC-6 using LANDSAT Thematic Mapper imagery as well as color and color infrared aerial photography. Soil samples were collected form the 0 to 7.5 cm layer at all transects in the wet season and at a subsample of the transects in the dry season and analyzed for pH, organic matter, conductivity, cation exchange capacity, exchangeable Ca, Mg, Na, K, and Al, available NH3-N, PO4-P, Cu, Fe, Mn, Zn, and TKN.

Schmaizer, Paul A.↗