We describe here the coherent formulation of electromagnetism in the nonrelativistic quantummechanical many-body theory. We use the mathematical frame of the field theory and its quantization in the spirit of the QED. This is necessary because of the manifold of misinterpretations emerging from the hystorical development of quantum mechanics, starting from the Schrödinger equation of a single particle in the presence of given electromagnetic fields, followed by the many-body theories of many charged identical particles having just Coulomb interactions inspired from the classical electromagnetic theory of point-like charges. However, this later is known to be inconsistent due to the self-interaction. This way could not be continued further to include properly the magnetic forces between the charged particles and lead to a lot of confusion about the interpretation of the magnetic field in the Hamiltonian, as well as about the gauge invariance. We emphasize the importance of the distinction between the applied (external fields) and the field in the matter. All these problems are length properly solved within the non-relativistic QED, nevertheless the confusion dominates in all the problems related to the magnetic properties of the solid state.

The relations between the kinetic energy spectrum and the second order longitudinal structure function in two dimensions are derived, and several examples are considered. The forward conversion (from spectrum to structure function) is illustrated first with idealized power law spectra, representing turbulent inertial ranges. The forward conversion is also applied to the zonal kinetic energy spectrum of Nastrom and Gage (1985) and the result agrees well with the longitudinal structure function of Lindborg (1999). The inverse conversion (from structure function to spectrum) is tested with data from 2D turbulence simulations. When applied to the theoretical structure function (derived from the forward conversion of the spectrum), the result closely resembles the original spectrum, except at the largest wavenumbers. However the inverse conversion is much less successful when applied to the structure function obtained from pairs of particles in the flow. This is because the inverse conversion favors large pair separations, which are typically noisy with particle data. Fitting the structure function to a polynomial improves the result, but not sufficiently to distinguish the correct inertial range dependencies. Furthermore the inversion of non-local spectra is largely unsuccessful. Thus it appears that focusing on structure functions with Lagrangian data is preferable to estimating spectra.

This paper gives both a general canonical quantum gravity theory and the general canonical quantum gravity theories of the Universe and general black hole, and discovers the relations reflecting symmetric properties of the standard nonlinear gravitational Lagrangian, which are not relevant to any concrete metric models. This paper concretely shows the general commutation relations of the general gravitational field operators and their zeroth, first, second and third style, respectively, of high order canonical momentum operators for the general nonlinear system of the standard gravitational Lagrangian, and then has finished all the four styles of the canonical quantization of the standard gravity.

The statistical bundle is the set of couples $(Q,W)$ of a probability density Q and a random variable W such that $\mathbb{E}$_Q [W] = 0. On a finite state space, we assume Q to be a probability density with respect to the uniform probability and give an affine atlas of charts such that the resulting manifold is a model for Information Geometry. Velocity and accelleration of a one-dimensional statistical model are computed in this set up. The Euler-Lagrange equations are derived from the Lagrange action integral. An example of Lagrangian using minus the entropy as potential energy is briefly discussed.

Although Eulerian approaches are standard in computational acoustics, they are less effective for certain classes of problems like bubble acoustics and combustion noise. A different approach for solving acoustic problems is to compute with individual particles following particle motion. In this paper, a Lagrangian approach to model sound propagation in moving fluid is presented and implemented numerically, using three meshfree methods to solve the Lagrangian acoustic perturbation equations (LAPE) in the time domain. The LAPE split the fluid dynamic equations into a set of hydrodynamic equations for the motion of fluid particles and perturbation equations for the acoustic quantities corresponding to each fluid particle. Then, three meshfree methods, the smoothed particle hydrodynamics (SPH) method, the corrective smoothed particle (CSP) method, and the generalized finite difference (GFD) method, are introduced to solve the LAPE and the linearized LAPE (LLAPE). The SPH and CSP methods are widely used meshfree methods, while the GFD method based on the Taylor series expansion can be easily extended to higher orders. Applications to modeling sound propagation in steady or unsteady fluids in motion are outlined, treating a number of different cases in one and two space dimensions. A comparison of the LAPE and the LLAPE using the three meshfree methods is also presented. The Lagrangian approach shows good agreement with exact solutions. The comparison indicates that the CSP and GFD method exhibit convergence in cases with different background flow. The GFD method is more accurate, while the CSP method can handle higher Courant numbers.

One open question remains regarding the theory of the generalized variational principle, that is, why the stress-strain relation still be derived from the generalized variational principle while the Lagrangian multiplier method is applied in vain? This study shows that the generalized variational principle can only be understood and implemented correctly within the framework of thermodynamics. As long as the functional has one of the combination $A(\epsilon_{ij})-\sigma_{ij}\epsilon_{ij}$ or $B(\sigma_{ij})-\sigma_{ij}\epsilon_{ij}$, its corresponding variational principle will produce the stress-strain relation without the need to introduce extra constraints by the Lagrangian multiplier method. It is proved herein that the Hu-Washizu functional $\Pi_{HW}[u_i,\epsilon_{ij},\sigma_{ij}]$ and Hu-Washizu variational principle comprise a real three-field functional.

By developing the previously proposed method of combining continuum mechanics with Einstein Field Equations, it has been shown that the classic relativistic description, curvilinear description, and quantum description of the physical system may be reconciled using the proposed Alena Tensor. For a system with an electromagnetic field, the Lagrangian density equal to the invariant of the electromagnetic field was obtained, vanishing four-divergence of canonical four-momentum appears to be consequence of the Poyinting theorem, and explicit form of one of gauges of the electromagnetic four-potential was introduced. The proposed method allows for further development with additional fields.

The Lagrangian formalism is developed for the population dynamics of interacting species that are described by several well-known models. The formalism is based on standard Lagrangians, which represent differences between the physical kinetic and potential energy-like terms. A method to derive these Lagrangians is presented and applied to selected theoretical models of the population dynamics. The role of the derived Lagrangians and the energy-like terms in the population dynamics is investigated, and it is shown that the obtained standard Lagrangians can be used to study oscillatory behavior of the models and the period of their oscillations, which may have ecological and environmental implications. These and other physical and biological insights gained from the derived standard Lagrangians are discussed.

In this work, we study the dynamical behaviors of the electromagnetic fields and material responses in the hyperbolic metamaterial consisting of periodically arranged metallic and dielectric layers. The thickness of each unit cell is assumed to be much smaller than the wavelength of the electromagnetic waves, so the effective medium concept can be applied. When electromagnetic (EM) fields are present, the responses of the medium in the directions parallel to and perpendicular to the layers are like that of Drude and Lorentz media, respectively. We derive the energy density of the EM fields and the power loss in the effective medium based on Poynting theorem and the dynamical equations of the polarization field. We also show that the Lagrangian density of the system can be constructed. The Euler-Lagrangian equations yield the correct dynamical equations of the electromagnetic fields and the polarization field in the medium. The canonical momentum conjugates to every dynamical field can be derived from the Lagrangian density via differentiation or variation with respect to that field. We apply Legendre transformation to this system, and find that the resultant Hamitonian density is identical to the energy density, up to an irrelevant divergence term.

The author demonstrates a stable Lagrangian solid modeling method, tracking the interactions of solid mass particles rather than using a meshed grid. This numerical method avoids the problem of tensile instability often seen with smooth particle applied mechanics by having the solid particles apply stresses expected with Hooke's law, as opposed to using a smoothing function for neighboring solid particles. This method has been tested successfully with a bar in tension, compression, and shear, as well as a disk compressed into a flat plate, and the numerical model consistently matched the analytical Hooke's law as well as Hertz contact theory for all examples. The solid modeling numerical method was then built into a 2-D model of a pressure vessel, which was tested with liquid water particles under pressure and simulated with smoothed particle hydrodynamics. This simulation was stable, and demonstrated the feasibility of Lagrangian specification modeling for fluid–solid interactions.

We have found a region in the northwestern Japan Sea near the Subpolar Front where mesoscale eddies regularly form and circulate. The anticyclonic eddy with a strong dynamic feature, vertically extended to the bottom (~3500 m), multilayer core structure and extreme values of temperature and salinity was sampled there in the cruise in May 2004. The altimetry-based eddy tracking showed that it had a life span of nine months with the size reaching 120 km. CTD observational data were used to explore the eddy’s features and vertical structure. The eddy had a multilayer core with different thermohaline characteristics which could be originated by eddy interaction with the surrounding waters during its life. The evolution of the eddy has been studied and analyzed with the help of altimetry-based Lagrangian indicators of water motion. Performing the particle-tracking numerical experiments, we computed the Lagrangian maps that contain information on the origin of water inside the eddy core and its ‘age’ on every day of the eddy’s life. Inspecting daily-computed Lagrangian maps, we documented the essential events in the evolution of the eddy including its formation, splitting, merger, entrainment and detrainment of water, erosion and eventual decay. All these observations have been verified with the infrared satellite images. The new Lagrangian technique has been applied to calculate day by day the content of water of various origin (coastal, subtropical and subarctic) entrained into the eddy core. It has been found that the surface core has been filled mainly with subtropical water originated in the southern flank of the Subpolar Front. The episodes with cardinal changes of ratio of different water types have been recorded and verified with the infrared satellite images and maps of the origin of water inside the eddy core.

A rapid-response Lagrangian model for simulating the transport of a chemical pollutant in the Arabian/Persian Guls is described. The model is well suited to provide a fast response after an emergency due to an accident or a deliberate spill. Baroclinic circulation was obtained from HYCOM ocean model and tides were calculated using a barotropic model. The interactions of pollutants with sediments (uptake/release processes) were described using a dynamic approach based on kinetic transfer coefficients and a stochastic numerical method. Some examples of model applications are shown.

The purpose of this paper is threefold: (i) to introduce and study the Euler-Lagrange prolongations of flatness PDEs solutions (best approximation of flatness) via associated least squares Lagrangian densities and integral functionals on Riemannian manifolds; (ii) to analyze some decomposable multivariate dynamics represented by Euler-Lagrange PDEs of least squares Lagrangians generated by flatness PDEs and Riemannian metrics; (iii) to give examples of explicit flat extremals and non-flat approximations.

The availability of high quality surface observations of precipitation and volume observations by polarimetric operational radars make it possible to constrain, evaluate and validate numerical models with a wide variety of microphysical schemes. In this article, a novel particle-based Monte-Carlo microphysical model (called “McSnow”) is used to simulate the outer rain bands of Hurricane Dorian which traversed the densely instrumented precipitation research facility operated by NASA at Wallops Island, Virginia. The rain bands showed steady stratiform vertical profiles with radar signature of dendritic growth layers near −15 °C and peak reflectivity in the bright band of 55 dBZ along with polarimetric signatures of wet snow with sizes inferred to exceed 15 mm. A 2D-video disdrometer measured frequent occurrences of large drops >5 mm and combined with an optical array probe the drop size distribution was well-documented in spite of uncertainty for drops <0.5 mm due to high wind gusts and turbulence. The 1D McSnow control run and four numerical “experiments” were conducted and compared with observations. One of the main findings is that even at the moderate rain rate of 10 mm/h collisional breakup is essential for the shape of the drop size distribution.

Abstract Using the Pauli algebra version of the Dirac equation six symmetric but non-equivalent copies of the QED Lagrangian are constructed. The group of symmetries is the group of outer automorphisms of the finite quaternion group, isomorphic to the group of permutations on three letters. Solutions to corresponding Dirac equations form six sectors of fermion fields. The intersections between different sectors can only contain massless fields, implying that the sectors are distinct. If the symmetry could be extended to the full Standard Model Lagrangian it would predict the ratio of dark matter to ordinary matter of 5 : 1 which is close to the ratio of 5.2 : 1 according to current observations.

Non-standard Lagrangians do not display any discernible energy-like terms, yet they give the same equations of motion as standard Lagrangians, which have easily identifiable energy-like terms. A new method to derive non-standard Lagrangians for second-order nonlinear differential equations with damping is developed and the limitations of this method are explored. It is shown that the limitations do not exist only for those nonlinear dynamical systems that can be converted into linear ones. The obtained results are applied to selected population dynamics models for which non-standard Lagrangians and their corresponding null Lagrangians and gauge functions are derived, and their roles in the population dynamics are discussed.

This article is devoted to Feller's diffusion equation which arises naturally in probabilities and physics (e.g. wave turbulence theory). If discretized naively, this equation may represent serious numerical difficulties since the diffusion coefficient is practically unbounded and most of its solutions are weakly divergent at the origin. In order to overcome these difficulties we reformulate this equation using some ideas from the Lagrangian fluid mechanics. This allows us to obtain a numerical scheme with a rather generous stability condition. Finally, the algorithm admits an elegant implementation and the corresponding Matlab code is provided with this article under an open source license.

Loop quantum gravity (LQG) quantizes gravity through the quantization of the space-time. It begins with general relativity (GR) and brings some concepts from quantum field theory to quantize the space-time. The conceptual background of GR is necessary to comprehend the LQG framework. Therefore, in this paper, prerequisite concepts such as special relativity, general relativity, covariant electromagnetism, tensor analysis, Lagrangian formulation, Hamiltonian formulation and basics of quantum mechanics are briefly introduced; that, are needed to study loop quantum gravity(LQG).

Self-organized criticality (SOC) is a universal mechanism for self-sustained critical behavior in large-scale systems evolving outside equilibrium. Our report explores a tentative link between SOC and Lagrangian field theory, with the long-term goal of bridging the gap between complex dynamics and the non-perturbative behavior of quantum fields.

The Lagrangian and Eulerian structure and dynamics of a strong wind event in the Tokar Gap region are described using a WRF model hindcast for 2008. Winds in the Tokar Gap reach 25 m s-1 and remain coherent as a jet far out over the Red Sea, whereas equally strong wind jets occurring in neighboring gaps are attenuated abruptly by a jump-like hydraulic transition that occur just offshore of the Sudan coast. The transition is made possible by the supercritical nature of the jets, which are fed by air that spills down from passes at relatively high elevation. By contrast, the spilling flow in the ravine-like Tokar Gap does not become substantially supercritical and therefore does not undergo a jump, and also carries more total horizontal momentum. The Tokar Wind Jet carries some air parcels across the Red Sea and into Saudi Arabia, whereas air parcel trajectories in the neighboring jets ascend as they cross through the jumps, then veer sharply to the southeast and do not cross the Red Sea. The mountain parameter Nh/U is estimated to lie in the rage 1.0-4.0 for the general region, a result roughly consistent with a primary gap jet having a long extension, and supercritical jets spilling down from higher elevation passes. The strong event is marked by the formation of a cyclonic cell near the upstream entrance to the Tokar Gap, a feature absent from the more moderate events that occur throughout the summer. The cell contains descending air parcels that are fed into the primary and secondary jets. An analysis of the Bernoulli function along air parcel trajectories reveals an approximate balance between the loss of potential energy and gain of internal energy and pressure, with surprisingly little contribution from kinetic energy, along the path of the descending flow. All jets attain the critical wind speed nominally required to loft dust into the atmosphere, though only the Tokar Gap has a broad, delta region with plentiful deposits of silt.

A cosmological model described by the perturbation metric is proposed, in which space is entrained by matter expanding as a result of the Big Bang. In the absence of an input of external energy and momentum, the solution of the Einstein equations assumes a pressure equal to zero and an accelerated expansion. The cosmological perturbation model explains the difference in measurements of the Hubble constant using gravitational lensing and a distance ladder.

In Lagrangian dynamics, the detection of coherent clusters can help understand the organization of transport by identifying regions with coherent trajectory patterns. Many clustering algorithms, however, rely on user-input parameters, requiring a priori knowledge about the flow and making the outcome subjective. Building on the conventional spectral clustering method of Hadjighasem et al (2016), a new parameter-free spectral clustering approach is developed that automatically identifies parameters and does not require any user-input choices. A noise-based metric for quantifying the coherence of the resulting coherent clusters is also introduced. The parameter-free spectral clustering is applied to two benchmark analytical flows, the Bickley Jet and the asymmetric Duffing oscillator, and to a realistic, numerically-generated oceanic coastal flow. In the latter case, the identified model-based clusters are tested using observed trajectories of real drifters. In all examples, our approach succeeded in performing the partition of the domain into coherent clusters with minimal inter-cluster similarity and maximum intra-cluster similarity. For the coastal flow, the resulting coherent clusters are qualitatively similar over the same phase of the tide on different days and even different years, whereas coherent clusters for the opposite tidal phase are qualitatively different.

A time integration scheme based on semi-Lagrangian advection and Laplace transform adjustment has been implemented in a baroclinic primitive equation model. The semi-Lagrangian scheme makes it possible to use large time steps. However, errors arising from the semi-implicit scheme increase with the time step size. In contrast, the errors using the Laplace transform adjustment remain relatively small for typical time steps used with semi-Lagrangian advection. Numerical experiments confirm the superior performance of the Laplace transform scheme relative to the semi-implicit reference model. The algorithmic complexity of the scheme is comparable to the reference model, making it computationally competitive, and indicating its potential for integrating weather and climate prediction models.

The Finite size Lyapunov exponent (FSLE) has been used extensively since the late 1990’s to diagnose turbulent regimes from Lagrangian experiments and to detect Lagrangian coherent structures in geophysical flows and two-dimensional turbulence. Historically, the FSLE was defined in terms of its computational method rather than via a mathematical formulation, and the behavior of the FSLE in the turbulent inertial ranges is based primarily on scaling arguments. Here we propose an exact definition of the FSLE based on conditional averaging of the finite amplitude growth rate (FAGR) of the particle pair separation. With this new definition, we show that the FSLE is a close proxy for the inverse structural time, a concept introduced a decade before the FSLE. The (in)dependence of the FSLE on initial conditions is also discussed, as well as the links between the FAGR and other relevant Lagrangian metrics, such as the finite time Lyapunov exponent and the second order velocity structure function.

Purpose: To improve numerical simulation of the non-contact tonometry test by using Arbitrary Eulerian-Lagrangian deforming mesh in the coupling between computational fluid dynamics model of an air jet and finite element model of the human eye. Methods: Computational fluid dynamics model simulated impingement of the air puff and consisted of 25920 wedge6 elements and employed Spallart-Allmaras model to simulate capture turbulence of the air jet. The time span of the jet wais 30 ms and maximum Reynolds number

Besides the Lagrangian and the Eulerian descriptions, the motion of a body can also be expressed relative to the present configuration of the body, known as the relative motion description. It is interesting to consider such a relative motion description in general to formulate the basic system of field equations for solid bodies. In doing so, when the time increment from the present state is small enough, the nonlinear constitutive equations can be linearized relative to the present state so that the resulting system becomes linear. This will be done for thermoelastic materials with a brief comment on the exploitation of entropy principle in general. Relative Lagrangian formulation is based on the well-known ``small-on-large'' idea, and can be implemented for solving problems with large deformation in successive incremental manner. Some applications of such a formulation in numerical simulations are briefly reviewed.

Measurements of temperature and salinity in the Amazon River plume in the open ocean of the western tropical North Atlantic are considered. The measurements were carried out using the AML CTD probe in the upper layer and a flow-through system that measures salinity and turbidity of seawater while the vessel is on the way. Additionally, archive oceanographic data from the WOD18 database, data of satellite altimetry measurements (AVISO), and satellite salinity data from Aquarius and SMOS were used. It is shown that the width of the Amazon plume is about 300–400 km, the depth of desalination is from 50 to 100 m. Surface salinity decreases compared to the background (36 PSU) by 0.5 PSU in February and more than by 3 PSU in September during the period of maximum development of the plume, which is determined from satellite measurements of surface salinity. Lagrangian modeling of the back-in-time advection of passive markers simulating freshwater particles was carried out. It was shown that the source of fresh water in the measurement area is the discharge of the Amazon.

In this paper, we prove that the sub convexlikeness introduced by V. Jeyakumar [1], and the subconvexlikeness defined in V. Jeyakumar [2] are equivalent in loccallly convex topological spaces. And then, we work with set-valued vector optimization problems and obtain some vector saddle-point theorems and vector Lagrangian theorems.

In remote sensing applications, one of the key points is the acquisition, real-time pre-processing and storage of information. Due to the large amount of information present in the form of images or videos, compression of this data is necessary. Compressed sensing (CS) is an efficient technique to meet this challenge. It consists in acquiring a signal, assuming that it can have a sparse representation, using a minimal number of non-adaptive linear measurements. After this CS process, a reconstruction of the original signal must be performed at the receiver. Reconstruction techniques are often unable to preserve the texture of the image and tend to smooth out its details. To overcome this problem, we propose in this work, a CS reconstruction method that combines the total variation regularization and the non-local self-similarity constraint. The optimization of this method is performed by the augmented Lagrangian which avoids the difficult problem of non-linearity and non-differentiability of the regularization terms. The proposed algorithm, called denoising compressed sensing by regularizations terms (DCSR), will not only perform image reconstruction but also denoising. To evaluate the performance of the proposed algorithm, we compare its performance with state-of-the-art methods, such as Nesterov's algorithm, group-based sparse representation and wavelet-based methods, in terms of denoising, and preservation of edges, texture and image details, as well as from the point of view of computational complexity. Our approach allows to gain up to 25% in terms of denoising efficiency, and visual quality using two metrics: PSNR and SSIM.

The design of thin-walled structures is commonly based on the solutions of linear boundary-valued problems, formulated within well-developed theories for elastic plates and shells. However, in modern appliances, especially in MEMS design, it is necessary to take into account non-linear mechanical effects that become decisive for flexible elements. Among substantial non-linear effects, that significantly change the deformation properties of thin plates, are the effects of residual stresses caused by incompatibility of deformations, which inevitably arise during the manufacture of ultra-thin elements. The development of new methods of mathematical modeling of residual stresses and incompatible finite deformations in plates is the subject of the paper. To this end the local unloading hypothesis is used. It makes possible to define smooth fields of local deformations (inverse implant field) for the mathematical formalization of incompatibility. The main outcomes are field equations, natural boundary conditions and conservation laws, derived from least action principle and variational symmetries with account of implant field. The derivations were carried out in the framework of elasticity theory for simple materials and, over and above, within the Cosserat’s theory of a two-dimensional continuum. As illustrative examples the distributions of incompatible deformations in a circular plate are considered.

The computational model that is able to estimate the temperature distribution inside a solid specimen during the film boiling phase of immersion quenching (water entry) process has been presented in paper, and it is based on the prescribed initial temperatures of the solid specimen and the liquid quenchant. In addition, the turbulence effects have to be considered using the assumed turbulence kinetic energy value, i.e., the “frozen turbulence” approach, that remains constant thorough the simulation. The studied material is nickel alloy, Inconel 600, for which the extensive experimental data is available. The work has been carried out using ANSYS Fluent computational fluid dynamics software and the methods for solution of Stefan problem by Eulerian two fluid VOF model.. A satisfactory agreement between the experimental and the calculated data has been achieved, yielding thereby the computationally obtained data that fits to a great extent to the prescribed error band of ± 10 % during the estimated film boiling phase of the immersion quenching process itself.

Meshfree particle method, which is always regarded as a pure Lagrangian approach, is easily represented complicated domain topologies, moving boundaries, and multiphase media. Solving acoustic problems with the mesfree particle method forms a branch of the acoustic wave modeling field, namely, particle-based computational acoustics (PCA). The aim of this paper is to improve the accuracy of using the PCA method to solve two-dimensional acoustic problems, and realize the particle representation with a hybrid meshfree and finite-difference time-domain (FDTD) method for acoustic boundary conditions at both the plane and curved surface. As a widely used Lagrangian meshfree method, the smoothed particle hydrodynamics (SPH) based on the support domain and the kernel function has developed rapidly in recent years. The traditional SPH method is easily implements parallel processing and has been applied in sound wave simulation. As a corrective method with higher accuracy than SPH, the acoustic propagation and scattering in the time domain is simulated with the corrective smoothed particle method (CSPM). Moreover, a hybrid meshfree-FDTD boundary treatment technique is utilized to represent different acoustic boundaries in the Lagrangian approach. In this boundary treatment technique, the parameter value of virtual particles is obtained with the FDTD method, which concerns truncation errors based on the Tayler series expansion. Soft, rigid, and Mur’s absorbing boundary conditions are developed to simulate sound waves in finite and infinite domain. Results of modeling acoustic propagation and scattering show that CSPM is accurate and convergence with exact solutions, and different acoustic boundaries are validated to be effective in the computation.

The development of computational acoustics allows simulation of sound generation and propagation in complex environment. In particular, meshfree methods are widely used to solve acoustics problems through arbitrarily distributed field points and approximation smoothness flexibility. As a Lagrangian meshfree method, smoothed particle hydrodynamics (SPH) method reduce the difficulty in solving problems with deformable boundaries, complex topologies, or multiphase medium. The traditional SPH method has been applied in acoustic simulation. This study presents the corrective smoothed particle method (CSPM), which is a combination of SPH kernel estimate and Taylor series expansion. The CSPM is introduced as a Lagrangian approach to improve accuracy in solving acoustic wave equations in the time domain. Moreover, a boundary treatment technique based on the hybrid meshfree and finite difference time domain (FDTD) method is proposed to represent different acoustic boundaries with particles. To model sound propagation in pipes with different boundaries, soft, rigid, and absorbing boundary conditions are built with this technique. Numerical results show that the CSPM algorithm is consistent and demonstrates convergence with exact solutions. Main computational parameters are discussed, and different boundary conditions are validated to be effective for benchmark problems in computational acoustics.

In recent years, there has been a growing interest in the preoperative modeling of fluid-structure interaction in the treatment of cerebral aneurysms. In this study, we investigate two cases involving laminar incompressible fluid flow interacting with hyperelastic materials (a vertical flap and aneurysm walls) under the effect of fluid flow. We present a finite element method (FEM) for these prototypical two-dimensional ($2D$) configurations, taking into account the complex flows and deformation of the variant's structure models. We use an Arbitrary Lagrangian-Eulerian (ALE) formulation in a continuum, fully monolithic coupled way, and discretize the fluid and solid domains using the quadratic LBB stable $P_2P_1$ finite element pair to approximate the displacement, velocity, and pressure spaces independently. The resulting discretized form of the nonlinear algebraic system is linearized using a variant of Newton's procedure, and the Jacobian matrices are approximated via the divided difference method. The resulting linear systems are solved using a direct solver MUMPS. We evaluate hydrodynamic forces such as drag and lift coefficients for nonlinear elastic material models, including Saint-Venant Kirchoff, Neo-Hookean, and Mooney-Rivlin (2 parameters) separately. To gain more physical insight into the problem, we verify the computed results by comparing the velocity, viscosity, and pressure fields. Our aim is to qualitatively analyze the changes in the mechanics of the wall behavior of the elasticity of the flap vs. the complexity of deformation, interface, details for prototypical flow situations, corresponding displacement, and deformation. We also provide a physical interpretation of blood flow magnitude and basic hemodynamic phenomena of wall shear stresses (WSS) in aneurysma and the range of deformation. Overall, this study presents a comprehensive approach to modeling fluid-structure interaction in cerebral aneurysms, which could help in the development of more effective treatment strategies.

In this paper, we study quantum gravity effect on the symmergent black hole which is derived from quadratic-curvature gravity. To do so, we use the Klein-Gordon equation which is modified by generalized uncertainty principle (GUP). After solving the field equations, we examine the symmergent black hole's tunneling and Hawking temperature. We explore the graphs of the temperature through the outer horizon to check the GUP influenced conditions of symmergent black hole stability. We also explain how symmergent black holes behave physically when influenced by quantum gravity. The impacts of thermal fluctuations on the thermodynamics of a symmergent black holes spacetime are examined. We first evaluate the model under consideration's thermodynamic properties, such as its Hawking temperature, angular velocity, entropy, and electric potential. We evaluate the logarithmic correction terms for entropy around the equilibrium state in order to examine the impacts of thermal fluctuations. In the presence of these correction terms, we also examine the viability of the first law of thermodynamics. Finally, we evaluate the system's stability using the Hessian matrix and heat capacity. It is determined that a stable model is generated by logarithmic corrections arising from thermal fluctuations.

In this paper, we investigate the accuracy and robustness of three classes of methods for solving two-phase incompressible flows on a staggered grid. Here, the unsteady two-phase flow equations are simulated by finite volumes and penalty methods using implicit and monolithic approaches (such as the augmented Lagrangian and the fully coupled methods), where all velocity components and pressure variables are solved simultaneously (as opposed to segregated methods). The interface tracking is performed with a Volume-of-Fluid (VOF) method, using the Piecewise Linear Interface Construction (PLIC) technique. The home code Fugu is used for implementing the various methods. Our target application is the simulation of two-phase flows at high density and viscosity ratios, which are known to be challenging to simulate. The resulting strategies of monolithic approaches will be proven to be considerably better suited for these two-phase cases, they also allow to use larger time step than segregated methods.

Mathematical modeling is the most important tool for constructing the theory of different kinds of two-phase flows. The review is devoted to the analysis of the advent of the approach of mathematical modeling of two-phase flows, where solid particles act mainly as the dispersed phase. The main problems and features of the study of gas-solid two-phase flows are excluded. The main characteristics of gas flows with solid particles are given and the classification of two-phase flows is developed on their basis. The Lagrangian and Euler approaches to modeling the motion of a dispersed phase (particles) are described. Much attention is paid to the consideration of numerical simulation methods that describe descriptions of turbulent gas flow at different hierarchical levels (RANS, LES, DNS), different levels of description of interphase interactions (one-way coupling (OWC), two-way coupling (TWC) and four-way coupling (FWC)), as well as at different levels of interface resolution (partial-point (PP) and particle-resolved (PR)). Excluded are examples of some studies carried out on the basis of the identified approaches, as well as their exclusion for mathematical modeling of various classes of gas-solid two-phase flows.

Micro-grids’ operations offer local reliability; in the event of faults or low voltage/frequency events on the utility side, micro-grids can disconnect from the main grid and operate autonomously while providing the continued supply of power to local customers. With the ever-increasing penetration of renewable generation, however, the operations of micro-grids become increasingly complicated because of the associated fluctuations of voltages. As a result, transformer taps are adjusted frequently, thereby leading to the fast degradation of expensive tap-changer transformers. In the islanding mode, the difficulties also come from the drop of voltage and frequency upon disconnecting from the main grid. To appropriately model the above, the nonlinear AC power flow constraints are necessary. Computationally, the discrete nature of tap-changer operations and the stochasticity caused by renewables add two layers of difficulty on top of a complicated AC-OPF problem. To resolve the above computational difficulties, the main principles of the recently-developed "l1-proximal" Surrogate Lagrangian Relaxation are extended. Testing results based on 9-bus system demonstrate the efficiency of the method to obtain the exact feasible solutions for micro-grid operations thereby avoiding approximations inherent to existing methods, while demonstrating that through the optimization, 1. the number of tap changes is drastically reduced, and 2. the method is capable of handling networks with meshed topologies.

Tissue engineering has the potential for repairing large bone defects, which impose a heavy financial burden on the public health. However, difficulties with O₂ delivery to the cells residing in the interior of tissue engineering scaffolds make it challenging to grow artificial tissues of clinically-relevant sizes. This study uses image-based simulation in order to provide insight into how to better optimize the scaffold manufacturing parameters, and the culturing conditions, in order to resolve the O₂ bottleneck. To do this, high resolution 3D X-ray images of two common scaffold types (salt leached foam and non-woven fiber mesh) are fed into a Lattice Boltzmann Method fluid dynamics and reactive Lagrangian Scalar Tracking mass transfer solvers. The obtained findings indicate that the scaffolds should have maximal surface area-to-solid volume ratios, for higher chances of the molecular collisions with the cells. Furthermore, the cell culture media should be flown through the scaffold pores as fast as practically possible (without detaching or killing the cells). Finally, we have provided a parametric sweep that maps how the molecular transport within the scaffolds is affected by variations in rates of O₂ consumption by the cells. Ultimately, the results of this study are expected to benefit the computer-assisted design of tissue engineering scaffolds and culturing experiments.

A controlled ammonia (NH₃) release experiment was performed at a grassland site to quantify the effect of dry deposition, at the field scale between the source and the receptors (NH₃ measurement locations), on the estimates of emission rates by means of inverse dispersion modelling. NH₃ was released for 3 hours at a constant rate of Q = 6.29 mg s⁻¹ from a grid of 36 orifices spread over an area of 250 m². The increase in line-integrated NH₃ concentration was measured with open-path optical miniDOAS devices at different locations downwind of the artificial source. Using a backward Lagrangian stochastic (bLS) dispersion model (bLSmodelR), the fraction of the modelled release rate to the emitted NH₃ (Q_bLS/Q) was calculated from the measurements of the individual instruments. Q_bLS/Q was found to be systematically lower than 1, on average between 0.69 and 0.91, depending on the location of the receptor. We hypothesized that NH₃ dry deposition to grass and soil surfaces was the main factor responsible for the observed depletion of NH₃ between source and receptor. A dry deposition algorithm based on a deposition velocity approach was included in the bLS modelling. Model deposition velocities were evaluated from a ‘big‑leaf’ canopy resistance analogy. Canopy resistances (generally termed Rc) that provided Q_bLS/Q = 1 ranged from 75 to 290 s m⁻¹, showing that surface removal of NH₃ by dry deposition can plausibly explain the original underestimation of Q_bLS/Q. The inclusion of a dry deposition process in dispersion modelling is crucial for emission estimates, which are based on concentration measurements of depositing tracers downwind of homogeneous area sources or heterogeneously distributed hot spots, such as e.g. urine patches on pastures in the case of NH₃.

Deterministic-stochastic empirical mode decomposition (EMD) is used to obtain low-frequency (non-diffusive, i.e., background velocity) and high-frequency (diffusive, i.e., eddies) components from a Lagrangian drifter‘s trajectory. Eddy characteristics are determined from the time series of eddy trajectories from individual Lagrangian drifter such as the eddy radial scale, eddy velocity scale, eddy Rossby number, and eddy-background kinetic energy ratio. A long-term dataset of the SOund Fixing And Ranging float time series obtained near the California coast by the Naval Postgraduate School from 1992 to 2004 at depth between 150 and 600 m (http://www.oc.nps.edu/npsRAFOS/) is used as an example to demonstrate the capability of the deterministic-stochastic EMD.

We delve into nonrelativistic quantum electrodynamics within a background magnetic field, ensuring gauge invariance via a vector potential. Our exploration extends to the Lagrangian, which incorporates electron self-interactions and electromagnetic field interactions. Through the path integral formalism, we elucidate the effective action, with a specific focus on the photon propagator and screening effects. Examining density and current components reveals Laurent expansions in the hydrodynamic limit. The equations of motion in real space are derived, leading to discussions on phenomena such as Poiseuille-like flow and the Navier-Stokes equation. To ground our theoretical framework at practical contexts, we consider applications such as the fluid dynamics at binary star systems. Depending on the velocity of the expelled fluid, our analysis allows for both nonrelativistic and relativistic treatments, providing a versatile tool for understanding the intricacies of fluid behavior at astrophysical scenarios. Additionally, we recognize that at the presence of a magnetic field, magnetohydrodynamics becomes crucial. While our current focus is at nonrelativistic quantum electrodynamics, the insights gained here contribute to a broader understanding, offering a comprehensive foundation for quantum electrodynamics analysis at diverse physical systems.