Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “predictive simulation”

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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

SU-F-T-471: Simulated External Beam Delivery Errors Detection with a Large Area Ion Chamber Transmission Detector

Purpose: The Integral Quality Monitor (IQM), developed by iRT Systems GmbH (Koblenz, Germany) is a large-area, linac-mounted ion chamber used to monitor photon fluence during patient treatment. Our previous work evaluated the change of the ion chamber’s response to deviations from static 1×1 cm2 and 10×10 cm2 photon beams and other characteristics integral to use in external beam detection. The aim of this work is to simulate two external beam radiation delivery errors, quantify the detection of simulated errors and evaluate the reduction in patient harm resulting from detection. Methods: Two well documented radiation oncology delivery errors were selected for simulation. The first error was recreated by modifying a wedged whole breast treatment, removing the physical wedge and calculating the planned dose with Pinnacle TPS (Philips Radiation Oncology Systems, Fitchburg, WI). The second error was recreated by modifying a static-gantry IMRT pharyngeal tonsil plan to be delivered in 3 unmodulated fractions. A radiation oncologist evaluated the dose for simulated errors and predicted morbidity and mortality commiserate with the original reported toxicity, indicating that reported errors were approximately simulated. The ion chamber signal of unmodified treatments was compared to the simulated error signal and evaluated in Pinnacle TPS again with radiation oncologist prediction of simulated patient harm. Results: Previous work established that transmission detector system measurements are stable within 0.5% standard deviation (SD). Errors causing signal change greater than 20 SD (10%) were considered detected. The whole breast and pharyngeal tonsil IMRT simulated error increased signal by 215% and 969%, respectively, indicating error detection after the first fraction and IMRT segment, respectively. Conclusion: The transmission detector system demonstrated utility in detecting clinically significant errors and reducing patient toxicity/harm in simulated external beam delivery. Future work will evaluate detection of other smaller magnitude delivery errors.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Predicting Mechanical Properties from Microstructure Images in Fiber-Reinforced Polymers Using Convolutional Neural Networks

Evaluating the mechanical response of fiber-reinforced composites can be extremely time-consuming and expensive. Machine learning (ML) techniques offer a means for faster predictions via models trained on existing input–output pairs and have exhibited success in composite research. This paper explores a fully convolutional neural network modified from StressNet, which was originally used for linear elastic materials, and extended here for a non-linear finite element (FE) simulation to predict the stress field in 2D slices of segmented tomography images of a fiber-reinforced polymer specimen. The network was trained and evaluated on data generated from the FE simulations of the exact microstructure. The testing results show that the trained network accurately captures the characteristics of the stress distribution, especially on fibers, solely from the segmented microstructure images. The trained model can make predictions within seconds in a single forward pass on an ordinary laptop, given the input microstructure, compared to 92.5 h to run the full FE simulation on a high-performance computing cluster. These results show promise in using ML techniques to conduct fast structural analysis for fiber-reinforced composites and suggest a corollary that the trained model can be used to identify the location of potential damage sites in fiber-reinforced polymers.

Sun, Yixuan (ORCID:0000000311093380)↗
DOE OSTI2021

Coupling 1D xRAGE simulations with machine learning for graded inner shell design optimization in double shell capsules

Advances in machine learning provide the ability to leverage data from expensive simulations of high-energy-density experiments to significantly cut down on computational time and costs associated with the search for optimal target designs. This study presents an application of cutting-edge Bayesian optimization methods to the one-dimensional (1D) design optimization of double shell graded layer targets for inertial confinement fusion experiments. This investigation attempts to reduce hydrodynamic instabilities while retaining high yields for future NIF experiments. Machine learning methods can use predictive physics simulations to identify graded layer designs from within the vast design space that demonstrate high predicted performance, including novel designs with high uncertainty in performance that may hold unexpected promise. By applying machine learning tools to the simulation design, we map the trade-off between 1D yield and instability, specifically isolating parameter ranges, which maintain high performance while showing significantly improved Rayleigh–Taylor stability over the point design. Furthermore, the groundwork laid in this study will be a useful design tool for future NIF experiments with graded layer targets.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

The non-linear power spectrum of the Lyman alpha forest

The Lyman alpha forest power spectrum has been measured on large scales by the BOSS survey in SDSS-III at z∼ 2.3, has been shown to agree well with linear theory predictions, and has provided the first measurement of Baryon Acoustic Oscillations at this redshift. However, the power at small scales, affected by non-linearities, has not been well examined so far. We present results from a variety of hydrodynamic simulations to predict the redshift space non-linear power spectrum of the Lyα transmission for several models, testing the dependence on resolution and box size. A new fitting formula is introduced to facilitate the comparison of our simulation results with observations and other simulations. The non-linear power spectrum has a generic shape determined by a transition scale from linear to non-linear anisotropy, and a Jeans scale below which the power drops rapidly. In addition, we predict the two linear bias factors of the Lyα forest and provide a better physical interpretation of their values and redshift evolution. The dependence of these bias factors and the non-linear power on the amplitude and slope of the primordial fluctuations power spectrum, the temperature-density relation of the intergalactic medium, and the mean Lyα transmission, as well as the redshift evolution, is investigated and discussed in detail. A preliminary comparison to the observations shows that the predicted redshift distortion parameter is in good agreement with the recent determination of Blomqvist et al., but the density bias factor is lower than observed. We make all our results publicly available in the form of tables of the non-linear power spectrum that is directly obtained from all our simulations, and parameters of our fitting formula.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Early Solar Wind and Dynamo Magnetic Field Topology Predictions for (16) Psyche and Other Asteroids

Abstract Asteroid (16) Psyche is a metal‐rich body that might record an ancient coherent magnetization if some relict crust or mantle is preserved. Herein, we use magnetohydrodynamic simulations to predict (16) Psyche's field topology for several distinct pathways: (i) an early solar wind‐induced magnetization imparted after a larger body was impacted, forming the present‐day asteroid, (ii) a core dynamo magnetization imparted in an asteroid that is either presently largely intact or was a rubble pile, and (iii) magnetization in the turbulent solar nebula disk. For pathway (i) we find the field to be predominantly dipolar and spin axis‐aligned. For pathway (ii) we find the field to be either dipolar and spin axis‐misaligned, or highly multipolar. We also find that (iii) a field produced earlier before the solar nebula cleared, would be highly multipolar. In cases (i) and (ii) we also place constraints on the field strength. Simple detection of a magnetic field without constraining its topology and temporal variability would be insufficient to confirm a remanent source, due to the influence of the present‐day solar wind, electromagnetic induction, and (16) Psyche's high obliquity. For sufficiently strong fields however, the field topology and orientation may reveal key observable consequences of the nature and history of (16) Psyche. Our framework is also broadly applicable to the study of magnetic fields from other asteroids.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Variability in the sensitivity among model simulations of permafrost and carbon dynamics in the permafrost region between 1960 and 2009

A significant portion of the large amount of carbon (C) currently stored in soils of the permafrost region in the Northern Hemisphere has the potential to be emitted as the greenhouse gases CO 2 and CH 4 under a warmer climate. In this study we evaluated the variability in the sensitivity of permafrost and C in recent decades among land surface model simulations over the permafrost region between 1960 and 2009. The 15 model simulations all predict a loss of near-surface permafrost (within 3 m) area over the region, but there are large differences in the magnitude of the simulated rates of loss among the models (0.2 to 58.8 × 10 3 km 2 yr –1 ). Sensitivity simulations indicated that changes in air temperature largely explained changes in permafrost area, although interactions among changes in other environmental variables also played a role. All of the models indicate that both vegetation and soil C storage together have increased by 156 to 954 Tg C yr –1 between 1960 and 2009 over the permafrost region even though model analyses indicate that warming alone would decrease soil C storage. Increases in gross primary production (GPP) largely explain the simulated increases in vegetation and soil C. The sensitivity of GPP to increases in atmospheric CO 2 was the dominant cause of increases in GPP across the models, but comparison of simulated GPP trends across the 1982–2009 period with that of a global GPP data set indicates that all of the models overestimate the trend in GPP. Disturbance also appears to be an important factor affecting C storage, as models that consider disturbance had lower increases in C storage than models that did not consider disturbance. Furthermore, to improve the modeling of C in the permafrost region, there is the need for the modeling community to standardize structural representation of permafrost and carbon dynamics among models that are used to evaluate the permafrost C feedback and for the modeling and observational communities to jointly develop data sets and methodologies to more effectively benchmark models.

58 GEOSCIENCES↗
DOE OSTI2022

Computational studies of hydraulic stressors for biological performance assessment in a hydropower plant with Kaplan turbine

We report Hydropower is currently one of the preeminent sources of renewable energy in the United States and globally. Hydropower plants also have detrimental impacts on the environment and ecology, including direct impacts to anadromous fish populations. The computational fluid dynamics (CFD) – based Biological Performance Assessment (BioPA) toolset is used for biological evaluations of fish passage through hydropower plants. The hydraulic environment of a hydropower plant was evaluated using CFD coupled with discrete element method (DEM) simulations. The predicted flow field and particle collision rate were validated against the experimental data in a water flume that has an idealized hydroturbine distributor geometry. Flow simulations were conducted to evaluate the hydraulic stressors, such as nadir pressure, fluid shear, runner collision, in a physical scale in a hydropower plant with Kaplan turbine which are responsible for injury and mortality of fish in a downstream migration. The cumulative exposure probability for the nadir pressure and collision with turbine runner was found to decrease with increased discharge rate. The lowest discharge rate shows the higher value of cumulative shear exposure probability. The maximum value of collision velocity increases with increased discharge rate. We offer the conclusions that will help in understanding various hydraulic stressors for biological assessment for environmentally sustainable hydroturbine passage.

13 HYDRO ENERGY↗
DOE OSTI2020

Recreating Giants Impacts in the Laboratory: Shock Compression of MgSiO 3 Bridgmanite to 14 Mbar

Understanding giant impacts requires accurate description of how extreme pressures and temperatures affect the physical properties of the constituent materials. Here, we report shock experiments on two polymorphs of MgSiO 3 : enstatite and bridgmanite (perovskite) crystals.We obtain pressure-density shock equation of state to 14 Mbar and more than 9 g/cm 3 , a 40% increase in density from previous data on MgSiO 3 . Density-functional-theory molecular dynamics (DFT-MD) simulations provide predictions for the shock Hugoniot curves for bridgmanite and enstatite and suggest that the Grüneisen parameter decreases with increasing density. The good agreement between the simulations and the experimental data, including for the shock temperature along the enstatite Hugoniot reveals that DFT-MD simulations reproduce well the behavior of dense fluid MgSiO 3 .We also reveal a high optical reflectance indicative of a metal-like electrical conductivity which supports the hypothesis that magma oceans may contribute to planetary magnetic field generation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2012

CoCoMANS Ocean App [Slides]

Thermohaline circulation is that part of ocean circulation which is driven by density gradients. It is widely believed that changes in the ocean thermohaline circulation on the 10 to 100 year time scale can instigate rapid climate change. Scientists world wide are currently using simulation tools based on semi-implicit and split methods that are first-order in time. From a numerical error perspective, the quality of such simulations is suspect. As a result of numerical stability issues resulting from semi-implicit and split methods they are forced to use low resolution grids; ie high resolution grids will require prohibitive numbers of timesteps. Having the ability to perform such simulation with reduced splitting errors, using high resolution grids, and using second-order accurate time integration methods will be a quantum leap forward in predictive ocean simulation.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Improved non-local electron thermal transport model for two-dimensional radiation hydrodynamics simulations

An implicit, non-local thermal conduction algorithm based on the algorithm developed by Schurtz, Nicolai, and Busquet (SNB) [Schurtz et al., Phys. Plasmas 7, 4238 (2000)] for non-local electron transport is presented and has been implemented in the radiation-hydrodynamics code DRACO. To study the model's effect on DRACO's predictive capability, simulations of shot 60 303 from OMEGA are completed using the iSNB model, and the computed shock speed vs. time is compared to experiment. Temperature outputs from the iSNB model are compared with the non-local transport model of Goncharov et al. [Phys. Plasmas 13, 012702 (2006)]. Effects on adiabat are also examined in a polar drive surrogate simulation. Results show that the iSNB model is not only capable of flux-limitation but also preheat prediction while remaining numerically robust and sacrificing little computational speed. Additionally, the results provide strong incentive to further modify key parameters within the SNB theory, namely, the newly introduced non-local mean free path. This research was supported by the Laboratory for Laser Energetics of the University of Rochester.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A multiphysics model of the versatile test reactor based on the MOOSE framework

The traditional modeling approach for sodium fast reactor cores relies on separate physics models, where the fuel performance, thermal–hydraulics, and neutronics calculations required to predict the core physics characteristics for nominal conditions are decoupled by relying on user-imposed boundary conditions. Here, this paper aims at evaluating the impact of multiphysics simulations for predicting the core characteristics of the Versatile Test Reactor, which is being designed as a 300-MWt sodium-cooled fast reactor. The purpose of the Versatile Test Reactor is to accelerate the testing of advanced nuclear materials in the United States. The proposed multiphysics model relies on the Griffin reactor physics code, the SAM thermal–hydraulic system code, the BISON fuel performance code, as well as generic Multiphysics Object-Oriented Simulation Environment capabilities implemented in the open-source tensor mechanics module. For k eff calculations, the introduction of a tight coupling between the neutronics, thermo-mechanical and thermal–hydraulics models induces a change of around 543 pcm in the eigenvalue, compared to the traditional standalone neutronics calculation where approximate temperature profiles are used. The multiphysics model is then employed for quantifying the impact of the thermal conductivity uncertainties on some of the key figures of merit, such as the fuel centerline temperature, assembly powers, and keff for nominal core conditions. As anticipated, uncertainties on fuel thermal conductivity mostly impact the fuel centerline temperature, and to a lesser extend the k eff .

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Viscoelastic necking dynamics between attractive microgels

Hypothesis: Microgels can deform and interpenetrate and display colloid/polymer duality. The effective interaction of microgels in the collapsed state is governed by the interplay of polymer–solvent interfacial tension and bulk elasticity. A connecting neck is shown to mediate microgel interaction, but its temporal evolution has not been addressed. Here we hypothesize that the necking dynamics of attractive microgels exhibits liquid-like or solid-like behavior over different time and length scales. Experiments: We simulate the merging and pinching of attractive microgels with different crosslinking densities in explicit solvent using dissipative particle dynamics. The temporal coalescence dynamics of microgels is investigated and compared with simple liquid and polymeric droplets. We model the neck growth on long time scales using Maxwell model of polymer relaxation and compare the theoretical prediction with simulation data. The mechanical strength of the neck is characterized systematically via simulated pinch-off of microgels by steered molecular dynamics. Findings: We evidence a crossover in the coalescence dynamics reflecting the viscoelastic signature of microgels. In contrast to the common knowledge that viscoelastic materials respond elastically on short time scales, the early expansion of the microgel neck exhibits a linear behavior, similar to the viscous coalescence of liquid droplets. However, the late regime with arrested dynamics resembles sintering of solid particles. Through an analytical model relating microgel dynamics to neck growth, we show that the long-term behavior is governed by stress relaxation of the polymers in the neck region and predict an exponential decay in the rate of growth, which agrees favorably with the simulation. Different from coalescence, the thread thinning in microgel breakup primarily highlights its polymeric characteristics.

42 ENGINEERING↗
DOE OSTI2019

Predictions for energy correlators probing substructure of groomed heavy quark jets

We develop an effective field theory (EFT) framework to perform an analytic calculation for energy correlator observables computed on groomed heavy-quark jets. A soft-drop grooming algorithm is applied to a jet initiated by a massive quark to minimize soft contamination effects such as pile-up and multi-parton interactions. We specifically consider the two-particle energy correlator as an initial application of this EFT framework to compute heavy quark jet substructure. We find that there are different regimes for the event shapes, depending on the size of the measured correlator observable, that require the use of different EFT formulations, in which the quark mass and grooming parameters may be relevant or not. We use the EFT to resum large logarithms in the energy correlator observable in terms of the momentum of a reconstructed heavy hadron to NLL' accuracy and subsequently match it to a full QCD O(α s ) cross section, which we also compute. We compare our predictions to simulations in Pythia for e + e – collisions. We find a good agreement with partonic simulations, as well as hadronic ones with an appropriate shape function used to describe nonperturbative effects and the heavy quark hadron decay turned off. We also predict the scaling behavior for the leading nonperturbative power correction due to hadronization. Consequently, we can give a prediction for the energy correlator distribution at the level of the reconstructed heavy hadron. This work provides a general framework for the analysis of heavy quark jet substructure observables.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

MB-Fit: Software infrastructure for data-driven many-body potential energy functions

Many-body potential energy functions (MB-PEFs), which integrate data-driven representations of many-body short-range quantum mechanical interactions with physics-based representations of many-body polarization and long-range interactions, have recently been shown to provide high accuracy in the description of molecular interactions from the gas to the condensed phase. Here, we present MB-Fit, a software infrastructure for the automated development of MB-PEFs for generic molecules within the TTM-nrg (Thole-type model energy) and MB-nrg (many-body energy) theoretical frameworks. Besides providing all the necessary computational tools for generating TTM-nrg and MB-nrg PEFs, MB-Fit provides a seamless interface with the MBX software, a many-body energy and force calculator for computer simulations. Given the demonstrated accuracy of the MB-PEFs, particularly within the MB-nrg framework, we believe that MB-Fit will enable routine predictive computer simulations of generic (small) molecules in the gas, liquid, and solid phases, including, but not limited to, the modeling of quantum isomeric equilibria in molecular clusters, solvation processes, molecular crystals, and phase diagrams.

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

Rapid Prediction of Adsorption Isotherms of a Diverse Range of Molecules in Hyper-Cross-Linked Polymers

Hyper-cross-linked polymers (HCPs) are materials with permanent porosity that have potential applications in adsorption-based chemical separations. We use molecular simulations to predict the largest set of adsorption isotherms available to date for HCPs. Specifically, we predict room-temperature isotherms for 24 distinct molecules in 9 different HCPs that vary in both the monomer composition and cross-link density. This information allows us to develop efficient approximate models that predict adsorption isotherms in HCPs based only on the Henry’s constant for adsorption, textural information about the HCP, and bulk phase properties of the adsorbing molecule. This approximate approach will be useful in considering the adsorption of diverse groups of molecules in HCPs in the future. Our data also enables us to describe the selectivity of adsorption for the 276 binary mixtures that can be formed from the 24 molecules we simulated. Furthermore, this information greatly expands the information that is available about adsorption selectivity in HCPs.

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

Exploratory Precipitation Metrics: Spatiotemporal Characteristics, Process-Oriented, and Phenomena-Based

Abstract Precipitation sustains life and supports human activities, making its prediction one of the most societally relevant challenges in weather and climate modeling. Limitations in modeling precipitation underscore the need for diagnostics and metrics to evaluate precipitation in simulations and predictions. While routine use of basic metrics is important for documenting model skill, more sophisticated diagnostics and metrics aimed at connecting model biases to their sources and revealing precipitation characteristics relevant to how model precipitation is used are critical for improving models and their uses. This paper illustrates examples of exploratory diagnostics and metrics including 1) spatiotemporal characteristics metrics such as diurnal variability, probability of extremes, duration of dry spells, spectral characteristics, and spatiotemporal coherence of precipitation; 2) process-oriented metrics based on the rainfall–moisture coupling and temperature–water vapor environments of precipitation; and 3) phenomena-based metrics focusing on precipitation associated with weather phenomena including low pressure systems, mesoscale convective systems, frontal systems, and atmospheric rivers. Together, these diagnostics and metrics delineate the multifaceted and multiscale nature of precipitation, its relations with the environments, and its generation mechanisms. The metrics are applied to historical simulations from phases 5 and 6 of the Coupled Model Intercomparison Project. Models exhibit diverse skill as measured by the suite of metrics, with very few models consistently ranked as top or bottom performers compared to other models in multiple metrics. Analysis of model skill across metrics and models suggests possible relationships among subsets of metrics, motivating the need for more systematic analysis to understand model biases for informing model development.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Development of a Comprehensive Two-Phase Flow Database for the Validation of NEK-2P

Three-dimensional (3-D) two-phase Computational Fluid Dynamics (CFD) codes are emerging as a powerful and potentially practical tool for applications in which detailed local flow information is needed. However, two-phase flow models and associated closure relations are not well established for CFD applications, which is partly due to the lack of high-quality validation data. The main objective of this work is to develop a comprehensive database of two-phase flows that can be used to validate two-phase CFD codes such as NEK-2P. In this project, four advanced local measurement systems, including Particle Image Velocimetry and Planar Laser-Induced Fluorescence (PIV-PLIF), high-speed imaging, x-ray densitometry, and multi-sensor conductivity probe are employed to measure the local two-phase flow parameters of both gas and liquid phases. By combining these techniques, the local void fraction, bubble velocity, interfacial area concentration, bubble frequency, liquid velocity, turbulence intensity, etc., in various two-phase flow regimes can be obtained. These local measurement techniques are first used in a 25.4 mm circular pipe test section. Seven air-water two-phase flow conditions spanning the bubbly, slug, churn-turbulent, and annular flow regimes are measured in this facility. The obtained database contains the radial profiles of both gas- and liquid-phase parameters at three axial locations along the test section. The second facility used in this work contains a 30 mm × 10 mm rectangular test section, in which three two-phase flow and two single-phase flow conditions are measured. Two-dimensional distributions of local two-phase flow parameters in the cross-sectional plane are measured in this facility at three axial locations as well. A facility featuring a 3×3 electrically heated rod bundle is also designed and being constructed in this project. This facility is specially designed for optical measurements and is expected to provide high-quality boiling data in the future. Preliminary analyses have been performed for the data taken in the 25.4 mm circular pipe. Both center-peaked and wall-peaked void fraction profiles have been observed in the data depending on the two-phase flow conditions and/or developing lengths. The 1-D drift-flux model was evaluated with the newly obtained datasets, in which both gas- and liquid-phases data were directly measured. The distribution parameter model has been optimized based on a new void-profile classification method proposed in this study. The optimized drift-flux model shows a significant improvement in predicting both gas velocity and void fraction. The measured liquid-phase turbulence was used to benchmark Sato’s turbulence model considering the bubble-induced shear stress for the three tested bubbly flows. The benchmark results showed good agreement between the PIV measurements and model predictions. In the bubbly flows tested that have low void fractions less than 3%, the effect of the bubble-induced turbulence was found not significant. However, the bubble-induced shear stress becomes important with the increase of the void fraction. The Conjugate Heat Transfer (CHT) model was developed and implemented in NEK-2P. This model allows the coupled simulation of the solid domain and two-phase fluid domain, allowing the specification of realistic boundary conditions. The CHT implementation was verified first with non-boiling simulations. The predicted temperatures in the fluid domain were shown to be identical in simulations with or without the CHT model. The CHT model was then validated through simulations of three Becker benchmark CHF tests performed under both Dryout (DO) and Departure from Nucleate Boiling (DNB) conditions. Reasonably good agreement was obtained between calculated wall temperatures and corresponding experimental data.

42 ENGINEERING↗
DOE OSTI2020

Validation of finite element analysis strategy to investigate acoustic levitation in a two-axis acoustic levitator

A two-axis acoustic levitator can be used to generate a standing pressure wave capable of levitating solid and liquid particles at appropriate input conditions. This work proposes a simulation framework to investigate the two-axis levitation particle stability using a commercial, computational fluid dynamics software based on the harmonic solution to the acoustic wave equation. The simulation produced predictions of the standing wave that include a strong "+" shaped pattern of nodes and anti-nodes that are aligned with the levitator axes. To verify the simulation, a levitator was built and used to generate the standing wave. The field was probed with a microphone and a motorized-scanning system. After scaling the simulated pressure to the measured pressure, the magnitudes of the sound pressure level at corresponding high-pressure locations were different by no more than 5%. This is the first time a measurement of a two-axis levitator standing pressure wave has been presented and shown to verify simulations. Furthermore, as an additional verification, the authors consulted high speed camera measurements of a reference-levitator transducer, which was found to have a maximum peak-to-peak displacement of 50 +/- 5 mu m. The reference-levitator is known to levitate water at 160 dB. The system for this work was simulated to match the operation of the reference-levitator so that it produced sound pressure levels of 160 dB. This pressure was achieved when the transducer maximum peak-to-peak displacement was 50.8 mu m. The agreement between the two levitators' displacements provides good justification that the modeling approach presented here produces reliable results.

2 axis levitator↗