Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “phase field method”

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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

The effects of cloud inhomogeneities upon radiative fluxes, and the supply of a cloud truth validation dataset

A series of cloud and sea ice retrieval algorithms are being developed in support of the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) Science Team objectives. These retrievals include the following: cloud fractional area, cloud optical thickness, cloud phase (water or ice), cloud particle effective radius, cloud top heights, cloud base height, cloud top temperature, cloud emissivity, cloud 3-D structure, cloud field scales of organization, sea ice fractional area, sea ice temperature, sea ice albedo, and sea surface temperature. Due to the problems of accurately retrieving cloud properties over bright surfaces, an advanced cloud classification method was developed which is based upon spectral and textural features and artificial intelligence classifiers.

Welch, Ronald M.↗
NASA NTRS2006

Tracking Energy Flow Using a Volumetric Acoustic Intensity Imager (VAIM)

A new measurement device has been invented at the Naval Research Laboratory which images instantaneously the intensity vector throughout a three-dimensional volume nearly a meter on a side. The measurement device consists of a nearly transparent spherical array of 50 inexpensive microphones optimally positioned on an imaginary spherical surface of radius 0.2m. Front-end signal processing uses coherence analysis to produce multiple, phase-coherent holograms in the frequency domain each related to references located on suspect sound sources in an aircraft cabin. The analysis uses either SVD or Cholesky decomposition methods using ensemble averages of the cross-spectral density with the fixed references. The holograms are mathematically processed using spherical NAH (nearfield acoustical holography) to convert the measured pressure field into a vector intensity field in the volume of maximum radius 0.4 m centered on the sphere origin. The utility of this probe is evaluated in a detailed analysis of a recent in-flight experiment in cooperation with Boeing and NASA on NASA s Aries 757 aircraft. In this experiment the trim panels and insulation were removed over a section of the aircraft and the bare panels and windows were instrumented with accelerometers to use as references for the VAIM. Results show excellent success at locating and identifying the sources of interior noise in-flight in the frequency range of 0 to 1400 Hz. This work was supported by NASA and the Office of Naval Research.

Klos, Jacob↗
DOE OSTI2023

Phase Field Modeling of Corrosion for Design of Next-Generation Magnesium-Aluminum Vehicle Joints

In this study, WPI and its partners sought to understand galvanic corrosion in Mg-Al friction stir weld (FSW) joints. Unlike fusion welds, FSW is done in the solid state. This limits (but does not eliminate) the formation of brittle Al-Mg intermetallics, resulting in better joint strength and corrosion resistance. This study aimed to understand corrosion behavior of Mg-Al FSW joints by using a Cahn-Hilliard phase field electrochemistry model based on the work of Pongsaksawad et al. Its context challenge problem was an ultra-light door designed by Magna with 6061-Al sheet outer and ZEK100-Mg sheet inner panels, joined by an adhesive and FSW. Unlike most FSW joints, which weld through the softer material into the harder one, the Al sheet’s hem joint around Mg required the FSW tool to through the harder Al into the softer Mg, which was itself a challenge. Despite this challenge, PNNL succeeded in consistently making very strong FSW lap joints, by using a triflute tool and power control. Characterization showed mostly separate Al and Mg regions in the weld, usually with a hook protruding from Al into Mg, which likely contributed to high joint strength. Nanohardness mapping showed higher hardness in thin Mg swirls into the Al nugget, likely due to the intermetallics. The study used the industry standard SAE J2334 Cyclic Corrosion Test (CCT) to simulate corrosion conditions over the life of a vehicle, as well as linear polarization testing and ASTM G71 pitting corrosion testing. Corrosion reactions were: H 2 O + ½O 2 + 2e - → 2OH - at the Al cathode, and either Mg → Mg 2+ + 2e - or Mg + 2OH - → Mg(OH) 2 + 2e - at the Mg anode, with electron transfer through the joint. Mass loss increased with Al section length, supporting the hypothesis that galvanic corrosion in Al-Mg joints is limited by the cathodic reaction. Both pitting corrosion (first anode reaction) and hydroxide film formation (second reaction) were observed. The J2334 CCT test method was slightly modified, adding a chromium solution cleaning step after each week of corrosion testing, in order to remove Mg oxides and hydroxides and accurately measure mass loss over time. That said, FSW joints treated in this way did not did not exhibit significant reduction in lap shear strength of the joints vs. newly welded samples which had not undergone corrosion. Corrosion appeared to consist mainly of pitting in Mg in a way which did not directly affect the bond between the materials where the joints failed. A Cahn-Hilliard phase field model coupled with electrical potential described the electronically mediated galvanic corrosion electrochemical reactions described above in the four-component Al-Mg-H 2 O-(H 2 O+½O 2 ) system. This 2-D model successfully predicted Mg(OH)₂ formation at the anode in some circumstances, and Mg pitting corrosion in others, and showed the correct electric field directions in both cases. The model predicted Mg pit corrosion depth within a factor of two of measured pits.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS1999

Sonic Boom Prediction and Minimization of the Douglas Reference OPT5 Configuration

Conventional CFD methods and grids do not yield adequate resolution of the complex shock flow pattern generated by a real aircraft geometry. As a result, a unique grid topology and supersonic flow solver was developed at Northrop Grumman based on the characteristic behavior of supersonic wave patterns emanating from the aircraft. Using this approach, it was possible to compute flow fields with adequate resolution several body lengths below the aircraft. In this region, three-dimensional effects are diminished and conventional two-dimensional modified linear theory (MLT) can be applied to estimate ground pressure signatures or sonic booms. To accommodate real aircraft geometries and alleviate the burdensome grid generation task, an implicit marching multi-block, multi-grid finite-volume Euler code was developed as the basis for the sonic boom prediction methodology. The Thomas two-dimensional extrapolation method is built into the Euler code so that ground signatures can be obtained quickly and efficiently with minimum computational effort suitable to the aircraft design environment. The loudness levels of these signatures can then be determined using a NASA generated noise code. Since the Euler code is a three-dimensional flow field solver, the complete circumferential region below the aircraft is computed. The extrapolation of all this field data from a cylinder of constant radius leads to the definition of the entire boom corridor occurring directly below and off to the side of the aircraft's flight path yielding an estimate for the entire noise "annoyance" corridor in miles as well as its magnitude. An automated multidisciplinary sonic boom design optimization software system was developed during the latter part of HSR Phase 1. Using this system, it was found that sonic boom signatures could be reduced through optimization of a variety of geometric aircraft parameters. This system uses a gradient based nonlinear optimizer as the driver in conjunction with a computationally efficient Euler CFD solver (NIIM3DSB) for computing the three-dimensional near-field characteristics of the aircraft. The intent of the design system is to identify and optimize geometric design variables that have a beneficial impact on the ground sonic boom. The system uses a simple wave drag data format to specify the aircraft geometry. The geometry is internally enhanced and analytic methods are used to generate marching grids suitable for the multi-block Euler solver. The Thomas extrapolation method is integrated into this system, and hence, the aircraft's centerline ground sonic boom signature is also automatically computed for a specified cruise altitude and yields the parameters necessary to evaluate the design function. The entire design system has been automated since the gradient based optimization software requires many flow analyses in order to obtain the required sensitivity derivatives for each design variable in order to converge on an optimal solution. Hence, once the problem is defined which includes defining the objective function and geometric and aerodynamic constraints, the system will automatically regenerate the perturbed geometry, the necessary grids, the Euler solution, and finally the ground sonic boom signature at the request of the optimizer.

Siclari, Michael J.↗
DOE OSTI2022

Bulk-like first-order magnetoelastic transition in FeRh particles

Near-equiatomic, chemically-ordered iron-rhodium (FeRh) alloy is a fundamentally interesting material that may become useful in niche applications making use of its unique magneto functional phenomena, for example, the giant inverse magnetocaloric effect near room temperature that is associated with a sharp first-order magnetic phase transition. The nearly discontinuous antiferromagnetic-ferromagnetic phase transformation in bulk FeRh is well-known; however, the transition broadens considerably in fine particles and films with thickness less than 50 nm, precluding their potential applications. Here, we report an abrupt, bulk-like first-order magnetoelastic transformation in powders consisting of sub-micron particles of nearly equiatomic FeRh compound synthesized via solid-state mechanochemical co-reduction of FeF 2 and RhCl3 and subsequent heat treatments. In this work, we demonstrate that annealing at temperatures ranging from 600 °C to 800 °C enables tailoring phase content, particle size, and magnetic properties of the powders. A maximum magnetic-field-induced entropy change of ~10 J/kg K at μ0ΔH = 1 T has been achieved in powders annealed at 800 °C. The retention of extraordinary responsiveness in sub-micron particles of FeRh is likely to open doors for system component fabrication using additive manufacturing methods, along with new opportunities to employ FeRh in theranostics.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Metastable Materials Discovery in the Age of Large-Scale Computation

Computational materials discovery has been successful in predicting novel, technologically relevant materials. However, it has remained focused almost exclusively on finding ground-state structures. Now that the lower-hanging fruit has been found in many fields of application, materials exploration is moving toward metastable materials: higher energy phases that are stable at practical time scales. Because of the challenges associated with predicting which phases are realistic, this class of materials has remained relatively unexplored, despite numerous examples of metastable structures with unmatched properties (e.g., diamond). This article highlights recent advances in developing computational and theoretical methods for predicting useful and realizable metastable materials. Topics discussed cover (1) the latest strategies for identifying potential metastable phases, (2) methodologies for assessing which phases can be realized experimentally, and (3) current approaches to estimate the lifetime of metastable materials.

amorphous materials↗
DOE OSTI2021

Eulerian simulation of complex suspensions and biolocomotion in three dimensions

Significance Fluid–structure interactions are ubiquitous in many natural and man-made environments. They are difficult to study analytically, and therefore accurate and flexible computational methods are an indispensable tool in the field. Typically, fluids are simulated with a fixed background computational mesh, whereas a solid is simulated with a mesh that moves with it, making it challenging to couple the two. Here we develop a three-dimensional computational method where both fluid and solid can be represented on a fixed computational grid, which simplifies the coupling between the two phases considerably. Our method can simulate scenarios that are difficult to do with other methods, such as complex suspensions containing soft particles that are both heavier and lighter than the fluid.

3D fluid-structure interaction↗
DOE OSTI2020

Magnetotransport and electronic structure of the antiferromagnetic semimetal YbAs

A number of rare-earth monopnictides have topologically nontrivial band structures together with magnetism and strong electronic correlations. In order to examine whether the antiferromagnetic (AFM) semimetal YbAs (TN=0.5 K) exhibits such a scenario, we have grown high-quality single crystals using a flux method, and characterized the magnetic properties and electronic structure using specific heat, magnetotransport, and angle-resolved photoemission spectroscopy (ARPES) measurements, together with density functional theory (DFT) calculations. Both ARPES and DFT calculations find no evidence for band inversions in YbAs, indicating a topologically trivial electronic structure. From low-temperature magnetotransport measurements, we map the field-temperature phase diagram, where we find the presence of a field stabilized phase distinct from the AFM phase at low temperatures. Additionally, an extremely large magnetoresistance (XMR) for both YbAs and the nonmagnetic counterpart LuAs is also observed, which can consistently be accounted for by the presence of electron-hole compensation. Moreover, an angle-dependent study of the Shubnikov-de Haas effect oscillations reveals very similar Fermi surfaces between YbAs and LuAs, with light effective masses down to at least 0.5 K, indicating that the Yb-4f electrons are well localized, and do not contribute to the Fermi surface. However, the influence of the localized Yb-4f electrons on the magnetotransport of YbAs can be discerned from the distinct temperature dependence of the XMR compared to that of LuAs, which we attribute to the influence of short-ranged spin correlations that appear well above TN.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

YbV 3 ⁢Sb 4 and EuV 3 ⁢Sb 4 vanadium-based kagome metals with Yb 2+ and Eu 2+ zigzag chains

Here, in this work, we present Yb⁢V 3 ⁢Sb 4 and Eu⁢V 3 ⁢Sb 4 , two compounds exhibiting slightly distorted vanadium-based kagome nets interleaved with zigzag chains of divalent Yb 2+ and Eu 2+ ions. Single crystal growth methods are reported alongside magnetic, electronic, and heat capacity measurements. Yb⁢V 3⁢ Sb 4 is a nonmagnetic metal with no collective phase transitions observed between 60 mK and 300 K. Conversely, Eu⁢V 3 ⁢Sb 4 is a magnetic kagome metal exhibiting easy-plane ferromagneticlike order below T C =32 K with hints of modulated spin texture under low field. With Yb⁢V 3 ⁢Sb 4 and Eu⁢V 3 ⁢Sb 4 we demonstrate another direction for the discovery and development of vanadium-based kagome metals while incorporating the chemical and magnetic degrees of freedom offered by a rare-earth sublattice.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Phase-contrast 3D tomography of HeLa cells grown in PLLA polymer electrospun scaffolds using synchrotron X-rays

Advanced imaging is useful for understanding the three-dimensional (3D) growth of cells. X-ray tomography serves as a powerful noninvasive, nondestructive technique that can fulfill these purposes by providing information about cell growth within 3D platforms. There are a limited number of studies taking advantage of synchrotron X-rays, which provides a large field of view and suitable resolution to image cells within specific biomaterials. In this study, X-ray synchrotron radiation microtomography at Diamond Light Source and advanced image processing were used to investigate cellular infiltration of HeLa cells within poly L-lactide (PLLA) scaffolds. This study demonstrates that synchrotron X-rays using phase contrast is a useful method to understand the 3D growth of cells in PLLA electrospun scaffolds. Two different fiber diameter (2 and 4 µm) scaffolds with different pore sizes, grown over 2, 5 and 8 days in vitro , were examined for infiltration and cell connectivity. After performing visualization by segmentation of the cells from the fibers, the results clearly show deeper cell growth and higher cellular interconnectivity in the 4 µm fiber diameter scaffold. Finally, this indicates the potential for using such 3D technology to study cell–scaffold interactions for future medical use.

3D culture↗
DOE OSTI2024

Photovoltaic Analysis and Response Support (PARS) Platform for Solar Situational Awareness and Resiliency Services

The project's primary objective is to develop a digital-twin based Photovoltaic (PV) Analysis and Response Support (PARS) platform, which aims to provide real-time situational awareness and optimal response plans. This platform is designed to enhance the performance of hybrid PV systems, making them competitive with or even superior to conventional generation resources. The PARS platform enabled the project team to develop and evaluate an extensive suite of grid support functionalities for the hybrid PV systems to enhance grid performance, across key areas including visibility, dispatchability, security, resilience, and reliability. Given the global push toward achieving 100% clean energy by 2035, there is a significant increase in the integration of inverter-based resources (IBRs) throughout the energy grid. Effectively managing the inherent variability and uncertainty associated with IBRs is crucial for ensuring cost-effectiveness, reliability, and security in both the main grid and islanded microgrids. Constrained to a limited array of IEEE test systems or standard feeder models, traditional IBR modeling struggles to assimilate new field data, accurately reflect system dynamics, and adapt to the evolving energy landscape. In our project, we embraced a Digital Twin (DT) strategy for crafting the PARS platform. A digital twin acts as a precise virtual counterpart of a physical system, built on historical data and continuously honed with real-time insights. This enables the high-fidelity DT to accurately mirror current system operations and forecast future scenarios. Consequently, the PARS platform becomes an ideal environment for testing and refining monitoring, control, power, and energy management algorithms designed to boost hybrid PV system performance. The defining feature of the PARS platform, distinguishing it from other advanced simulation tools, is its exceptional adaptability. This is achieved by employing actual network topologies and utilizing real-time field data for fine-tuning and calibration, ensuring a close emulation of real-world conditions. The project deliverables include: 1) High-fidelity IBR models and tools for real-time parameterization, utilizing real-time field measurements to refine IBR models for enhanced accuracy and performance; 2) Grid-forming and Grid-following capabilities to deliver resilience services, including blackstart, voltage and frequency support, cold-load pick-up, power reserves, and three-phase load balancing across grid-connected and microgrid settings; 3) Machine learning-based forecasting tools and methods for generating synthetic data and topologies, creating diverse and realistic simulation environments for evaluating varied operational scenarios; 4) Advanced microgrid power and energy management algorithms for optimizing the integration and operation of PV, storage, and demand response resources within both feeder and community scales. The power grid data sets are provided by four utility companies in North Carolina and the New York Power Administration. Acting as industry advisors, our industry partners communicated stakeholder needs and regulatory standards to the research teams, aiding technology transfer by incorporating the developed methodologies into their daily operations. This collaboration ensures that the PARS platform, functioning as a power system digital twin, enhances our understanding of IBR dynamic behaviors and enables the development and evaluation of IBR control functions that match or exceed the capabilities of conventional synchronous generators.

14 SOLAR ENERGY↗
DOE OSTI2024

GENERATING CORCOS-COHERENT SIGNAL SERIES USING PHASE PERTURBATION​

Phase perturbations are used to simulate the time-domain behavior of pressure loads from turbulent boundary layers while faithfully representing the coherence decay behaviors of turbulence-generated pressure fields. Turbulence "evolves" as it proceeds downstream, and coherence decays as the decaying exponential of the product of a flow-related constant with frequency and distance. Time domain methods for a line, a plane, a closed cyclic surface, and most of a cone were given.

Swift, Stephen Hales↗
NASA NTRS1993

Particle tracing in the magnetosphere - New algorithms and results

A fast and efficient method is employed to trace charged particles through realistic magnetospheric electric and magnetic fields, greatly reducing computer simulation times. The method is applicable for particles having arbitrary charge, energy, or pitch angle but which conserve the first two adiabatic invariants. An efficient method is used to classify drift orbits, which greatly facilitates identification of the last closed drift path or other drift boundaries. The time-independent evolution of the bounce-averaged phase space density along convective drift orbits is calculated. These three tools can be used to quantitatively describe convective evolution of the particle distribution from the tail, an essential step in understanding the production of unstable distributions in the magnetosphere.

Sheldon, R. B.↗
NASA NTRS1996

Numerical Modeling of Spray Combustion with an Unstructured-Grid Method

The present unstructured-grid method follows strictly the basic finite volume forms of the conservation laws of the governing equations for the entire flow domain. High-order spatially accurate formulation has been employed for the numerical solutions of the Navier-Stokes equations. A two-equation k-epsilon turbulence model is also incorporated in the unstructured-grid solver. The convergence of the resulted linear algebraic equation is accelerated with preconditioned Conjugate Gradient method. A statistical spray combustion model has been incorporated into the present unstructured-grid solver. In this model, spray is represented by discrete particles, rather than by continuous distributions. A finite number of computational particles are used to predict a sample of total population of particles. Particle trajectories are integrated using their momentum and motion equations and particles exchange mass, momentum and energy with the gas within the computational cell in which they are located. The interaction calculations are performed simultaneously and eliminate global iteration for the two-phase momentum exchange. A transient spray flame in a high pressure combustion chamber is predicted and then the solution of liquid-fuel combusting flow with a rotating cup atomizer is presented and compared with the experimental data. The major conclusion of this investigation is that the unstructured-grid method can be employed to study very complicated flow fields of turbulent spray combustion. Grid adaptation can be easily achieved in any flow domain such as droplet evaporation and combustion zone. Future applications of the present model can be found in the full three-dimensional study of flow fields of gas turbine and liquid propulsion engine combustion chambers with multi-injectors.

Shang, H. M.↗
NASA NTRS1958

Three-Dimensional Orbits of Earth Satellites, Including Effects of Earth Oblateness and Atmospheric Rotation

The principal purpose of the present paper is to present sets of equations which may be used for calculating complete trajectories of earth satellites from outer space to the ground under the influence of air drag and gravity, including oblateness effects, and to apply these to several examples of entry trajectories starting from a circular orbit. Equations of motion, based on an "instantaneous ellipse" technique, with polar angle as independent variable, were found suitable for automatic computation of orbits in which the trajectory consists of a number of revolutions. This method is suitable as long as the trajectory does not become nearly vertical. In the terminal phase of the trajectories, which are nearly vertical, equations of motion in spherical polar coordinates with time as the independent variable were found to be more suitable. In the first illustrative example the effects of the oblateness component of the earth's gravitational field and of atmospheric rotation were studied for equatorial orbits. The satellites were launched into circular orbits at a height of 120 miles, an altitude sufficiently high that a number of revolutions could be studied. The importance of the oblateness component of the earth's gravitational field is shown by the fact that a satellite launched at circular orbital speed, neglecting oblateness, has a perigee some 67,000 feet lower when oblateness forces are included in the equations of motion than when they are not included. Also, the loss in altitude per revolution is double that of a satellite following an orbit not subject to oblateness. The effect of atmospheric rotation on the loss of altitude per revolution was small. As might be surmised, the regression of the line of nodes as predicted by celestial mechanics is unchanged when drag is included. It is clear that the inclination of the orbital plane to the equator will be relatively unaffected by drag for no atmospheric rotation since the drag lies in the orbital plane in this case. With the inclusion of atmospheric rotation it was found that the inclination of the plane changed about one-millionth of a radian per revolution. Thus the prediction of the position of the orbital plane of an earth satellite is not complicated by the introduction of drag. The line of apsides, which without drag but with oblateness moves slowly in space, tends to move with the satellite when drag is included in the calculations. As a results, the usual linearized solutions based on oblateness alone must be basically altered when drag is included to take into account the rapid movement of the line of apsides. In the second illustrative example the final revolution was calculated to impact for a number of trajectories in an orbital plane inclined at 650 to the equator. Of particular interest is the large effect the oblateness gravitational field and atmospheric rotation can have on the impact point. For a value of CDA/m of unity, and for an initial downward angle at 80 miles altitude of 0.01 radian, such as might be utilized for manned re-entry, oblateness had an influence of about 300 miles in the impact point, and atmospheric rotation had about a 150-mile influence.

Nielsen, Jack N.↗
NASA NTRS2003

Application of Powder Diffraction Methods to the Analysis of Short- and Long-Range Atomic Order in Nanocrystalline Diamond and SiC: The Concept of the Apparent Lattice Parameter (alp)

Two methods of the analysis of powder diffraction patterns of diamond and SiC nanocrystals are presented: (a) examination of changes of the lattice parameters with diffraction vector Q ('apparent lattice parameter', alp) which refers to Bragg scattering, and (b), examination of changes of inter-atomic distances based on the analysis of the atomic Pair Distribution Function, PDF. Application of these methods was studied based on the theoretical diffraction patterns computed for models of nanocrystals having (i) a perfect crystal lattice, and (ii), a core-shell structure, i.e. constituting a two-phase system. The models are defined by the lattice parameter of the grain core, thickness of the surface shell, and the magnitude and distribution of the strain field in the shell. X-ray and neutron experimental diffraction data of nanocrystalline SiC and diamond powders of the grain diameter from 4 nm up to micrometers were used. The effects of the internal pressure and strain at the grain surface on the structure are discussed based on the experimentally determined dependence of the alp values on the Q-vector, and changes of the interatomic distances with the grain size determined experimentally by the atomic Pair Distribution Function (PDF) analysis. The experimental results lend a strong support to the concept of a two-phase, core and the surface shell structure of nanocrystalline diamond and SiC.

Palosz, B.↗
NASA NTRS2011

Application of Laser Mass Spectrometry to Art and Archaeology

REMPI laser mass spectrometry is a combination of resonance enhanced multiphoton ionization spectroscopy and time of flight mass spectrometry, This technique enables the collection of mass specific optical spectra as well as of optically selected mass spectra. Analytes are jet-cooled by entrainment in a molecular beam, and this low temperature gas phase analysis has the benefit of excellent vibronic resolution. Utilizing this method, mass spectrometric analysis of historically relevant samples can be simplified and improved; Optical selection of targets eliminates the need for chromatography while knowledge of a target's gas phase spectroscopy allows for facile differentiation of molecules that are in the aqueous phase considered spectroscopically indistinguishable. These two factors allow smaller sample sizes than commercial MS instruments, which in turn will require less damage to objects of antiquity. We have explored methods to optimize REMPI laser mass spectrometry as an analytical tool to archaeology using theobromine and caffeine as molecular markers in Mesoamerican pottery, and are expanding this approach to the field of art to examine laccaic acid in shellacs.

Gulian, Lase Lisa E.↗