Bulletin of the American Physical Society
73rd Annual Meeting of the APS Division of Fluid Dynamics
Sunday–Tuesday, November 22–24, 2020; Virtual, CT (Chicago time)
Session J12: Multiphase Flows: Computational Methods (8:00am  8:45am CST)Interactive On Demand

Hide Abstracts 

J12.00001: Conditioned quadrature moment methods for cavitating bubble dispersions Spencer Bryngelson, Rodney Fox, Tim Colonius Bubbles dispersed in a fluid flow oscillate in response to pressure changes. These volume changes are sensitive to properties of the surrounding fluid and the bubble contents. The effect of the oscillations on the carrier fluid flow can be significant for even small void fractions. These interactions can be represented in a simulation framework via phaseaveraging. However, this model requires statistical moments of the oscillating bubble dispersion. This can be problematic: for polydisperse bubble dispersions, the computational expense of classbased simulations is dominated by these moment computations. We show that quadraturebased moment methods can alleviate this cost and provide a more general description of the evolving bubble quantities. Conditioned hyperbolic QMOM (CHyQMOM) represents and inverts the internal coordinates of the bubble dynamics model. Comparisons to Monte Carlo simulations show that this approach can represent the evolving moment system, despite extrapolation out of the moment space. Fullycoupled simulations of acousticallyexcited bubble screens are compared to classbased solutions. This comparison shows that a QMOM approach can reach a broader range of physical conditions, and thus is expected to be useful for matching experimental configurations. [Preview Abstract] 

J12.00002: A new Ghost Fluid Method for eliminating overheating errors in compressible multimedium flows Pedram Bigdelou, Praveen Ramaprabhu, Chen Liu, Prashant Tarey We describe a new version of the Ghost Fluid Method (GFM) called the Efficient GFM (EGFM), capable of eliminating overheating errors observed at fluid interfaces in compressible multimedium flows. Earlier approaches, while mitigating overheating errors to some extent, have not been able to eliminate them completely. The proposed algorithm makes use of a lemma in [1] to apply the exact boundary conditions at the interface. In addition, the isentropic and shock relations proposed in [2] are applied to fix the density values near the interface, and are also extended here to address velocity corrections. The EGFM algorithm is validated using a wide array of 1D single as well as multimedium shock tube problems and shockinterface interactions, shockbubble interaction problems and the RichtmyerMeshkov instability. When compared with the original GFM approach and its subsequent variants, the EGFM scheme proposed here is robust, and has been demonstrated to result in highly accurate solutions. $^{\mathrm{1}}$L. Xu {\&} T. Liu, J. Comput. Phys.,~\textbf{230}, 4975 (2011). $^{\mathrm{2}}$L. Xu, C. Feng {\&} T. Liu, Commun. Comput. Phys., \textbf{20(3)}, 619 (2016). [Preview Abstract] 

J12.00003: An arbitrary Lagrangian Eulerian approach with exact mass conservation for the numerical simulation of a Taylor bubble problem Mehmet Sahin, Cagatay Guventurk An arbitrary Lagrangian Eulerian approach is developed for incompressible multiphase fluid flows with exact mass conservation. The incompressible NavierStokes equations are discretized over unstructured moving quadrilateral/hexahedral meshes using the divstable sidecentered finite volume formulation. The pressure field is treated to be discontinuous across the interface. The surface tension term is considered as a force tangent to the interface. Interface kinematic boundary condition in normal direction is applied by satisfying both local and global discrete geometric conservation law [Guventurk and M. Sahin, 2017]. The resulting algebraic equations are solved in a fully coupled manner and preconditioning is performed by using a new matrix factorization proposed by Cetin and Sahin [2019]. The numerical algorithm is validated by simulating a rising Taylor bubble in a tube in two and three dimensions. The numerical results are consistent with the results of an experimental work of Bugg and Saad [2017]. [Preview Abstract] 

J12.00004: An interfacepreserving and conservative phasefield method for the surface tension dynamics Xiaoyu Mao, Vaibhav Joshi, Rajeev Jaiman We present an interfacepreserving and conservative AllenCahn phasefield formulation for the modeling of incompressible twophase flows with surface tension dynamics. To preserve the hyperbolic tangent interface profile, the mobility coefficient is adjusted adaptively as a function of gradients of the velocity and the order parameter in the diffuse interface region in such a way that the free energy minimization properly opposes the convective distortion. In the phasefield formulation, the mass conservation is achieved by enforcing a Lagrange multiplier with both temporal and spatial dependence on the phasefield function. We integrate the interfacepreserving and conservative phasefield formulation with the incompressible NavierStokes equations and the continuum surface tension force model for the simulation of incompressible twophase flows. A positivity preserving scheme designed for the boundedness and stability of the solution is employed for the variational discretization using unstructured finite elements. We examine the convergence and accuracy of the AllenCahn phasefield solver through a generic onedimensional bistable diffusionreaction system in a stretching flow. We quantify and systematically assess the relative interface thickness error and the relative surface tension force error with respect to the convective distortion parameter. Two and threedimensional rising bubble cases are further simulated to examine the effectiveness of the proposed model on the volumepreserving mean curvature flow and the interfacepreserving capability. [Preview Abstract] 

J12.00005: An approach to address numerical stiffness in the simulation of cavitating flows Mrugank Bhatt, Krishnan Mahesh A numerical approach based on the preconditioning and the dualtime stepping methodologies is proposed for LES of cavitating flows at low Mach numbers. The methodology is based on a fully compressible homogeneous mixture model as discussed in Gnanaskandan and Mahesh (IJMF. vol. 70, April 2016) which has shown promising results for capturing sheet to cloud cavitation. The goal of the present work is to allow the method to simulate wetted conditions, incipient cavitation, and cavitation inception regimes that typically require low freestream nuclei, which imply predominantly low Mach number regions in the water. The numerical formulation of the preconditioning equations is discussed in detail. The method is used to simulate (i) cavitation inception over the unsteady flow over a cylinder at Re$=$200, and (ii) the wetted conditions over a marine propeller at Re$=$894000 using LES. Overall, a significant saving in the total runtime as compared to the original solver is obtained, without compromising accuracy. [Preview Abstract] 

J12.00006: Modeling and simulation of NphaseMcomponent incompressible flows Ziyang Huang, Guang Lin, Arezoo Ardekani A consistent and conservative model is proposed for NphaseMcomponent incompressible flows, where N (N≥1) and M (M≥0) are the numbers of phases and components, respectively. Phases are immiscible with each other, while components are dissolvable in given phases. Several consistency conditions are proposed for multiphase and multicomponent flows to avoid unphysical behaviors, e.g., incompatibility between the mass and momentum equations, generating fictitious phases, local void, or overfilling, and leakage of components at phase boundaries. The mass of individual phases, the amount of each component, and the momentum of the fluid mixture are conserved by the proposed model. The model also satisfies the 2nd law of thermodynamics and Galilean invariance. A 2ndorder scheme is developed, which preserves the consistency and conservation properties of the model, and the numerical solution is demonstrated to preserve the Galilean invariance and energy law. Various challenging multiphase and multicomponent flows, including large density and/or viscosity ratios, can be effectively studied by using the proposed model and scheme. The proposed model and scheme are also applicable for multiphase problems where different phases have different miscibility with each other. [Preview Abstract] 

J12.00007: A vortexparticle mesh method for mixed rigid and soft body fluidstructure interaction Tejaswin Parthasarathy, Yashraj Bhosale, Mattia Gazzola We outline a 2D algorithm for solving incompressible flowstructure interaction problems for rigid and soft bodies, within the framework of remeshedvortex methods. We adopt a twofluid model to represent the solid and fluid phases on an Eulerian grid, separated by a diffused interface. Rigid solids are treated using the Brinkman penalization approach while a reference map technique is used to account for elastic stresses in hyperelastic solids. We then test our approach against a variety of increasingly challenging benchmarks, and further illustrate the versatility of our solver by studying elastic effects in the context of selfpropelled active swimmers, multibody contact, heat transfer and rectified viscous streaming. Through these illustrations, we showcase the ability of the solver to deal with different constitutive models and boundary conditions, solve multiphysics problems and achieve fast timetosolution by sidestepping CFL timestep restrictions. [Preview Abstract] 

J12.00008: A Numerical Formulation to Study Interactions Between Fluids and Deformable Solid in Extension to Thin Layer Geometries Jiazhen Qiao, Amir Riaz, Akash Dhruv, Elias Balaras In the present study, level set formulations are used to track solidfluid interface as well as to track a dynamic grid which captures solid deformation. The deformable solid is governed by a NeoHookean model. A unified framework of equation of motion is used to solve for both fluid and solid dynamics. FluidStructure Interaction is accounted for by adding a volumetric body force term in the solid region. This solid force is diffused into the fluid by a Heaviside function. Linear extrapolation has been used to reconstruct the dynamic grid in prevention of the distortion in the grid by artificial advection outside the solid region. Validation simulations have been performed in comparison with literatures. Implementation to thin layer geometries has also been explored by defining a thin region of solid for a few computational cells where solid shear modulus quickly goes from solid to zero. Numerical results are shown to present the robustness of this method. [Preview Abstract] 

J12.00009: A kinetic energy and entropy preserving scheme for compressible twophase turbulence simulations Parviz Moin, Suhas S Jain Compressible twophase turbulent flows have wide range of applications such as highspeed atomization, droplet combustion, bubble cavitation, cloud dynamics, RayleighTaylor and RichtmyerMeshkov instability flows. Accurate numerical modeling of such flows requires a method that is nondissipative and is stable for highdensity ratios and at high Reynolds (\textit{Re}) numbers. For a compressible flow, it is known that discrete conservation of kinetic energy is not a sufficient condition for numerical stability unlike in incompressible flows (Honein {\&} Moin, \textit{JCP}, 2004; Chandrashekar, \textit{CCP}, 2013; Kuya et al., \textit{JCP}, 2018). In this study, we adopt the recently developed diffuseinterface method (Jain, Mani {\&} Moin, \textit{JCP}, 2020) that is nondissipative. We propose discrete consistency conditions between the numerical fluxes of mass, momentum, kinetic energy, and internal energy, such that an exact conservation of kinetic energy and approximate conservation of entropy of each phase are achieved in the absence of pressure work, viscosity, and thermal diffusion. To this end, we present numerical simulations of compressible twophase turbulent flows at finite and infinite \textit{Re}, to illustrate the stability and accuracy of the method in isotropic turbulence and TaylorGreen vortex flows. [Preview Abstract] 

J12.00010: Modeling scalar transport in twophase flows with a diffuseinterface method Suhas S Jain, Ali Mani Transport of scalar quantities in twophase flows is an important problem that finds applications in wide range of natural phenomena and industrial processes. In a wide range of applications transport coefficients on the two sides of the interface can be drastically different resulting in an effective confinement of the scalar in one of the phases in the time scales of interest. This results in the formation of sharp gradients in the scalar concentration at the interface, which is a numerically challenging problem to simulate. To overcome this challenge, we developed a novel transport model (Jain {\&} Mani, \textit{CTR Annual Research Briefs}, 2019) for the simulation of scalar quantities in twophase flows with a conservative diffuseinterface method. The provable strengths of the model are that: (a) it maintains the positivity property of the scalar concentration field despite using the centraldifference scheme, (b) the transport of the scalar field is consistent with the transport of the volume fraction field, and therefore, prevents the artificial leakage of the scalar across the interface. We present numerical simulations using the proposed model in a wide range of twophase flow regimes, spanning laminar to turbulent flows. [Preview Abstract] 

J12.00011: A Deep Neural NetworkBased Approach for the Force Predictions of Particulate Multiphase Flows Majid Allahyari, W. C. Moore, Bhargav Sriram Siddani, S. Balachandar The hydrodynamics forces within random distributions of fixed monodisperse spherical particles are predicted using deep neural networkbased technique. The test data are obtained from particle resolveddirect numerical simulations. The fully resolved simulations are based on direct forcing immersed boundary method. The methodology is examined for wide ranges of Reynolds number and mean volume fraction. A data augmentation strategy is implemented to achieve 3D rotationinvariant trained network and to improve the accuracy of the predictions. To evaluate the performance of the model, direct force predictions are compared with the corresponding fully resolved solutions. [Preview Abstract] 

J12.00012: Large eddy simulation of sprays using the spectralelement method. Juan Colmenares Fernandez, Muhsin Ameen, Saumil Patel A novel approach for modeling sprays from highspeed liquid injection is presented here. The liquid phase is modeled stochastically by Lagrangian parcels, while the gas phase is modeled in an Eulerian reference frame using the spectralelement method. The coupling between gas and liquid phase is described. The gasphase solution is stabilized using a highpass filter relaxation term, which is suitable for LES using spectral methods and bypasses the need for Smagorinskytype subgrid scale models. Sensitivity of the solution to grid resolution, temporal resolution, and spray parameters is analyzed. Simulation results are compared quantitatively and qualitatively against experimental data from the Engine Combustion Network. Some of the challenges and advantages of the current approach for modeling sprays are discussed. The aim of this work is to enable highfidelity internal combustion engine simulations. [Preview Abstract] 

J12.00013: Verification of a DualScale Approach to Modeling SubFilter ShearInduced Instabilities Austin Goodrich, Marcus Herrmann A method to predict subfilter shearinduced velocities on a liquidgas material interface for use in a dualscale LESDNS method is verified with Direct Numerical Simulations. The dualscale method maintains both an LES flow solver grid and an additional Refined Local Surface Grid (RLSG) in flow solver cells that contain a material interface. The RLSG is tasked with maintaining a fully resolved realization of the material interface and transporting the interface with subgrid velocities that model subgrid physical effects. The subgrid model reconstructs velocities due the shearinduced instabilities using the OrrSommerfeld equations and appropriate boundary and interface conditions. The OrrSommerfeld equations are then solved by a Chebyshev collocation method and the results are compared against DNS solutions of the same flow conditions. The subgrid model is shown to accurately predict the dynamics of the filtered material interface in the presence of unstable subfilter corrugations within a linear regime of growth, and the model is shown to be well suited in predicting the dynamics of an interface with small random fluctuations. [Preview Abstract] 

J12.00014: Simulating Primary Atomization of a Liquid Jet in a Supersonic Crossflow Karthik Kannan, Fabian Fritz, Carlos Ballesteros, Nico Fleischmann, Marcus Herrmann Numerical simulations can provide valuable insight to understanding the atomization process of fuel within the combustor of a scramjet engine. A highfidelity extension to a diffuseinterface approach for simulating compressible multiphase flows (Garrick et al., 2017) was recently developed (Fritz et al., 2019) targeting the accurate simulation of primary atomization. To this end, a novel, unstructured, cellbased adaptive mesh refinement (AMR) framework (Ballesteros, 2019) is used to achieve higher mesh resolution in regions of liquid and shock discontinuities. Furthermore, low dissipation, high spatial accuracy is obtained by using a WENOZ (Borges et al., 2008) reconstruction, while the numerical smearing of the material interface is controlled using a THINC reconstruction scheme (Shyue and Xiao, 2014). Curvature is computed using a stretched variant of the standard height function method (Cummins et al., 2005) to account for the smearing of the material interface. Simulations of liquid jet in supersonic crossflow are performed to study the jet penetration height, droplet size and velocities in the near field of the injected plume. [Preview Abstract] 

J12.00015: Multicomponent Nearcritical Flow Simulations: Reducing Spurious Pressure Oscillations Bradley Boyd, Dorrin Jarrahbashi Spurious pressure oscillations are the most common numerical instabilities observed in multiphase and multicomponent flow problems near the critical point. A diffuseinterface model is developed to simulate transcritical mixing in multispecies and multiphase systems where spurious pressure oscillations are problematic. To reduce the spurious pressure oscillations, three methods have been proposed and implemented: (1) artificially thickening the interface between different species, (2) reconstruction of the primitive variables in the characteristics space, and (3) developing a hybrid method that switches between quasiconservative doubleflux and the classical fullyconservative numerical procedures based on the changes in the effective specific heat ratio and the effective reference internal energy. Characteristicwise reconstruction tends to reduce the spurious pressure oscillations compared to primitivewise reconstruction. The hybrid model is found to effectively reduce the magnitude of the spurious pressure oscillations, the loss of energy conservation, and provides more accurate results for all tested cases. [Preview Abstract] 

J12.00016: Eulerian Label Advection method for bubble tracking in twofluid NavierStokes simulations Declan Gaylo, Kelli Hendrickson, Dick Yue The dynamics of bubbles near the airsea interface have an important impact on many freesurface entrainment driven problems including airsea gas exchange, waves in the surfzone, and wake around a vessel. The dynamics are driven by breakup, coalescence, degassing, entrainment, and dissolution, and understanding these processes individually in numerical simulations requires a reliable method to track individual bubbles through time. Tracking methods that solve an inverse problem based Lagrangian properties can be problematic in the presence of multiple processes, which is inherent near the airsea interface. We present a forward method that uses an Eulerian advection step built upon the conservative Volume of Fluid (cVOF) method to identify the sources of fluid in a grid cell based on labels previously assigned using informed connectedcomponent labeling (ICL). Eulerian Label Advection (ELA) allows us to identify a bubble's history from its air sources, making Lagrangian tracking achievable using a forward Eulerian method. We demonstrate using 3D numerical simulations of twofluid turbulent flow that ELA is able to track bubbles in the presence of multiple processes, ultimately enabling accurate and independent quantification of the processes involved bubble dynamics. [Preview Abstract] 

J12.00017: Methodology for DNS Datadriven Machine Learning Bubble Drag Model and Its Integration to OpenFOAM ChengKai Tai, Ilya Evdokimov, Fabian Schlegel, Dirk Lucas, Igor Bolotnov This work aims to develop a twophase DNS datadriven bubble drag model and to implement it into a multiphase CFD simulation. To achieve the goal, a Tensorflow (TF)OpenFOAM(OF) integration interface has been established. The interface can call and make machine learning model predict quantities of interest on the fly. A benchmark case for the bubble drag coefficient is proposed to validate the interface. A feedforward neural network was utilized to approximate the drag correlation (Tomiyama et al., 1998) using artificially generated data. Results of the integration show good consistency in radial void fraction and velocity profiles. Next, actual DNS bubble tracking datasets are used as a data source (Fang et al., 2017, Cambareri et al., 2019). The data segments where bubbles have quasistable mainstream velocity were filtered out for drag coefficient calculation. The DNSinformed model predicts drag coefficient by taking bubble Reynolds and Eötvös number as input. The model is applied in an EulerEuler twophase flow simulation of a bubbly pipe flow in OF. The closure terms, except the drag model, utilize the baseline model of Liao et al. (2020). The results of radial void fraction and velocity profiles are compared to reference by baseline. [Preview Abstract] 

J12.00018: A multiscale numerical model for the multiphase flow in porous media Zhipeng Qin, Soheil Esmaeilzadeh, Amir Riaz, Hamdi Tchelepi Accurate prediction of multiphase flow dynamics in porous media is difficult due to the challenges associated with capturing the interfacial deformations, resolving the subgrid thin films, and accounting for the complex solid geometries. Here, we propose an accurate and efficient multiscale numerical framework for modeling the pore scale multiphase flow in the presence of complex solid geometries and thin films. Within this framework, we solve the NavierStokes equations using a projection method with approximate factorization and fractional time stepping on a staggered Cartesian grid, and employ a topologically preserved level set method to capture the evolution of immiscible fluidfluid interfaces. In order to account for the effect of complex solid boundaries and capture the evolution of thin films, we couple an immersed boundary method based on a direct forcing approach with a subgrid thinfilm model. We evaluate the accuracy of our multiscale framework by studying the immiscible oilwater displacement in a convergingdiverging capillary tube for both drainage and imbibition scenarios, and compare the results with a fully resolved solution on a much finer grid resolution. [Preview Abstract] 

J12.00019: A PhaseField Method for Computing ElastoCapillary Flows of Liquid Crystals Mingfeng Qiu, James Feng, JeanChristophe Loudet We propose a phasefield model to compute elastocapillary flows of nematic liquid crystals. The formulation provides a consistent description of nematic microstructure, in particular topological defects, and recovers macroscopic surface tension and liquid crystal anchoring stress. This is made possible by incorporating the BerisEdwards theory for nematic hydrodynamics based on a tensor order parameter in a phasefield formalism approximating the sharpinterface limit. We apply the method to the problem of drop retraction in the presence of a nematicisotropic interface. We characterize a variety of different cases and examine their dynamics. In this regime, our calculation uncovers a proportional relationship between steadystate drop deformation and the elastocapillary length, signifying competition between bulk distortional elasticity and surface tension, mediated by the anchoring energy. The new computational framework opens doors to a large class of fundamental problems concerning colloidal interaction in coupled elastocapillary fields. [Preview Abstract] 

J12.00020: Solving compressible, multiphase flows with pressure projection and relaxation Michael Kuhn, Olivier Desjardins Simulating compressible, multiphase flows presents unique challenges in terms of stability, with a wide range of scales and discontinuities from shocks and material interfaces. A wellrounded solver must address how to represent interfaces, transport discontinuous variables, and limit numerical dissipation. Also, since the sound speed and Mach number can vary significantly, pressure projection schemes are useful in providing stability and computational efficiency by removing the acoustic limit on the time step. However, the application of projection methods to these flows is limited in the literature, leaving open questions as to how the pressure in each phase should contribute to the pressure projection equation. Pressure relaxation frameworks described by Saurel et al. (JCP 2009) and Pelanti and Shyue (JCP 2014) provide a straightforward approach to reconciling the pressures in each phase, but these have been exclusively applied in Riemannbased solvers. In this work, we integrate pressure projection and pressure relaxation in the same framework, enabling improvements in stability via consistency. We demonstrate the capabilities of our scheme by simulating several test cases, including underwater bubblewave interactions and a liquid jet in supersonic crossflow. [Preview Abstract] 

J12.00021: Simulation and Modeling of Droplet Collisions using a FilmCapturing VOF Method Robert Chiodi, Olivier Desjardins Due to the presence of the liquidgas interface and corresponding discontinuities in physical properties, simulating multiphase flows can be difficult. Furthermore, interface instabilities and their interaction with the surrounding flow field can occur over a range of temporal and spatial scales spanning multiple orders of magnitude. Even with improvements in computing power, it is not possible to directly resolve all relevant scales in the flow and subgrid scale models are needed to enable the simulation of realistic multiphase systems. In this talk, we discuss the modeling of droplet coalescence and rebound during droplet collisions. These models are developed in the context of a novel interface reconstruction algorithm for geometric volume of fluid methods, R2P, which is capable of capturing subgrid scale fluid films. First, we implement existing droplet collision models for use with R2P and compare their results to experiments. We then go on to discuss opportunities for physicsbased modeling informed by resolved flowfeatures and filmmetrics extracted from the R2P interface reconstruction. [Preview Abstract] 

J12.00022: Surface Tension Advancements to a DualScale Approach for Turbulent Phase Interfaces Dominic Kedelty, Marcus Herrmann Direct Numerical Simulation remains an expensive task for atomization simulations. To decrease the burden of DNS, a dualscale modeling approach (Gorokhovski and Herrmann, 2008) that describes turbulent phase interface dynamics in a Large Eddy Simulation spatial filtering context is proposed. Spatial filtering of the equations of fluid motion introduce several subfilter terms that require modeling. Instead of developing individual closure models for the interface associated terms, the dualscale approach uses an exact closure by explicitly filtering a fully resolved realization of the phase interface. This resolved realization is maintained using a Refined Local Surface Grid approach (Herrmann, 2008) employing an unsplit geometric VolumeofFluid method (Owkes and Desjardins, 2014). Advection of the phase interface on this DNS scale requires a reconstruction of the fully resolved interface velocity. In this work, results from the dualscale LES model employing subfilter turbulent eddy reconstruction by combined approximate deconvolution and nonlinear spectral enrichment (Bassenne et al. 2019) and subfilter surface tension model (Herrmann 2013) are compared to DNS results for a phase interface in a homogeneous isotropic turbulent flow at two different Weber numbers. [Preview Abstract] 

J12.00023: Paraboloidbased models of liquidgas interfaces from volume fraction data Austin Han, Robert Chiodi, Fabien Evrard, Olivier Desjardins In this talk, we investigate the use of semianalytical volume integration for extracting a paraboloid model of a captured interface in the context of volume of fluid (VOF). In contrast with the heightfunction (HF) method, which also represents the interface as a paraboloid, we allow for the paraboloid to be arbitrarily rotated, making possible the application of the volume integration to unstructured meshes. We evaluate the performance of the paraboloid model on a range of canonical problems and discuss the calculation of error metrics such as interface normal, curvature, and neighborcell volume fraction. Finally, we investigate the robustness of the model to volume fraction error and the reliability of piecewiselinear interface reconstruction (PLIC) information used in the formation of the paraboloid model. [Preview Abstract] 
Follow Us 
Engage
Become an APS Member 
My APS
Renew Membership 
Information for 
About APSThe American Physical Society (APS) is a nonprofit membership organization working to advance the knowledge of physics. 
© 2021 American Physical Society
 All rights reserved  Terms of Use
 Contact Us
Headquarters
1 Physics Ellipse, College Park, MD 207403844
(301) 2093200
Editorial Office
1 Research Road, Ridge, NY 119612701
(631) 5914000
Office of Public Affairs
529 14th St NW, Suite 1050, Washington, D.C. 200452001
(202) 6628700