Projectile motion with Mathematica

This note describes how to use the computer algebra system (CAS) Mathematica to analyse projectile motion with and without air resistance. For a projectile fired from ground level with an initial velocity ν ft/s at an angle θ degrees from the horizontal (0 < θ < 90°), it is well known that in the absence of air resistance, the projectile follows a parabolic path. However, this is not true if air resistance is taken into account. In the presence of air resistance, the equations of motion become complicated, thus making traditional handcalculation methods quite ineffective, and a powerful CAS such as Mathematica becomes an invaluable tool to better understand projectile motion. The note discusses how Mathematica can be used to create simulated experiments of projectiles with and without air resistance. These experiments result in several conjectures, leading to theorems.     

Zur Berechnung der Geschossabweichung unter dem Einfluss eines Seitenwindes

Summary The first order correction which should be made to thez-coordinate of the standard trajectory of a projectile to account for a cross-wind is given by (1). It is shown that this formula is a very good approximation also for angles of departure near 90°.     

Analysis of projectile motion in view of fractional calculus

The projectile motion is examined by means of the fractional calculus. The fractional differential equations of the projectile motion are introduced by generalizing Newton’s second law and Caputo’s fractional derivative is considered to use the physical initial conditions. In the absence of air resistance it is found that at certain conditions, the range and the maximum height of the projectile obtained by using the fractional calculus give the same results of the classical calculus. It is also found that, unlike the classical projectile motion, the launching angle that maximizes the horizontal range is a function of the arbitrary order of the fractional derivative α. Moreover, in a resistant medium, the parametric equations are expressed in terms of Mittag-Leffler function and the results agree with those of the classical projectile as α → 2. Moreover, the trajectories of the projectile are discussed in graphs and compared with those of the classical calculus. In order to explore the validity of modelling the projectile motion by the fractional approach, we compared our results with the experimental data of mortar.     

Motion of a circular viscoplastic plate subject to projectile impact

A method is presented for the solution of problems involving the impact on circular plates of projectiles moving with velocities high enough to produce plastic deformations. The plate material is assumed to obey the Mises-Huber yield condition and its associated flow rule for static deformations and behave as a viscoplastic solid for dynamic deformations. The nonlinear constitutive equations of the material are partially linearized in such a way as to allow the required result to be obtained by superposition of the static rigid perfectly-plastic solution and the dynamical elastic solution of equivalent problems. The method is illustrated by application to the case of a clamped plate struck at the center by a projectile of negligible radius.In addition a method for the rapid estimation of the motion of the projectile after impact is given and the result for the time required to bring the projectile to rest and the consequent central deflection of the plate are obtained.

---

Zusammenfassung Es wird eine Methode zur Lösung von Problemen beschrieben, die sich im Zusammenhang mit dem Aufschlag von Geschossen auf kreisförmige Platten ergeben; dabei soll die Geschossgeschwindigkeit so hoch sein, dass plastische Deformationen auftreten. Das Plattenmaterial soll der Mises-Huberschen Fliessbedingung und dem zugehörigen Fliessgesetz für statische Deformationen genügen und sich für dynamische Deformation wie ein viskoplastischer Festkörper verhalten. Die nicht linearen Stoffgleichungen des Materials werden teilweise linearisiert, so dass das gesuchte Resultat als Superposition der statischen starr-ideal-plastischen Lösung und der dynamisch-elastischen Lösung erhalten werden kann. Die Methode wird durch Anwendung auf den Fall einer eingespannten Platte veranschaulicht, auf deren Mittelpunkt ein Geschoss von vernachlässigbarem Radius auftrifft.Ausserdem wird eine Methode zur Abschätzung der Geschossbewegung nach dem Einschlag gegeben und die Zeit vom Aufprall bis zum Stillstand des Geschosses sowie die resultierende Auslenkung des Scheibenmittelpunktes berechnet.

     

An iteration method for the symmetric solutions and the optimal approximation solution of the matrix equation AXB=C

An iteration method is constructed to solve the linear matrix equation AXB=C over symmetric X. By this iteration method, the solvability of the equation AXB=C over symmetric X can be determined automatically, when the equation AXB=C is consistent over symmetric X, its solution can be obtained within finite iteration steps, and its least-norm symmetric solution can be obtained by choosing a special kind of initial iteration matrix, furthermore, its optimal approximation solution to a given matrix can be derived by finding the least-norm symmetric solution of a new matrix equation . Finally, numerical examples are given for finding the symmetric solution and the optimal approximation symmetric solution of the matrix equation AXB=C.     

Minimum norm solutions for cooperative games

A solution f for cooperative games is a minimum norm solution, if the space of games has a norm such that f(v) minimizes the distance (induced by the norm) between the game v and the set of additive games. We show that each linear solution having the inessential game property is a minimum norm solution. Conversely, if the space of games has a norm, then the minimum norm solution w.r.t. this norm is linear and has the inessential game property. Both claims remain valid also if solutions are required to be efficient. A minimum norm solution, the least square solution, is given an axiomatic characterization.      

On solutions of matrix equations and

With the help of the Kronecker map, a complete, general and explicit solution to the Yakubovich matrix equation V−AVF=BW, with F in an arbitrary form, is proposed. The solution is neatly expressed by the controllability matrix of the matrix pair (A,B), a symmetric operator matrix and an observability matrix. Some equivalent forms of this solution are also presented. Based on these results, explicit solutions to the so-called Kalman–Yakubovich equation and Stein equation are also established. In addition, based on the proposed solution of the Yakubovich matrix equation, a complete, general and explicit solution to the so-called Yakubovich-conjugate matrix is also established by means of real representation. Several equivalent forms are also provided. One of these solutions is neatly expressed by two controllability matrices, two observability matrices and a symmetric operator matrix.     

A conflict between sequential rationality and consistency principles

It is shown that no solution concept that selects sequentially rational (perfect, proper, persistent, or members of some stable set of) equilibria satisfies the following consistency property. Suppose that in every solution of the game G, player i's action is a, and denote by G ^a the game in which player i is restricted to choose a. Then some player j≠i has an action c that is used with positive probability in both some solution of G and some solution of G ^a. This result illustrates a conflict between a mild consistency condition and sequential rationality. Received: January 2001/Final version: April 2002     

On Interacting Systems of Hilbert-Space-Valued Diffusions

A nonlinear Hilbert-space-valued stochastic differential equation where L ^-1 (L being the generator of the evolution semigroup) is not nuclear is investigated in this paper. Under the assumption of nuclearity of L ^-1 , the existence of a unique solution lying in the Hilbert space H has been shown by Dawson in an early paper. When L ^-1 is not nuclear, a solution in most cases lies not in H but in a larger Hilbert, Banach, or nuclear space. Part of the motivation of this paper is to prove under suitable conditions that a unique strong solution can still be found to lie in the space H itself. Uniqueness of the weak solution is proved without moment assumptions on the initial random variable. A second problem considered is the asymptotic behavior of the sequence of empirical measures determined by the solutions of an interacting system of H -valued diffusions. It is shown that the sequence converges in probability to the unique solution Λ ₀ of the martingale problem posed by the corresponding McKean—Vlasov equation. Accepted 4 April 1996     

Construction of Solutions and L^1-error Estimates of Viscous Methods for Scalar Conservation Laws with Boundary

This paper is concerned with an initial boundary value problem for strictly convex conservation laws whose weak entropy solution is in the piecewise smooth solution class consisting of finitely many discontinuities. By the structure of the weak entropy solution of the corresponding initial value problem and the boundary entropy condition developed by Bardos-Leroux Nedelec, we give a construction method to the weak entropy solution of the initial boundary value problem. Compared with the initial value problem, the weak entropy solution of the initial boundary value problem includes the following new interaction type： an expansion wave collides with the boundary and the boundary reflects a new shock wave which is tangent to the boundary. According to the structure and some global estimates of the weak entropy solution, we derive the global L^1-error estimate for viscous methods to this initial boundary value problem by using the matching travelling wave solutions method. If the inviscid solution includes the interaction that an expansion wave collides with the boundary and the boundary reflects a new shock wave which is tangent to the boundary, or the inviscid solution includes some shock wave which is tangent to the boundary, then the error of the viscosity solution to the inviscid solution is bounded by O（ε^1/2） in L^1-norm; otherwise, as in the initial value problem, the L^1-error bound is O（ε｜ In ε｜）.     

On the unique solvability of fuzzy relational equations

The solution set of a consistent system of fuzzy relational equations with max-min composition can be characterized by one maximum solution and a finite number of minimal solutions. A polynomial-time method of O(mn) complexity is proposed to determine whether such a system has a unique minimal solution and/or a unique solution, where m, n are the dimensions of the input data. The proposed method can be extended to examining a system of fuzzy relational equations with max-T composition where T is a continuous triangular norm.     

Reproducing Kernel method for the solution of nonlinear hyperbolic telegraph equation with an integral condition

In this article, an iterative method is proposed for solving nonlinear hyperbolic telegraph equation with an integral condition. Its exact solution is represented in the form of series in the reproducing kernel space. In the mean time, the n‐term approximation u_n(x, t) of the exact solution u(x, t) is obtained and is proved to converge to the exact solution. Moreover, the partial derivatives of u_n(x, t) are also convergent to the partial derivatives of u(x, t). Some numerical examples have been studied to demonstrate the accuracy of the present method. Results obtained by the method have been compared with the exact solution of each example and are found to be in good agreement with each other. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 27: 867–886, 2011     

A two-space dimensional semilinear heat equation perturbed by (Gaussian) white noise

A two-space dimensional heat equation perturbed by a white noise in a bounded volume is considered. The equation is perturbed by a non-linearity of the type λ : f(AU) :, where :: means Wick (re)ordering with respect to the free solution;λ, A are small parameters, U denotes a solution, f is the Fourier transform of a complex measure with compact support. Existence and uniqueness of the solution in a class of Colombeau-Oberguggenberger generalized functions is proven. An explicit construction of the solution is given and it is shown that each term of the expansion in a power series in λ is associated with an L ²-valued measure when A is a small enough. Received: 20 July 1997 / Revised version: 1 February 2001 / Published online: 9 October 2001     

On the Logistic Equation in the Complex Plane

The famous logistic differential equation is studied in the complex plane. The method used is based on a functional analytic technique which provides a unique solution of the ordinary differential equation (ODE) under consideration in H ₂(𝔻) or H ₁(𝔻) and gives rise to an equivalent difference equation for which a unique solution is established in ?₂ or ?₁. For the derivation of the solution of the logistic differential equation this discrete equivalent equation is used. The obtained solution is analytic in {z ∈ ?: |z| <T}, T > 0. Numerical experiments were also performed using the classical 4th order Runge–Kutta method. The obtained results were compared for real solutions as well as for solutions of the form y(t) = u(t) + iv(t), t ∈ ?. For t ∈ ? the solution derived by the present method, seems to have singularities, that is, points where it ceases to be analytic, in certain sectors of the complex plane. These sectors, depending on the values of the involved parameters, can move at different directions, join forming common sectors, or pass through each other and continue moving independently. Moreover, the real and imaginary part of the solution seem to exhibit oscillatory behavior near these sectors.     

An efficient algorithm for the generalized centro-symmetric solution of matrix equation <Emphasis Type="Italic">A X B</Emphasis> = <Emphasis Type="Italic">C</Emphasis>

In this paper, an iterative algorithm is constructed for solving linear matrix equation AXB = C over generalized centro-symmetric matrix X. We show that, by this algorithm, a solution or the least-norm solution of the matrix equation AXB = C can be obtained within finite iteration steps in the absence of roundoff errors; we also obtain the optimal approximation solution to a given matrix X ₀ in the solution set of which. In addition, given numerical examples show that the iterative method is efficient.     

On the sensitivity matrix of the Nash bargaining solution

In this note we provide a characterization of a subclass of bargaining problems for which the Nash solution has the property of disagreement point monotonicity. While the original d-monotonicity axiom and its stronger notion, strong d-monotonicity, were introduced and discussed by Thomson (J Econ Theory, 42: 50–58, 1987), this paper introduces local strong d-monotonicity and derives a necessary and sufficient condition for the Nash solution to be locally strongly d-monotonic. This characterization is given by using the sensitivity matrix of the Nash bargaining solution w.r.t. the disagreement point d. Moverover, we present a sufficient condition for the Nash solution to be strong d-monotonic.     

Global existence to a reaction–diffusion equation through a self‐similar‐like upper solution

A demonstration method is presented, which will ensure the existence of positive global solutions in time to the reaction–diffusion equation ?u_t+Δu+u^p=0 in ?ⁿ×[0, ∞), for exponents p?3 and space dimensions n?3. This method does not require the initial value to have a specific uniform smallness condition, but rather to satisfy a bell‐like form. The method is based on a specific upper solution, which models the diffusion process of the heat equation. The upper solution is not self‐similar, but does have a self‐similar‐like form. After transforming the reaction–diffusion problem into an equivalent one, whose initial value is uniformly very small, a local solution is obtained in the time interval [0, 1] by the use of this upper solution. This local solution is then extended to [0, ∞) through an infinite sequence of extensions. At each step, an appropriate change of variables will transform the extension into a problem nearly identical to the local problem in [0, 1]. These transformations exploit the diffusive and self‐similar‐like nature of the upper solution. Copyright © 2009 John Wiley & Sons, Ltd.     

Numerical evaluation of gauge invariants for <Emphasis Type="Italic">a</Emphasis>-gauge solutions in open string field theory

We explore the vacuum structure in the bosonic open string field theory expanded near an identity-based solution parameterized by a (≥ −1/2). Analyzing the expanded theory using the level-truncation approximation up to the level 20, we find that the theory has the tachyon vacuum solution for a ≥ −1/2. We also find that at a = −1/2, there exists an unstable vacuum solution in the expanded theory and the solution is expected to be the perturbative open string vacuum. These results reasonably support the hypothesis that the identity-based solution is a trivial pure gauge configuration for a > −1/2, but it can be regarded as the tachyon vacuum solution at a = −1/2.     

Well-posedness and stability of solutions for set optimization problems

In this paper, some characterizations for the generalized l-B-well-posedness and the generalized u-B-well-posedness of set optimization problems are given. Moreover, the Hausdorff upper semi-continuity of l-minimal solution mapping and u-minimal solution mapping are established by assuming that the set optimization problem is l-H-well-posed and u-H-well-posed, respectively. Finally, the upper semi-continuity and the lower semi-continuity of solution mappings to parametric set optimization problems are investigated under some suitable conditions.