Modeling material responses by arbitrary Lagrangian Eulerian formulation and adaptive mesh refinement method

In this paper we report an efficient numerical method combining a staggered arbitrary Lagrangian Eulerian (ALE) formulation with the adaptive mesh refinement (AMR) method for materials modeling including elastic–plastic flows, material failure, and fragmentation predictions. Unlike traditional AMR applied on fixed domains, our investigation focuses on the application to moving and deforming meshes resulting from Lagrangian motion. We give details of this numerical method with a capability to simulate elastic–plastic flows and predict material failure and fragmentation, and our main focus of this paper is to create an efficient method which combines ALE and AMR methods to simulate the dynamics of material responses with deformation and failure mechanisms. The interlevel operators and boundary conditions for these problems in AMR meshes have been investigated, and error indicators to locate material deformation and failure regions are studied. The method has been applied on several test problems, and the solutions of the problems obtained with the ALE–AMR method are reported. Parallel performance and software design for the ALE–AMR method are also discussed.     

Arbitrary Lagrangian Eulerian formulation for two-dimensional flows using dynamic meshes with edge swapping

The dynamic modification of the computational grid due to element displacement, deformation and edge swapping is described here in terms of suitably-defined continuous (in time) alterations of the geometry of the elements of the dual mesh. This new interpretation allows one to describe all mesh modifications within the arbitrary Lagrangian Eulerian framework, thus removing the need to interpolate the solution across computational meshes with different connectivity. The resulting scheme is by construction conservative and it is applied here to the solution of the Euler equations for compressible flows in two spatial dimensions. Preliminary two dimensional numerical simulations are presented to demonstrate the soundness of the approach. Numerical experiments show that this method allows for large time steps without causing element invalidation or tangling and at the same time guarantees high quality of the mesh elements without resorting to global re-meshing techniques, resulting in a very efficient solver for the analysis of e.g. fluid-structure interaction problems, even for those cases that require large mesh deformations or changes in the domain topology.     

Discretization of hyperelasticity on unstructured mesh with a cell-centered Lagrangian scheme

We present a new cell-centered Lagrangian scheme on unstructured mesh for hyperelasticity. It is based on the recently proposed Glace scheme [11] for compressible gas dynamics. We show how to use the multiplicative decomposition of the gradient of deformation and the entropy property to derive the new scheme. We also prove the compatibility of this discretization with usual calculations of mass. Our motivation is to use hyperelasticity models for the study of finite plasticity, which is an extension of hypoelasticity to finite deformations. Hyperelasticity is a natural choice for extended models in solid mechanics, because of its mathematical structure which is a system of conservation laws with full rotational invariance. We study these properties for the Lagrangian system, and detail the various Eulerian formulations.     

On a Lagrangian formulation of gravitoelectromagnetism

We examine a Lagrangian formulation of gravity based on an approach analogous to electromagnetism, called Gravitoelectromagnetism (GEM). The gravitational analogue of the electromagnetic field tensor is a three-index tensor, \({\mathcal {F}_{\mu\nu\lambda}}\), defined in terms of a two-index gravitoelectromagnetic potential, \({\mathcal {A}_{\mu\nu}}\). The energy-momentum tensor is derived and is symmetric. We construct a Lagrangian which allows us to describe interactions between fermions, photons and gravitons. We calculate transition amplitudes of various processes involving gravitons: gravitational Møller scattering, gravitational Compton scattering, and the graviton photoproduction.     

Finite volume TVD formulation of lattice Boltzmann simulation on unstructured mesh

A numerical scheme is presented for accurate simulation of fluid flow using the lattice Boltzmann equation (LBE) on unstructured mesh. A finite volume approach is adopted to discretize the LBE on a cell-centered, arbitrary shaped, triangular tessellation. The formulation includes a formal, second order discretization using a Total Variation Diminishing (TVD) scheme for the terms representing advection of the distribution function in physical space, due to microscopic particle motion. The advantage of the LBE approach is exploited by implementing the scheme in a new computer code to run on a parallel computing system. Performance of the new formulation is systematically investigated by simulating four benchmark flows of increasing complexity, namely (1) flow in a plane channel, (2) unsteady Couette flow, (3) flow caused by a moving lid over a 2D square cavity and (4) flow over a circular cylinder. For each of these flows, the present scheme is validated with the results from Navier–Stokes computations as well as lattice Boltzmann simulations on regular mesh. It is shown that the scheme is robust and accurate for the different test problems studied.     

Lagrangian formulation of a general Hamiltonian theory with constraints

A Lagrangian formulation is constructed for a general Hamiltonian theory with constraints. A modification is proposed of the standard procedure of the Hamiltonianization of a Lagrangian theory in the case when the Lagrangian theory has primary constraints. The obtained results are used to establish the Lagrangian form of infinitesimally small canonical transformations in the Hamiltonian formalism.Translated from Izvestiya Vysshikh Uchebnykh Zavedenii, Fizika, No. 3, pp. 104–112, March, 1986.     

Innovative investigation on ion flows by an ad hoc Lagrangian formulation

A specialised Lagrangian formulation is here recommended to describe the main performances of unipolar ion flows driven by partial discharges. The subject fits into a number of electrostatic applications also governed by fluid-dynamic laws, which is why it needs careful examination by a solid and comprehensive approach. A reasoned selection of referential case studies are proposed in order to better appreciate the value of this predictive framework. The interdisciplinary character of the general problem under examination makes the role of the Lagrangian approach pivotal for the physical interpretation of theoretical results. These are given by resolving the Euler-Lagrange differential equation associated to individual examples.     

Radiation diffusion for multi-fluid Eulerian hydrodynamics with adaptive mesh refinement

Block-structured meshes provide the ability to concentrate grid points and computational effort in interesting regions of a flow field, without sacrificing the efficiency and low memory requirements of a regular grid. We describe an algorithm for simulating radiation diffusion on such a mesh, coupled to multi-fluid gasdynamics. Conservation laws are enforced by using locally conservative difference schemes along with explicit synchronization operations between different levels of refinement. In unsteady calculations each refinement level is advanced at its own optimal timestep. Particular attention is given to the appropriate coupling between the fluid energy and the radiation field, the behavior of the discretization at sharp interfaces, and the form of synchronization between levels required for energy conservation in the diffusion process. Two- and three-dimensional examples are presented, including parallel calculations performed on an IBM SP-2.     

Volume consistency in a staggered grid Lagrangian hydrodynamics scheme

Staggered grid Lagrangian schemes for compressible hydrodynamics involve a choice of how internal energy is advanced in time. The options depend on two ways of defining cell volumes: an indirect one, that guarantees total energy conservation, and a direct one that computes the volume from its definition as a function of the cell vertices. It is shown that the motion of the vertices can be defined so that the two volume definitions are identical. A so modified total energy conserving staggered scheme is applied to the Coggeshall adiabatic compression problem, and now also entropy is basically exactly conserved for each Lagrangian cell, and there is increased accuracy for internal energy. The overall improvement as the grid is refined is less than what might be expected.     

Development of a 3D Eulerian/Lagrangian Aircraft Icing Simulation Solver Based on OpenFOAM

A 3D icing simulation code is developed in the open-source CFD toolbox OpenFOAM. A hybrid Cartesian/body-fitted meshing method is used to generate high-quality meshes around complex ice shapes. Steady-state 3D Reynolds-averaged Navier-Stokes (RANS) equations are solved to provide the ensemble-averaged flow around the airfoil. Considering the multi-scale nature of droplet size distribution, and more importantly, to represent the less uniform nature of the Super-cooled Large Droplets (SLD), two droplet tracking methods are realized: the Eulerian method is used to track the small-size droplets (below 50 μm) for the sake of efficiency; the Lagrangian method with random sampling is used to track the large droplets (above 50 μm); the heat transfer of the surface overflow is solved on a virtual surface mesh; the ice accumulation is estimated via the Myers model; finally, the final ice shape is predicted by time marching. Limited by the availability of experimental data, validations are performed on 3D simulations of 2D geometries using the Eulerian and Lagrangian methods, respectively. The code proves to be feasible and accurate enough in predicting ice shapes. Finally, an icing simulation result of the M6 wing is presented to illustrate the full 3D capability.     

Lagrangian formulation of a geometrical scalar-tensor theory of gravitation

The author's geometrical theory of the scalar-tensor gravitational field is extended by formulating it in terms of a Lagrangian. An exact solution of the coupled nonlinear field equations for a static point mass is also presented. This theory which is conformally equivalent to the empty spaceEinstein equations predicts the same results for experiments as the usual theory of Brans and Dicke which has a non-zero energy momentum tensor.     

用电子秤测液体的表面张力系数

介绍了用电子秤测液体表面张力系数的原理和方法。     

耦合界面张力的三维流体界面不稳定性的格子Boltzmann模拟

采用介观格子Boltzmann方法模拟界面张力作用下三维流体界面的Rayleigh-Taylor (RT)不稳定性的增长过程,主要分析表面张力对流体界面动力学行为及尖钉和气泡后期增长的影响机制.首先发现三维RT不稳定性的发生存在临界表面张力(σ_c),其值随着流体Atwood数的增大而增大,且数值预测值与理论分析结果σ_c=(ρ_h-ρ₁)g/k~2一致.另外,随着表面张力的增大,不稳定性演化过程中界面卷吸程度和结构复杂性逐渐减弱,系统中界面破裂形成离散液滴的数目也显著减少.相界面的后期动力学行为也从非对称发展转向始终保持关于中轴线对称.尖钉与气泡振幅在表面张力较小时对其变化不显著,当表面张力增大到一定值后,可以有效地抑制尖钉与气泡振幅的增长.进一步发现,高雷诺数三维RT不稳定性在不同表面张力下均经历4个不同的发展阶段:线性阶段、饱和速度阶段、重加速和混沌混合阶段.尖钉与气泡在饱和速度阶段以近似恒定的速度增长,其渐进速度的值与修正的势流理论模型结果一致.受非线性Kelvin-Helmholtz旋涡的剪切作...     

Lagrangian formulation of a generalized Lane-Emden equation and double reduction

Abstract

We classify the Noether point symmetries of the generalized Lane-Emden equation y″+ ny′/x+ f(y)?=?0 with respect to the standard Lagrangian L = xⁿy′²/2 — xⁿ ∫f(y)dy for various functions f(y). We obtain first integrals of the various cases which admit Noether point symmetry and find reduction to quadratures for these cases. Three new cases are found for the function f(y). One of them is f(y) = αy^r , where r ≠ 0,1. The case r?=?5 was considered previously and only a one-parameter family of solutions was presented. Here we provide a complete integration not only for r?= 5 but for other r values. We also give the Lie point symmetries for each case. In two of the new cases, the single Noether symmetry is also the only Lie point symmetry.     

Classical Lorentz-invariant theory of systems with self-action: Lagrangian formulation

A Lorentz-invariant theory of singular Lagrangian particle systems with selfaction, treated as a dependence of the Lagrangian upon acceleration, is developed. The Lagrangian equations and their exact special solutions are found. The first integrals of 4-momentum and angular momentum are calculated. Particles possessing weak self-action are treated as classical analogs of particles with half-integral spin.     

Dispersion simulation of a wind tunnel experiment with Lagrangian particle models

Summary A particle model LAMDA suitable for dealing with the atmospheric dispersion is presented. The reliability of the model is tested comparing the results of its simulations to the wind tunnel measurements by Khurshudyanet al. (1981). Two versions of the model, both based on the Langevin equation and askewed distribution of the vertical wind velocity fluctuations, are considered. To develop the second version of the model we derived a proper scheme to produce skewed distributions of particle velocities consistent with the observed first and second moments of turbulent fluctuating velocities, with their horizontal and vertical derivatives, and their cross-correlations. The 2D geometry of the obstacle (in the wind tunnel experiment considered) allowed some simplifications. Due to the lack of some input data (such as the vertical profiles of crosswind standard deviation of wind velocity fluctuations and of the Lagrangian time scales) we looked for three different parametrizations. It was found that particular combinations of measured and parametrized data could give rise to critical vertical regions in which the derived scheme for the generation of random vertical velocity fluctuations cannot be applied. The best vertical distributions of the Lagrangian time scales (third parametrization) were estimated by fitting simple formulations for the average plume height and lateral variances to the measured data. The main results of this work were the model performance in simulating dispersion in shear flow over flat terrain, and its sensivity to the shape of the crosscorrelation term .     

Scattering on plasmonic nanostructures arrays modeled with a surface integral formulation

The surface integral formulation is a flexible, multiscale and accurate tool to simulate light scattering on nanostructures. Its generalization to periodic arrays is introduced in this paper. The general electromagnetic scattering problem is reduced to a discretizated model using the Method of Moments on the surface of the scatterers in the unit cell. The study of the resonances of an array of bowtie antennas illustrates the main features of the method. When placed into an array, the bowtie antennas show additional resonances compared to those of an individual antenna. Using the surface integral formulation, we are able to investigate both near-field and far-field properties of these resonances, with a high level of accuracy.