Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userNumerical simulation of knotted solutions for Maxwell equations
57
0
0.0
(
0
)
Added by
Juan Jos\\'e Omiste
Publication date
2020
fields
Physics
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
In this work, we use the finite differences in time domain (FDTD) numerical method to compute and assess the validity of Hopf solutions, or hopfions, for the electromagnetic field equations. In these solutions, field lines form closed loops characterized by different knot topologies which are preserved during their time evolution. Hopfions have been studied extensively in the past from an analytical perspective but never, to the best of our knowledge, from a numerical approach. The implementation and validation of this technique eases the study of more complex cases of this phenomena; e.g. how these fields could interact with materials (e.g. anisotropic or non-linear), their coupling with other physical systems (e.g. plasmas), and also opens the path on their artificial generation by different means (e.g. antenna arrays or lasers).
rate research
Read More
We pursue here the development of models for complex (viscoelastic) fluids in shallow free-surface gravity flows which was initiated by [Bouchut-Boyaval, M3AS (23) 2013] for 1D (translation invariant) cases. The models we propose are hyperbolic quasilinear systems that generalize Saint-Venant shallow-water equations to incompressible Maxwell fluids. The models are compatible with a formulation of the thermo-dynamics second principle. In comparison with Saint-Venant standard shallow-water model, the momentum balance includes extra-stresses associated with an elastic potential energy in addition to a hydrostatic pressure. The extra-stresses are determined by an additional tensor variable solution to a differential equation with various possible time rates. For the numerical evaluation of solutions to Cauchy problems, we also propose explicit schemes discretizing our generalized Saint-Venant systems with Finite-Volume approximations that are entropy-consistent (under a CFL constraint) in addition to satisfy exact (discrete) mass and momentum conservation laws. In comparison with most standard viscoelastic numerical models, our discrete models can be used for any retardation-time values (i.e. in the vanishing solvent-viscosity limit). We finally illustrate our hyperbolic viscoelastic flow models numerically using computer simulations in benchmark test cases. On extending to Maxwell fluids some free-shear flow testcases that are standard benchmarks for Newtonian fluids, we first show that our (numerical) models reproduce well the viscoelastic physics, phenomenologically at least, with zero retardation-time. Moreover, with a view to quantitative evaluations, numerical results in the lid-driven cavity testcase show that, in fact, our models can be compared with standard viscoelastic flow models in sheared-flow benchmarks on adequately choosing the physical parameters of our models. Analyzing our models asymptotics should therefore shed new light on the famous High-Weissenberg Number Problem (HWNP), which is a limit for all the existing viscoelastic numerical models.
It has been shown in the authors companion paper that solutions of Maxwell-Klein-Gordon equations in $mathbb{R}^{3+1}$ possess some form of global strong decay properties with data bounded in some weighted energy space. In this paper, we prove pointwise decay estimates for the solutions for the case when the initial data are merely small on the scalar field but can be arbitrarily large on the Maxwell field. This extends the previous result of Lindblad-Sterbenz cite{LindbladMKG}, in which smallness was assumed both for the scalar field and the Maxwell field.
A simple model of random Brownian walk of a spherical mesoscopic particle in viscous liquids is proposed. The model can be both solved analytically and simulated numerically. The analytic solution gives the known Eistein-Smoluchowski diffusion law $<r^2> = Dt$ where the diffusion constant $D$ is expressed by the mass and geometry of a particle, the viscosity of a liquid and the average effective time between consecutive collisions of the tracked particle with liquid molecules. The latter allows to make a simulation of the Perrin experiment and verify in detailed study the influence of the statistics on the expected theoretical results. To avoid the problem of small statistics causing departures from the diffusion law we introduce in the second part of the paper the idea of so called Artificially Increased Statistics (AIS) and prove that within this method of experimental data analysis one can confirm the diffusion law and get a good prediction for the diffusion constant even if trajectories of just few particles immersed in a liquid are considered.
We present a computer simulation of entangled polymer solutions at equilibrium. The chains repel each other via a soft Gaussian potential, appropriate for semi-dilute solutions at the scale of a correlation blob. The key innovation to suppress chain crossings is to use a pseudo-continuous model of a backbone which effectively leaves no gaps between consecutive points on the chain, unlike the usual bead-and-spring model. Our algorithm is sufficiently fast to observe the entangled regime using a standard desktop computer. The simulated structural and mechanical correlations are in fair agreement with the expected predictions for a semi-dilute solution of entangled chains.
In this article we develop a numerical scheme to deal with interfaces between touching numerical grids when solving the second-order wave equation. We show that it is possible to implement an interface scheme of penalty type for the second-order wave equation, similar to the ones used for first-order hyperbolic and parabolic equations, and the second-order scheme used by Mattsson et al (2008). These schemes, known as SAT schemes for finite difference approximations and penalties for spectral ones, and ours share similar properties but in our case one needs to pass at the interface a smaller amount of data than previously known schemes. This is important for multi-block parallelizations in several dimensions, for it implies that one obtains the same solution quality while sharing among different computational grids only a fraction of the data one would need for a comparable (in accuracy) SAT or Mattsson et al.s scheme. The semi-discrete approximation used here preserves the norm and uses standard finite-difference operators satisfying summation by parts. For the time integrator we use a semi-implicit IMEX Runge-Kutta method. This is crucial, since the explicit Runge-Kutta method would be impractical given the severe restrictions that arise from the stiff parts of the equations.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
