Search NASA⌕ Search

SEARCH · Search NASA

Results for “generalized finite differences”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,405 records · Page 78
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Assembling Multiphysics Nuclear Reactor Simulations Using the MOOSE Framework

The Multiphysics Object Oriented Simulation Environment (MOOSE) [1] is an open-source, parallel finite element framework which provides the foundation for many advanced modeling and simulation tools developed under the Department of Energy (DOE) Nuclear Energy Advanced Modeling and Simulation (NEAMS) Program [2] for the analysis of advanced reactors. The MOOSE framework provides the common foundational capability on which many NEAMS codes for reactor analysis are built. The MOOSE framework also includes several systems to assemble unique workflows and couplingamong MOOSE-based applications. In particular, the MultiApp and Transfer Systems are widely used to assemble different MOOSE-based or MOOSE-wrapped physics applications together to perform loosely or tightly coupled multiphysics simulations. The National Reactor Innovation Center (NRIC) Virtual Test Bed (VTB) [3] hosts publicly available nuclear reactor multiphysics simulation examples which leverage MOOSE’s MultiApp System to meet the modeling needs of different reactor types. The flexibility and robustness of coupling provided by MOOSE permits rapid development of coupled physics models for a wide range of reactor types and events

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2018

A quantum hydrodynamical description for scrambling and many-body chaos

Recent studies of out-of-time ordered thermal correlation functions (OTOC) in holographic systems and in solvable models such as the Sachdev-Ye-Kitaev (SYK) model have yielded new insights into manifestations of many-body chaos. So far the chaotic behavior has been obtained through explicit calculations in specific models. In this paper we propose a unified description of the exponential growth and ballistic butterfly spreading of OTOCs across different systems using a newly formulated “quantum hydrodynamics,” which is valid at finite ℏ and to all orders in derivatives. The scrambling of a generic few-body operator in a chaotic system is described as building up a “hydrodynamic cloud,” and the exponential growth of the cloud arises from a shift symmetry of the hydrodynamic action. The shift symmetry also shields correlation functions of the energy density and flux, and time ordered correlation functions of generic operators from exponential growth, while leads to chaotic behavior in OTOCs. The theory also predicts an interesting phenomenon of the skipping of a pole at special values of complex frequency and momentum in two-point functions of energy density and flux. This pole-skipping phenomenon may be considered as a “smoking gun” for the hydrodynamic origin of the chaotic mode. Here, we also discuss the possibility that such a hydrodynamic description could be a hallmark of maximally chaotic systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

The Transient Multi-Level method for Monte Carlo reactor statics calculations

The Transient Multi-Level (TML) method is applied to a time-dependent Monte Carlo transport solver to offload some of the computational burden of the expensive Monte Carlo solve to lower-order Coarse Mesh Finite Difference (CMFD) and Exact Point Kinetics Equations (EPKE) solvers via factorization of the neutron flux at the transport and CMFD levels using the Predictor Corrector Quasi-Static Method (PCQM). The Monte Carlo transient is solved by a modified fission source iteration scheme that introduces a single transient source bank. The method is implemented in the production-level Monte Carlo code, Shift, and verified with prescribed reactivity ramps from the two-dimensional version of the C5G7-TD reactor benchmark. The results show that, as compared to other quasi-static methods, the TML reduces the stochastic noise inherent to the transient Monte Carlo solver by factors of ~2 to 6 for various norm comparisons of the reactor power amplitude. Finally, the TML additionally reduces the number of Monte Carlo evaluations needed to simulate the transient, leading to roughly an order of magnitude improvement in CPU time relative to the standard PCQM for the problems tested.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2016

Strain induced novel quantum magnetotransport properties of topological insulators

Recent theoretical and experimental researches have revealed that the strained bulk HgTe can be regarded as a three-dimensional topological insulator (TI). Motivated by this, we explore the strain effects on the transport properties of the HgTe surface states, which are modulated by a weak 1D in-plane electrostatic periodic potential in the presence of a perpendicular magnetic field. We analytically derive the zero frequency (dc) diffusion conductivity for the case of quasielastic scattering in the Kubo formalism, and find that, in strong magnetic field regime, the Shubnikov–de Haas oscillations are superimposed on top of the Weiss oscillations due to the electric modulation for null and finite strain. Furthermore, the strain is shown to remove the degeneracy in inversion symmetric Dirac cones on the top and bottom surfaces. This accordingly gives rise to the splitting and mixture of Landau levels, and the asymmetric spectrum of the dc conductivity. These phenomena, not known in a conventional 2D electron gas and even in a strainless TI and graphene, are a consequence of the anomalous spectrum of surface states in a fully stained TI. These results should be valuable for electronic and spintronic applications of TIs, and thus we fully expect to see them in the further experiment. - Highlights: • The strain removes the degeneracy in inversion symmetric Dirac cones. • The strain gives rise to the splitting and mixture of the Landau levels. • The strain leads to the asymmetric spectrum of the dc conductivity. • Shubnikov de Haas oscillations are shown to be superimposed on Weiss oscillations. • Interplay between strain and electric field causes different occupancy of TI states.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

3D numerical simulation of the long range propagation of acoustical shock waves through a heterogeneous and moving medium

Many situations involve the propagation of acoustical shock waves through flows. Natural sources such as lightning, volcano explosions, or meteoroid atmospheric entries, emit loud, low frequency, and impulsive sound that is influenced by atmospheric wind and turbulence. The sonic boom produced by a supersonic aircraft and explosion noises are examples of intense anthropogenic sources in the atmosphere. The Buzz-Saw-Noise produced by turbo-engine fan blades rotating at supersonic speed also propagates in a fast flow within the engine nacelle. Simulating these situations is challenging, given the 3D nature of the problem, the long range propagation distances relative to the central wavelength, the strongly nonlinear behavior of shocks associated to a wide-band spectrum, and finally the key role of the flow motion. With this in view, the so-called FLHOWARD (acronym for FLow and Heterogeneous One-Way Approximation for Resolution of Diffraction) method is presented with three-dimensional applications. A scalar nonlinear wave equation is established in the framework of atmospheric applications, assuming weak heterogeneities and a slow wind. It takes into account diffraction, absorption and relaxation properties of the atmosphere, quadratic nonlinearities including weak shock waves, heterogeneities of the medium in sound speed and density, and presence of a flow (assuming a mean stratified wind and 3D turbulent ? flow fluctuations of smaller amplitude). This equation is solved in the framework of the one-way method. A split-step technique allows the splitting of the non-linear wave equation into simpler equations, each corresponding to a physical effect. Each sub-equation is solved using an analytical method if possible, and finite-differences otherwise. Nonlinear effects are solved in the time domain, and others in the frequency domain. Homogeneous diffraction is handled by means of the angular spectrum method. Ground is assumed perfectly flat and rigid. Due to the 3D aspect, the code was massively parallelized using the single program, multiple data paradigm with the Message Passing Interfaces (MPI) for distributed memory architectures. This allows us to handle problems in the order of a thousand billion mesh points in the four dimensions (3 dimensions of space plus time). The validity of the method has been thoroughly evaluated on many cases with known solutions: linear piston, scattering of plane wave by a heterogeneous sphere, propagation in a waveguide with a shear flow, scattering by a finite amplitude vortex and nonlinear propagation in a thermoviscous medium. This validation process allows for a detailed assessment of the advantages and limitations of the method. Finally, applications to atmospheric propagation of shock waves will be presented.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Structure-dependent localized surface plasmon resonance characteristics and surface enhanced Raman scattering performances of quasi-periodic nanoarrays: Measurements and analysis

A set of periodic and quasi-periodic Au nanoarrays with different morphologies have been fabricated by using electron beam lithography technique, and their optical properties have been examined experimentally and analyzed theoretically by scanning near-field optical microscope and finite element method, respectively. Results present that the localized surface plasmon resonance of the as-prepared Au nanoarrays exhibit the structure-depended characteristics. Comparing with the periodic nanoarrays, the quasi-periodic ones demonstrate stronger electric field enhancement, especially for Thue-Morse nanoarray. Meanwhile, the surface enhanced Raman scattering (SERS) spectra of 4-mercaptobenzoic acid molecular labeled nanoarrays show that the quasi-periodic nanoarrays exhibit distinct SERS enhancement, for example, a higher enhancement factor of ∼10{sup 7} is obtained for the Thue-Morse nanoarray consisted of square pillars of 100 nm size. Therefore, it is significant to optimally design and fabricate the chip-scale quasi-periodic nanoarrays with high localized electric field enhancement for SERS applications in biosensing field.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Far-field coupling in nanobeam photonic crystal cavities

We optimized the far-field emission pattern of one-dimensional photonic crystal nanobeams by modulating the nanobeam width, forming a sidewall Bragg cross-grating far-field coupler. By setting the period of the cross-grating to twice the photonic crystal period, we showed using three-dimensional finite-difference time-domain simulations that the intensity extracted to the far-field could be improved by more than three orders of magnitude compared to the unmodified ideal cavity geometry. We then experimentally studied the evolution of the quality factor and far-field intensity as a function of cross-grating coupler amplitude. High quality factor (>4000) blue (λ = 455 nm) nanobeam photonic crystals were fabricated out of GaN thin films on silicon incorporating a single InGaN quantum well gain medium. Micro-photoluminescence spectroscopy of sets of twelve identical nanobeams revealed a nine-fold average increase in integrated far-field emission intensity and no change in average quality factor for the optimized structure compared to the unmodulated reference. These results are useful for research environments and future nanophotonic light-emitting applications where vertical in- and out-coupling of light to nanocavities is required.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Review of CTF's Fuel Rod Modeling Using FRAPCON-4.0's Centerline Temperature Predictions

Coolant Boiling in Rod Arrays-Two Fluid (COBRA-TF), or CTF1, is a nuclear thermal hydraulic subchannel code used throughout academia and industry. CTF's fuel rod modeling is originally developed for VIPRE code. Its methodology is based on GAPCON and FRAP fuel performance codes, and material properties are included from MATPRO handbook. This work focuses on review of CTF's fuel rod modeling to address shortcomings in CTF's temperature predictions. CTF is compared to FRAPCON which is U.S. NRC's steady-state fuel performance code for light-water reactor fuel rods. FRAPCON calculates the changes in fuel rod variables and temperatures including the effects of cladding hoop strain, cladding oxidation, hydriding, fuel irradiation swelling, densification, fission gas release and rod internal gas pressure. It uses fuel, clad and gap material properties from MATPRO. Additionally, it has its own models for fission gas release, cladding corrosion and cladding hydrogen pickup. It allows finite difference or finite element approaches for its mechanical model. In this study, FRAPCON-4.0 is used as a reference fuel performance code. In comparison, Halden Reactor Data for IFA432 Rod 1 and Rod 3. CTF simulations are performed in two ways; informing CTF with gap conductance value from FRAPCON, and using CTF's dynamic gap conductance model. First case is chosen to show temperature is predicted correctly with CTF's models for thermal and cladding conductivities once gap conductance is provided. Latter is to review CTF's dynamic gap conductance model. These Halden test cases are selected to be representative of cases with and without any physical contact between fuel-pellet and clad while reviewing functionality of CTF's dynamic gap conductance model. Improving the CTF's dynamic gap conductance model will allow prediction of fuel and cladding thermo-mechanical behavior under irradiation, and better temperature feedbacks from CTF in transient calculations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Spectral properties of three-dimensional Anderson model

The three-dimensional Anderson model represents a paradigmatic model to understand the Anderson localization transition. In this work we first review some key results obtained for this model in the past 50 years, and then study its properties from the perspective of modern numerical approaches. Our main focus is on the quantitative comparison between the level sensitivity statistics and the level statistics. While the former studies the sensitivity of Hamiltonian eigenlevels upon inserting a magnetic flux, the latter studies the properties of unperturbed eigenlevels. We define two versions of dimensionless conductance, the first corresponding to the width of the level curvature distribution relative to the mean level spacing, and the second corresponding to the ratio of the Heisenberg time and the Thouless time obtained from the spectral form factor. We show that both conductances look remarkably similar around the localization transition, in particular, they predict a nearly identical critical point consistent with other well-established measures of the transition. We then study some further properties of those quantities: for level curvatures, we discuss particular similarities and differences between the width of the level curvature distribution and the characteristic energy studied by Edwards and Thouless (1972) in their pioneering work, in which the hopping at one lattice edge is changed from periodic to antiperiodic boundary conditions. In the context of the spectral form factor, we show that at the critical point it enters a broad time-independent regime, in which its value is consistent with the level compressibility obtained from the level variance. Finally, we test the scaling solution of the average level spacing ratio in the crossover regime using the cost function minimization approach introduced recently in Šuntajs et al. (2020). The latter approach seeks for the optimal scaling solution in the vicinity of the crossing point, while at the same time it allows for the crossing point to drift due to finite-size corrections. We find that the extracted transition point and the scaling coefficient agree with those from the literature to high numerical accuracy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Quantum efficiency enhancement in simulated nanostructured negative electron affinity GaAs photocathodes

We report nanostructured negative electron affinity GaAs photocathodes for a polarized electron source are studied using finite difference time domain optical simulation. The structures studied are nanosquare columns, truncated nanocones, and truncated nanopyramids. Mie-type resonances in the 700–800 nm waveband, suitable for generation of polarized electrons, are identified. At resonance wavelengths, the nanostructures can absorb up to 99% of the incident light. For nanosquare columns and truncated nanocones, the maximum quantum efficiency (QE) at 780 nm obtained from simulation is 27%, whereas for simulated nanopyramids, the QE is ~21%. The high photocathode quantum efficiency is due to the shift of Mie resonance toward the longer wavelength, leading to increased light absorption. The field profile distribution shows the excitation of dipole and quadrupole modes within the nanostructures at resonant frequencies. This leads to enhanced photoabsorption and photoelectron generation closer to emission surfaces than for a flat photocathode. The enhanced photoabsorption and reduced electron transport distance for the nanostructured photocathode enhance its QE compared to that for the flat surface wafer.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Development of High-Performance Cast Crankshafts. Final Technical Report

The objective of this project was to develop technologies that would enable the production of cast crankshafts that can replace high performance forged steel crankshafts. To achieve this, the Ultimate Tensile Strength (UTS) of the new material needs to be 850 MPa with a desired minimum Yield Strength (YS; 0.2% offset) of 615 MPa and at least 10% elongation. Perhaps more challenging, the cast material needs to be able to achieve sufficient local fatigue properties to satisfy the durability requirements in today’s high performance gasoline and diesel engine applications. The project team focused on the development of cast steel alloys for application in crankshafts to take advantage of the higher stiffness over other potential material choices. The material and process developed should be able to produce high-performance crankshafts at no more than 110% of the cost of current production cast units, perhaps the most difficult objective to achieve. To minimize costs, the primary alloy design strategy was to design compositions that can achieve the required properties with minimal alloying and post-casting heat treatments. An Integrated Computational Materials Engineering (ICME) based approach was utilized, rather than relying only on traditional trial-and-error methods, which has been proven to accelerate alloy development time. Prototype melt chemistries designed using ICME were cast as test specimens and characterized iteratively to develop an alloy design within a stage-gate process. Standard characterization and material testing was done to validate the alloy performance against design targets and provide feedback to material design and manufacturing process models. Finally, the project called for Caterpillar and General Motors (GM) to develop optimized crankshaft designs using the final material and manufacturing processing path developed. A multi-disciplinary effort was to integrate finite element analyses by engine designers and geometry-specific casting simulations with existing materials models to optimize crankshaft cost and performance. Prototype crankshafts of the final design were to be produced and validated using laboratory bench testing and on-engine durability testing. ICME process simulation tools were used to investigate a broad range of processing concepts. These concepts included casting orientation, various mold and core materials, and various filling and feeding strategies. Each crankshaft was first simulated without gating and risers, which is termed natural solidification. The natural solidification results were used as a baseline for strategy development of each concept. Casting process simulations and ICME tools were proven to be reasonable predictors of real world results. Potential alloys were developed that could meet the project material property goals with appropriate normalization and temper treatments. For the alloys considered, post-normalization temper treatments proved to be necessary to achieve the desired yield strengths and elongations and appropriate heat treatments were designed using ICME tools. The experimental data of all the alloys were analyzed in combination with ICME tools to establish chemistry-process-structure relations. Several GM small gas engine (SGE) crankshafts were successfully cast in sand molds using two different sprue, runner, gate, riser, chill designs. These crankshafts were cast in two different steel alloys developed during the project, but casting finishing (e.g. riser removal) remains a cost challenge. A long list of future work was left unfinished when this project was unexpectedly terminated.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Vaidya spacetime in the diagonal coordinates

We have analyzed the transformation from initial coordinates (v, r) of the Vaidya metric with light coordinate v to the most physical diagonal coordinates (t, r). An exact solution has been obtained for the corresponding metric tensor in the case of a linear dependence of the mass function of the Vaidya metric on light coordinate v. In the diagonal coordinates, a narrow region (with a width proportional to the mass growth rate of a black hole) has been detected near the visibility horizon of the Vaidya accreting black hole, in which the metric differs qualitatively from the Schwarzschild metric and cannot be represented as a small perturbation. It has been shown that, in this case, a single set of diagonal coordinates (t, r) is insufficient to cover the entire range of initial coordinates (v, r) outside the visibility horizon; at least three sets of diagonal coordinates are required, the domains of which are separated by singular surfaces on which the metric components have singularities (either g{sub 00} = 0 or g{sub 00} = ∞). The energy–momentum tensor diverges on these surfaces; however, the tidal forces turn out to be finite, which follows from an analysis of the deviation equations for geodesics. Therefore, these singular surfaces are exclusively coordinate singularities that can be referred to as false fire-walls because there are no physical singularities on them. We have also considered the transformation from the initial coordinates to other diagonal coordinates (η, y), in which the solution is obtained in explicit form, and there is no energy–momentum tensor divergence.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Evaluation of the ring tension test (RTT) for robust determination of material strengths

The ring tension test (RTT) is an experimental method for determining mechanical behavior in a material’s circumferential or hoop direction. It is a crucial test for testing anisotropic materials with tube geometry, such as nuclear fuel cladding or irradiated pipes. Several RTT configurations exist, each with their own advantages and disadvantages. However, this test is significantly more complex than traditional tensile testing and can be especially sensitive to small differences and inconsistencies in the test setup and geometry, ultimately affecting the derived mechanical properties. Previous research has focused on method development, and little work has been done on understanding the subtle differences between an ideal test and experiments, specifically when the tests are performed on highly irradiated materials in hot cells. In this work, a finite element-based investigation of the RTT is conducted to fill this gap. The two most used test configurations are investigated, comparing their ability to determine accurate material strengths through plastic deformation. Several non-ideal conditions and uncontrollable effects which are likely to occur during experimental testing such as machining tolerances, variations of specimen geometry from nominal dimensions, rotation of specimens and fixturing, and other test setup discrepancies are studied. The sensitivity of measured strengths to these conditions is presented. A mechanics-based approach to describing and correcting raw data to determine actual strengths is also included for one of the configurations, resulting in a robust correction method with highly accurate material strength measurements. Finally, based on these analyses, the hemicylindrical mandrel configuration is recommended with a gauge region oriented at a 45° angle.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Modelling the brittle failure of graphite induced by the controlled impact of runaway electrons in DIII-D

The thermo-mechanical response of an ATJ graphite sample to controlled runaway electron (RE) dissipation, realized in DIII-D, is modelled with a novel work-flow that features the RE orbit code KORC, the Monte Carlo particle transport code Geant4 and the finite element multiphysics software COMSOL. KORC provides the RE striking positions and momenta, Geant4 calculates the volumetric energy deposition and COMSOL simulates the thermoelastic response. Brittle failure is predicted according to the maximum normal stress criterion, which is suitable for ATJ graphite owing to its linear elastic behavior up to fracture and its isotropic mechanical properties. Measurements of the conducted energy, damage topology, explosion timing and blown-off material volume, impose a number of empirical constraints that suffice to distinguish between different RE impact scenarios and to identify RE parameters which provide the best match to the observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Reverse Monte Carlo modeling for low-dimensional systems

Reverse Monte Carlo (RMC) is one of the commonly used approaches for modeling total scattering data. However, to extend the capability of the RMC method for refining the structure of nanomaterials, the dimensionality and finite size need to be considered when calculating the pair distribution function (PDF). To achieve this, the simulation box must be set up to remove the periodic boundary condition in one, two or three of the dimensions. This then requires a correction to be applied for the difference in number density between the real system and the simulation box. In certain circumstances an analytical correction for the uncorrelated pairings of atoms is also applied. The validity and applicability of our methodology is demonstrated by applying the algorithms to simulate the PDF patterns of carbon systems with various dimensions, and also by using them to fit experimental data of CuO nanoparticles. This alternative approach for characterizing the local structure of nano-systems with the total scattering technique will be made available via the RMCProfile package. The theoretical formulation and detailed explanation of the analytical corrections for low-dimensional systems – 2D nanosheets, 1D nanowires and 0D nanoparticles – is also given.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Loop corrections in Minkowski spacetime away from equilibrium. Part II. Finite-time results

Loop corrections to finite-time correlation functions in quantum field theories away from equilibrium can be calculated using the in-in path integral approach. In this paper, we calculate the unequal-time two-point correlator for different massless self-interacting scalar quantum field theories on a Minkowski background, starting the field evolution at an arbitrary initial time. We find the counterterms that need to be added to UV-renormalize the result, including usual in-out counterterms in the dynamics and additional initial state counterterms that are required to cancel all UV divergences. We find that the late-time limit of the renormalized correlation function exhibits a linear or logarithmic growth in time, depending on whether the interaction strength is dimension-one or dimensionless, respectively. The late-time correlations match those obtained in our companion paper and, as shown there, the divergences do not indicate a real IR issue, consistent with what one would expect in Minkowski.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

On the questions of asymptotic recoverability of information and subsystems in quantum gravity

A longstanding question in quantum gravity regards the localization of quantum information; one way to formulate this question is to ask how subsystems can be defined in quantum-gravitational systems. The gauge symmetry and necessity of solving the gravitational constraints appear to imply that the answers to this question here are different than in finite quantum systems, or in local quantum field theory. Specifically, the constraints can be solved by providing a “gravitational dressing” for the underlying field-theory operators, but this modifies their locality properties. It has been argued that holography itself may be explained through this role of the gauge symmetry and constraints, at the nonperturbative level, but there are also subtleties in constructing a holographic map in this approach. There are also claims that holography is implied even by perturbative solution of the constraints. This short note provides further examination of these questions, and in particular investigates to what extent perturbative or nonperturbative solution of the constraints implies that information naïvely thought to be localized can be recovered by asymptotic measurements, and the relevance of this in defining subsystems. In the leading perturbative case, the relevant effects are seen to be exponentially suppressed and asymptotically vanishing, for massive fields. These questions are, for example, important in sharply characterizing the unitarity problem for black holes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Axial acoustic radiation force on a sphere in Gaussian field

Based on the finite series method, the acoustical radiation force resulting from a Gaussian beam incident on a spherical object is investigated analytically. When the position of the particles deviating from the center of the beam, the Gaussian beam is expanded as a spherical function at the center of the particles and the expanded coefficients of the Gaussian beam is calculated. The analytical expression of the acoustic radiation force on spherical particles deviating from the Gaussian beam center is deduced. The acoustic radiation force affected by the acoustic frequency and the offset distance from the Gaussian beam center is investigated. Results have been presented for Gaussian beams with different wavelengths and it has been shown that the interaction of a Gaussian beam with a sphere can result in attractive axial force under specific operational conditions. Results indicate the capability of manipulating and separating spherical spheres based on their mechanical and acoustical properties, the results provided here may provide a theoretical basis for development of single-beam acoustical tweezers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗