Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “field validation”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Sensor Distortion Effects in Photon Monte Carlo Simulations

We present a detailed method to simulating sensor distortions using a photon and electron Monte Carlo method. We use three-dimensional electrostatic simulations to parameterize the perturbed electric field profile for nonideal sensor details. We follow the conversion of simulated photons, and the subsequent response of the converted electrons to the electric field pattern. These nonideal sensor details can be implemented efficiently in a Monte Carlo approach. We demonstrate that the nonideal sensor distortions have a variety of observable consequence including the modification of the astrometric pattern; the distortion of the electron diffusion size and shape; and the distortion of flats. We show analytic validation of the diffusion physics, reproduce two kinds of edge distortion, and show qualitative validation of field-free regions, lithography errors, and fringing. We also demonstrate that there are two related effects of doping variation having different observable consequences. We show that field distortions from accumulated electrons lead to intensity-dependent point-spread functions and the sublinear variance in flats. The method is implemented in the Photon Simulator and the code is publicly available.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Prediction of Climate-Specific Degradation Rate for Photovoltaic Encapsulant Discoloration

Encapsulant discoloration is a well-known field degradation mode of crystalline-silicon photovoltaic modules, particularly in the hot climate zones. The discoloration rate is influenced by several weathering factors, such as UV light, module temperature, and humidity, as well as the permeability of oxygen into the module. Here, a rate dependence model employing the modified Arrhenius equations to predict the degradation rate for encapsulant discoloration in different climates is presented. Two modeling approaches are introduced, which utilize the field and accelerated UV testing degradation data in conjunction with the field meteorological data to determine the acceleration factor for encapsulant browning. A novel method of accelerated UV stress testing at three simultaneous module temperatures in a single environmental chamber test run is implemented to estimate the activation energy for browning. The test was performed on three field-retrieved modules to capture the wear-out failure mechanism. The degradation in short-circuit current Isc rather than maximum power is used as a decisive parameter for the discoloration analysis. Furthermore, the developed model has been used to predict the Isc degradation rate for the Arizona field characterized by a hot and dry climate and is validated against the field-measured value. It has also been applied to other climate types, e.g., the cold and dry climate of New York.

14 SOLAR ENERGY↗
DOE OSTI2022

2.3.3.404 - National Lab and University Collaboration for MHK Instrumentation and Data Processing Tools

Field and laboratory validation, testing, demonstration, and operation are critical steps for increasing the technology readiness level of marine energy (ME) converters because they provide high-quality testing and performance data that are critical information used to feed all aspects of technology development. This project, in partnership with industry, enables the marine and hydrokinetic energy (MHK) community to reliably and efficiently collect, process, manage, and share quality data by facilitating access to and development of instrumentation, guidelines and data processing/QA tools. Under this project, open-source data processing code (MHKiT) and tools (ME Data Pipeline, MRE Code Hub, PRIMRE Code Catalog), instrumentation (loads measurements), data acquisition systems (miniDAQ), and measurement guidance tools (Telesto, high EMI guidance) were developed to facilitate the collection and processing of quality laboratory and field data. Overall, this project is intended to improve the quality of the data collected during laboratory and field demonstration projects by standardizing the collection and processing techniques, as well as by improving access to instrumentation, code, and measurement guidance. Quality data will, in turn, lead to improved knowledge capture following ME device testing.

data processing↗
DOE OSTI2025

Diagnosing electric and magnetic fields in laser-driven coil targets

Laser-driven capacitor coils are widely used to generate intense magnetic fields for various applications in high-energy-density (HED) physics research. Accurate measurement of the magnetic fields is essential but challenging, due to the overlapping contributions from magnetic and electric fields in proton radiography, which is the primary tool diagnosing the field generation around the coils. In this study, we systematically analyze proton radiographs obtained from laser-driven capacitor-coil targets along two orthogonal axes under various electromagnetic field conditions, including magnetic field only, electric field only, and combined electromagnetic fields. By analyzing key features in the radiographs, we distinguish and characterize the respective contributions from magnetic and electric fields. Using detailed simulations validated by experimental benchmarks, methods to isolate and quantify the magnetic field and electric field are given. The methods are successfully applied to determine the electric current and charge distribution in a double coil configuration. Our findings provide insights into improving the diagnostic capability of proton radiography, potentially leading to more accurate measurements of electromagnetic fields and enhancing the utility of laser-driven capacitor coils in HED experiments.

electromagnetic field diagnostics↗
DOE OSTI2025

Validation of a comprehensive first-principles-based framework for predicting the performance of future stellarators

This paper presents the validation of the GENE-KNOSOS-Tango framework for recovering both the steady-state plasma profiles in the considered radial domain and selected turbulence trends in a stellarator. This framework couples the gyrokinetic turbulence code GENE, the neoclassical transport code KNOSOS, and the transport solver Tango in a multi-timescale simulation feedback loop. Ion-scale kinetic-electron and electron-scale adiabatic-ion flux-tube simulations were performed to evolve the density and temperature profiles for four OP1.2b W7-X scenarios. The simulated density and temperature profiles showed good agreement with the experimental data using a reasonable set of boundary conditions. Equally important was the reproduction of observed trends for several turbulence properties, such as density fluctuations and turbulent heat diffusivities. Key effects were also touched upon, such as electron-scale turbulence and the neoclassical radial electric field shear. The validation of the GENE-KNOSOS-Tango framework enables credible predictions of physical phenomena in stellarators and reactor performance based on a given set of edge parameters.

Gyrokinetic simulations↗
DOE OSTI2023

Viewfactor and Raytracing for AgriPV Modeling

View factor models are used in due diligence software to calculate rear irradiance for bifacial modules. An intermediate step in this calculation is the irradiance at the ground level, which can be leveraged for evaluating the Photo Active Radiation available for crops in Agrivoltaic setups. This paper presents the metrics and modifications to the model for ground irradiance study with the view factor approach, compares it to the raytracing method, and validates it with field measurements of ground irradiance. It is found that for the clearances, row-to-row setups, and tilts studied, the view factor method matches with raytracing results within 2% MBD. The comparison is performed for the nine most common agriPV configurations using high-performance computing and the NSRDB database for the whole US, with results and data made available open-source on the InSPIRE AgriPV website.

AgriPV↗
DOE OSTI2021

A Heterogeneous System for Eagle Detection, Deterrent, and Wildlife Collision Detection for Wind Turbines (Final Technical Report)

This report summarizes the design, implementation, and test of an integrated system for automated detection and deterrence of eagles, with included wind turbine blade strike detection and imaging functionality. A machine learning approach was used in conjunction with a 360° camera system for automated detection and classification of golden eagles. This was developed using footage obtained from trained golden eagles and other raptors, in collaboration with wildlife biologists and professional bird handlers. Oregon State University developed a visual deterrent system, which uses inflatable anthropomorphic sculptures with random, kinetic motion to deter eagles, and conducted limited field testing on live eagles; the deterrent can be triggered by the visual detection of eagles using the vision system. Finally, a multi-sensor module was developed that is mounted at the turbine blade root. This module measures vibration and other motions to detect blade strikes, and an integrated on-blade camera captures an image of any impacting objects. Long-term, this blade strike detection system is intended to support an automatic monitoring and certification system for the eagle detection and deterent system. Independent field testing of each system component is described. Testing of the integrated system on an operational wind turbine was conducted across three separate field tests. This includes multi-day fields tests on a General Electric 1.5MW wind turbine at the National Renewable Energy Laboratory (NREL) National Wind Technology Center (NWTC) in Boulder, CO in October 2018 and July 2019; installation procedures, test procedures, and a summary of collected data are presented. A third multi-day on-turbine field test is also presented, which was performed using a General Electric 1.5MW wind turbine at the North American Wind Research and Training Center (NAWRTC) at Mesalands Community College, Tucumcari, NM in April 2019. Across these field tests, the vision system was demonstrated using unmanned aerial vehicles (UAV), and the eagle classification algorithm was not tested; the visual deterrent system was demonstrated, including automatic, remote deployment following surrogate visual detections; and, multi-sensor on-blade data was recorded across multiple wind turbine operational conditions and through more than 100 surrogate blade strikes using soft projectiles, including the successful demonstration of automatic image capture of striking objects. This data set was also used for offline development and validation of enhanced collision detection algorithms. As summarized in this report, the development and field validation of an integrated detection, deterrent, and blade collision detection system represents a critical proof of concept for future technology development of related detection and deterrent technologies, where both deterrent as well as collision detection recording devices are needed for future siting, monitoring, and operation of wind turbine installations, both onshore and offshore.

17 WIND ENERGY↗
DOE OSTI2023

Quantifying optical loss of high-voltage degradation modes in photovoltaic modules using spectral analysis

The direct current bias for photovoltaic (PV) modules interconnected in series-strings may include both high voltage negative (“HV-”) and positive (“HV+”) polarity with respect to the electrical ground. Multiple degradation modes, resulting in quantifiable optical loss, were found to occur during HV-/HV+ sequential stress, including corrosion of the external glass surface, encapsulant delamination (at its interfaces with the glass and the PV cell), internal haze formation (resulting from a chemical interaction between the glass and the encapsulant), corrosion and migration of the gridlines, and corrosion of the silicon nitride (Si x N y ) antireflective coating on the cell. The effects of these separate modes were examined using monolithic (e.g., glass or PV cell) and laminated-coupon (glass/encapsulant/glass or glass/encapsulant/cell/encapsulant/backsheet) specimens. Characterizations during and after unbiased accelerated testing at 85°C/85% relative humidity included spectrophotometry, optical microscopy, electron microscopy, and ellipsometry. For some module components (i.e., the glass and the Si x N y coating), the optical performance was determined through iterative analysis of empirical measurements. Concentrating on just their spectral effect, a novel model was then developed to estimate the transfer of light to the PV cell and the return of light from the PV module with simultaneous degradation mechanisms, which was compared with a mini-module previously subjected to HV-/HV+ stress. Here the model suggests that one third of the current loss observed for the mini-module can be attributed to the optical degradation of the packaging materials. The dominant degradation modes include encapsulant delamination and corrosion of the Si x N y coating. Recommendations are given so that the optical model may be improved relative to accelerated testing and validated relative to field aging.

14 SOLAR ENERGY↗
DOE OSTI2025

Energy performance of an operational government building retrofitted with ceiling phase change material tiles in a mixed-humid climate

The aging U.S. building stock requires various retrofit measures to enhance their energy efficiency. Here, this study explores the integration of thermal energy storage and advanced building controls as viable retrofit solutions for load flexibility and peak demand response while maintaining the occupants' comfort. A detailed assessment is conducted on the energy use of an administrative building in Sumner County, Kansas, focusing on the implementation of phase change materials (PCMs) in the ceiling of occupied zones. First, a time-resolved, whole-building energy model is developed in EnergyPlus, incorporating complex thermal behaviors such as air exchange between the plenum space and occupied zones, envelope leakage, and operational schedules. The model is then validated using experimental field test data, and subsequently a parametric assessment of key PCM properties and application strategies is performed to evaluate cooling electricity demand benefits. The parametric study shows that the optimal retrofit strategy, comprising a PCM with 23°C peak melting temperature, 0.125 in. (3.17 mm) thickness, and 150 kJ/kg latent heat, combined with active controls that include 8 h of precooling, forced convection under the ceiling, and a 2°C thermostat setback during peak hours, can result in a maximum load shift during the peak period of 99.6 % and the total electricity savings during the peak period of 98.9 % for the optimum case and thus provide significant cost savings under time-of-use pricing scenarios.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Mismatch losses in a PV system due to shortened strings

Numerous events may require intentional removal of one or more photovoltaic modules from a string, shortening the length of the string relative to others within the array, resulting in a string length mismatch. The impact of such a mismatch is not well understood either in measurable operational effects (voltage, current, power) or in the potential effects on long-term module health. It is impractical to solely approach this problem experimentally due to the size and complexity of arrays that may experience string length shortening. Here this work presents simulations, validated through limited field experiments on a two string array, providing a basis from which more complex arrays and scenarios may be explored. Refinement of the simulation achieves an overall error in IMPP for the nominal (S 1 ) and test strings (S 2 ) between the simulation and experimental values, through all test conditions, of +0.35 ± 1.46% and -0.36 ± 1.58% respectively. Shortening one of two strings by one module results in a power loss greater than the power contribution of the module alone (1.29 module equivalents); the impact increases through the maximum test case of a six module mismatch with a power loss equivalent to more than 11 modules. The impact of using string-end blocking diodes is presented with an emphasis at the array maximum power point and at open circuit. Implications are discussed for arrays of higher complexity.

14 SOLAR ENERGY↗
DOE OSTI2020

Hydrogen wide area monitoring of LH2 releases

The characterization of liquid hydrogen (LH2) releases has been identified as an international research priority to expand the safe use of hydrogen as an energy carrier. The elucidation of LH2 release behavior will require the development of dispersion and other models, guided and validated by empirical field measurements such as those afforded by Hydrogen Wide Area Monitoring (HyWAM). HyWAM can be defined as the quantitative spatial and temporal three-dimensional monitoring of planned or unintentional hydrogen releases. With support provided through the FCH JU Prenormative Research for the Safe Use of Liquid Hydrogen (PRESLHY) program, HSE performed a series of LH2 releases to characterize the dispersion and pooling behavior of cold hydrogen releases. The NREL Sensor Laboratory developed a HyWAM system based upon a distributed array of point sensors that is amenable for profiling cold hydrogen plumes. The NREL Sensor Laboratory and HSE formally committed to collaborate on profiling the LH2 releases. This collaboration included the integration of the NREL HyWAM into the HSE LH2 release hardware. This was achieved through a deployment plan jointly developed by the NREL and HSE personnel. Under this plan, the NREL Sensor Laboratory provided multiple HyWAM modules that accommodated 32 sampling points for near-field hydrogen profiling during the HSE PRESLHY LH2 releases. The NREL HyWAM would be utilized throughout the LH2 release study performed under PRESLHY by HSE, including Work Package 3 (WP3—Release and Mixing--Rainout) and subsequent work packages (WP4—Ignition and WP5—Combustion). Under the auspices of the PRESLHY WP6 (Implementation), data and findings from the HSE LH2 Releases are to be made available to stakeholders in the hydrogen community. Comprehensive data analysis and dissemination is ongoing, but the integration of the NREL HyWAM into the HSE LH2 Release Apparatus and its performance as well as some key outcomes of the LH2 releases in WP3 are presented.

08 HYDROGEN↗
DOE OSTI2022

Development and application of ReaxFF methodology for understanding the chemical dynamics of metal carbonates in aqueous solutions

A new ReaxFF reactive force field has been developed for metal carbonate systems including Na + , Ca 2+ , and Mg 2+ cations and the CO 3 2– anion. This force field is fully transferable with previous ReaxFF water and water/electrolyte descriptions. The Me–O–C (Me = metal) three-body valence angle parameters and Me–C non-reactive parameters of the force field have been optimized against quantum mechanical calculations including equations of state, heats of formation, heats of reaction, angle distortions and vibrational frequencies. The new metal carbonate force field has been validated using molecular dynamics simulations to study the solvation and reactivity of metal and carbonate ions in water at 300 K and 700 K. The coordination radius and self-diffusion coefficient show good consistency with existing experimental and simulation results. The angular distribution analysis explains the structural preference of carbonate ions to form carbonates and bicarbonates, where Na + predominantly forms carbonates due to weaker angular strain, while Ca 2+ and Mg 2+ prefer to form bicarbonate monodentate in nature. Residence time distribution analyses on different systems reveal the role of ions in accelerating and decelerating the dynamics of water and carbonate ions under different thermodynamic conditions. The formation and dissolution of bicarbonates and carbonates in solution were explored on the basis of the protonation capability in different systems. Here, the nucleation phenomenon of metal carbonates at ambient and supercritical conditions is explained from the perspective of cluster formation over time: Ca 2+ ions can form prenucleation clusters at ambient temperature but show saturation with increasing temperature, whereas Na + and Mg 2+ ions show a rapid increase in cluster size and amount upon increasing time and temperature.

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

Evaluation of nonlinear control performance of air handling units under variable operation conditions using root distribution approach

The control systems at Air handling units (AHUs) have nonlinear nature, resulting potentially in sluggish, aggressive, or oscillatory responses under variable operation conditions. This paper is to develop an approach for evaluating the nonlinear control performance of AHUs under variable operation conditions and for identifying the impact of varying multiple operation parameters which guides controller design. To do so, an operation parameter-based root distribution technique, intended to sketch the paths traced by closed-loop poles as multiple operation parameters are varied, is first introduced. The technique—which may be viewed as an extension of conventional, single-variable root locus—is then applied to a fan speed control system that includes a motor with a variable frequency drive and a supply fan with an air duct under variable fan speed and duct static pressure. Using the technique, the system’s stability region is graphically identified. The system is also found to be aggressive and gradually lose its stability as the fan speed decreases and the duct static pressure increases, which is validated by the field tests. Here, the root distribution approach is an effective tool to evaluate the performance of any nonlinear control systems at AHUs with multiple variable operation parameters without requiring time-domain simulations.

42 ENGINEERING↗
DOE OSTI2024

Magic Angles and Fractional Chern Insulators in Twisted Homobilayer Transition Metal Dichalcogenides

We explain the appearance of magic angles and fractional Chern insulators in twisted K-valley homobilayer transition metal dichalcogenides by mapping their continuum model to a Landau level problem. Here, our approach relies on an adiabatic approximation for the quantum mechanics of valence band holes in a layer-pseudospin field that is valid for sufficiently small twist angles and on a lowest Landau level approximation that is valid for sufficiently large twist angles. It provides a simple qualitative explanation for the nearly ideal quantum geometry of the lowest moir´e miniband at particular twist angles, predicts that topological flat bands occur only when the valley-dependent moir´e potential is sufficiently strong compared to the interlayer tunneling amplitude, and provides a convenient starting point for the study of interactions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Holo-Omics disentangle drought response and biotic interactions among plant, endophyte and pathogen

Holo-omics provide a novel opportunity to study the interactions among fungi from different functional guilds in host plants in field conditions. We address the entangled responses of plant pathogenic and endophytic fungi associated with sorghum when droughted through the assembly of the most abundant fungal, endophyte genome from rhizospheric metagenomic sequences followed by a comparison of its metatranscriptome with the host plant metabolome and transcriptome. The rise in relative abundance of endophytic Acremonium persicinum (operational taxonomic unit 5 (OTU5)) in drought co-occurs with a rise in fungal membrane dynamics and plant metabolites, led by ethanolamine, a key phospholipid membrane component. The negative association between endophytic A. persicinum (OTU5) and plant pathogenic fungi co-occurs with a rise in expression of the endophyte's biosynthetic gene clusters coding for secondary compounds. Endophytic A. persicinum (OTU5) and plant pathogenic fungi are negatively associated under preflowering drought but not under postflowering drought, likely a consequence of variation in fungal fitness responses to changes in the availability of water and niche space caused by plant maturation over the growing season. Our findings suggest that the dynamic biotic interactions among host, beneficial and harmful microbiota in a changing environment can be disentangled by a blending of field observation, laboratory validation, holo-omics and ecological modelling.

Chen, Peilin↗
DOE OSTI2023

AmeriFlux FLUXNET-1F US-Bar Bartlett Experimental Forest

This is the AmeriFlux Management Project (AMP) created FLUXNET-1F version of the carbon flux data for the site US-Bar Bartlett Experimental Forest. This is the FLUXNET version of the carbon flux data for the site US-Bar Bartlett Experimental Forest produced by applying the standard ONEFlux (1F) software. Site Description - The Bartlett Experimental Forest (448170 N, 71830 W) is located within the White Mountains National Forest in north-central New Hampshire, USA. The 1050 ha forest extends across an elevational range from 200 to 900 m a.s.l. It was established in 1931 and is managed by the USDA Forest Service Northeastern Research Station in Durham, NH. The climate is humid continental with short, cool summers (mean July temperature, 19.8C) and long, cold winters (mean January temperature, 9.8C). Annual precipitation averages 130 cm and is distributed evenly throughout the year. Soils are developed from glacial till and are predominantly shallow, well-drained spodosols. At lowto mid-elevation, vegetation is dominated by northern hardwoods (American beech, Fagus grandifolia; sugar maple, Acer saccharum; yellow birch, Betula alleghaniensis; with some red maple, Acer rubrum and paper birch, Betula papyrifera). Conifers (eastern hemlock, Tsuga canadensis; eastern white pine, Pinus strobus; red spruce, Picea rubens) are occasionally found intermixed with the more abundant deciduous species but are generally confined to the highest (red spruce) and lowest (hemlock and pine) elevations. In 2003, the site was adopted as a NASA North American Carbon Program (NACP) Tier-2 field research and validation site. A 26.5 m high tower was installed in a low-elevation northern hardwood stand in November, 2003, for the purpose of making eddy covariance measurements of the forest–atmosphere exchange of CO2, H2O and radiant energy. Continuous flux and meteorological measurements began in January, 2004, and are ongoing. Average canopy height in the vicinity of the tower is approximately 20–22 m. In the tower footprint, the forest is predominantly classified into red maple, sugar maple, and American beech forest types. Leaf area index in the vicinity of the tower is 3.6 as measured by seasonal litterfall collection, and 4.5 as measured by the optically based Li-Cor LAI-2000 instrument. Further site information: http://www.fs.fed.us/ne/durham/4155/bartlett.htm

Richardson, Andrew↗
DOE OSTI2024

Heart Shape to Fracture Distance: Characterizing Hydraulic Fracture Propagation before Hits

Estimating the distance from the hydraulic fracture tip to the monitor well can be useful for fracture characterization, well spacing optimization, and preventing parent-child well interference. A heart-shaped signal is referred to as the extensional precursor of a fracture hit recorded by crosswell strain measurements and can serve as a vital tool for such estimation. This study incorporates the 3D displacement discontinuity method (DDM) to understand the impact of fracture geometry and monitor well offset on the heart-shaped signal’s characteristics. Results from numerical simulation and analytical solutions reveal a strong linear correlation between the spatial extent of the heart-shaped signal and the fracture tip distance. This relationship was further developed to predict tip distance using field data from the Hydraulic Fracture Test Site 2 (HFTS2). A reasonable approximation result from field data further validates the methodology. In addition, it is worth noting that the estimation accuracy depends on the ratio between fracture dimension and tip distance. The findings of this study offer a novel approach for real-time monitoring and characterizing hydraulic fracture propagation, which can be further used for well spacing optimization in unconventional and enhanced geothermal system reservoir development, as well as caprock integrity monitoring for carbon sequestration projects.

58 GEOSCIENCES↗
DOE OSTI2020

VTO FY20Q4 quarterly report

Project Objective.This project has three objectives: (1) tuning microstructures of 3D-printed SSE separators,(2) determining material and processing compatibilities with cathode printing,and (3)3D printing of sintering-free SSE separators.Project Impact.All-solid-state lithium batteries are difficult to process due to the brittleness of ceramic materials, poor solid-solid contact, and electrolyte-electrode stability issues. As a result, the energy and power density and also cycling stability are far from satisfying. This project will address fabrication difficulties by using state-of-the-art 3D-printing techniques that can introduce 3D interfaces and architectures to enhance solid-solid contact and reduce charge transfer resistance. Success will benefit the DOE by establishing the best manufacturing methods for ASSBsto achieve VTO goals on the performance of beyond Li-ion batteries. Approach.The project employs 3D-printing techniques to manufacture SSEs and related components for ASSBs. The team starts with direct ink writing (DIW) to develop ink recipes with desired rheological properties and explore post-sintering approaches to achieve high densification. In parallel with DIW3D printing, the team also explores other 3D printing options, such as projection microstereolithography (PuSL). Both sintering and sintering-free approaches can gain benefit from 3D printing and therefore will be investigated.Out-Year Goals.The team will determine particle sizes and morphologies as well as post-processing conditions that deliver good ionic conductivity and charge transfer properties of LLZTO separators. They will examine co-sintering of LLZTO/cathode and determine the effects of material, particle size, surface coating, and sintering conditions.Collaborations.Microstructures, ionic conductivities, and mechanical properties will be provided to the LLNL simulation team, led by PI B. Wood, for establishing and validating phase-field modeling methods.

36 MATERIALS SCIENCE↗