Formulation Of Boundary Conditions Research Articles

Thermodynamics of Irreversible Processes: Fundamental Constraints, Representations, and Formulation of Boundary Conditions

Starting from the analysis of the lack of positivity of the Cattaneo heat equation, this work addresses the thermodynamic relevance of the positivity constraint in irreversible thermodynamics, that is at least as significant as the entropic constraints. The fulfillment of this condition in hyperbolic models leads to the parametrization of the concentration fields with respect to internal variables associated with the microscopic dynamics. Using Brownian motion theory as a landmark example for deriving macroscopic transport equations from the equations of motion at the particle/molecular level, we discuss two typical problems involving hydrodynamic interactions at the microscale: surface chemical reactions at a solid interface of a diffusing reactant, and mass-balance equations in a complex viscoelastic fluid, in which the physics of the interaction leads either to overcoming the parabolic diffusion model or to considering the parametrization of the concentration with respect to the degrees of freedom associated with the relaxation dynamics of the solvent fluid.

Comparison of Computational Methods and Approaches Applied in Formulation of Boundary Conditions in Lagrange’s Ballistic Problem

This paper compares algorithms for constructing a differential solution to the gas dynamics equations in the surrounding of a moving boundary. In order to compare the algorithms, the Lagrange problem, also known as the piston problem, was chosen as a test problem. A comparison was made between the author’s algorithm based on the method of characteristics and the classical approach with a fictitious cell using the solution of the Riemann problem. Calculations were carried out for the case of subsonic and supersonic gas motion. Comparisons were made on a fixed grid with a dynamically expanding cell (Euler grid) and on a uniformly expanding grid (Arbitrary Lagrangian-Eulerian – ALE grid). The problem was solved using numerical schemes of second order in time and space Lax - Wendroff type: the Richtmyer scheme and the MacCormack scheme. The results of calculations using the characteristics method were used as a benchmark solution. The comparative analysis carried out allows conclusions to be drawn regarding the accuracy of the individual approaches, as well as the scope of their applicability.

Force relaxation response in the poroelastic axisymmetric Boussinesq problem for an indenter of arbitrary profiles I: Permeable and impermeable surface drainage conditions

In this study, we have derived a theoretical solution for the poroelastic axisymmetric Boussinesq problem when an indenter with arbitrary profiles is subjected to step displacement loading. The analysis is conducted within the framework of Biot’s theory, employing the McNamee–Gibson displacement function method. We explore two distinct scenarios for surface drainage boundary conditions: one where the surface allows full permeability and another where it is fully impermeable. The formulation of the mechanical boundary condition at the surface relies on an equation that establishes the relationship between the depth of indentation and the contact radius. Our solution encompasses two aspects: force asymptotes at the undrained and drained limits, and a normalized force relaxation that accounts for transient responses. Specific findings are presented for three indenter shapes: paraboloidal, conical, and cylindrical. Of particular significance, we demonstrate that for these specific indenter shapes, the normalized force relaxation at ω=0 can be expressed in closed form in the Laplace domain. In addition to this analytical breakthrough, we conducted numerical simulations employing actual indenters. Remarkably, the normalized force relaxation curves from actual indenters show negligible dependence on material properties and closely align with our closed-form solution at ω=0. This finding suggests that the closed-form solution could act as a universal master curve capable of characterizing the normalized force relaxation response in instances of a real indenter, irrespective of material properties.

The Einstein equation solution inside a ball with uniform density

A great number of solutions of the Einstein field equation are known. They describe the gravitational field in the empty space-time, in the space-time with electromagnetic field and for a ball filled with a liquid under pressure. The present work is devoted to gravitational field generated by some mass distribution. One of the simplest cases is considered, when mass is uniformly distributed inside a ball and is not moving. The boundary problem for the Einstein equation is formulated. Solution outside the ball is the Schwartzschild solution in vacuum. The coordinates at which the Schwartzschild solution is written are different from the coordinates used in equations for components of the metric tensor inside the ball. Relations between internal and external coordinates are found on the ball surface. They allow to use the Schwartzschild solution for formulation of boundary conditions for internal solution. The solution of the boundary problem is found for the case of weak field. This solution can be used as an example in the analysis of laws of conservation for the gravitational field, in which interaction of mass with field generated by the mass gives a contribution to momentum and energy of the gravitational field.

A highly robust gas network simulation approach through an inherently solvable problem formulation for network states far from intended design points

Pipeline networks are an efficient and widespread transportation system for supplying natural gas as well as increasingly green gas to Europe. However, they are exposed to risks arising from environmental factors, accidents, crime and political issues. This work contributes to the assessment of gas transmission networks and enables a targeted improvement of their resilience. Therefore, the objective of this work is to develop a simulation tool that enables gas transmission system operators (TSOs) to identify weak spots, e.g. potential bottle necks, in their pipeline networks. To this end, a steady-state gas network simulation approach is developed that provides fast results for gas supply deficiencies for complete sets of multi-event disruption scenarios, enabling the identification of the most important network elements. This paper describes how to implement a computationally fast, numerically robust, and physically accurate gas network simulation and assessment tool. By solving mass and momentum balance equations, it accounts for key network elements such as flow control valves, pressure control valves, compressor stations, gas sources, and gas consumers. The central feature that distinguishes this approach from classical steady state solvers is the robustness of the gas flow calculation through the choice of a problem boundary condition formulation, which results in an inherently solvable system of equations. This feature is key in ensuring mathematical convergence of predictions for gas network states that strongly deviate from their original design point, where boundary conditions must not become mathematically invalid but need to reflect physically reasonable behavior. The capabilities of the system are demonstrated using a fictitious gas network and a representative national gas transmission network. Finally, the limits of applicability are outlined based on the size of the network under consideration and the types of analyses.

Study on axial response of semi-parallel wire suspender with spherical plain bearings

Semi-parallel wire suspenders with spherical plain bearings have untwisting under axial force, and the collision in the untwist process bring adverse effects on the structure. To investigate the axial mechanical response of semi-parallel wire suspenders with spherical plain bearings, the equilibrium equations of suspender are deduced by assuming the distribution of steel wire after deformation, the forms of boundary conditions are distinguished according to the use of spherical plain bearing, and the linear equation expressed by equivalent stiffness coefficients is deduced. The Costello theory and finite element modeling are used to verify the proposed calculation method in this paper, and the influence of suspender parameters on the stiffness equivalent coefficient is analyzed. The results show that the proposed calculation method is suitable for multi-layer semi-parallel wire suspenders; the twist angle and the total number of layers of wires have the greatest effect on the stiffness equivalence factor; when the twisting angle is close to 0, the torsional stiffness of suspender is 0.010–0.21 of that of the steel rod with the same section; the axial tensile stiffness of the conventional suspender is basically the same as that of the steel rod with the same section, but the torsional stiffness is only 0.012–0.023 of the steel rod with the same section, and the tension-torsion coupling effect cannot be ignored. The research results provide a concise method to calculate the axial response of semi-parallel wire suspenders, which lays a foundation for the engineering treatment of the untwisting problem of semi-parallel wire suspender using spherical plain bearings.

A CPC fractional model of the heat and mass transport mechanism in Carbon nanotubes with slip effects on velocity

A fractional technique is used to evaluate the temperature, mass, and velocity flow of single and double wall CNTs over a vertical plate. Slip boundary conditions and applied magnetic force are addressed. Human blood is used to examine how base fluid behaves. Applying the proper dimensionless variables results in the dimensionless formulation of initial and boundary conditions related to the governed dimensional concentration, momentum, and energy equations. The Laplace transform technique is used to resolve the dimensionless governing partial differential equations and get the solutions. The constant proportional Caputo (CPC) time-fractional derivative is a unique class of fractional model used in the simulation technique. The fundamental definitions are used to support the said model first. Using MWCNTs and SWCNTs in comparison to the flow characteristics, a thermal and mass study is given. The heat and mass transfer processes for single-walled carbon nanotubes (SWCNTs) have been shown to typically be progressive. The momentum profile decreases as the fractional variables rise. Multi-walled carbon nanotubes (MWCNTs) show more progressive velocity control as a result of the magnetic parameter. Graphs demonstrate the influence of embedding factors on the velocity, energy, and concentration profiles.

A 2D Hydraulic Simulation Model Including Dynamic Piping and Overtopping Dambreach

Numerical simulation of unsteady free surface flows using depth averaged equations that consider the presence of initial discontinuities has been often reported for situations dealing with dam break flow. The usual approach is to assume a sudden removal of the gate at the dam location. Additionally, in order to prevent any kind of dam risk in earthen dams, it is very interesting to include the possibility of a progressive dam breach leading to dam overtopping or dam piping so that predictive hydraulic models benefit the global analysis of the water flow. On the other hand, when considering a realistic large domain with complex topography, a fine spatial discretization is mandatory. Hence, the number of grid cells is usually very large and, therefore, it is necessary to use parallelization techniques for the calculation, with the use of Graphic Processing Units (GPU) being one of the most efficient, due to the leveraging of thousands of processors within a single device. The aim of the present work is to describe an efficient GPU-based 2D shallow water flow solver (RiverFlow2D-GPU) supplied with the formulation of internal boundary conditions to represent dynamic dam failure processes. The results obtained indicate that it is able to develop a transient flow analysis taking into account several scenarios. The efficiency of the model is proven in two complex domains, leading to >76× faster simulations compared with the traditional CPU computation.

Dirichlet and Neumann boundary conditions for a lattice Boltzmann scheme for linear elastic solids on arbitrary domains

We propose novel, second-order accurate boundary formulations of Dirichlet and Neumann boundary conditions for arbitrary curved boundaries, within the scope of our recently introduced lattice Boltzmann method for linear elasticity. The proposed methodology systematically constructs and analyzes the boundary formulations on the basis of the asymptotic expansion technique; thus, it is expected to be of general applicability to boundary conditions for lattice Boltzmann formulations, well beyond the focus of this paper. For Dirichlet boundary conditions, a modified version of the bounce-back method is developed, whereas the Neumann boundary conditions are based on a novel generalized ansatz proposed in this work. The Neumann boundary formulation requires information from one additional neighbor node, for which a finite difference algorithm is introduced that enables the handling of very general boundary geometries. The theoretical derivations are verified by convergence studies using manufactured solutions, as well as by comparison with the analytical solution of the classical plate with hole problem.

Evaluation of Accelerator Pedal Strength under Critical Loads Using the Finite Element Method

The core idea of the research consists in a formulation of boundary conditions of a mechanical accelerator pedal’s strength in an Ansys environment, whose conditions are equivalent to full-scale tests under the critical loads defined by the UNECE’s Regulation No. 13. The lack of regulatory requirements for the strength of pedal types other than brake pedals is a major gap in vehicle certification, especially when it comes to agricultural machinery. In such cases, the authors suggest being guided by UNECE R 13 regarding the strength of the accelerator and other types of pedals and checking their behavior under loads of at least 1000 N. The real value of the yield strength of the material (Silumin 4000) is very important, both in the physical real-life experiments and in FEA simulation. The critical case of a short-term shock loading of the pedal in its extreme position has been considered separately. With the help of the Ansys Explicit Dynamics module, results of a pedal’s behavior were obtained; it lost its integrity and suffered destruction. It is also necessary to check the intermediate stress values depending on the loads for direct and hybrid tasks using the Transient Structural module in order to estimate other critical cases of the pedal behavior.

Numerical Analysis of Unsteady Heat Transfer in the Chamber in the Piston Hybrid Compressor with Regenerative Heat Exchange

The adjoint thermodynamic and heat exchange processes in the new class piston compressor with regenerative heat exchange are considered in the paper. The implicit tridiagonal matrix algorithm is implemented to study the unsteady thermal conductivity in the cylinder–piston group. After the formulation of proper initial and boundary conditions, the optimal relationship between temporal and spatial discretization steps is determined. Two different time steps are used in the numerical solution of the two-way coupled fast thermodynamic and slow heat exchange models. The relationship between those time steps is determined as well. The conducted numerical experiments allow the analysis of the dynamics of heat exchange in the cylinder–piston group, temperature variation in different parts of the cylinder–piston group, the impact of the heat transfer processes on isothermal and adiabatic efficiency, the impact of heat exchange dynamics on the thermodynamic cycle, as well as other thermodynamic and energetic effects.

Effect of boundary conditions on particle spurious movement in Smoothed Particle Hydrodynamics method

SPH method due to its Lagrangian particle nature suffers from non-physical particle motion resulting from its numerical properties. We study these effects associated with solid walls realization.We used stationary test case to investigate behavior of multiple formulations of boundary conditions for solid walls. Multi-layer boundaries with ghost particles are considered (Dynamic Boundary Conditions, Modified Dynamic Boundary Conditions) as well as one-layer boundary conditions based on discretization of boundary terms and usage of boundary integrals. Some aspects of each realization with respect to the stationary case and particle spurious motion are discussed. Moreover we examined the levels of spurious motion caused by different sources as the free surface, straight walls, sharp corners, … The results show that for the boundary conditions good enough the most significant portion of spurious motion arises in place where the water surface connects to the wall and the spurious motion caused by this interaction is several orders of magnitude higher, then artificial motion caused by other sources.

LaBCof: Lattice Boltzmann boundary condition framework

The boundary condition is one of the most challenging aspects of lattice Boltzmann simulations. One of the reasons that make it difficult to implement a new boundary condition is the dependency of formulation on the orientation of the boundary nodes. So, any boundary condition formulation should be derived, implemented, and verified in terms of the boundary orientation. In this study a new framework for the implementation of boundary conditions in the lattice Boltzmann method has been developed. The main goal of this framework is to implement a boundary condition for one side/corner of a domain which can then be applied to any other side/corner. Since implementing a boundary condition depends on the orientation of the boundary nodes, using this framework makes it easier to implement any set of boundary conditions. The framework is developed in C++ using object-oriented programming. It is then verified by implementing some common boundary conditions in LBM. This framework can be used for any kind of data structure of lattice Boltzmann code.

A case study on the migration of nuclides in the coupled vadose zone-groundwater system at a proposed spent fuel reprocessing site in Gansu, China

Nuclear energy represents a relatively novel clean energy source, and potential nuclear contamination must be considered while developing the nuclear industry. In the context of a hypothetical nuclear accident, a nuclide migration model was established based on regional geological and hydrological data to investigate the migration characteristics of the radionuclides 3H, 99Tc(VII), 137Cs, and 79Se(IV) in a proposed spent fuel reprocessing site (SFRS) in Gansu, China. The temporal and spatial evolution of these four nuclides as well as their effects on the environment and local population were discussed in detail. The model of nuclide migration in the vadose zone was coupled to a large-scale regional nuclide migration model of groundwater in the form of boundary conditions. The results under the least favorable conditions showed that nuclides with high adsorption rates (137Cs and 79Se(IV)) remained in the vadose zone and decayed out completely. In contrast, nuclides with low adsorption rates (3H and 99Tc(VII)) penetrated the vadose zone and polluted the groundwater. The vadose zone delayed the arrival of 3H and 99Tc(VII) to the phreatic surface and reduced their concentrations. It was predicted that the concentration of 3H in the groundwater would fall below the threshold after 40 years, and the distance from the leakage area to the leading edge of the contamination plume would be 374 m. 99Tc(VII), which has a slightly higher adsorption rate and a lower decay rate than 3H, spread 833 m after 3000 years when its concentration was below the threshold. 3H and 99Tc(VII) did not migrate to downstream farms to threaten the safety of residential drinking water. This research presents a safety assessment of a nuclear area in Gansu, China, and provides a significant reference for the future construction of nuclear facilities and nuclear contamination control.

Bilateral XVA pricing under stochastic default intensity: PDE modelling and computation

Adjusting derivative prices to take into account default risk has attracted the attention of several researchers and practitioners, especially after the 2007-2008 financial crisis. We derive a novel partial differential equation (PDE) model for derivative pricing including the adjustment for default risk, assuming that the default risk of one of the counterparties (the buyer) follows a Cox-Ingersoll-Ross (CIR) process, while the other party has constant default risk. The time-dependent PDE derived is of Black-Scholes type and involves two “space” variables, namely the asset price and the buyer default intensity, as well as a nonlinear source term. We formulate boundary conditions appropriate for the default intensity variable. The numerical solution of the PDE is based on standard finite differences, and a penalty-like iteration for handling the nonlinearity. We also develop and analyze a novel asymptotic approximation formula for the adjusted price of derivatives, resulting in a very efficient, accurate, and convenient for practitioners formula. We present numerical results that indicate stable second order convergence for the 2D PDE solution in terms of the discretization size. We compare the effectiveness of the 2D PDE and asymptotic approximations. We study the effect of various numerical and market parameters to the values of the adjusted prices and to the accuracy of the computed solutions.

To E or Not to E : Numerical Nuances of Global Coronal Models

In recent years, global coronal models have experienced an ongoing increase in popularity as tools for forecasting solar weather. Within the domain of up to 21.5 R ⊙, magnetohydrodynamics (MHD) is used to resolve the coronal structure using magnetograms as inputs at the solar surface. Ideally, these computations would be repeated with every update of the solar magnetogram so that they could be used in the ESA Modeling and Data Analysis Working Group magnetic connectivity tool (http://connect-tool.irap.omp.eu/). Thus, it is crucial that these results are both accurate and efficient. While much work has been published showing the results of these models in comparison with observations, not much of it discusses the intricate numerical adjustments required to achieve these results. These range from the details of boundary condition formulations to adjustments as large as enforcing parallelism between the magnetic field and velocity. By omitting the electric field in ideal MHD, the description of the physics can be insufficient and may lead to excessive diffusion and incorrect profiles. We formulate inner boundary conditions that, along with other techniques, reduce artificial electric field generation. Moreover, we investigate how different outer boundary condition formulations and grid design affect the results and convergence, with special focus on the density and radial component of the B -field. The significant improvement in accuracy of real magnetic map–driven simulations is illustrated for an example of the 2008 eclipse.

A local meshless method for transient nonlinear problems: Preliminary investigation and application to phase-field models

Transient nonlinear problems play an important role in many engineering problems. Phase-field equations, including the well-known Allen-Cahn and Cahn-Hilliard equations, fall in this category, and have applications in cutting-edge technologies such as modeling the diffusion of lithium (Li) ions in two-phase electrode particles of Li-ion batteries. In this paper, a local meshless method for solving this category of partial differential equations (PDEs) is proposed. The Newton-Kantorovich scheme is employed to transform the nonlinear PDEs to an iterative series of linear ones which can be solved with the proposed method. The accuracy and performance of the method are examined in various linear and nonlinear problems, such as Laplace equation, three dimensional elasticity as well as some abstract mathematical equations with linear or nonlinear boundary conditions. The main focus of the work is on applying the proposed method in solution of the phase-field equations, including the Allen-Cahn and Cahn-Hilliard equations. In addition to homogeneous Neumann boundary condition which has been widely examined in the literature, we also employ a practical nonlinear, inhomogeneous Neumann boundary condition formulation specialized for modeling the diffusion of lithium ions in electrode particles of Li-ion batteries. The generalized- α method is used for time integration of diffusion-type equations to overcome the intrinsic stiffness of the phase-field equations. It is shown that the method is capable of capturing the main features of the phase-field models i.e. phase separation, coarsening and energy decay in closed systems.

A time domain boundary element method for acoustic scattering by lined surfaces in a subsonic uniform mean flow

A time domain boundary integral equation with Burton-Miller reformulation is presented for acoustic scattering by surfaces with liners in a uniform mean flow. The Ingard-Myers impedance boundary condition is implemented using a broadband multipole impedance model which is in turn converted into time domain differential equations to augment the boundary integral equation. The coupled integral-differential equations are solved numerically by a March-On-in-Time (MOT) scheme. While the Ingard-Myers condition is known to support Kelvin-Helmholtz instability due to its use of a vortex sheet interface between the flow and the lined surface, it is found that by neglecting a second-order derivative term in the current time domain impedance boundary condition formulation, the instability can be effectively avoided in computation. Proposed formulation and implementation are validated with numerical examples. Moreover, a minimization procedure for finding the poles and coefficients of the broadband multiple impedance expansion is formulated by which, unlike the commonly used vector-fitting method, passivity of the model is ensured. Numerical tests show the proposed minimization approach is effective for modeling liners that are commonly used in aeroacoustic applications.

An Empirical Model of the Equatorial Electron Pitch Angle Distributions in Earth's Outer Radiation Belt

AbstractIn this study, we present an empirical model of the equatorial electron pitch angle distributions (PADs) in the outer radiation belt based on the full data set collected by the Magnetic Electron Ion Spectrometer (MagEIS) instrument onboard the Van Allen Probes in 2012–2019. The PADs are fitted with a combination of the first, third and fifth sine harmonics. The resulting equation resolves all PAD types found in the outer radiation belt (pancake, flat‐top, butterfly and cap PADs) and can be analytically integrated to derive omnidirectional flux. We introduce a two‐step modeling procedure that for the first time ensures a continuous dependence on L, magnetic local time and activity, parametrized by the solar wind dynamic pressure. We propose two methods to reconstruct equatorial electron flux using the model. The first approach requires two uni‐directional flux observations and is applicable to low‐PA data. The second method can be used to reconstruct the full equatorial PADs from a single uni‐ or omnidirectional measurement at off‐equatorial latitudes. The model can be used for converting the long‐term data sets of electron fluxes to phase space density in terms of adiabatic invariants, for physics‐based modeling in the form of boundary conditions, and for data assimilation purposes.

Effective property calculation and its numerical implementation of spatially graded plate structures based on asymptotic homogenization

Spatially graded heterogeneous plate structures that are mapped from periodic plates based on mapping functions are attracting increasing attention due to their excellent mechanical performances and high tailorability, and homogenization approaches are a powerful tool for their efficient numerical analysis. In this work, based on asymptotic homogenization, a unified framework of the effective properties prediction and the corresponding unit cell problems is first theoretically established for spatially-varying plates with arbitrarily-shaped microstructures, which are distinct from those of periodic plates due to the Jacobian of mapping functions. Moreover, the FE formulation for efficient numerical implementation are also proposed for not only solid elements, but also shell and beam elements, where numerical treatment details in stiffness matrix formulation and periodic boundary conditions are elaborated, so that highly efficient effective stiffness computation can be achieved for thin-walled heterogeneous plate structures. At last, unit cells under different Jacobian matrices, which are related to different kinds of geometric deformation patterns, and two spatially graded plates are analyzed to corroborate the correctness and efficiency of the proposed method.