Curvilinear Geometry Research Articles

Master curves for unidirectional flows of FENE-P fluids in rectilinear and curvilinear geometries

We demonstrate that velocity profiles for steady, unidirectional shear flows of the FENE-P (Finitely-Extensible Nonlinear Elastic, with Peterlin closure) fluid, undergoing canonical rectilinear (pressure-driven flow in a rectangular channel or a circular pipe) or curvilinear (in Taylor–Couette or Dean configurations) flows, obey universal master curves that are a function only of the ratio Wi/L , for a fixed solvent to solution viscosity parameter β. Here, Wi is the Weissenberg number defined as the product of the dumbbell relaxation time and an appropriate shear rate, while L is the ratio of the maximum extension of the polymer to its equilibrium root-mean-square end-to-end distance. The data collapse and the resulting master curves for the velocity profile is a generalization of the recent demonstration of master curves for polymer viscosity and first normal stress coefficient for a FENE-P fluid under steady simple shear flow (Yamani and McKinley, 2023). For pressure-driven channel and pipe flows, we derive simple analytical expressions for the velocity profiles, in the high shear-rate regime of Wi/L≫1, that readily elucidate the role of finite extensibility of the polymer on the velocity profiles. In the Wi/L≫1 regime, for all the flows considered, the limit of zero solvent (β=0) is shown to be singular, owing to the absence of a high-shear plateau in the total solution viscosity, resulting in very different velocity profiles for β=0 and β→0.

Moving Morphable Components Using Strain-Based Beam Geometry Description for Topology Optimization

In the moving-morphable-component topology optimization, morphable components are introduced as a geometrical model mapped onto the background finite elements, and their shape parameters are utilized as design variables for topology optimization. Whereas a complex curved geometry ensuring C1 continuity can be generated using existing curved components, the component curvatures cannot be selected as design variables in the existing methods; thus geometric constraints associated with curvatures cannot also be directly imposed. To address this issue, this study proposes a curvature-based morphable component by introducing the curvilinear geometry representation in the strain-based beam formulation. Since the proposed component is parameterized by curvatures using the curvilinear equation, the component curvatures can be utilized as the design variables. This allows for directly imposing curvature constraints on structural members, thereby accounting for the manufacturability of an optimal topology. It is demonstrated that a symmetric placement of the design variables using the midpoint curvilinear coordinate system is critical in ensuring convergence of the proposed curvature-based component optimization. The symmetric curvature component is further extended to account for multiple curvatures within a single component while ensuring C1 continuity. Several examples are presented to demonstrate the benefits of the proposed multicurvature component for topology optimization.

Irrational choices via a curvilinear representational geometry for value

We make decisions by comparing values, but it is not yet clear how value is represented in the brain. Many models assume, if only implicitly, that the representational geometry of value is linear. However, in part due to a historical focus on noisy single neurons, rather than neuronal populations, this hypothesis has not been rigorously tested. Here, we examine the representational geometry of value in the ventromedial prefrontal cortex (vmPFC), a part of the brain linked to economic decision-making, in two male rhesus macaques. We find that values are encoded along a curved manifold in vmPFC. This curvilinear geometry predicts a specific pattern of irrational decision-making: that decision-makers will make worse choices when an irrelevant, decoy option is worse in value, compared to when it is better. We observe this type of irrational choices in behavior. Together, these results not only suggest that the representational geometry of value is nonlinear, but that this nonlinearity could impose bounds on rational decision-making.

Numerical analysis of MHD free convection in curvilinear porous enclosure with two active cylinders filled with NE-phase change material

In nuclear power industry, the energy conservation of heat exchangers is a crucial demand. The current study addresses the influence of nano-encapsulated phase change material (NEPCM) on free convection under magnetic force in a porous enclosure involving two active cylinders. The curvilinear geometry of the enclosure is set with undulated top wall. The top, horizontal and the inclined walls are set cold while the two cylinders presented the hot surface. All other walls are assumed to be adiabatically isolated. The considered parameters are Rayleigh number, Ra (103–105), Hartmann number, Ha (0–64), fusion temperature, θf (0.1–0.9), Stefan number, Ste (0.1–0.9), volume concentration, ϕ (0–0.05), Darcy number, Da (10−5–10−2). Finite element method has been used in the numerical analysis. Some of the most important results stated that fusion temperature has a significant influence on heat transfer only with high Ra number and low Ste number. Furthermore, the melting solidification region notably changed with the change of Da and θf with less influec under the magnetic force. Also, the heat transfer is enhanced with the increase of Ra, Da and ϕ while it diminished with the increase of Stephan number and MHD. It found that a promotion in the percentage of Nuavg when the volume concentration is raised from 0 to 0.05 is about 65 %, while the Nuavg dwindles by 14 % when Hartman number is raised from 0 to 64. It can be concluded that using NEPCMs is crucial for enhancing heat transfer and hence the improving the related engineering applications.

Flexural behavior of singly curved GFRP sandwich panel with polyurethane foam core

Curvilinear architecture and structures are becoming a worldwide trend for its streamlined appearance and synchronicity with the surroundings. Nonetheless, the conventional construction materials might not be able to achieve the desired curvilinear geometries in an effective and economic manner. In this regard, glass fiber-reinforced polymer (GFRP) curved sandwich panel with polyurethane foam core is able to meet the requirement of the construction efficiency and economy as well as satisfy the need for a lower energy consumption. However, few studies are available on the mechanical behavior of such panels. In this regard, this work was to conduct an experimental and numerical investigation on the singly curved GFRP sandwich panels to evaluate the flexural performance. In particular, the effect of the foam density, inclined angle, core thickness and stacking sequence of face sheet on the load-carrying capacity and failure mode was addressed. Four failure modes were identified. The singly curved GFRP sandwich panel showed an improved stiffness when compared to the flat panel, but showed more complicated failure mechanism. The stiffness of the curved panel with inclined angle 30° was 250 % higher than that of the flat panel in the test results. It was shown that the strength and stiffness increased as the density and thickness of the foam core increased. In addition, a finite element analysis was conducted via MSC.Marc and validated by the experimental results. A parametric study was performed to address the design parameters. It was found that the specimens with inclined angles of 30° and 45° showed the greatest stiffness.

Magnetopause MHD surface wave theory: progress & challenges

Sharp boundaries are a key feature of space plasma environments universally, with their wave-like motion (driven by pressure variations or flow shears) playing a key role in mass, momentum, and energy transfer. This review summarises magnetohydrodynamic surface wave theory with particular reference to Earth’s magnetopause, due to its mediation of the solar-terrestrial interaction. Basic analytic theory of propagating and standing surface waves within simple models are presented, highlighting many of the typically-used assumptions. We raise several conceptual challenges to understanding the nature of surface waves within a complex environment such as a magnetosphere, including the effects of magnetic topology and curvilinear geometry, plasma inhomogeneity, finite boundary width, the presence of multiple boundaries, turbulent driving, and wave nonlinearity. Approaches to gain physical insight into these challenges are suggested. We also discuss how global simulations have proven a fruitful tool in studying surface waves in more representative environments than analytic theory allows. Finally, we highlight strong interdisciplinary links with solar physics which might help the magnetospheric community. Ultimately several upcoming missions provide motivation for advancing magnetopause surface wave theory towards understanding their global role in filtering, accumulating, and guiding turbulent solar wind driving.

Mesoscopic lattice Boltzmann modeling of dense gas flows in curvilinear geometries

Mesoscopic lattice Boltzmann modeling of dense gas flows in curvilinear geometries

Implicit particle-in-cell development for ion source plasmas

Particle-in-Cell (PIC) codes used to study plasma dynamics within ion sources typically use an explicit scheme. These methods can be slow when simulating regions of high electron density in ion sources, which require resolving the Debye length in space and the plasma frequency in time. Recent developments on fully-implicit PIC models in curvilinear geometries have shown that these spatial/time scales can be significantly decreased/increased respectively, allowing for notable speed-ups in simulation time, and thus making it a potential tool for studying the physics of ion sources. For this purpose, a charge and energy conserving implicit PIC code has been developed in 1D to determine its potential for simulating bounded plasmas. In this paper, we use this model to simulate a 1D benchmark of a bounded plasma with fixed plasma density and electron/ion temperatures. The results are shown to compare well to the benchmark and to the results using an explicit PIC code. It is shown that the total amount of macro-particles used in the simulation is a better figure of merit for accurate results than the standard particles per cell used in literature. Significant speed-ups in computation time can be achieved for high plasma densities if the accuracy requirements are relaxed. In this case, we demonstrate the ability of the implicit PIC code to speed-up simulation time by nearly a factor of 12 compared to explicit PIC.

Elasto-inertial instability in torsional flows of shear-thinning viscoelastic fluids

It is well known that inertia-free shearing flows of a viscoelastic fluid with curved streamlines, such as the torsional flow between a rotating cone and plate or the flow in a Taylor–Couette geometry, can become unstable to a three-dimensional time-dependent instability at conditions exceeding a critical Weissenberg ( $Wi$ ) number. However, the combined effects of fluid elasticity, shear thinning and finite inertia (as quantified by the Reynolds number $Re$ ) on the onset of elasto-inertial instabilities are not fully understood. Using a set of cone–plate geometries, we experimentally explore the entire $Wi$ – $Re$ phase space for a series of nonlinear viscoelastic fluids (with the dependence on shear rate $\dot{\gamma}$ quantified using a shear-thinning parameter $\beta _P(\dot {\gamma })$ ). We tune $\beta _P(\dot {\gamma })$ by varying the dissolved polymer concentration in solution. This progressively reduces shear thinning but leads to finite inertial effects before the onset of elastic instability, and thus naturally results in elasto-inertial coupling. Time-resolved rheometric measurements and flow visualization experiments allow us to investigate the effects of flow geometry, and document the combined effects of varying $Wi, Re$ and $\beta _P(\dot {\gamma })$ on the emergence of secondary motions at the onset of instability. The resulting critical state diagram quantitatively depicts the competition between the stabilizing effects of shear thinning and the destabilizing effects of inertia. We extend the curved streamline instability criterion of Pakdel & McKinley (Phys. Rev. Lett., vol. 77, no. 12, 1996, p. 2459) for the onset of purely elastic instability in curvilinear geometries by using scaling arguments to incorporate shear thinning and finite inertial effects. The augmented condition facilitates predictions of the onset of instability over a broader range of flow conditions, and thus bridges the gap between purely elastic and elasto-inertial curved streamline instabilities.

Exact domain truncation for the Morse-ingard equations

Morse and Ingard [1] give a coupled system of time-harmonic equations for the temperature and pressure of an excited gas. These equations form a critical aspect of modeling trace gas sensors. Like other wave propagation problems, the computational problem must be closed with suitable far-field boundary conditions. Working in a scattered-field formulation, we adapt a nonlocal boundary condition proposed in [2] for the Helmholtz equation to this coupled system. This boundary condition uses a Green's formula for the true solution on the boundary, giving rise to a nonlocal perturbation of standard transmission boundary conditions. However, the boundary condition is exact and so Galerkin discretization of the resulting problem converges to the restriction of the exact solution to the computational domain. Numerical results demonstrate that accuracy can be obtained on relatively coarse meshes on small computational domains, and the resulting algebraic systems may be solved by GMRES using the local part of the operator as an effective preconditioner. These numerical results taken together combine several advanced techniques, including higher-order finite elements, geometric multigrid in curvilinear geometry, native use of complex arithmetic, and incorporation of nonlocal operators. These are tied together in a high-level simulation using the Firedrake library.

Crystal Self-Assembly under Confinement: Bridging Nanomaterials to Integrated Devices.

ConspectusSelf-assembly, a spontaneous process that organizes disordered constituents into ordered structures, has revolutionized our fundamental understanding of living matter, nanotechnology, and molecular science. From the perspective of nanomaterials, self-assembly serves as a bottom-up method for creating long-range-ordered materials. This is accomplished by tailoring the geometry, chemistry, and interactions of the components, thereby facilitating the efficient fabrication of high-quality materials and high-performance functional devices. Over the past few decades, we have seen controllable organization and diverse phases in self-assembled materials, such as organic crystals, biomolecular structures, and colloidal nanoparticle supercrystals. However, most self-assembled ordered materials and their assembly mechanisms are derived from constituents in a liquid bulk medium, where the effects of boundaries and interfaces are negligible. In the context of nanostructure patterning, self-assembly occurs in confined spaces, with feature sizes ranging from a few to hundreds of nanometers. In such settings, ubiquitous boundaries and interfaces can trap the system in a kinetically favored but metastable state, devoid of long-range order. This makes it extremely difficult to achieve ordered structures in micro/nano-patterning techniques that rely on sessile microdroplets, such as inkjet printing, dip-pen lithography, and contact printing.In stark contrast to sessile droplets, capillary bridges─formed by liquids confined between two solid surfaces─provide unique opportunities for understanding the long-range-ordered self-assembly of crystalline materials under spatial confinement. Because capillary bridges are stabilized by Laplace pressure, which is inversely proportional to the feature size, the confinement and manipulation of solutions or suspensions of functional materials at the nanoscale become accessible through the rational design of surface chemistry and geometry. Although global thermodynamic equilibrium is unattainable in evaporative systems, ordered nucleation and packing of constituent components can be locally realized at the contact line of capillary bridges. This enables the unprecedented fabrication of long-range-ordered micro/nanostructures with deterministic patterns.In this Account, we review the advancements in long-range-ordered self-assembly of crystalline micro/nanostructures under confinement. First, we briefly introduce crystalline materials characterized by strong intramolecular interactions and relatively weak intermolecular forces, analyzing both the opportunities and challenges inherent to self-assembled nanomaterials. Next, we delve into the construction and manipulation of confined liquids, focusing especially on capillary bridges controlled by engineered chemistry and geometry to regulate Laplace pressure. Through this approach, we have achieved capillary bridges with thicknesses on the order of a few nanometers and wafer-scale homogeneity, facilitating the self-assembly of ordered structures. Supported by factors such as local free-volume entropy, electrostatic interactions, curvilinear geometry, directional microfluidics, and nanoconfinement, we have achieved long-range-ordered, deterministic patterning of organic semiconductors, metal-halide perovskites, and colloidal nanocrystal superlattices using this capillary-bridge platform. These long-range microstructures serve as a bridge between nanomaterials and integrated devices, enabling emergent functionalities like intrinsic stretchability, giant photoconductivity, propagating and interacting exciton polaritons, and spin-valley-locked lasing, which are otherwise unattainable in disordered materials. Finally, we discuss potential directions for both the fundamental understanding and practical applications of confined self-assembly.

Homogenization of composite materials reinforced with unidirectional fibres with complex curvilinear cross section: a virtual element approach

<p>The paper presents an augmented curvilinear virtual element method to determine homogenized in-plane shear material moduli of long-fibre reinforced composites in the framework of asymptotic homogenization method. The new virtual element combine an exact representation of the curvilinear computational geometry for complex fibre cross section shapes through an innovative two-dimensional cubature suite for NURBS-like polygonal domains. A selection of representative numerical tests supports the accuracy and efficiency of the proposed approach for both doubly periodic and random fibre arrangement with matrix domain.</p>

Massively parallel nodal discontinous Galerkin finite element method simulator for room acoustics

We present a massively parallel and scalable nodal discontinuous Galerkin finite element method (DGFEM) solver for the time-domain linearized acoustic wave equations. The solver is implemented using the libParanumal finite element framework with extensions to handle curvilinear geometries and frequency dependent boundary conditions of relevance in practical room acoustics. The implementation is benchmarked on heterogeneous multi-device many-core computing architectures, and high performance and scalability are demonstrated for a problem that is considered expensive to solve in practical applications. In a benchmark study, scaling tests show that multi-GPU support gives the ability to simulate large rooms, over a broad frequency range, with realistic boundary conditions, both in terms of computing time and memory requirements. Furthermore, numerical simulations on two non-trivial geometries are presented, a star-shaped room with a dome and an auditorium. Overall, this shows the viability of using a multi-device accelerated DGFEM solver to enable realistic large-scale wave-based room acoustics simulations.

Relativistic superconductivity in 2D: An example

Recently, we formulated an effective scalar electrodynamics model of relativistic type-II superconductivity. We discussed possible electromagnetic behaviors of superconductivity along the lines of the nonrelativistic Ginzburg–Landau theory. This brief work considers a [Formula: see text]D condensate subjected not only to a homogeneous magnetic field but also to a static radial electric field (polar coordinates). The electric field is regarded to have two components: a central constant field and a Coulomb-like field with “charge” g, where g can take only evenly spaced real numbers. We provide the single-vortex supercurrent density together with its associated magnetic flux density for the ground state of the system. Feasible extensions of this example in [Formula: see text]D curvilinear geometries could explain some superconductivity phenomena which do not occur in planar [Formula: see text]D geometries.

Melon and rostral muscle morphology of Gervais' beaked whale (Mesoplodon europaeus): Alternating patterns of bilateral asymmetry.

Toothed whales utilize specialized nasal structures such as the lipid-rich melon to produce sound and propagate it into the aquatic environment. Very little nasal morphology of mesoplodont beaked whales has been described in the literature, and the anatomy of the melon and associated musculature of Gervais' beaked whale (Mesoplodon europaeus) remains undescribed. Heads of three (n = 3) Gervais' beaked whales were examined in detail via dissection as well as computed tomography (CT) and magnetic resonance imaging (MRI). Two additional Gervais' beaked whale individuals (n = 2) were studied via archived CT and MRI scans. Representative transverse dissection sections of the melon were processed for polarized light imaging to verify the presence of tendons inserting into the melon tissue. Three-dimensional (3D) CT reconstructions of the melon, rostral muscles, and associated structures were performed to assess morphology and spatial relationships. In all individuals, the melon's main body demonstrated a bilaterally asymmetrical, curvilinear geometry. This curvilinear shape was defined by a pattern of alternating asymmetry in the medial rostral muscles that projected into the melon's tissue. In transverse polarized light imaging, a network of tendons originating from these asymmetrical rostral muscle projections was observed permeating the melon's lipid tissue. This curvilinear melon morphology and associated asymmetrical musculature suggest a means of lengthening the lipid pathway within a relatively short dimensional footprint. In addition, the species-specific arrangement of muscular projections suggests complex fine-tuning of the melon's geometry during echolocation. Further studies may lend additional insight into the function of this unusual melon morphology.

Properties of superconducting MgB2 spherical shells deposited on 2 mm and 1 mm diameter Si3N4 spheres

Superconducting spherical shells may not only provide an appealing platform for exploring superconductivity in three-dimensional (3D) geometries with curved surfaces, but also be practically applied in important fields, such as in gyroscopes and gravimeters. In inertial confinement fusion (ICF), the perturbation on the target capsule caused by traditional contact support methods is recognized as one of the major contributors to performance degradation in ICF implosions. It was proposed recently that coating the capsule with a thin superconducting MgB2 spherical shell to realize non-contact support by means of magnetic levitation at the required temperatures K might offer a promising solution to this longstanding problem. To simulate the coating of ICF capsules, we deposited MgB2 spherical shells onto polycrystalline Si3N4 spheres (diameter d = 2 mm and 1 mm) via a hybrid physical–chemical vapour deposition technique. For both diameters, the spherical shells, of about 1 µm in thickness and covering the whole spheres completely, were polycrystalline with thin plate-shaped MgB2 crystallites oriented randomly and closely connected. The spherical shells exhibited superconducting transition at a zero resistance temperature of 38.5–39.4 K in four-probe resistance measurements and ideal diamagnetism at low temperatures in magnetization measurements, suggesting the good crystallinity and homogeneity of the shells. The upper critical field , the irreversibility field and the lower critical field were characterized and showed similar magnitudes and temperature dependencies for the d = 2 mm and 1 mm shells, yielding of 7–8 T and of 7–12 mT at 20 K. The superconducting critical current density values, evaluated from the magnetization hysteresis loops, were A cm−2 and A cm−2 at 20 K under zero and 2 T applied field, respectively, for the d = 1 mm shell. The occurrence of small flux jumps was observed at low temperatures up to 12 K and in low fields below 0.3 T. The results demonstrate the feasibility of fabricating MgB2 spherical shells with properties applicable to potential areas employing magnetic levitation such as in ICF and that MgB2 spherical shells may be exploited as an experimental system to study superconductivity in 3D curvilinear geometry.

Lattice Boltzmann model in general curvilinear coordinates applied to exactly solvable 2D flow problems

Numerical simulation results of basic exactly solvable fluid flows using the previously proposed by H. Chen Lattice Boltzmann Method (LBM) formulated on a general curvilinear coordinate system are presented. As was noted in the theoretical work of H. Chen, such curvilinear Lattice Boltzmann Method preserves a fundamental one-to-one exact advection feature in producing minimal numerical diffusion, as the Cartesian lattice Boltzmann model. As we numerically show, the new model converges to exact solutions of basic fluid flows with the increase of grid resolution in the presence of both natural curvilinear geometry and/or grid non-uniform contraction, both for near equilibrium and non-equilibrium LBM parameter choices.

Stadiums based on curvilinear geometry: Approximation of the ellipsoid offset surface

Abstract This article provides an effective method of determining the coordinates of points and the angles for various surfaces approximating the ellipsoid surface. The stands of the stadium (with a capacity of up to 82,000 seats) with rows shaped by ellipse arcs were designed (in accordance with the European Standard EN 13200-1). An algorithm was proposed to determine the optimal height of the ring steps (with the guarantee of unobstructed vision of the entire playing field from each seat of the stands). The second algorithm ensures a comfortable space in each individual seating place in the stands (deviation of the depth of the row generated by ellipse arcs did not exceed 0.000009 m at any place of the stands). This article encourages readers to program DXF files (read in the AutoCAD system) to streamline design work.

A discontinuous Galerkin method for sequences of earthquakes and aseismic slip on multiple faults using unstructured curvilinear grids

SUMMARYPhysics-based simulations provide a path to overcome the lack of observational data hampering a holistic understanding of earthquake faulting and crustal deformation across the vastly varying space–time scales governing the seismic cycle. However, simulations of sequences of earthquakes and aseismic slip (SEAS) including the complex geometries and heterogeneities of the subsurface are challenging. We present a symmetric interior penalty discontinuous Galerkin (SIPG) method to perform SEAS simulations accounting for the aforementioned challenges. Due to the discontinuous nature of the approximation, the spatial discretization natively provides a means to impose boundary and interface conditions. The method accommodates 2-D and 3-D domains, is of arbitrary order, handles subelement variations in material properties and supports isoparametric elements, that is, high-order representations of the exterior boundaries, interior material interfaces and embedded faults. We provide an open-source reference implementation, Tandem, that utilizes highly efficient kernels for evaluating the SIPG linear and bilinear forms, is inherently parallel and well suited to perform high-resolution simulations on large-scale distributed memory architectures. Additional flexibility and efficiency is provided by optionally defining the displacement evaluation via a discrete Green’s function approach, exploiting advantages of both the boundary integral and volumetric methods. The optional discrete Green’s functions are evaluated once in a pre-computation stage using algorithmically optimal and scalable sparse parallel solvers and pre-conditioners. We illustrate the characteristics of the SIPG formulation via an extensive suite of verification problems (analytic, manufactured and code comparison) for elastostatic and quasi-dynamic problems. Our verification suite demonstrates that high-order convergence of the discrete solution can be achieved in space and time and highlights the benefits of using a high-order representation of the displacement, material properties and geometries. We apply Tandem to realistic demonstration models consisting of a 2-D SEAS multifault scenario on a shallowly dipping normal fault with four curved splay faults, and a 3-D intersecting multifault scenario of elastostatic instantaneous displacement of the 2019 Ridgecrest, CA, earthquake sequence. We exploit the curvilinear geometry representation in both application examples and elucidate the importance of accurate stress (or displacement gradient) representation on-fault. This study entails several methodological novelties. We derive a sharp bound on the smallest value of the SIPG penalty ensuring stability for isotropic, elastic materials; define a new flux to incorporate embedded faults in a standard SIPG scheme; employ a hybrid multilevel pre-conditioner for the discrete elasticity problem; and demonstrate that curvilinear elements are specifically beneficial for volumetric SEAS simulations. We show that our method can be applied for solving interesting geophysical problems using massively parallel computing. Finally, this is the first time a discontinuous Galerkin method is published for the numerical simulations of SEAS, opening new avenues to pursue extreme scale 3-D SEAS simulations in the future.

PAX (Passive-Active Crossing) Method for Sub-Millimeter Coregistration of Passive Acoustic Mapping and B-Mode Images.

Nonlinear ultrasonic emissions produced during a therapeutic ultrasound procedure can be detected, localized, and quantified through a class of methods that can be referred to as passive acoustic mapping (PAM). While a variety of PAM beamforming algorithms may be employed, they share a common limitation that a single sound speed is specified for both phase steering of array elements and for calculation of source power or energy. The specified value may be inadequate whether derived from B-mode-based metrics or literature values for constituent materials. This study employed experiments and simulations with linear and curvilinear array geometries to investigate the impact of in situ sound speed uncertainties on source localization in layered media. The data were also used to evaluate a new method for optimizing coregistration of PAM and B-mode images. Coregistration errors as large as 10 mm were observed with the curvilinear array, which also showed much greater sound speed sensitivity than the linear array. Errors with both array geometries were typically reduced to the order of 0.1 mm using the proposed optimization method regardless of beamformer choice or whether the array was calibrated. In a further step toward reliable implementation of PAM, the current work provides an approach that can help ensure that therapeutic ultrasound procedures are accurately guided by cavitation emissions.