Least Squares Fitting of Chacón-Gielis Curves by the Particle Swarm Method of Optimization

Ricardo Chacón generalized Johan Gielis's superformula by introducing elliptic functions in place of trigonometric functions. In this paper an attempt has been made to fit the Chacón-Gielis curves (modified by various functions) to simulated data by the least squares principle. Estimation has been done by the Particle Swarm (PS) methods of global optimization. The Repulsive Particle Swarm optimization algorithm has been used. It has been found that although the curve-fitting exercise may be satisfactory, a lack of uniqueness of Chacón-Gielis parameters to data (from which they are estimated) poses an insurmountable difficulty to interpretation of findings.Least squares multimodal nonlinear curve-fitting; Ricardo Chacón; Jacobian Elliptic functions; Weierstrass ; Gielis super-formula; supershapes; Particle Swarm method; Repulsive Particle Swarm method of Global optimization; nonlinear programming; multiple sub-optima; global; local optima; fit; empirical; estimation; cellular automata; fractals

A database of genus 2 curves over the rational numbers

We describe the construction of a database of genus 2 curves of small discriminant that includes geometric and arithmetic invariants of each curve, its Jacobian, and the associated L-function. This data has been incorporated into the L-Functions and Modular Forms Database (LMFDB).Comment: 15 pages, 7 tables; bibliography formatting and typos fixe

Proving the triviality of rational points on Atkin-Lehner quotients of Shimura curves

In this paper we give a method for studying global rational points on certain quotients of Shimura curves by Atkin-Lehner involutions. We obtain explicit conditions on such quotients for rational points to be ``trivial'' (coming from CM points only) and exhibit an explicit infinite family of such quotients satisfying these conditions.Comment: 25 pages. To appear in Mathematische Annale

Stationary solutions of the one-dimensional nonlinear Schroedinger equation: I. Case of repulsive nonlinearity

All stationary solutions to the one-dimensional nonlinear Schroedinger equation under box and periodic boundary conditions are presented in analytic form. We consider the case of repulsive nonlinearity; in a companion paper we treat the attractive case. Our solutions take the form of stationary trains of dark or grey density-notch solitons. Real stationary states are in one-to-one correspondence with those of the linear Schr\"odinger equation. Complex stationary states are uniquely nonlinear, nodeless, and symmetry-breaking. Our solutions apply to many physical contexts, including the Bose-Einstein condensate and optical pulses in fibers.Comment: 11 pages, 7 figures -- revised versio

Horizon Pretracking

We introduce horizon pretracking as a method for analysing numerically generated spacetimes of merging black holes. Pretracking consists of following certain modified constant expansion surfaces during a simulation before a common apparent horizon has formed. The tracked surfaces exist at all times, and are defined so as to include the common apparent horizon if it exists. The method provides a way for finding this common apparent horizon in an efficient and reliable manner at the earliest possible time. We can distinguish inner and outer horizons by examining the distortion of the surface. Properties of the pretracking surface such as its expansion, location, shape, area, and angular momentum can also be used to predict when a common apparent horizon will appear, and its characteristics. The latter could also be used to feed back into the simulation by adapting e.g. boundary or gauge conditions even before the common apparent horizon has formed.Comment: 14 pages, 8 figures, minor change

Boundary-Conforming Finite Element Methods for Twin-Screw Extruders using Spline-Based Parameterization Techniques

This paper presents a novel spline-based meshing technique that allows for usage of boundary-conforming meshes for unsteady flow and temperature simulations in co-rotating twin-screw extruders. Spline-based descriptions of arbitrary screw geometries are generated using Elliptic Grid Generation. They are evaluated in a number of discrete points to yield a coarse classical mesh. The use of a special control mapping allows to fine-tune properties of the coarse mesh like orthogonality at the boundaries. The coarse mesh is used as a 'scaffolding' to generate a boundary-conforming mesh out of a fine background mesh at run-time. Storing only a coarse mesh makes the method cheap in terms of memory storage. Additionally, the adaptation at run-time is extremely cheap compared to computing the flow solution. Furthermore, this method circumvents the need for expensive re-meshing and projections of solutions making it efficient and accurate. It is incorporated into a space-time finite element framework. We present time-dependent test cases of non-Newtonian fluids in 2D and 3D for complex screw designs. They demonstrate the potential of the method also for arbitrarily complex industrial applications

Edge-promoting reconstruction of absorption and diffusivity in optical tomography

In optical tomography a physical body is illuminated with near-infrared light and the resulting outward photon flux is measured at the object boundary. The goal is to reconstruct internal optical properties of the body, such as absorption and diffusivity. In this work, it is assumed that the imaged object is composed of an approximately homogeneous background with clearly distinguishable embedded inhomogeneities. An algorithm for finding the maximum a posteriori estimate for the absorption and diffusion coefficients is introduced assuming an edge-preferring prior and an additive Gaussian measurement noise model. The method is based on iteratively combining a lagged diffusivity step and a linearization of the measurement model of diffuse optical tomography with priorconditioned LSQR. The performance of the reconstruction technique is tested via three-dimensional numerical experiments with simulated measurement data.Comment: 18 pages, 6 figure