scholarly journals SELCIE: a tool for investigating the chameleon field of arbitrary sources

2021 ◽  
Vol 2021 (12) ◽  
pp. 043
Author(s):  
Chad Briddon ◽  
Clare Burrage ◽  
Adam Moss ◽  
Andrius Tamosiunas

Abstract The chameleon model is a modified gravity theory that introduces an additional scalar field that couples to matter through a conformal coupling. This `chameleon field' possesses a screening mechanism through a nonlinear self-interaction term which allows the field to affect cosmological observables in diffuse environments whilst still being consistent with current local experimental constraints. Due to the self-interaction term the equations of motion of the field are nonlinear and therefore difficult to solve analytically. The analytic solutions that do exist in the literature are either approximate solutions and or only apply to highly symmetric systems. In this work we introduce the software package SELCIE (https://github.com/C-Briddon/SELCIE.git). This package equips the user with tools to construct an arbitrary system of mass distributions and then to calculate the corresponding solution to the chameleon field equation. It accomplishes this by using the finite element method and either the Picard or Newton nonlinear solving methods. We compared the results produced by SELCIE with analytic results from the literature including discrete and continuous density distributions. We found strong (sub-percentage) agreement between the solutions calculated by SELCIE and the analytic solutions.

2020 ◽  
Vol 25 (2) ◽  
pp. 29
Author(s):  
Desmond Adair ◽  
Aigul Nagimova ◽  
Martin Jaeger

The vibration characteristics of a nonuniform, flexible and free-flying slender rocket experiencing constant thrust is investigated. The rocket is idealized as a classic nonuniform beam with a constant one-dimensional follower force and with free-free boundary conditions. The equations of motion are derived by applying the extended Hamilton’s principle for non-conservative systems. Natural frequencies and associated mode shapes of the rocket are determined using the relatively efficient and accurate Adomian modified decomposition method (AMDM) with the solutions obtained by solving a set of algebraic equations with only three unknown parameters. The method can easily be extended to obtain approximate solutions to vibration problems for any type of nonuniform beam.


1972 ◽  
Vol 94 (1) ◽  
pp. 1-7 ◽  
Author(s):  
O. B. Dale ◽  
R. Cohen

A method is presented for obtaining and optimizing the frequency response of one-dimensional damped linear continuous systems. The systems considered are assumed to contain unknown constant parameters in the boundary conditions and equations of motion which the designer can vary to obtain a minimum resonant response in some selected frequency interval. The unknown parameters need not be strictly dissipative nor unconstrained. No analytic solutions, either exact or approximate, are required for the system response and only initial value numerical integrations of the state and adjoint differential equations are required to obtain the optimal parameter set. The combinations of state variables comprising the response and the response locations are arbitrary.


Mathematics ◽  
2021 ◽  
Vol 9 (23) ◽  
pp. 3038
Author(s):  
Zi-Liang Li ◽  
Jin-Qing Liu

The horizontal equations of motion for an inviscid homogeneous fluid under the influence of pressure disturbance and waves are applied to investigate the nonlinear process of solitary waves and cyclone genesis forced by a moving pressure disturbance in atmosphere. Based on the reductive perturbation analysis, it is shown that the nonlinear evolution equation for the wave amplitude satisfies the Korteweg–de Vries equation with a forcing term (fKdV equation for short), which describes the physics of a shallow layer of fluid subject to external pressure forcing. Then, with the help of Hirota’s direct method, the analytic solutions of the fKdV equation are studied and some exact vortex solutions are given as examples, from which one can see that the solitary waves and vortex multi-pole structures can be excited by external pressure forcing in atmosphere, such as pressure perturbation and waves. It is worthy to point out that cyclone and waves can be excited by different type of moving atmospheric pressure forcing source.


2021 ◽  
pp. 106-115
Author(s):  
A.D. Matveev

It is important to know the error or the approximate solution when calculating the strength of elastic composite structures (bodies) using the finite element method (FEM). To construct a sequence of solutions according to the FEM is necessary for the evaluation of the approximation error. It involves the grinding procedure for discrete models. The implementation of the grinding procedure for basic models that take into account the inhomogeneous, micro-homogeneous structures of bodies within the microapproach requires ample computer resources. This paper proposes a method of equivalent strength conditions (MESC) to calculate the static strength of elastic bodies with a non-uniform, microuniform regular structure. The calculation of composite bodies strength according to the MESC is reduced to the calculation of isotropic homogeneous bodies strength using equivalent strength conditions. Adjusted equivalent strength conditions are used in the numerical implementation of the MESC. They take into account the error of the approximate solutions. If a set of loads is specified for a composite body, then generalized equivalent strength conditions are applied. The FEM-based calculation of composite bodies strength that follows the MESC using multigrid finite elements requires 103 ÷ 105 times less computer memory than a similar calculation using ground basic models of composite bodies. The provided example of strength calculation for a beam with an inhomogeneous regular fiber structure using the MESC shows its high efficiency. The main MESC implementation procedures are outlined.


Open Physics ◽  
2018 ◽  
Vol 16 (1) ◽  
pp. 143-148 ◽  
Author(s):  
Emran Khoshrouye Ghiasi ◽  
Reza Saleh

AbstractIn this paper, homotopy analysis method (HAM) and variational iteration method (VIM) are utilized to derive the approximate solutions of the Tricomi equation. Afterwards, the HAM is optimized to accelerate the convergence of the series solution by minimizing its square residual error at any order of the approximation. It is found that effect of the optimal values of auxiliary parameter on the convergence of the series solution is not negligible. Furthermore, the present results are found to agree well with those obtained through a closed-form equation available in the literature. To conclude, it is seen that the two are effective to achieve the solution of the partial differential equations.


Author(s):  
Shiro Kobayashi ◽  
Soo-Ik Oh ◽  
Taylan Altan

The concept of the finite-element procedure may be dated back to 1943 when Courant approximated the warping function linearly in each of an assemblage of triangular elements to the St. Venant torsion problem and proceeded to formulate the problem using the principle of minimum potential energy. Similar ideas were used later by several investigators to obtain the approximate solutions to certain boundary-value problems. It was Clough who first introduced the term “finite elements” in the study of plane elasticity problems. The equivalence of this method with the well-known Ritz method was established at a later date, which made it possible to extend the applications to a broad spectrum of problems for which a variational formulation is possible. Since then numerous studies have been reported on the theory and applications of the finite-element method. In this and next chapters the finite-element formulations necessary for the deformation analysis of metal-forming processes are presented. For hot forming processes, heat transfer analysis should also be carried out as well as deformation analysis. Discretization for temperature calculations and coupling of heat transfer and deformation are discussed in Chap. 12. More detailed descriptions of the method in general and the solution techniques can be found in References [3-5], in addition to the books on the finite-element method listed in Chap. 1. The path to the solution of a problem formulated in finite-element form is described in Chap. 1 (Section 1.2). Discretization of a problem consists of the following steps: (1) describing the element, (2) setting up the element equation, and (3) assembling the element equations. Numerical analysis techniques are then applied for obtaining the solution of the global equations. The basis of the element equations and the assembling into global equations is derived in Chap. 5. The solution satisfying eq. (5.20) is obtained from the admissible velocity fields that are constructed by introducing the shape function in such a way that a continuous velocity field over each element can be denned uniquely in terms of velocities of associated nodal points.


1970 ◽  
Vol 21 (1) ◽  
pp. 18-48 ◽  
Author(s):  
C. W. Stammers

SummaryThe nature of flapping torsion flutter of a helicopter blade in forward flight is discussed. The essential complication in the analysis is the presence of periodic coefficients in the equations of motion; approximate solutions are obtained by use of a perturbation procedure. An unusual behaviour of the flutter equations which occurs when the fundamental frequency of flutter is a half-integer multiple of rotational speed is studied. Two different instability mechanisms can be distinguished and are related to the two energy sources in the system, namely the rotation of the rotor and the forward speed of the helicopter. It is found that forward flight can have a significant stabilising influence on flutter and that, as the tip speed ratio increases, flutter occurs predominantly at half-integer frequencies. The results are confirmed by the use of a numerical method.


1982 ◽  
Vol 37 (8) ◽  
pp. 848-858 ◽  
Author(s):  
D. Correa-Restrepo

Resistive ballooning modes in general three-dimensional configurations are studied on the basis of the equations of motion of resistive MHD. Assuming small, constant resistivity and perturbations localized transversally to the magnetic field, a stability criterion is derived in the form of a coupled system of two second-order differential equations. This criterion contains several limiting cases, in particular the ideal ballooning mode criterion and criteria for the stability of symmetric systems. Assuming small growth rates, analytical results are derived by multiple-length-scale expansion techniques. Instabilities are found, their growth rates scaling as fractional powers of the resistivity


1975 ◽  
Vol 69 ◽  
pp. 237-244
Author(s):  
P. O. Vandervoort

This paper reviews a series of investigations of the orbits of stars in the regions of the Lindblad resonances of a spiral galaxy. The analysis is formulated in an epicyclic approximation. Analytic solutions of the epicyclic equations of motion are obtained by the method of harmonic balance of Bogoliubov and Mitropolsky. These solutions represent the resonance phenomena exhibited by the orbits in generally excellent agreement with numerical solutions.


Sign in / Sign up

Export Citation Format

Share Document