Closed-form Analytical Solution Research Articles

An analytical study on the influence of magnetic field-dependent viscosity on viscous and ohmic dissipative second-grade fluid boundary layers

This article aims to investigate the magnetohydrodynamic (MHD) boundary layer flow of a second-grade non-Newtonian fluid over a horizontally stretching sheet, considering magnetic field-dependent (MFD) viscosity, as well as viscous and Ohmic dissipations. Analytical solutions for previously unexplored momentum and heat transfer equations involving MFD viscosity are derived. Two boundary conditions, namely Prescribed Surface Temperature (PST) and Prescribed Heat Flux (PHF), are taken into account. Governing dimensional partial differential equations (PDEs) are transformed into non-dimensional ordinary differential equations (ODEs) using similarity transformations. Closed-form analytical solutions for flow are derived, considering stretching velocity, MFD viscosity, second-grade fluid properties, and suction impacts. Heat transfer equations are transformed into Gauss-hypergeometric form, yielding solutions in terms of confluent hypergeometric functions. Analytical expressions for skin friction, local Nusselt number, and non-dimensional wall temperature are derived. A unique solution is obtained for the flow equation. It is found that both second-grade and magnetic viscosity parameters expand the momentum boundary layer while the magnetic parameter reduces thickness. Thermal boundary layer thickness increases with a higher second-grade parameter, magnetic viscosity, and Eckert number. Moreover, analytical solutions are validated against published results for a special case.

Closed-form asymptotic solution for the transport of chlorine concentration in composite pipes

Recently, some researchers have revisited the analysis of chlorine transportation in cylindrical pipes by deploying a coupling between the Laplace transform method and the complex analysis’ residue approach for inverting complex integrals. This method yielded interesting results after the incorporation of root-finding numerical schemes. Thus, away from incorporating numerical tools, the present study makes consideration of the same formulation of chlorine transport in a single-layered pipe and further extends it to the case of a bi-layered pipe using the hybrid of the Laplace transform method and the asymptotic approximations method. The need for asymptotic approximations for the modified Bessel functions, which arise in the reduced ordinary differential equations, necessitates the quest for closed-form analytical solutions, which are largely considered benchmark solutions for numerical investigations. Moreover, the obtained closed-form asymptotic solutions have been examined graphically; where it was observed that both the radial diffusion coefficient η and the spatial radial variable are contributory in the transport of chorine concentration in the media.

Estimating Position-Dependent and Anisotropic Diffusivity Tensors from Molecular Dynamics Trajectories: Existing Methods and Future Outlook.

Confinement can substantially alter the physicochemical properties of materials by breaking translational isotropy and rendering all physical properties position-dependent. Molecular dynamics (MD) simulations have proven instrumental in characterizing such spatial heterogeneities and probing the impact of confinement on materials' properties. For static properties, this is a straightforward task and can be achieved via simple spatial binning. Such an approach, however, cannot be readily applied to transport coefficients due to lack of natural extensions of autocorrelations used for their calculation in the bulk. The prime example of this challenge is diffusivity, which, in the bulk, can be readily estimated from the particles' mobility statistics, which satisfy the Fokker-Planck equation. Under confinement, however, such statistics will follow the Smoluchowski equation, which lacks a closed-form analytical solution. This brief review explores the rich history of estimating profiles of the diffusivity tensor from MD simulations and discusses various approximate methods and algorithms developed for this purpose. Besides discussing heuristic extensions of bulk methods, we overview more rigorous algorithms, including kernel-based methods, Bayesian approaches, and operator discretization techniques. Additionally, we outline methods based on applying biasing potentials or imposing constraints on tracer particles. Finally, we discuss approaches that estimate diffusivity from mean first passage time or committor probability profiles, a conceptual framework originally developed in the context of collective variable spaces describing rare events in computational chemistry and biology. In summary, this paper offers a concise survey of diverse approaches for estimating diffusivity from MD trajectories, highlighting challenges and opportunities in this area.

A neural network inverse problem solution for a MIMO magnetic levitation actuator

A multiple-input and multiple-output planar actuator is proposed, which utilizes mechanically driven stator magnet arrays to levitate a permanent magnet mover. A state of levitation and actuation is obtained by mechanically altering the orientation of the stator magnets to control the forces and torques on the mover. A challenge for the design and control of the actuator is inverting the relationship between the force and stator magnet rotation angles, as there is no closed-form analytical solution. In this study, a feed-forward neural network is applied to model the forward relation between stator magnet angle input and a force and torque output to reduce the forward computation time for the design process and for error estimation in real-time applications. Additionally, the neural network is considered for inverting the solution for a motion profile sampled at 1000 Hz. The developed forward model is able to calculate the forces and torques on the mover a factor 10 faster than the equivalent charge or Fourier model with an absolute error of 3 mN and 0.1 mNm for the forces and torques, respectively, and a feed-forward neural network is able to accurately learn an inverse solution for small motion profiles.

Fermion-antifermion pairs in a magnetized space-time with non-zero cosmological constant

We analyze a fermion-antifermion pair within the (2+1)-dimensional Bonnor-Melvin magnetic (BMM) spacetime featuring a nonzero cosmological constant (Λ>0). To do this, we rigorously investigate analytical solutions to the corresponding form of a fully-covariant two-body Dirac equation. We derive a non-perturbative wave equation and obtain an exact closed-form analytical solution for the system under consideration. In this magnetic background spacetime, described by a surface with positive Gaussian curvature, we noticed that fermion-antifermion pairs display a collective behavior resembling that of a single quantum oscillator, carrying a total rest mass of the particles and oscillating with a frequency (ωo∝Λ), even in the absence of mutual interaction between the particles. This indicates a compelling correlation between the cosmological constant and the energy spectrum of fermion-antifermion systems, implying connections between quantum realms and cosmology.

Anisotropic functionally graded nano-beam models and closed-form solutions in plane gradient elasticity

This study delves into the investigation of exact analytical solutions for the plane stress and displacement fields within linear homogeneous anisotropic nano-beam models of gradient elasticity. It focuses on solving the Helmholtz equation, which encompasses a second-order non-homogeneous linear partial differential equation in plane gradient elasticity theory, utilizing polynomial series-type solutions. The analysis centers on the utilization of gradient Airy stress functions to derive stress fields in the gradient theory. The research yields closed-form analytical solutions for Airy stress functions, stress, and strain fields in both classical and gradient theories. The study considers five distinct types of two-dimensional functionally graded cantilever beams with various boundary conditions: a cantilever anisotropic nano-beam subjected to a concentrated force at the free end, a cantilever anisotropic nano-beam under a uniform load, a simple anisotropic nano-beam under a uniform load, a propped cantilever anisotropic nano-beam under a uniform load, and a fixed end cantilever anisotropic nano-beam under a uniform load. General analytical solutions for the gradient stress and displacement fields of two-dimensional and one-dimensional anisotropic nano-beams under different boundary conditions are provided. The study showcases significant strain gradient size effects at the nano-scale through the derived analytical solutions for anisotropic beams. Additionally, it demonstrates that the strain gradient theory results for the limit case of the gradient coefficient c precisely align with results for isotropic and anisotropic materials in elasticity theory and classical theory. Furthermore, real-world applications are discussed, considering stress and displacement fields in real anisotropic materials such as TiSi2 single crystals and orthotropic materials, which are special cases of anisotropic materials like wood and epoxy as documented in the literature.

Large displacement model for the nonlinear behaviour of the face plate component

This paper deals with the face plate component in steel joints and addresses the full characterization of the nonlinear behaviour of the face plate component. An equivalent beam strip model is proposed that tackles the connection of a beam to the column web of open I-sections in a minor axis joint or the face of tubular columns, covering endplate or fin plate joint typologies. Closed-form analytical solutions are obtained for the elastic large displacement and the elastic-plastic large displacement behaviour of the equivalent beam strip. Criteria for the establishment of the design resistance using the continuous strength method are also proposed. It was concluded that the model is easy to apply, was validated against a large parametric study using finite element beam models, leads to accurate solutions and demonstrates the need to consider membrane effects in design.

The Pathologically Evolving Aggregation-State of Cells in Cancerous Tissues as Interpreted by Fractal and Multi-Fractal Dispersion Theory in Saturated Porous Formations.

A recent author's fractal fluid-dynamic dispersion theory in porous media has focused on the derivation of the associated nonergodic (or effective) macrodispersion coefficients by a 3-D stochastic Lagrangian approach. As shown by the present study, the Fickian (i.e., the asymptotic constant) component of a properly normalized version of these coefficients exhibits a clearly detectable minimum in correspondence with the same fractal dimension (d ≅ 1.7) that seems to characterize the diffusion-limited aggregation state of cells in advanced stages of cancerous lesion progression. That circumstance suggests that such a critical fractal dimension, which is also reminiscent of the colloidal state of solutions (and may therefore identify the microscale architecture of both living and non-living two-phase systems in state transition conditions) may actually represent a sort of universal nature imprint. Additionally, it suggests that the closed-form analytical solution that was provided for the effective macrodispersion coefficients in fractal porous media may be a reliable candidate as a physically-based descriptor of blood perfusion dynamics in healthy as well as cancerous tissues. In order to evaluate the biological meaningfulness of this specific fluid-dynamic parameter, a preliminary validation is performed by comparison with the results of imaging-based clinical surveys. Moreover, a multifractal extension of the theory is proposed and discussed in view of a perspective interpretative diagnostic utilization.

Approximating accelerator impedances with resonator networks

It is common in the accelerator community to use the impedance of accelerator components to describe wake interactions in the frequency domain. However, it is often desirable to understand such wake interactions in the time domain in a general manner for excitations that are not necessarily Gaussian in nature. The conventional method for doing this involves taking the inverse Fourier Transform of the component impedance, obtaining the Green's Function, and then convolving it with the desired excitation distribution. This method can prove numerically cumbersome, for a convolution integral must be evaluated for each individual point in time when the wake function is desired. An alternative to this method would be to compute the wake function analytically, which would sidestep the need for repetitive integration. Only a handful of cases, however, are simple enough for this method to be tenable. One of these cases is the case where the component in question is an RLC resonator, which has a closed-form analytical wake function solution. This means that a component which can be represented in terms of resonators can leverage this solution. As it happens, common network synthesis techniques may be used to map arbitrary impedance profiles to RLC resonator networks in a manner the accelerator community has yet to take advantage of. In this work, we will use Foster Canonical Resonator Networks and partial derivative descent optimization to develop a technique for synthesizing resonator networks that well approximate the impedances of real-world accelerator components. We will link this synthesis to the closed-form resonator wake function solution, giving rise to a powerful workflow that may be used to streamline beam dynamics simulations.

Analytical solution for determining wall closure associated with stope excavation underneath sill mat constructed by cemented backfill

In the underhand cut-and-fill mining method, a sill mat (i.e. an artificial horizontal pillar) constructed by cemented backfill is essential to prevent mine workers from being directly exposed under problematic rock roofs. A critical issue is to determine the minimum required strength of the sill mat to ensure a safe and cost-effective design. Until now, Mitchell's analytical solution is the only available option, considering two stiff and immobile rock walls. Unavoidable rock wall closure associated with stope excavation below the sill mat was neglected. This, along with other undefined parameters, explains why Mitchell's solution is rarely used in sill mat design. A new analytical solution for determining the minimum required strength of the sill mat accounting for wall closure is necessary. In this study, a closed-form analytical solution for estimating rock wall closure generated by stope excavation below a sill mat is developed by using Salamon's and Flamant's models. The proposed analytical solution does not contain any coefficients of correction or calibration. Despite several assumptions (or somewhat of oversimplifications) necessary to render a simple analytical solution possible, good agreements are obtained between the rock wall closures predicted by applying the proposed analytical solution and those obtained numerically with FLAC3D for many cases with arbitrarily chosen geometrical and material parameters. The proposed analytical solution is therefore validated and can be used to evaluate the rock wall closure generated by stope excavation below a sill mat.

Prediction of horizontal and vertical turning maneuvers of a generic submarine model by closed-form solutions and CFD

Turning maneuvering performance of a generic submarine in both the horizontal and vertical planes is investigated using closed-form (analytical) solutions and the Computational Fluid Dynamics (CFD) method. The aim is to provide a simple yet improvable method throughout the design phase of an underwater vehicle. Both the turning circle maneuver in the horizontal plane and the rising maneuver in the vertical plane are performed. Vehicle trajectories are obtained using closed-form solutions and a linearized maneuvering model with a limited set of hydrodynamic coefficients computed from the CFD model. The derivation of the closed-form solution is based on the basic mechanical principles of steady circular motion, while the hydrodynamic coefficients are determined via steady CFD simulations using a single grid topology. The steady turning diameter serves as the figure of merit in a turning circle maneuver. For the rising maneuver, trajectories predicted by both methods are directly compared. The methods are found to be in good agreement for both maneuvers. Maneuvering predictions made with the linear coefficient set are also open to development. As the design of the submarine evolve, nonlinear and coupled effects can be included in the mathematical model by incorporating related hydrodynamic coefficients to obtain more accurate results.

Asymptotically consistent and computationally efficient modeling of short-ranged molecular interactions between curved slender fibers undergoing large 3D deformations

This article proposes a novel computational modeling approach for short-ranged molecular interactions between curved slender fibers undergoing large 3D deformations, and gives a detailed overview how it fits into the framework of existing fiber or beam interaction models, either considering microscale molecular or macroscale contact effects. The direct evaluation of a molecular interaction potential between two general bodies in 3D space would require to integrate molecule densities over two 3D volumes, leading to a sixfold integral to be solved numerically. By exploiting the short-range nature of the considered class of interaction potentials as well as the fundamental kinematic assumption of undeformable fiber cross-sections, as typically applied in mechanical beam theories, a recently derived, closed-form analytical solution is applied for the interaction potential between a given section of the first fiber (slave beam) and the entire second fiber (master beam), whose geometry is linearly expanded at the point with smallest distance to the given slave beam section. This novel approach based on a pre-defined section–beam interaction potential (SBIP) requires only one single integration step along the slave beam length to be performed numerically. In addition to significant gains in computational efficiency, the total beam–beam interaction potential resulting from this approach is shown to exhibit an asymptotically consistent angular and distance scaling behavior. Critically for the numerical solution scheme, a regularization of the interaction potential in the zero-separation limit as well as the finite element discretization of the interacting fibers, modeled by the geometrically exact beam theory, are presented. In addition to elementary two-fiber systems, carefully chosen to verify accuracy and asymptotic consistence of the proposed SBIP approach, a potential practical application in form of adhesive nanofiber-grafted surfaces is studied. Involving a large number of helicoidal fibers undergoing large 3D deformations, arbitrary mutual fiber orientations as well as frequent local fiber pull-off and snap-into-contact events, this example demonstrates the robustness and computational efficiency of the new approach.

Solutions for the unsteady motion of porous elastic solids within the context of an implicit constitutive theory

In this work we investigate well-posedness for the partial differential equation stemming from the balance of linear momentum for an implicit constitutive relation that can describe the response of porous elastic solids whose material moduli depend on the density. We study heteroclinic travelling waves and obtain closed form analytic solutions for the resulting ordinary differential equation. Finally, we consider some special solutions for the partial differential equation and solve the resulting equations numerically.

Minkowskian Approach to the Pseudorange Navigation Equations

Our starting point is the covariant coordinate transformation equation of a relativistic positioning system in Minkowski space–time that maps the receiver’s emission coordinates (proper times broadcast by the emitters) to its coordinates in some inertial reference frame. Bancroft’s analytical (closed-form) solution to the basic pseudorange navigation equations with four emitters is recovered, and the subjacent elements are geometrically interpreted. The case of four static beacons is analysed as a clarifying situation.

Influence of Chemical and Radiation on an Unsteady MHD Oscillatory Flow using Artificial Neural Network (ANN)

This paper delves into the intricate interplay between chemical and thermal radiation in the context of an unstable magnetohydrodynamic(MHD) oscillatory flow through a porous medium. The fluid under investigation is presumed to be incompressible, electrically conductive, and radiating with the additional influence of a homogeneous magnetic field applied perpendicular to the channel’s plane. Analytical closedform solutions are derived for the momentum, energy, and concentration equations providing a comprehensive understanding of the system’s behavior. The investigation systematically explores the impact of various flow factors, presenting their effects through graphical representations. The governing partial differential equations (PDE) of the boundary layer are transformed into a set of coupled nonlinear ordinary differential equations (ODE) using a closed-form method. Subsequently, an artificial neural network (ANN) is applied to these ODEs, and the obtained results are validated against numerical simulations. The temperature profiles exhibit oscillatory behavior with changes in the radiation parameter (N), revealing insights into the system’s dynamic response. Furthermore, the paper uncovers that higher heat sources lead to increased temperature profiles. Additionally, concentration profiles demonstrate a decrease with escalating chemical reaction parameters, with a reversal observed as the Schmidt number (Sc) increases. This study highlights the efficacy of an ANN model in providing highly efficient estimates for heat transfer rates from an engineering standpoint. This innovative approach leverages the power of artificial intelligence to enhance our understanding of complex fluid magnetohydrodynamics and porous media flows.

Higher-order composition of short- and long-period effects for satellite analytical ephemeris computation

The construction of an analytic orbit theory in closed form of the eccentricity that takes into account the main effects of the Geopotential is notably simplified when splitting the removal of periodic effects in several stages. Conversely, this splitting of the closed-form analytical solution into several transformations reduces the evaluation efficiency for dense ephemeris output. However, the advantage is twofold when the different parts of the mean-to-osculating transformation are composed into a single transformation. To show that, Brouwer’s solution is extended to the second order of the zonal harmonic of the second degree by the sequential elimination of short and long period terms. Then, the generating functions of the different transformations are composed into a single one, from which a single mean-to-osculating transformation is derived. The new, unique transformation notably speeds up the evaluation process, commonly improving evaluation efficiency by at least one third with respect to the customary decomposition of the analytical solution into three different parts.

Optical second harmonic generation in anisotropic multilayers with complete multireflection of linear and nonlinear waves using ♯SHAARP.ml package

Optical second harmonic generation (SHG) is a nonlinear optical effect widely used for nonlinear optical microscopy and laser frequency conversion. Closed-form analytical solution of the nonlinear optical responses is essential for evaluating materials whose optical properties are unknown a priori. A recent open-source code, ♯SHAARP.si, can provide such closed form solutions for crystals with arbitrary symmetries, orientations, and anisotropic properties at a single interface. However, optical components are often in the form of slabs, thin films on substrates, and multilayer heterostructures with multiple reflections of both the fundamental and up to ten different SHG waves at each interface, adding significant complexity. Many approximations have therefore been employed in the existing analytical approaches, such as slowly varying approximation, weak reflection of the nonlinear polarization, transparent medium, high crystallographic symmetry, Kleinman symmetry, easy crystal orientation along a high-symmetry direction, phase matching conditions and negligible interference among nonlinear waves, which may lead to large errors in the reported material properties. To avoid these approximations, we have developed an open-source package named Second Harmonic Analysis of Anisotropic Rotational Polarimetry in Multilayers (♯SHAARP.ml). The reliability and accuracy are established by experimentally benchmarking with both the SHG polarimetry and Maker fringes using standard and commonly used nonlinear optical materials as well as twisted 2-dimensional heterostructures.

The formation of invariant optical soliton structures to electric-signal in the telegraph lines on basis of the tunnel diode and chaos visualization, conserved quantities: Lie point symmetry approach

This study critically examines the nonlinear Lonngren wave model, a mathematical representation used to describe waves in shallow water environments such as beaches, lakes, and rivers, as well as electrical signals in semiconductor materials like tunnel diodes. The symmetry generators take into account the Lie invariance criteria. A three-dimensional Lie algebra can be obtained by using the suggested methodology. It links translation symmetries in space and time, respectively, with conservation of mass and energy. The third symmetry pertains to scaling or dilation. Recognizing the constraints of the inverse scattering transform in addressing the Cauchy problem for this equation, the investigation adopts new auxiliary equation approach to ordinary differential equations. This approach yields closed-form analytical solutions for the wave model. The obtained solutions include a plane solution, mixed hyperbolic solution, periodic and mixed periodic solutions, mixed trigonometric solution, smooth soliton, trigonometric solution, mixed singular solution, complex solitary and singular solution. Furthermore, the study develops the Hamiltonian function of the dynamical system to discuss the total energy of the system in terms of its position and momentum. Chaos analysis is performed to describe periodic, quasi-periodic, and chaotic behaviours, and a sensitivity analysis highlights the model’s sensitivity to different wave numbers. Graphical representations of selected solutions are presented in both 3-D and 2-D, illustrating the behaviour under specific physical conditions. The study concludes by identifying a comprehensive set of local conservation laws for the nonlinear Lonngren wave equation, applicable to arbitrary constant coefficients, using conservation laws multipliers.

A Numerical and Closed-Form Analytical Solution of the Global Buckling Critical Load of Tall Buildings Including Soil Flexibility: Timoshenko Beam

A Numerical and Closed-Form Analytical Solution of the Global Buckling Critical Load of Tall Buildings Including Soil Flexibility: Timoshenko Beam

The Allen–Cahn model with a time-dependent parameter for motion by mean curvature up to the singularity

In this article, we investigate the temporal evolution of arbitrary, simple, closed two-dimensional (2D) and three-dimensional (3D) interfaces under motion driven by mean curvature up to a singularity. To facilitate this investigation, we propose a novel Allen–Cahn (AC) model with a time-dependent interfacial thickness parameter. The original AC equation was developed to model the phase separation of a binary mixture. It is well known that a level set or interface of the solution of the AC equation obeys the dynamics of motion by curvature as the interfacial thickness parameter approaches zero. Generally, it is difficult to find a closed-form analytic solution of the AC equation with any initial condition. Therefore, we need to estimate the solution of the AC equation through computational approaches such as finite difference method (FDM), finite element method (FEM), Fourier-spectral method (FSM), and finite volume method (FVM). Any simple, closed curves and convex surfaces eventually shrink to a point due to motion by mean curvature. Therefore, it becomes necessary to use adaptive mesh techniques to resolve this small size problem. However, even though we use adaptive mesh techniques, we still confront the relatively thick interfacial transition layer when the curves or interfaces become very small. To avoid this problem, we can start with a very small mesh size for a small value of the interfacial parameter, which results in an extremely high computational cost even when using adaptive mesh techniques. To resolve these issues, we present the AC equation with a time-dependent interfacial parameter and develop an adaptive mesh refinement system. To show the superior performance of the proposed mathematical equation and its computational algorithm, we present various numerical experiments and investigate the motion by mean curvature up to the singularity of curves in 2D space and interfaces in 3D space.