-
Frequency Control and Disturbance Containment Using Grid-Forming Embedded Storage Networks
Authors:
Kaustav Chatterjee,
Ramij Raja Hossain,
Sai Pushpak Nandanoori,
Soumya Kundu,
Diane Baldwin,
Ronald Melton
Abstract:
The paper discusses fast frequency control in bulk power systems using embedded networks of grid-forming energy storage resources. Differing from their traditional roles of regulating reserves, the storage resources in this work operate as fast-acting grid assets shaping transient dynamics. The storage resources in the network are autonomously controlled using local measurements for distributed fr…
▽ More
The paper discusses fast frequency control in bulk power systems using embedded networks of grid-forming energy storage resources. Differing from their traditional roles of regulating reserves, the storage resources in this work operate as fast-acting grid assets shaping transient dynamics. The storage resources in the network are autonomously controlled using local measurements for distributed frequency support during disturbance events. Further, the grid-forming inverter systems interfacing with the storage resources, are augmented with fast-acting safety controls designed to contain frequency transients within a prescribed tolerance band. The control action, derived from the storage network, improves the frequency nadirs in the system and prevents the severity of a disturbance from propagating far from the source. The paper also presents sensitivity studies to evaluate the impacts of storage capacity and inverter controller parameters on the dynamic performance of frequency control and disturbance localization. The performance of the safety-constrained grid-forming control is also compared with the more common grid-following control. The results are illustrated through case studies on an IEEE test system.
△ Less
Submitted 18 October, 2024;
originally announced October 2024.
-
Coordinated Frequency Regulation in Grid-Forming Storage Network via Safety-Consensus
Authors:
Ramij Raja Hossain,
Kaustav Chatterjee,
Sai Pushpak Nandanoori,
Soumya Kundu,
Laurentiu Marinovici,
Karan Kalsi,
Diane Baldwin
Abstract:
Inverter-based storages are poised to play a prominent role in future power grids with massive renewable generation. Grid-forming inverters (GFMs) are emerging as a dominant technology with synchronous generators (SG)-like characteristics through primary control loops. Advanced secondary control schemes, e.g., consensus algorithms, allow GFM-interfaced storage units to participate in frequency reg…
▽ More
Inverter-based storages are poised to play a prominent role in future power grids with massive renewable generation. Grid-forming inverters (GFMs) are emerging as a dominant technology with synchronous generators (SG)-like characteristics through primary control loops. Advanced secondary control schemes, e.g., consensus algorithms, allow GFM-interfaced storage units to participate in frequency regulations and restore nominal frequency following grid disturbances. However, it is imperative to ensure transient frequency excursions do not violate critical safety limits while the grid transitions from pre- to post-disturbance operating point. This paper presents a hierarchical safety-enforced consensus method -- combining a device-layer (decentralized) transient safety filter with a secondary-layer (distributed) consensus coordination -- to achieve three distinct objectives: limiting transient frequency excursions to safe limits, minimizing frequency deviations from nominal, and ensuring coordinated power sharing among GFM-storage units. The proposed hierarchical (two-layered) safety-consensus technique is illustrated using a GFM-interfaced storage network on an IEEE 68-bus system under multiple grid transient scenarios.
△ Less
Submitted 18 October, 2024;
originally announced October 2024.
-
Grid-Forming Storage Networks: Analytical Characterization of Damping and Design Insights
Authors:
Kaustav Chatterjee,
Ramij Raja Hossain,
Sai Pushpak Nandanoori,
Soumya Kundu,
Subhrajit Sinha,
Diane Baldwin,
Ronald Melton
Abstract:
The paper presents a theoretical study on small-signal stability and damping in bulk power systems with multiple grid-forming inverter-based storage resources. A detailed analysis is presented, characterizing the impacts of inverter droop gains and storage size on the slower eigenvalues, particularly those concerning inter-area oscillation modes. From these parametric sensitivity studies, a set of…
▽ More
The paper presents a theoretical study on small-signal stability and damping in bulk power systems with multiple grid-forming inverter-based storage resources. A detailed analysis is presented, characterizing the impacts of inverter droop gains and storage size on the slower eigenvalues, particularly those concerning inter-area oscillation modes. From these parametric sensitivity studies, a set of necessary conditions are derived that the design of droop gain must satisfy to enhance damping performance. The analytical findings are structured into propositions highlighting potential design considerations for improving system stability. The findings are illustrated via numerical studies on an IEEE 68-bus grid-forming storage network.
△ Less
Submitted 5 September, 2024;
originally announced September 2024.
-
Off-axis Hartmann wavefront sensing for the GMT-Consortium Large Earth Finder (G-CLEF) red camera optics
Authors:
Matthew C. H. Leung,
Colby A. Jurgenson,
Andrew Szentgyorgyi,
Brian McLeod,
Cem Onyuksel,
Joseph Zajac,
David Charbonneau,
William Podgorski,
Abigail Unger,
Mark Mueller,
Matthew Smith,
Daniel Baldwin,
V. Ashley Villar
Abstract:
The Hartmann test is a method used to measure the wavefront error in a focal optical system, wherein a mask with a pattern of small holes is placed at the system's aperture stop. By taking an image at a defocused plane, the differences between the ideal and real positions of the reimaged holes (called the transverse ray aberrations) can be measured, which can then be used to estimate the wavefront…
▽ More
The Hartmann test is a method used to measure the wavefront error in a focal optical system, wherein a mask with a pattern of small holes is placed at the system's aperture stop. By taking an image at a defocused plane, the differences between the ideal and real positions of the reimaged holes (called the transverse ray aberrations) can be measured, which can then be used to estimate the wavefront error. However, the Hartmann test is usually used with an on-axis field. In this paper, we present a wavefront sensing method which generalizes the classical Hartmann test for off-axis field angles and arbitrary reference wavefronts. Our method involves taking images at two defocused planes, and then using the real reimaged hole positions on both planes to estimate the trajectories of rays from the system's exit pupil, at which the reference wavefront is situated. We then propagate the rays forward from the reference wavefront to one of the two defocused planes, in order to find the ideal reimaged hole positions, from which we can compute transverse ray aberrations. We derive and solve a pair of nonlinear partial differential equations relating transverse ray aberrations to wavefront error, using Zernike decomposition and nonlinear least squares. Our method has been verified on simulated data from the 7-lens f/2.25 red camera system of the GMT-Consortium Large Earth Finder (G-CLEF), a high resolution optical echelle spectrograph which will be a first light instrument for the Giant Magellan Telescope (GMT).
△ Less
Submitted 29 July, 2024;
originally announced July 2024.
-
Answering real-world clinical questions using large language model based systems
Authors:
Yen Sia Low,
Michael L. Jackson,
Rebecca J. Hyde,
Robert E. Brown,
Neil M. Sanghavi,
Julian D. Baldwin,
C. William Pike,
Jananee Muralidharan,
Gavin Hui,
Natasha Alexander,
Hadeel Hassan,
Rahul V. Nene,
Morgan Pike,
Courtney J. Pokrzywa,
Shivam Vedak,
Adam Paul Yan,
Dong-han Yao,
Amy R. Zipursky,
Christina Dinh,
Philip Ballentine,
Dan C. Derieg,
Vladimir Polony,
Rehan N. Chawdry,
Jordan Davies,
Brigham B. Hyde
, et al. (2 additional authors not shown)
Abstract:
Evidence to guide healthcare decisions is often limited by a lack of relevant and trustworthy literature as well as difficulty in contextualizing existing research for a specific patient. Large language models (LLMs) could potentially address both challenges by either summarizing published literature or generating new studies based on real-world data (RWD). We evaluated the ability of five LLM-bas…
▽ More
Evidence to guide healthcare decisions is often limited by a lack of relevant and trustworthy literature as well as difficulty in contextualizing existing research for a specific patient. Large language models (LLMs) could potentially address both challenges by either summarizing published literature or generating new studies based on real-world data (RWD). We evaluated the ability of five LLM-based systems in answering 50 clinical questions and had nine independent physicians review the responses for relevance, reliability, and actionability. As it stands, general-purpose LLMs (ChatGPT-4, Claude 3 Opus, Gemini Pro 1.5) rarely produced answers that were deemed relevant and evidence-based (2% - 10%). In contrast, retrieval augmented generation (RAG)-based and agentic LLM systems produced relevant and evidence-based answers for 24% (OpenEvidence) to 58% (ChatRWD) of questions. Only the agentic ChatRWD was able to answer novel questions compared to other LLMs (65% vs. 0-9%). These results suggest that while general-purpose LLMs should not be used as-is, a purpose-built system for evidence summarization based on RAG and one for generating novel evidence working synergistically would improve availability of pertinent evidence for patient care.
△ Less
Submitted 29 June, 2024;
originally announced July 2024.
-
Ricci Flat Metrics, Flat Connections and $G_{2}$-Manifolds
Authors:
Bobby Samir Acharya,
Daniel Andrew Baldwin
Abstract:
Inspired by considerations in $M$-theory, we prove the equivalence between the moduli spaces of (suitably complexified) torsion free $G_{2}$-structures on 7-manifolds which are families of hyperKähler ALE 4-manifolds fibered over compact flat 3-manifolds $T^{3}/K$ and certain moduli spaces of flat $G^{\mathbb{C}}_{ADE}$ connections on $T^{3}/K$.
Inspired by considerations in $M$-theory, we prove the equivalence between the moduli spaces of (suitably complexified) torsion free $G_{2}$-structures on 7-manifolds which are families of hyperKähler ALE 4-manifolds fibered over compact flat 3-manifolds $T^{3}/K$ and certain moduli spaces of flat $G^{\mathbb{C}}_{ADE}$ connections on $T^{3}/K$.
△ Less
Submitted 19 December, 2023;
originally announced December 2023.
-
Coulomb and Higgs Phases of $G_2$-manifolds
Authors:
Bobby Samir Acharya,
Daniel Andrew Baldwin
Abstract:
Ricci flat manifolds of special holonomy are a rich framework as models of the extra dimensions in string/$M$-theory. At special points in vacuum moduli space, special kinds of singularities occur and demand a physical interpretation. In this paper we show that the topologically distinct $G_2$-holonomy manifolds arising from desingularisations of codimension four orbifold singularities due to Joyc…
▽ More
Ricci flat manifolds of special holonomy are a rich framework as models of the extra dimensions in string/$M$-theory. At special points in vacuum moduli space, special kinds of singularities occur and demand a physical interpretation. In this paper we show that the topologically distinct $G_2$-holonomy manifolds arising from desingularisations of codimension four orbifold singularities due to Joyce and Karigiannis correspond physically to Coulomb and Higgs phases of four dimensional gauge theories. The results suggest generalisations of the Joyce-Karigiannis construction to arbitrary ADE-singularities and higher order twists which we explore in detail in explicitly solvable local models. These models allow us to derive an isomorphism between moduli spaces of Ricci flat metrics on these non-compact $G_2$-manifolds and flat ADE-connections on compact flat 3-manifolds which we establish explicitly for $\operatorname{SU}(n)$.
△ Less
Submitted 22 September, 2023;
originally announced September 2023.
-
Tabulating Knot Mosaics: Crossing Number 10 or Less
Authors:
Aaron Heap,
Douglas Baldwin,
James Canning,
Greg Vinal
Abstract:
The study of knot mosaics is based upon representing knot diagrams using a set of tiles on a square grid. This branch of knot theory has many unanswered questions, especially regarding the efficiency with which we draw knots as mosaics. While any knot or link can be displayed as a mosaic, for most of them it is still unknown what size of mosaic (mosaic number) is necessary and how many non-blank t…
▽ More
The study of knot mosaics is based upon representing knot diagrams using a set of tiles on a square grid. This branch of knot theory has many unanswered questions, especially regarding the efficiency with which we draw knots as mosaics. While any knot or link can be displayed as a mosaic, for most of them it is still unknown what size of mosaic (mosaic number) is necessary and how many non-blank tiles (tile number) are necessary to depict a given knot or link. We implement an algorithmic programming approach to find the mosaic number and tile number of all prime knots with crossing number 10 or less. We also introduce an online repository which includes a table of knot mosaics and a tool that allows users can create and identify their own knot mosaics.
△ Less
Submitted 21 March, 2023;
originally announced March 2023.
-
A Malmquist-like bias in the inferred areas of diamond caustics and consequences for inferred time delays of gravitationally lensed quasars
Authors:
Derek Baldwin,
Paul L. Schechter
Abstract:
Quasars are quadruply lensed only when they lie within the diamond caustic of a lensing galaxy. This precondition produces a Malmquist-like selection effect in observed populations of quadruply lensed quasars, overestimating the true caustic area. The bias toward high values of the inferred logarithmic area, $\ln A_{inf}$, is proportional to the square of the error in that area, $σ^2_{\ln{A}}$. In…
▽ More
Quasars are quadruply lensed only when they lie within the diamond caustic of a lensing galaxy. This precondition produces a Malmquist-like selection effect in observed populations of quadruply lensed quasars, overestimating the true caustic area. The bias toward high values of the inferred logarithmic area, $\ln A_{inf}$, is proportional to the square of the error in that area, $σ^2_{\ln{A}}$. In effect, Malmquist's correction compensates post-hoc for a failure to incorporate a prior into parameter optimization. Inferred time delays are proportional to the square root of the inferred caustic area of the lensing galaxy. Model time delays are biased long, leading to overestimates of the Hubble constant. Crude estimates of $σ_{\ln A}$ for a sample of 13 quadruple systems give a median value of 0.16.
We identify a second effect, "inferred magnification bias,'' resulting from the combination of selection by apparent magnitude and errors in model magnification. It is strongly anti-correlated with caustic area bias, and almost always leads to underestimates of the Hubble constant. Malmquist's scheme can be adapted to priors on multiple parameters, but for quad lenses, the negative covariances between caustic area and absolute magnitude are poorly known. Inferred magnification bias may even cancel out caustic area bias, depending upon (among other things) the slope of the number magnitude relation for the sample.
Proper correction for these combined effects can, in principle, be built into Bayesian modeling schemes as priors, eliminating the need for Malmquist-style approximation, but is likely to be challenging in practice.
△ Less
Submitted 19 December, 2023; v1 submitted 12 October, 2021;
originally announced October 2021.
-
Infrared Photometric Properties of 709 Candidate Stellar Bowshock Nebulae
Authors:
Henry A. Kobulnicky,
Danielle P. Schurhammer,
Daniel J. Baldwin,
William T. Chick,
Don M. Dixon,
Daniel Lee,
Matthew S. Povich
Abstract:
Arcuate infrared nebulae are ubiquitous throughout the Galactic Plane and are candidates for partial shells, bubbles, or bowshocks produced by massive runaway stars. We tabulate infrared photometry for 709 such objects using images from the Spitzer Space Telescope (SST), Wide-Field Infrared Explorer (WISE), and Herschel Space Observatory (HSO). Of the 709 objects identified at 24 or 22 microns, 42…
▽ More
Arcuate infrared nebulae are ubiquitous throughout the Galactic Plane and are candidates for partial shells, bubbles, or bowshocks produced by massive runaway stars. We tabulate infrared photometry for 709 such objects using images from the Spitzer Space Telescope (SST), Wide-Field Infrared Explorer (WISE), and Herschel Space Observatory (HSO). Of the 709 objects identified at 24 or 22 microns, 422 are detected at the HSO 70 micron bandpass. Of these, only 39 are detected at HSO 160 microns. The 70 micron peak surface brightnesses are 0.5 to 2.5 Jy/square arcminute. Color temperatures calculated from the 24 micron to 70 micron ratios range from 80 K to 400 K. Color temperatures from 70 micron to 160 micron ratios are systematically lower, 40 K to 200 K. Both of these temperature are, on average, 75% higher than the nominal temperatures derived by assuming that dust is in steady-state radiative equilibrium. This may be evidence of stellar wind bowshocks sweeping up and heating --- possibly fragmenting but not destroying --- interstellar dust. Infrared luminosity correlates with standoff distance, R_0, as predicted by published hydrodynamical models. Infrared spectral energy distributions are consistent with interstellar dust exposed to a either single radiant energy density, U=10^3 to 10^5 (in more than half of the objects) or a range of radiant energy densities U_min=25 to U_max=10^3 to 10^5 times the mean interstellar value for the remainder. Hence, the central OB stars dominate the energetics, making these enticing laboratories for testing dust models in constrained radiation environments. SEDs are consistent with PAH fractions q_PAH <1% in most objects.
△ Less
Submitted 22 October, 2017;
originally announced October 2017.
-
Interactions and Asymptotics of Dispersive Shock Waves -- Korteweg-de Vries Equation
Authors:
Mark J. Ablowitz,
Douglas E. Baldwin
Abstract:
The long-time asymptotic solution of the Korteweg-de Vries equation for general, step-like initial data is analyzed. Each sub-step in well-separated, multi-step data forms its own single dispersive shock wave (DSW); at intermediate times these DSWs interact and develop multiphase dynamics. Using the inverse scattering transform and matched-asymptotic analysis it is shown that the DSWs merge to for…
▽ More
The long-time asymptotic solution of the Korteweg-de Vries equation for general, step-like initial data is analyzed. Each sub-step in well-separated, multi-step data forms its own single dispersive shock wave (DSW); at intermediate times these DSWs interact and develop multiphase dynamics. Using the inverse scattering transform and matched-asymptotic analysis it is shown that the DSWs merge to form a single-phase DSW, which is the `largest' one possible for the boundary data. This is similar to interacting viscous shock waves (VSW) that are modeled with Burgers' equation, where only the single, largest-possible VSW remains after a long time.
△ Less
Submitted 6 January, 2013;
originally announced January 2013.
-
Nonlinear shallow ocean wave soliton interactions on flat beaches
Authors:
Mark J. Ablowitz,
Douglas E. Baldwin
Abstract:
Ocean waves are complex and often turbulent. While most ocean wave interactions are essentially linear, sometimes two or more waves interact in a nonlinear way. For example, two or more waves can interact and yield waves that are much taller than the sum of the original wave heights. Most of these nonlinear interactions look like an X or a Y or two connected Ys; at other times, several lines appea…
▽ More
Ocean waves are complex and often turbulent. While most ocean wave interactions are essentially linear, sometimes two or more waves interact in a nonlinear way. For example, two or more waves can interact and yield waves that are much taller than the sum of the original wave heights. Most of these nonlinear interactions look like an X or a Y or two connected Ys; at other times, several lines appear on each side of the interaction region. It was thought that such nonlinear interactions are rare events: they are not. Here we report that such nonlinear interactions occur every day, close to low tide, on two flat beaches that are about 2,000 km apart. These interactions are closely related to the analytic, soliton solutions of a widely studied multi-dimensional nonlinear wave equation. On a much larger scale, tsunami waves can merge in similar ways.
△ Less
Submitted 14 August, 2012;
originally announced August 2012.
-
Migration of nuclear shell gaps studied in the d(24Ne,p gamma)25Ne reaction
Authors:
W. N. Catford,
C. N. Timis,
R. C. Lemmon,
M. Labiche,
N. A. Orr,
B. Fernandez-Dominguez,
R. Chapman,
M. Freer,
M. Chartier,
H. Savajols,
M. Rejmund,
N. L. Achouri,
N. Amzal,
N. I. Ashwood,
T. D. Baldwin,
M. Burns,
L. Caballero,
J. M. Casadjian,
N. Curtis,
G. de France,
W. Gelletly,
X. Liang,
S. D. Pain,
V. P. E. Pucknell,
B. Rubio
, et al. (4 additional authors not shown)
Abstract:
The transfer of neutrons onto 24Ne has been measured using a reaccelerated radioactive beam of 24Ne to study the (d,p) reaction in inverse kinematics. The unusual raising of the first 3/2+ level in 25Ne and its significance in terms of the migration of the neutron magic number from N=20 to N=16 is put on a firm footing by confirmation of this state's identity. The raised 3/2+ level is observed sim…
▽ More
The transfer of neutrons onto 24Ne has been measured using a reaccelerated radioactive beam of 24Ne to study the (d,p) reaction in inverse kinematics. The unusual raising of the first 3/2+ level in 25Ne and its significance in terms of the migration of the neutron magic number from N=20 to N=16 is put on a firm footing by confirmation of this state's identity. The raised 3/2+ level is observed simultaneously with the intruder negative parity 7/2- and 3/2- levels, providing evidence for the reduction in the N=20 gap. The coincident gamma-ray decays allowed the assignment of spins as well as the transferred orbital angular momentum. The excitation energy of the 3/2+ state shows that the established USD shell model breaks down well within the sd model space and requires a revised treatment of the proton-neutron monopole interaction.
△ Less
Submitted 9 April, 2010;
originally announced April 2010.
-
TIARA: a large solid angle silicon array for direct reaction studies with radioactive beams
Authors:
M. Labiche,
W. N. Catford,
R. C. Lemmon,
C. N. Timis,
R. Chapman,
N. A. Orr,
B. Fernandez-Dominguez,
G. Moores,
N. L. Achouri,
N. Amzal,
S. Appleton,
N. I. Ashwood,
T. D. Baldwin,
M. Burns,
L. Caballero,
J. Cacitti,
J. M. Casadjian,
M. Chartier,
N. Curtis,
K. Faiz,
G. de France,
M. Freer,
J. M. Gautier,
W. Gelletly,
G. Iltis
, et al. (17 additional authors not shown)
Abstract:
A compact, quasi-4pi position sensitive silicon array, TIARA, designed to study direct reactions induced by radioactive beams in inverse kinematics is described here. The Transfer and Inelastic All-angle Reaction Array (TIARA) consists of 8 resistive charge division detectors forming an octagonal barrel around the target and a set of double-sided silicon-strip annular detectors positioned at eac…
▽ More
A compact, quasi-4pi position sensitive silicon array, TIARA, designed to study direct reactions induced by radioactive beams in inverse kinematics is described here. The Transfer and Inelastic All-angle Reaction Array (TIARA) consists of 8 resistive charge division detectors forming an octagonal barrel around the target and a set of double-sided silicon-strip annular detectors positioned at each end of the barrel. The detector was coupled to the -ray array EXOGAM and the spectrometer VAMOS at the GANIL Laboratory to demonstrate the potential of such an apparatus with radioactive beams. The 14N(d,p)15N reaction, well known in direct kinematics, has been carried out in inverse kinematics for that purpose. The observation of the 15N ground state and excited states at 7.16 and 7.86 MeV is presented here as well as the comparison of the measured proton angular distributions with DWBA calculations. Transferred l-values are in very good agreement with both theoretical calculations and previous experimental results obtained in direct kinematics.
△ Less
Submitted 13 January, 2010; v1 submitted 20 August, 2009;
originally announced August 2009.
-
A Symbolic Algorithm for Computing Recursion Operators of Nonlinear PDEs
Authors:
D. E. Baldwin,
W. Hereman
Abstract:
A recursion operator is an integro-differential operator which maps a generalized symmetry of a nonlinear PDE to a new symmetry. Therefore, the existence of a recursion operator guarantees that the PDE has infinitely many higher-order symmetries, which is a key feature of complete integrability. Completely integrable nonlinear PDEs have a bi-Hamiltonian structure and a Lax pair; they can also be…
▽ More
A recursion operator is an integro-differential operator which maps a generalized symmetry of a nonlinear PDE to a new symmetry. Therefore, the existence of a recursion operator guarantees that the PDE has infinitely many higher-order symmetries, which is a key feature of complete integrability. Completely integrable nonlinear PDEs have a bi-Hamiltonian structure and a Lax pair; they can also be solved with the inverse scattering transform and admit soliton solutions of any order.
A straightforward method for the symbolic computation of polynomial recursion operators of nonlinear PDEs in (1+1) dimensions is presented. Based on conserved densities and generalized symmetries, a candidate recursion operator is built from a linear combination of scaling invariant terms with undetermined coefficients. The candidate recursion operator is substituted into its defining equation and the resulting linear system for the undetermined coefficients is solved.
The method is algorithmic and is implemented in Mathematica. The resulting symbolic package PDERecursionOperator.m can be used to test the complete integrability of polynomial PDEs that can be written as nonlinear evolution equations. With PDERecursionOperator.m, recursion operators were obtained for several well-known nonlinear PDEs from mathematical physics and soliton theory.
△ Less
Submitted 15 May, 2009;
originally announced May 2009.
-
Soliton Generation and Multiple Phases in Dispersive Shock and Rarefaction Wave Interaction
Authors:
M. J. Ablowitz,
D. E. Baldwin,
M. A. Hoefer
Abstract:
Interactions of dispersive shock (DSWs) and rarefaction waves (RWs) associated with the Korteweg-de Vries equation are shown to exhibit multiphase dynamics and isolated solitons. There are six canonical cases: one is the interaction of two DSWs which exhibit a transient two-phase solution, but evolve to a single phase DSW for large time; two tend to a DSW with either a small amplitude wave train…
▽ More
Interactions of dispersive shock (DSWs) and rarefaction waves (RWs) associated with the Korteweg-de Vries equation are shown to exhibit multiphase dynamics and isolated solitons. There are six canonical cases: one is the interaction of two DSWs which exhibit a transient two-phase solution, but evolve to a single phase DSW for large time; two tend to a DSW with either a small amplitude wave train or a finite number of solitons, which can be determined analytically; two tend to a RW with either a small wave train or a finite number of solitons; finally, one tends to a pure RW.
△ Less
Submitted 21 April, 2009;
originally announced April 2009.
-
Paradoxical popups: Why are they hard to catch?
Authors:
Michael K. McBeath,
Alan M. Nathan,
A. Terry Bahill,
David G. Baldwin
Abstract:
Even professional baseball players occasionally find it difficult to gracefully approach seemingly routine pop-ups. This paper describes a set of towering pop-ups with trajectories that exhibit cusps and loops near the apex. For a normal fly ball, the horizontal velocity is continuously decreasing due to drag caused by air resistance. But for pop-ups, the Magnus force (the force due to the ball…
▽ More
Even professional baseball players occasionally find it difficult to gracefully approach seemingly routine pop-ups. This paper describes a set of towering pop-ups with trajectories that exhibit cusps and loops near the apex. For a normal fly ball, the horizontal velocity is continuously decreasing due to drag caused by air resistance. But for pop-ups, the Magnus force (the force due to the ball spinning in a moving airflow) is larger than the drag force. In these cases the horizontal velocity decreases in the beginning, like a normal fly ball, but after the apex, the Magnus force accelerates the horizontal motion. We refer to this class of pop-ups as paradoxical because they appear to misinform the typically robust optical control strategies used by fielders and lead to systematic vacillation in running paths, especially when a trajectory terminates near the fielder. In short, some of the dancing around when infielders pursue pop-ups can be well explained as a combination of bizarre trajectories and misguidance by the normally reliable optical control strategy, rather than apparent fielder error. Former major league infielders confirm that our model agrees with their experiences.
△ Less
Submitted 30 March, 2008;
originally announced March 2008.
-
Symbolic Software for the Painleve Test of Nonlinear Ordinary and Partial Differential Equations
Authors:
Douglas Baldwin,
Willy Hereman
Abstract:
The automation of the traditional Painleve test in Mathematica is discussed. The package PainleveTest.m allows for the testing of polynomial systems of ordinary and partial differential equations which may be parameterized by arbitrary functions (or constants). Except where limited by memory, there is no restriction on the number of independent or dependent variables. The package is quite robust…
▽ More
The automation of the traditional Painleve test in Mathematica is discussed. The package PainleveTest.m allows for the testing of polynomial systems of ordinary and partial differential equations which may be parameterized by arbitrary functions (or constants). Except where limited by memory, there is no restriction on the number of independent or dependent variables. The package is quite robust in determining all the possible dominant behaviors of the Laurent series solutions of the differential equation. The omission of valid dominant behaviors is a common problem in many implementations of the Painleve test, and these omissions often lead to erroneous results. Finally, our package is compared with the other available implementations of the Painleve test.
△ Less
Submitted 2 February, 2006; v1 submitted 1 May, 2005;
originally announced May 2005.
-
Symbolic computation of hyperbolic tangent solutions for nonlinear differential-difference equations
Authors:
Douglas Baldwin,
Unal Goktas,
Willy Hereman
Abstract:
A new algorithm is presented to find exact traveling wave solutions of differential-difference equations in terms of tanh functions. For systems with parameters, the algorithm determines the conditions on the parameters so that the equations might admit polynomial solutions in tanh.
Examples illustrate the key steps of the algorithm. Parallels are drawn through discussion and example to the ta…
▽ More
A new algorithm is presented to find exact traveling wave solutions of differential-difference equations in terms of tanh functions. For systems with parameters, the algorithm determines the conditions on the parameters so that the equations might admit polynomial solutions in tanh.
Examples illustrate the key steps of the algorithm. Parallels are drawn through discussion and example to the tanh-method for partial differential equations.
The new algorithm is implemented in Mathematica. The package DDESpecialSolutions.m can be used to automatically compute traveling wave solutions of nonlinear polynomial differential-difference equations. Use of the package, implementation issues, scope, and limitations of the software are addressed.
△ Less
Submitted 1 December, 2003;
originally announced December 2003.
-
Symbolic algorithms for the Painleve test, special solutions, and recursion operators for nonlinear PDEs
Authors:
Douglas Baldwin,
Willy Hereman,
Jack Sayers
Abstract:
This paper discusses the algorithms and implementations of three Mathematica packages for the study of integrability and the computation of closed-form solutions of nonlinear polynomial PDEs.
The first package, PainleveTest.m, symbolically performs the Painleve integrability test. The second package, PDESpecialSolutions.m, computes exact solutions expressible in hyperbolic or elliptic function…
▽ More
This paper discusses the algorithms and implementations of three Mathematica packages for the study of integrability and the computation of closed-form solutions of nonlinear polynomial PDEs.
The first package, PainleveTest.m, symbolically performs the Painleve integrability test. The second package, PDESpecialSolutions.m, computes exact solutions expressible in hyperbolic or elliptic functions. The third package, PDERecursionOperator.m, generates and tests recursion operators.
△ Less
Submitted 16 November, 2003;
originally announced November 2003.
-
Symbolic computation of exact solutions expressible in hyperbolic and elliptic functions for nonlinear PDEs
Authors:
D. Baldwin,
U. Goktas,
W. Hereman,
L. Hong,
R. S. Martino,
J. Miller
Abstract:
Algorithms are presented for the tanh- and sech-methods, which lead to closed-form solutions of nonlinear ordinary and partial differential equations (ODEs and PDEs). New algorithms are given to find exact polynomial solutions of ODEs and PDEs in terms of Jacobi's elliptic functions.
For systems with parameters, the algorithms determine the conditions on the parameters so that the differential…
▽ More
Algorithms are presented for the tanh- and sech-methods, which lead to closed-form solutions of nonlinear ordinary and partial differential equations (ODEs and PDEs). New algorithms are given to find exact polynomial solutions of ODEs and PDEs in terms of Jacobi's elliptic functions.
For systems with parameters, the algorithms determine the conditions on the parameters so that the differential equations admit polynomial solutions in tanh, sech, combinations thereof, Jacobi's sn or cn functions. Examples illustrate key steps of the algorithms.
The new algorithms are implemented in Mathematica. The package DDESpecialSolutions.m can be used to automatically compute new special solutions of nonlinear PDEs. Use of the package, implementation issues, scope, limitations, and future extensions of the software are addressed.
A survey is given of related algorithms and symbolic software to compute exact solutions of nonlinear differential equations.
△ Less
Submitted 16 November, 2003; v1 submitted 7 January, 2002;
originally announced January 2002.