Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “time and motion study”

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 271 records · Page 15

Jupiter's Great Red Spot as a shallow water system

Voyager cloud-top velocity data for Jupiter's Great Red Spot (GRS) is used to derive the bottom topography up to a constant that depends on the unknown radius of deformation. The bottom topography is inferred from the Bernoulli streamfunction, kinetic energy per unit mass, and absolute vorticity values derived from the velocity data. The results are used to calculate potential vorticity versus latitude far from the vortex. It is found that the deep atmosphere is in differential motion and that the far-field potential vorticity gradient changes sign at several latitudes. Numerical experiments are conducted to study the time-dependent behavior of the shallow water analog of Jupiter's analog.

Dowling, Timothy E.↗

Insight into the Physical and Dynamical Processes that Control Rapid Increases in Total Flash Rate

Rapid increases in total lightning (also termed "lightning jumps") have been observed for many decades. Lightning jumps have been well correlated to severe and hazardous weather occurrence. The main focus of lightning jump work has been on the development of lightning algorithms to be used in real‐time assessment of storm intensity. However, in these studies it is typically assumed that the updraft "increases" without direct measurements of the vertical motion, or specification of which updraft characteristic actually increases (e.g., average speed, maximum speed, or convective updraft volume). Therefore, an end‐to‐end physical and dynamical basis for coupling rapid increases in total flash rate to increases in updraft speed and volume must be understood in order to ultimately relate lightning occurrence to severe storm metrics. Herein, we use polarimetric, multi‐Doppler, and lightning mapping array measurements to provide physical context as to why rapid increases in total lightning are closely tied to severe and hazardous weather.

Schultz, Christopher J.↗

Numerical simulations of three-dimensional ion crystal dynamics in a Penning trap using the fast multipole method

We simulate the dynamics, including laser cooling, of three-dimensional (3-D) ion crystals confined in a Penning trap using a newly developed molecular dynamics-like code. The numerical integration of the ions’ equations of motion is accelerated using the fast multipole method to calculate the Coulomb interaction between ions, which allows us to efficiently study large ion crystals with thousands of ions. In particular, we show that the simulation time scales linearly with ion number, rather than with the square of the ion number. By treating the ions’ absorption of photons as a Poisson process, we simulate individual photon scattering events to study laser cooling of 3-D ellipsoidal ion crystals. Initial simulations suggest that these crystals can be efficiently cooled to ultracold temperatures, aided by the mixing of the easily cooled axial motional modes with the low frequency planar modes. In our simulations of a spherical crystal of 1000 ions, the planar kinetic energy is cooled to several millikelvin in a few milliseconds while the axial kinetic energy and total potential energy are cooled even further. This suggests that 3-D ion crystals could be well suited as platforms for future quantum science experiments.

Zaris, John (ORCID:0009000196476323)↗

Growth of organized flow coherent motions within a single-stream shear layer: 4D-PTV measurements

Abstract This study investigates the evolution of a single-stream shear layer (SSSL) originating from a wall boundary layer past a backward-facing step. Utilizing a time-resolved 3D-Particle Tracking Velocimetry (4D-PTV) technique, we track the trajectories of fluorescent particles to gain insight into the flow characteristics of the SSSL. A compact water tunnel facility ( $$\textrm{Re}_\tau =1\,240$$ Re τ = 1 240 ) is fabricated to obtain an SSSL with a perpendicular slow entrainment stream past the separation edge. A hybrid interpolation approach that combines ensemble binning and Gaussian weighting is implemented to derive minimally filtered mean and instantaneous lower- and higher-order flow field parameters. Spanwise-dominant coherent motion accompanied by finer flow scales is observed to grow due to flow entrainment through “nibbling” actions of small-scale vortices, “engulfing” by large-scale vortices, and vortex pairing events. Furthermore, the non-zero-speed stream edge grows relatively faster than the zero-speed stream edge, showing a strong asymmetry in mixing composition across a mixing layer. The SSSL reaches self-similarity at a streamwise distance of $$\approx 55\,\theta _{0}$$ ≈ 55 θ 0 , where $$\theta _0$$ θ 0 is the initial momentum thickness from the separation edge, i.e., considerably shorter than reported in previous studies. A literature comparison of growth rate parameters raises intriguing questions regarding a potential inclusive growth scaling unifying the free shear layers. A turbulent kinetic energy (TKE) budget analysis reveals a negative production region immediately downstream of the separation edge attributed to a large positive streamwise gradient of streamwise velocity. In the self-similar region, the phase-averaged flow mapping demonstrates a larger concentration of turbulence production rate around the outer edges of spanwise vortices, specifically at the intersection of braids and vortices. Furthermore, a spatial separation exists in the regions of peak production and dissipation rates within the vortex core region favoring dissipation. The braids exhibit a larger concentration of turbulence diffusion rates, indicating their function as a conduit for exchanging turbulence between neighboring coherent motions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Pairwise‐Parallel Entangling Gates on Orthogonal Modes in a Trapped‐Ion Chain

Abstract Parallel operations are important for both near‐term quantum computers and larger‐scale fault‐tolerant machines because they reduce execution time and qubit idling. This study proposes and implements a pairwise‐parallel gate scheme on a trapped‐ion quantum computer. The gates are driven simultaneously on different sets of orthogonal motional modes of a trapped‐ion chain. This work demonstrates the utility of this scheme by creating a Greenberger‐Horne‐Zeilinger (GHZ) state in one step using parallel gates with one overlapping qubit. It also shows its advantage for circuits by implementing a digital quantum simulation of the dynamics of an interacting spin system, the transverse‐field Ising model. This method effectively extends the available gate depth by up to two times with no overhead when no overlapping qubit is involved, apart from additional initial cooling. This scheme can be easily applied to different trapped‐ion qubits and gate schemes, broadly enhancing the capabilities of trapped‐ion quantum computers.

Optics↗

Dry electrodes with a printed cellulose–graphene ink for low-profile strain sensors in electromyography

Dihydrolevoglucosenone, commonly known as Cyrene, is a renewable and fully biodegradable cellulose-waste derived, environmentally friendly solvent, presenting a non-toxic alternative to N-methyl-2-pyrrolidone (NMP). Currently, solution-based processing of graphene and other similar van der Waals solids favor toxic solvents such as NMP, limiting their use for biosensing. However, with the use of Cyrene, bio-compatible printable devices are possible, and studies have already demonstrated its use in temperature and other biosensing methods through screen-printing. Screen-printing unfortunately often requires masks that constrain the minimum acquirable feature size to be above hundreds of centimeters and wastes material, adding to process complexity and cost. Conversely, inkjet-printing is an attractive alternative for the maskless patterning of hierarchically assembled structures, with micron length scales attainable. Graphene's high conductivity positions it ideally for long-wear sensors such as dry electrodes or respiration monitors. Here, we demonstrate the potential of Cyrene-based graphene inks through few-layer inkjet printing on flexible substrates for the first time, to produce non-toxic conductors toward a strain-mediated mechanism for biosensing, used to detect bodily motion for wearable electronics. The challenges overcome in this study include engineering ink chemistry and printing parameters such that Cyrene's relatively high viscosity compared to typical inkjet solvents, still allows for droplet ejection in a conventional material printer, yielding well-resolved clean line-edges in contrast to other solvents that exhibit diffuse line-edges possibly from stray droplets and ink-splashing. Temperature-dependent transport measurements on the inkjet-printed Cyrene-based graphene films showed the conductivity to be largely temperature-invariant but at lower temperatures below 100 K, conductivity decreased, likely as a result of increased inter-membrane separation arising from thermal contraction. Additionally, temperature-dependent Raman spectroscopy showed the red-shift in the G-band, 2D-band and D-band peaks, as temperature increased. As a result, by validating flexion motion detection of the proximal interphalangeal joint demonstrated in this study, our work is the first of its kind to successfully additively manufacture inkjet-printed Cyrene-based graphene strain sensors on flexible substrates for bio-sensing and wearables.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Diagnosis of vertical motion from VAS retrievals

Satellite-derived temperature profiles are used to determine if reliable estimates of synoptic-scale vertical motion can be obtained from the adiabatic, vorticity, and omega equation techniques. The period of study contains a short-wave trough over the Midwest and a convective outbreak over the middle Mississippi River Valley. Satellite soundings are available at 1-3 h intervals at five times. The emphasis is on assessing the strengths and weaknesses of the three vertical motion procedures, and determining the effects of short-interval observations on the calculations. Results show that the quasi-geostrophic omega equation provided patterns and magnitudes most consistent with observed weather events and 12 h radiosonde-derived motions. The vorticity method produced less satisfactory results, while adiabatic motions were unacceptable. The time derivative term dominated adiabatic motions and was a major influence in the vorticity method. Unrealistic temperature tendencies resulted from the retrieval algorithm; i.e., a diurnal temperature bias extended upwards to 500 mb, and there was a compensating effect at higher levels.

Fuelberg, Henry E.↗

Correlation studies of Pioneer Venus imagery obtained from PV experiments with near IR imagery obtained from ground-based observations during Venus inferior conjunction

The purpose of this study is to attempt to find correlations between data taken by experiments aboard the Pioneer-Venus Orbiter (PVO) and those obtained from Earth-based near-infrared (NIR) measurements of Venus during periods near inferior conjunction. Since the NIR measurements have been found to provide data on the middle atmosphere cloud morphology and motion, it is assumed that any correlations will also indicate that the PVO experiments are also documenting cloud behavior. If such correlations are found, then a further task is to attempt to study the long term behavior of the cloud features implied by the correlations. Many PVO data have been obtained over an extended period extending from 1978 until the PV demise in 1992. There exists a long, somewhat ill-conditioned time series of data that may contain valuable information on the long time, as well as short term behavior of the clouds, and, derivatively from cloud motion, atmospheric dynamics and wave activity in the Venus atmosphere. For example, determination of the zonal velocities of any OCPP (Cloud Photopolarimeter) 0.935 micron features could then be used for comparisons with data from other sources to attempt to fix the altitude region in which such features existed. A further task of this study is to attempt to correlate any features found in simultaneously obtained data, for example, the OCPP 0.365 and 0.935 micron data. The existence of such correlations may imply that data was obtained in overlapping altitude regions of the atmosphere.

Ragent, Boris↗

Lunar motion analysis and laser data management

Work completed in lunar motion analysis and laser data management during the period July 1, 1971 - September 30, 1975 was reported. In this context, analysis refers to theoretical or numerical studies involving real or potential applications of such observations to improvement of the physical model, and data management refers to the process by which observed photon events are turned into observations and are made available to potential users. The data analysis work included: (1) bringing to operational status of computer programs for the numerical integration of the lunar orbit motion and for the application of lunar laser time delays for the improvement of the parameters of the physical model, (2) program improvement and program integrity, (3) three-dimensional ephemeris, and (4) miscellaneous independent studies. The data management work included: (1) data identification, (2) observatory interfaces, and (3) data distribution.

Mulholland, J. D.↗

Dynamic Nanoscale Spatial Heterogeneity in a Perovskite-to-Brownmillerite Topotactic Phase Transformation

Phase transitions are omnipresent in modern condensed matter physics and its applications. In solids, first-order phase transformations typically occur by nucleation and growth under nonequilibrium conditions. Under constant external conditions, e.g., constant annealing temperature and pressure, the nucleation and growth dynamics are often thought of as spatially and temporally independent. Here, in situ Bragg X-ray photon correlation spectroscopy (XPCS) reveals nanoscale spatial and dynamical heterogeneity in the perovskite-to-brownmillerite topotactic phase transformation in La 0.7 Sr 0.3 CoO 3 thin films annealed under constant reducing conditions over a time span of multiple hours. Specifically, a time scale associated with domain growth remains stable, with a corresponding domain wall speed of v d = 6 ± 0.5 × 10 –4 nm/s (2 ± 0.2 nm/h), while a slower time scale, associated with temperature-driven depinning of domains, leads to accelerating dynamics with time scales following an aging power law with exponent −2.2 ± 0.5. This experiment demonstrates that Bragg XPCS is a powerful tool to study nanoscale dynamics in structural phase transformations, with the ability to extract quantitative average values related to nanodomain motion in situ. Furthermore, the results are relevant for phase engineering of phase-change devices, as they show that nanoscale dynamics, linked to domain and domain-wall motion, can continuously evolve and speed up with time, even hours after the initiation of the phase transformation, with potential repercussions on electrical performance.

X-ray photon correlation spectroscopy↗

Unsteady Transonics of a Wing With Tip Store

The presence of tip stores influences both the aerodynamic and aeroelastic performances of wings. Such effects are more pronounced in the transonic regime. In this study, a theoretical method is developed, for the first time, to compute unsteady transonics of oscillating wings with tip stores. The method is based on the small-disturbance aerodynamic equations or motion from the potential-flow theory. To validate the method, subsonic and transonic aerodynamic computations are made for a wing of low aspect ratio, and they are compared with the available experimental data. The comparisons are favorable. The strong effects of the tip store on the transonic aerodynamics on the wing are also illustrated. The method developed in this study can be used for transonic aeroelastic computations of wings with tip stores.

Guruswamy, Guru P.↗

Intermediate scattering functions of a rigid body monoclonal antibody protein in solution studied by dissipative particle dynamic simulation

In the past decade, there was increased research interest in studying internal motions of flexible proteins in solution using Neutron Spin Echo (NSE) as NSE can simultaneously probe the dynamics at the length and time scales comparable to protein domain motions. However, the collective intermediate scattering function (ISF) measured by NSE has the contributions from translational, rotational, and internal motions, which are rather complicated to be separated. Widely used NSE theories to interpret experimental data usually assume that the translational and rotational motions of a rigid particle are decoupled and independent to each other. To evaluate the accuracy of this approximation for monoclonal antibody (mAb) proteins in solution, dissipative particle dynamic computer simulation is used here to simulate a rigid-body mAb for up to about 200 ns. The total ISF together with the ISFs due to only the translational and rotational motions as well as their corresponding effective diffusion coefficients is calculated. The aforementioned approximation introduces appreciable errors to the calculated effective diffusion coefficients and the ISFs. For the effective diffusion coefficient, the error introduced by this approximation can be as large as about 10% even though the overall agreement is considered reasonable. Thus, we need to be cautious when interpreting the data with a small signal change. In addition, the accuracy of the calculated ISFs due to the finite computer simulation time is also discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Acquisition and Processing Considerations for Infrared Images of Rotating Turbine Blades

Abstract As designers aim to increase efficiency in gas turbines for aircraft propulsion and power generation, spatially resolved experimental measurements are needed to validate computational models and compare improvement gains of new cooling designs. Infrared (IR) thermography is one such method for obtaining spatially resolved temperature measurements. As technological advances in thermal detectors enable faster integration times, surface temperature measurements of rotating turbine blades become possible to capture including the smallest features. This paper outlines opportunities enabled by the latest IR detector technologies for capturing spatially resolved rotating blade temperatures, while also addressing some of the challenges of implementing IR for turbine rigs such as the one in the Steady Thermal Aero Research Turbine (START) Laboratory. This paper documents critical steps in achieving accurate measurements including calibration, integration times, spatial noise, and motion blur. From these results, recommendations are provided for achieving accurate IR measurements collected in a rotating turbine facility to study film cooling.

Engineering↗

Low-frequency creep in CoNiFe films.

Domain wall motion excited by slow rise-time, bipolar, hard-axis pulses in vacuum deposited CoNiFe films from 1500 to 2000 A thick is studied. The results are consistent with those of comparable NiFe films. Furthermore, the wall coercivity is found to be the most significant sample property correlated to the low-frequency creep properties of all the samples.

Bartran, D. S.↗

Transonic aeroelasticity of wings with tip stores

The presence of tip stores influences both aerodynamic and aeroelastic performances of wings. Such effects are more pronounced in the transonic regime. In this study, transonic aeroelasticity of wings with tip stores is studied for the first time by a theoretical method using the unsteady-small disturbance transonic aerodynamic equations coupled with modal structural equations of motion. The aerodynamic and structural equations of motion are simultaneously integrated by a time-accurate numerical scheme. To validate the tip store simulation, aeroelastic computations are made for a typical rectangular wing with a tip store and results are correlated with available wind tunnel data for the corresponding wing without a tip store at various flight conditions. Aeroelastic computations are also made for a typical fighter wing with a tip store. Present computations show that it is important to account for the aerodynamics of the tip store, particularly in the transonic regime where the tip store can make the wing aeroelastically unstable.

Guruswamy, G. P.↗

Dynamics of Single Chains of Suspended Ferrofluid Particles

We present an experimental study of the dynamics of isolated chains made of super-paramagnetic particles under the influence of a magnetic field. The motivation of this work is to understand if the chain fluctuations exist and, if it does, how does the fluctuation affect chain aggregation. We find that single chains strongly fluctuate and that the characteristic frequency of their fluctuations is inversely proportional to the magnetic field strength. The higher the field the lower the characteristic frequency of the chain fluctuations. In the high magnetic field limit, chains behave like rigid rods without any internal motions. In this work, we used ferrofluid particles suspended in water. These particles do not have any intrinsic magnetization. Once a magnetic field is applied, a dipole moment is induced in each particle, proportional to the magnetic field. A dipolar magnetic interaction then occurs between particles. If dipole-dipole magnetic energy is higher than the thermal energy, the result is a structure change inside the dipolar fluid. The ratio of these two energies is expressed by a coupling constant lambda as: lambda = (pi(a(exp 3))(chi(exp 2))(mu(sub 0))(H(sub 0))(exp 2))/18kT Where a is the particle radius, mu(sub 0) is the vacuum magnetic permeability, H(sub 0) the applied magnetic field, k the Boltzmann constant and T the absolute temperature. If lambda > 1, magnetic particles form chains along the field direction. The lateral coalescence of several chains may form bigger aggregates especially if the particle volume fraction is high. While many studies and applications deal with the rheological properties and the structural changes of these dipolar fluids, this work focuses on the understanding of the chain dynamics. In order to probe the chain dynamics, we used dynamic light scattering (DLS) in self-beating mode as our experimental technique. The experimental geometry is such that the scattering plane is perpendicular to the magnetic field. Therefore, only motions in this plane are probed. A very dilute sample of a ferrofluid emulsion with a particle volume fraction of 10(exp -5) is used in this experiment. We chose such a low volume fraction to avoid multiple light scattering as well as lateral chain-chain aggregation. DLS measures the dynamic structure factor S(q,t) of the sample (q is the scattering wave vector, t is the time). In the absence of the magnetic field, identical particles of ferrofluid droplets are randomly distributed and S(q,t) reduces to exp(-q(exp 2)2D(sub 0)t). D(sub 0)=(kT/(6(pi)(eta)(a)) is the diffusion coefficient of Brownian particles (where Xi = (6(pi)(eta)(a)) is the Stokes frictional coefficient of a spherical particle in a fluid of viscosity eta). If interactions or polydispersity can not be ignored, an effective diffusion coefficient is introduced. Formally, D(sub eff) is defined as: D(sub eff) = - q(exp -2) partial derivative of (ln(S(q,t)) with respect to time, as t goes to 0. D(sub eff) reduces to D(sub 0) if no interactions and only a few particles size are present. Therefore, we can use DLS to measure particle size. The particle radius was found to be a=0.23 mu m with 7% of polydispersity. In this case, if we vary the scattering angle theta (and so q) we do not have any change in the measured diffusion coefficient: it is q-independent. When a magnetic field is applied, particles aggregate into chains if lambda > 1. We first studied the kinetics of the chain formation when lambda = 406. At a fixed scattering angle, we measured diffusion coefficient D(sub eff) as a function of time. Experimentally, we find that D(sub eff) decreases monotonously with time. Physically, this means that chains are becoming longer and longer. Since we are only sensitive to motions in the scattering plane and since chains have their main axis perpendicular to this plane, the measured diffusion coefficient is the trans-verse diffusion coefficient. We can relate D(sub eff) to the mean number of particles per chain N(t) at a given time and to the diffusion coefficient of an isolated particle D(sub 0) as D(sub eff)=f(N(t))D(sub 0). Since f(N) is known from other recent work, N can be expressed as a function of the time. We found a square root dependency: N(t) proportional to the square root of t. As expected for very low volume fraction, this behavior is characteristic of a diffusion-limited aggregation as suggested by several authors and by our previous work. In this study, we focus on the dependence of the effective diffusion coefficient on the scattering angle and the magnetic field strength. After the magnetic field is applied (lambda = 406) for a long time, typically 6 hours, kinetics of chain formation becomes very slow. Chain size does not vary much over the next hour period. Thus, we can perform different interesting experiments. First, at a fixed magnetic field, we measure the effective diffusion coefficient as a function of the scattering angle (from 5 to 130 deg). Our results show that the measured diffusion coefficient increases linearly with the scattering angle: D(sub eff) proportional to q. If we do the same experiment for different lambda values, D(sub eff) depends on lambda as D(sub eff) proportional to lambda(exp -1/2). We also find for different lambda values that the same asymptotic D(sub eff) value is obtained when q approaches zero. The angle dependency of D(sub eff) suggests that an additional motion exists besides chain drifting. Chain size is constant during experiment, which was verified by measuring the same diffusion coefficient at the beginning and at the end of the angle switching. If chains are rigid, D(sub eff) is independent of q. Therefore, we found that D(sub eff) not only measures the motion of the entire chain but also its internal fluctuations. These internal motions are the fluctuations of the particles in the chain. To understand the q dependency of D(sub eff), let us look at the probing length used. In our study, the characteristic length scale probed is l=2pi/q which is in the range of 0.9<l/a<20. When l is much larger than the particle radius (l/a going to infinity) we are mainly sensitive to the center of mass diffusion of the chain. D(sub eff) is then the diffusion of the entire chain and depends only on N but not on q and lambda. The value D(sub eff) allows us to obtain the number of particles per chain N. In the opposite limit (i.e. l/a < 1), we are sensitive to motions on the size of individual particles. The main contributions to the measured diffusion coefficient come from internal motions of the chains (i.e. particles' fluctuations inside the chain). We investigated also the effect of the chain size on D(sub eff). For the same value of the magnetic interaction lambda and different chain sizes, we found that the slope of D(sub eff) versus q decreases when N increases as a power law D(sub eff) proportional to N(exp -0.7). This is an important behavior because this means that the longer the chain the slower its internal fluctuations besides the whole chain motion.

Cutillas, S.↗

Effects of human presence on African mammal waterhole attendance and temporal activity patterns

Abstract Human impacts on the environment and wildlife populations are increasing globally, threatening thousands of species with extinction. While wildlife‐based tourism is beneficial for educating tourists, generating income for conservation efforts, and providing local employment, more information is needed to understand how this industry may impact wildlife. In this study, we used motion‐activated cameras at 12 waterholes on a private game reserve in northern Namibia to determine if the presence of humans and permanent infrastructure affected mammal visits by examining their (1) number of visits, (2) time spent, and (3) diel activity patterns. Our results revealed no differences in the number of visits based on human presence for any of the 17 mammal species studied. However, giraffes ( Giraffe camelopardalis ) spent more time at waterholes before observer presence compared to during. Additionally, several species changed diel activity patterns when human observers were present. Notably, several carnivore and ungulate species increased overlap in their activity patterns during periods while humans were present relative to when humans were absent. These modifications of mammal temporal activity patterns due to human presence could eventually lead to changes in community structure and trophic dynamics because of altered predator–prey interactions. As humans continue to expand into wildlife habitats, and wildlife‐based tourism increases globally, it is imperative that we fully understand the effects of anthropogenic pressures on mammal behavior. Monitoring of wildlife behavioral changes in response to human activity is crucial to further develop wildlife tourism opportunities in a way that optimizes the impact of conservation goals.

Patterson, J. R. [Savannah River Ecology Lab, Warn↗

Time-reversibility and particle sedimentation

This paper studies an ODE model, called the Stokeslet model, and describes sedimentation of small clusters of particles in a highly viscous fluid. This model has a trivial solution in which the n particles arrange themselves at the vertices of a regular n-sided polygon. When n = 3, Hocking and Caflisch et al. (1988) proved the existence of periodic motion (in the frame moving with the center of gravity in the cluster) in which the particles form an isosceles triangle. Here, the study of periodic and quasi-periodic solutions of the Stokeslet model is continued, with emphasis on the spatial and time-reversal symmetry of the model. For three particles, the existence of a second family of periodic solutions and a family of quasi-periodic solutions is proved. It is also indicated how the methods generalize to the case of n particles.

Golubitsky, Martin↗