This paper aims to investigate the multi-soliton solutions of the coupled Lakshmanan–Porsezian–Daniel equations with variable coefficients under nonzero boundary conditions.These equations are utilized to model the ...This paper aims to investigate the multi-soliton solutions of the coupled Lakshmanan–Porsezian–Daniel equations with variable coefficients under nonzero boundary conditions.These equations are utilized to model the phenomenon of nonlinear waves propagating simultaneously in non-uniform optical fibers.By analyzing the Lax pair and the Riemann–Hilbert problem,we aim to provide a comprehensive understanding of the dynamics and interactions of solitons of this system.Furthermore,we study the impacts of group velocity dispersion or the fourth-order dispersion on soliton behaviors.Through appropriate parameter selections,we observe various nonlinear phenomena,including the disappearance of solitons after interaction and their transformation into breather-like solitons,as well as the propagation of breathers with variable periodicity and interactions between solitons with variable periodicities.展开更多
The Riemann–Hilbert approach is demonstrated to investigate the defocusing Lakshmanan–Porsezian–Daniel equation under fully asymmetric nonzero boundary conditions.In contrast to the symmetry case,this paper focuses...The Riemann–Hilbert approach is demonstrated to investigate the defocusing Lakshmanan–Porsezian–Daniel equation under fully asymmetric nonzero boundary conditions.In contrast to the symmetry case,this paper focuses on the branch points related to the scattering problem rather than using the Riemann surfaces.For the direct problem,we analyze the Jost solution of lax pairs and some properties of scattering matrix,including two kinds of symmetries.The inverse problem at branch points can be presented,corresponding to the associated Riemann–Hilbert.Moreover,we investigate the time evolution problem and estimate the value of solving the solutions by Jost function.For the inverse problem,we construct it as a Riemann–Hilbert problem and formulate the reconstruction formula for the defocusing Lakshmanan–Porsezian–Daniel equation.The solutions of the Riemann–Hilbert problem can be constructed by estimating the solutions.Finally,we work out the solutions under fully asymmetric nonzero boundary conditions precisely via utilizing the Sokhotski–Plemelj formula and the square of the negative column transformation with the assistance of Riemann surfaces.These results are valuable for understanding physical phenomena and developing further applications of optical problems.展开更多
The determination of the ultimate load-bearing capacity of structures made of elastoplastic heterogeneous materials under varying loads is of great importance for engineering analysis and design. Therefore, it is nece...The determination of the ultimate load-bearing capacity of structures made of elastoplastic heterogeneous materials under varying loads is of great importance for engineering analysis and design. Therefore, it is necessary to accurately predict the shakedown domains of these materials. The static shakedown theorem, also known as Melan's theorem, is a fundamental method used to predict the shakedown domains of structures and materials. Within this method, a key aspect lies in the construction and application of an appropriate self-equilibrium stress field(SSF). In the structural shakedown analysis, the SSF is typically constructed by governing equations that satisfy no external force(NEF) boundary conditions. However, we discover that directly applying these governing equations is not suitable for the shakedown analysis of heterogeneous materials. Researchers must consider the requirements imposed by the Hill-Mandel condition for boundary conditions and the physical significance of representative volume elements(RVEs). This paper addresses this issue and demonstrates that the sizes of SSFs vary under different boundary conditions, such as uniform displacement boundary conditions(DBCs), uniform traction boundary conditions(TBCs), and periodic boundary conditions(PBCs). As a result, significant discrepancies arise in the predicted shakedown domain sizes of heterogeneous materials. Built on the demonstrated relationship between SSFs under different boundary conditions, this study explores the conservative relationships among different shakedown domains, and provides proof of the relationship between the elastic limit(EL) factors and the shakedown loading factors under the loading domain of two load vertices. By utilizing numerical examples, we highlight the conservatism present in certain results reported in the existing literature. Among the investigated boundary conditions, the obtained shakedown domain is the most conservative under TBCs.Conversely, utilizing PBCs to construct an SSF for the shakedown analysis leads to less conservative lower bounds, indicating that PBCs should be employed as the preferred boundary conditions for the shakedown analysis of heterogeneous materials.展开更多
Laminated composites are widely used in many engineering industries such as aircraft, spacecraft, boat hulls, racing car bodies, and storage tanks. We analyze the 3D deformations of a multilayered, linear elastic, ani...Laminated composites are widely used in many engineering industries such as aircraft, spacecraft, boat hulls, racing car bodies, and storage tanks. We analyze the 3D deformations of a multilayered, linear elastic, anisotropic rectangular plate subjected to arbitrary boundary conditions on one edge and simply supported on other edge. The rectangular laminate consists of anisotropic and homogeneous laminae of arbitrary thicknesses. This study presents the elastic analysis of laminated composite plates subjected to sinusoidal mechanical loading under arbitrary boundary conditions. Least square finite element solutions for displacements and stresses are investigated using a mathematical model, called a state-space model, which allows us to simultaneously solve for these field variables in the composite structure’s domain and ensure that continuity conditions are satisfied at layer interfaces. The governing equations are derived from this model using a numerical technique called the least-squares finite element method (LSFEM). These LSFEMs seek to minimize the squares of the governing equations and the associated side conditions residuals over the computational domain. The model is comprised of layerwise variables such as displacements, out-of-plane stresses, and in- plane strains, treated as independent variables. Numerical results are presented to demonstrate the response of the laminated composite plates under various arbitrary boundary conditions using LSFEM and compared with the 3D elasticity solution available in the literature.展开更多
Boundary conditions for momentum and vorticity have been precisely derived, paying attention to the physical meaning of each mathematical expression of terms rigorously obtained from the basic equations: Navier-Stokes...Boundary conditions for momentum and vorticity have been precisely derived, paying attention to the physical meaning of each mathematical expression of terms rigorously obtained from the basic equations: Navier-Stokes equation and the equation of vorticity transport. It has been shown first that a contribution of fluid molecules crossing over a conceptual surface moving with fluid velocity due to their fluctuating motion is essentially important to understanding transport phenomena of momentum and vorticity. A notion of surface layers, which are thin layers at both sides of an interface, has been introduced next to elucidate the transporting mechanism of momentum and vorticity from one phase to the other at an interface through which no fluid molecules are crossing over. A fact that a size of δV, in which reliable values of density, momentum, and velocity of fluid are respectively defined as a volume-averaged mass of fluid molecules, a volume-averaged momentum of fluid molecules and a mass-averaged velocity of fluid molecules, is not infinitesimal but finite has been one of the key factors leading to the boundary conditions for vorticity at an interface between two fluids. The most distinguished characteristics of the boundary conditions derived here are the zero-value conditions for a normal component of momentum flux and tangential components of vorticity flux, at an interface.展开更多
The objective of this paper is to present a new method for designing absorbing or non-reflective boundary conditions (ABC) or (NRBC), illustrated by the case study of the modelling of a solid body in water, specifical...The objective of this paper is to present a new method for designing absorbing or non-reflective boundary conditions (ABC) or (NRBC), illustrated by the case study of the modelling of a solid body in water, specifically the capillary gravity waves generated by its motion at the surface. The study analyses the flow of an inviscid, barotropic, and compressible fluid around the stationary solid body. The dynamic behaviour of the fluid is analysed using a two-dimensional coupled Neumann-Kelvin model extended with capillarity and inertia terms. For computational purposes, it is necessary to truncate the unbounded spatial domain with artificial boundaries and then introduce appropriate absorbing boundary conditions. The propagation of short wavelength waves in a convective fluid medium with significant differences in properties between the interior and the surface of the fluid presents a number of difficulties in the design of these conditions. The results are illustrated numerically and commented upon.展开更多
The seismic response analysis of a tailing dam is studied using a fully coupled effective stress approach in conjunction with an advanced multi yield surface plastic constitutive model for tailing material.Strain cont...The seismic response analysis of a tailing dam is studied using a fully coupled effective stress approach in conjunction with an advanced multi yield surface plastic constitutive model for tailing material.Strain controlled static and cyclic triaxial tests were carried out to obtain the constitutive model for the tailing material.The tailing materials were collected from the Rampura Agucha tailing dam(Rajasthan State,India).A 2D nonlinear finite element(FE)model was then developed using different boundary conditions from the tailing embankment constructed using the downstream and upstream method of rising using OpenSees software.In first case,the model boundary was fixed in both the X and Y directions,and in the second case,viscous dashpots were introduced for both side and horizontal boundaries.The model was validated with experimental results on tailing material.Analyses were carried out considering five different earthquake motions,which were applied at the base.Comparisons of the different boundary conditions in terms of displacement flow vectors,pore pressure and stress-strain curves during shaking are presented.From the analysis,it was observed that the viscous boundary condition replicates the actual field conditions more accurately than the fixed boundary condition.In addition,it was found that the tailing embankment constructed by the downstream and upstream method of rising is not susceptible to liquefaction and lateral spreading for earthquake motions,even for a magnitude>5.5.展开更多
Terahertz(THz) radiation can be generated due to the instability of THz plasma waves in field-effect transistors(FETs). In this work, we discuss the instability of THz plasma waves in the channel of FETs with spin and...Terahertz(THz) radiation can be generated due to the instability of THz plasma waves in field-effect transistors(FETs). In this work, we discuss the instability of THz plasma waves in the channel of FETs with spin and quantum effects under non-ideal boundary conditions. We obtain a linear dispersion relation by using the hydrodynamic equation, Maxwell equation and spin equation. The influence of source capacitance, drain capacitance, spin effects, quantum effects and channel width on the instability of THz plasma waves under the non-ideal boundary conditions is investigated in great detail. The results of numerical simulation show that the THz plasma wave is unstable when the drain capacitance is smaller than the source capacitance;the oscillation frequency with asymmetric boundary conditions is smaller than that under non-ideal boundary conditions;the instability gain of THz plasma waves becomes lower under non-ideal boundary conditions. This finding provides a new idea for finding efficient THz radiation sources and opens up a new mechanism for the development of THz technology.展开更多
Based on Kirchhoff plate theory and the Rayleigh-Ritz method,the model for free vibration of rectangular plate with rectangular cutouts under arbitrary elastic boundary conditions is established by using the improved ...Based on Kirchhoff plate theory and the Rayleigh-Ritz method,the model for free vibration of rectangular plate with rectangular cutouts under arbitrary elastic boundary conditions is established by using the improved Fourier series in combination with the independent coordinate coupling method(ICCM).The effect of the cutout is taken into account by subtracting the energies of the cutouts from the total energies of the whole plate.The vibration displacement function of the hole domain is based on the coordinate system of the hole domain in this method.From the continuity condition of the vibration displacement function at the cutout,the transition matrix between the two coordinate systems is constructed,and the mass and stiffness matrices are completely obtained.As a result,the calculation is simplified and the computational efficiency of the solution is improved.In this paper,numerical examples and modal experiments are presented to validate the effectiveness of the modeling methods,and parameters related to influencing factors of the rectangular plate are analyzed to study the vibration characteristics.展开更多
We address the issue of point value reconstructions from cell averages in the context of third-order finite volume schemes,focusing in particular on the cells close to the boundaries of the domain.In fact,most techniq...We address the issue of point value reconstructions from cell averages in the context of third-order finite volume schemes,focusing in particular on the cells close to the boundaries of the domain.In fact,most techniques in the literature rely on the creation of ghost cells outside the boundary and on some form of extrapolation from the inside that,taking into account the boundary conditions,fills the ghost cells with appropriate values,so that a standard reconstruction can be applied also in the boundary cells.In Naumann et al.(Appl.Math.Comput.325:252–270.https://doi.org/10.1016/j.amc.2017.12.041,2018),motivated by the difficulty of choosing appropriate boundary conditions at the internal nodes of a network,a different technique was explored that avoids the use of ghost cells,but instead employs for the boundary cells a different stencil,biased towards the interior of the domain.In this paper,extending that approach,which does not make use of ghost cells,we propose a more accurate reconstruction for the one-dimensional case and a two-dimensional one for Cartesian grids.In several numerical tests,we compare the novel reconstruction with the standard approach using ghost cells.展开更多
By using cone theory and the MSnch fixed theorem combined with a monotone iterative technique, we investigate the existence of positive solutions for systems of second- order nonlinear singular differential equations ...By using cone theory and the MSnch fixed theorem combined with a monotone iterative technique, we investigate the existence of positive solutions for systems of second- order nonlinear singular differential equations with integral boundary conditions on infinite interval and establish the existence theorem of positive solutions and iterative sequence for approximating the positive solutions. The results in this paper improve some known results.展开更多
In this paper, we study a boundary value problem of nonlinear fractional dif- ferential equations of order q （1 〈 q 〈 2） with non-separated integral boundary conditions. Some new existence and uniqueness results a...In this paper, we study a boundary value problem of nonlinear fractional dif- ferential equations of order q （1 〈 q 〈 2） with non-separated integral boundary conditions. Some new existence and uniqueness results are obtained by using some standard fixed point theorems and Leray-Schauder degree theory. Some illustrative examples are also presented. We extend previous results even in the integer case q = 2.展开更多
Inverse method was used in single crystal superalloy DD6 processing simulation during solidification. Numerical modeling coupled with experiments has been used to estimate the interface heat transfer coefficient （IHT...Inverse method was used in single crystal superalloy DD6 processing simulation during solidification. Numerical modeling coupled with experiments has been used to estimate the interface heat transfer coefficient （IHTC） between the surface of slab casting and inner mold. Calculated temperature dependent values of IHTC were obtained from a numerical solution. The calculated temperatures agreed well with the measurement of cooling profile.展开更多
We developed a Global Ocean Circulation and Tide Model （GOCTM） with coarse grids in the open deep ocean degrading ‘smoothly’ into the highly resolved China Seas （CS） of refined grids to study the tides and circu...We developed a Global Ocean Circulation and Tide Model （GOCTM） with coarse grids in the open deep ocean degrading ‘smoothly’ into the highly resolved China Seas （CS） of refined grids to study the tides and circulation there.GOCTM is based on the framework of the Finite Volume approach for better mass conservation through improved transports across the discrete individual control volume.It also takes a full advantage of the geometric flexibility of unstructured mesh using a realistic global topography including the Arctic Ocean.The CS are given a special focus by refining the unstructured grids,but they are embedded into global domain naturally.Furthermore,GOCTM not only successfully avoids the treatment of the open boundaries,but also optimizes the trade-off between computational cost and model accuracy.Meanwhile,GOCTM is driven by the astronomical tide-generating potential and the secondary tide-generating potential directly,together with the wind stress and heat flux.GOCTM succeeds in reproducing the global eight principal tidal harmonic constants.Particularly,the simulated tidal results in the CS are improved compared to some other regional models with the discrepancy of 3.9 cm for M 2 tide.This idea of GOCTM can also be referred for other regional ocean study.展开更多
The numerical dynamic model of Chinese mainland lithosphere's stress and strain field was constructed withelasto-viscous creep constitive relation. The most recent data of the stress field of Chinese mainland and ...The numerical dynamic model of Chinese mainland lithosphere's stress and strain field was constructed withelasto-viscous creep constitive relation. The most recent data of the stress field of Chinese mainland and thehorizontal movement velocity of the crust blocks were used as constraint conditions. The values of the boundalsforce were computed by trial-and-error method. The effect of the Qinghai-Xizang Plateau's excess potential energyto the movement of Chinese mainland was studied also in this model. The results of the numerical computing showthat, recent rapid raising of the Qinghai-Xizang Plateau and the generation of normal faults in the southem part ofthe plateau resulted from the convergence of Indian Plate to the Eurasian Plate, and also from the rapid convectivethinning of the lower lithosphere. Horizontal extension of the Qinghai-Xizang Plateau is the main dynamic factorto torm the present tectonic framework of the Chinese mainland. The compressive loads on the eastern boundaryof the model were mainly applied by the subduction of the Pacific Plate. The compression from the PhilippinePlate maybe slight.展开更多
A time-dependent finite element method (FEM) is developed to analyze the transient hydroelastic responses of very large floating structures (VLFS) subjected to dynamic loads. The hydrodynamic problem is formulated bas...A time-dependent finite element method (FEM) is developed to analyze the transient hydroelastic responses of very large floating structures (VLFS) subjected to dynamic loads. The hydrodynamic problem is formulated based on the linear theory of fluid and the structural response is analyzed based on the thin plate theory. The FEM truncates the unbounded fluid domain by introducing an artificial boundary surface, thus defining a finite computational domain. At this boundary surface an impedance boundary conditions are applied so that no wave reflections occur. In the proposed scheme, all of the procedures are processed directly in time domain, which is efficient for nonlinear analyses of structure floating on unbounded fluid. Numerical results indicate acceptable accuracy of the proposed method.展开更多
A variable coefficient viscoelastic equation with a time-varying delay in the boundary feedback and acoustic boundary conditions and nonlinear source term is considered.Under suitable assumptions, general decay result...A variable coefficient viscoelastic equation with a time-varying delay in the boundary feedback and acoustic boundary conditions and nonlinear source term is considered.Under suitable assumptions, general decay results of the energy are established via suitable Lyapunov functionals and some properties of the convex functions. Our result is obtained without imposing any restrictive growth assumption on the damping term and the elements of the matrix A and the kernel function g.展开更多
This article concentrates on the steady magnetohydrodynamic （MHD） flow of viscous nanofluid. The flow is caused by a permeable exponentially stretching surface. An incompressible fluid fills the porous space. A comp...This article concentrates on the steady magnetohydrodynamic （MHD） flow of viscous nanofluid. The flow is caused by a permeable exponentially stretching surface. An incompressible fluid fills the porous space. A comparative study is made for the nanoparticles namely Copper （Cu）, Silver （Ag）, Alumina （A1203） and Titanium Oxide （TiO2）. Water is treated as a base fluid. Convective type boundary conditions are employed in modeling the heat transfer process. The non-linear partial differential equations governing the flow are reduced to an ordinary differential equation by similarity transformations. The obtained equations are then solved for the development of series solutions. Convergence of the obtained series solutions is explicitly discussed. The effects of different parameters on the velocity and temperature profiles are shown and analyzed through graphs.展开更多
Recent years have witnessed growing interests in solving partial differential equations by deep neural networks,especially in the high-dimensional case.Unlike classical numerical methods,such as finite difference meth...Recent years have witnessed growing interests in solving partial differential equations by deep neural networks,especially in the high-dimensional case.Unlike classical numerical methods,such as finite difference method and finite element method,the enforcement of boundary conditions in deep neural networks is highly nontrivial.One general strategy is to use the penalty method.In the work,we conduct a comparison study for elliptic problems with four different boundary conditions,i.e.,Dirichlet,Neumann,Robin,and periodic boundary conditions,using two representative methods:deep Galerkin method and deep Ritz method.In the former,the PDE residual is minimized in the least-squares sense while the corresponding variational problem is minimized in the latter.Therefore,it is reasonably expected that deep Galerkin method works better for smooth solutions while deep Ritz method works better for low-regularity solutions.However,by a number of examples,we observe that deep Ritz method can outperform deep Galerkin method with a clear dependence of dimensionality even for smooth solutions and deep Galerkin method can also outperform deep Ritz method for low-regularity solutions.Besides,in some cases,when the boundary condition can be implemented in an exact manner,we find that such a strategy not only provides a better approximate solution but also facilitates the training process.展开更多
Anti-ram bollards used in perimeter protection are tested to meet performance requirements of established standards such as the US Department of State Specification SD-STD-02.01. Under these standards, tests are condu...Anti-ram bollards used in perimeter protection are tested to meet performance requirements of established standards such as the US Department of State Specification SD-STD-02.01. Under these standards, tests are conducted in prescribed conditions that should be representative of the service installation. In actual project, conditions encountered on site may vary from the test environment and it would be expensive and time consuming to validate each deviation with a physical test. High-fidelity physics-based (HFPB) finite element modeling can provide precise simulations of the behavior of anti-ram bollards. This paper presents the use of HFPB finite element modeling, using LS-DYNA, in an actual project to evaluate the performance of an anti-ram bollard design subjected to various boundary conditions representing the physical conditions encountered on site. The study shows that boundary conditions can have a significant influence on the performance of the anti-ram bollards. This suggests that anti-ram bollards must be designed and engineered according to actual conditions that are found on site. It also shows that HFPB modeling can be an effective tool that supplements physical testing of anti-ram bollards.展开更多
基金supported by the Natural Science Foundation of Hebei Province,China (Grant No.A2021502004)the Fundamental Research Funds for the Central Universities (Grant No.2024MS126).
文摘This paper aims to investigate the multi-soliton solutions of the coupled Lakshmanan–Porsezian–Daniel equations with variable coefficients under nonzero boundary conditions.These equations are utilized to model the phenomenon of nonlinear waves propagating simultaneously in non-uniform optical fibers.By analyzing the Lax pair and the Riemann–Hilbert problem,we aim to provide a comprehensive understanding of the dynamics and interactions of solitons of this system.Furthermore,we study the impacts of group velocity dispersion or the fourth-order dispersion on soliton behaviors.Through appropriate parameter selections,we observe various nonlinear phenomena,including the disappearance of solitons after interaction and their transformation into breather-like solitons,as well as the propagation of breathers with variable periodicity and interactions between solitons with variable periodicities.
基金the Fundamental Research Funds for the Central Universities(Grant No.2024MS126).
文摘The Riemann–Hilbert approach is demonstrated to investigate the defocusing Lakshmanan–Porsezian–Daniel equation under fully asymmetric nonzero boundary conditions.In contrast to the symmetry case,this paper focuses on the branch points related to the scattering problem rather than using the Riemann surfaces.For the direct problem,we analyze the Jost solution of lax pairs and some properties of scattering matrix,including two kinds of symmetries.The inverse problem at branch points can be presented,corresponding to the associated Riemann–Hilbert.Moreover,we investigate the time evolution problem and estimate the value of solving the solutions by Jost function.For the inverse problem,we construct it as a Riemann–Hilbert problem and formulate the reconstruction formula for the defocusing Lakshmanan–Porsezian–Daniel equation.The solutions of the Riemann–Hilbert problem can be constructed by estimating the solutions.Finally,we work out the solutions under fully asymmetric nonzero boundary conditions precisely via utilizing the Sokhotski–Plemelj formula and the square of the negative column transformation with the assistance of Riemann surfaces.These results are valuable for understanding physical phenomena and developing further applications of optical problems.
基金Project supported by the National Natural Science Foundation of China (Nos. 52075070 and12302254)the Dalian City Supports Innovation and Entrepreneurship Projects for High-Level Talents (No. 2021RD16)the Liaoning Revitalization Talents Program (No. XLYC2002108)。
文摘The determination of the ultimate load-bearing capacity of structures made of elastoplastic heterogeneous materials under varying loads is of great importance for engineering analysis and design. Therefore, it is necessary to accurately predict the shakedown domains of these materials. The static shakedown theorem, also known as Melan's theorem, is a fundamental method used to predict the shakedown domains of structures and materials. Within this method, a key aspect lies in the construction and application of an appropriate self-equilibrium stress field(SSF). In the structural shakedown analysis, the SSF is typically constructed by governing equations that satisfy no external force(NEF) boundary conditions. However, we discover that directly applying these governing equations is not suitable for the shakedown analysis of heterogeneous materials. Researchers must consider the requirements imposed by the Hill-Mandel condition for boundary conditions and the physical significance of representative volume elements(RVEs). This paper addresses this issue and demonstrates that the sizes of SSFs vary under different boundary conditions, such as uniform displacement boundary conditions(DBCs), uniform traction boundary conditions(TBCs), and periodic boundary conditions(PBCs). As a result, significant discrepancies arise in the predicted shakedown domain sizes of heterogeneous materials. Built on the demonstrated relationship between SSFs under different boundary conditions, this study explores the conservative relationships among different shakedown domains, and provides proof of the relationship between the elastic limit(EL) factors and the shakedown loading factors under the loading domain of two load vertices. By utilizing numerical examples, we highlight the conservatism present in certain results reported in the existing literature. Among the investigated boundary conditions, the obtained shakedown domain is the most conservative under TBCs.Conversely, utilizing PBCs to construct an SSF for the shakedown analysis leads to less conservative lower bounds, indicating that PBCs should be employed as the preferred boundary conditions for the shakedown analysis of heterogeneous materials.
文摘Laminated composites are widely used in many engineering industries such as aircraft, spacecraft, boat hulls, racing car bodies, and storage tanks. We analyze the 3D deformations of a multilayered, linear elastic, anisotropic rectangular plate subjected to arbitrary boundary conditions on one edge and simply supported on other edge. The rectangular laminate consists of anisotropic and homogeneous laminae of arbitrary thicknesses. This study presents the elastic analysis of laminated composite plates subjected to sinusoidal mechanical loading under arbitrary boundary conditions. Least square finite element solutions for displacements and stresses are investigated using a mathematical model, called a state-space model, which allows us to simultaneously solve for these field variables in the composite structure’s domain and ensure that continuity conditions are satisfied at layer interfaces. The governing equations are derived from this model using a numerical technique called the least-squares finite element method (LSFEM). These LSFEMs seek to minimize the squares of the governing equations and the associated side conditions residuals over the computational domain. The model is comprised of layerwise variables such as displacements, out-of-plane stresses, and in- plane strains, treated as independent variables. Numerical results are presented to demonstrate the response of the laminated composite plates under various arbitrary boundary conditions using LSFEM and compared with the 3D elasticity solution available in the literature.
文摘Boundary conditions for momentum and vorticity have been precisely derived, paying attention to the physical meaning of each mathematical expression of terms rigorously obtained from the basic equations: Navier-Stokes equation and the equation of vorticity transport. It has been shown first that a contribution of fluid molecules crossing over a conceptual surface moving with fluid velocity due to their fluctuating motion is essentially important to understanding transport phenomena of momentum and vorticity. A notion of surface layers, which are thin layers at both sides of an interface, has been introduced next to elucidate the transporting mechanism of momentum and vorticity from one phase to the other at an interface through which no fluid molecules are crossing over. A fact that a size of δV, in which reliable values of density, momentum, and velocity of fluid are respectively defined as a volume-averaged mass of fluid molecules, a volume-averaged momentum of fluid molecules and a mass-averaged velocity of fluid molecules, is not infinitesimal but finite has been one of the key factors leading to the boundary conditions for vorticity at an interface between two fluids. The most distinguished characteristics of the boundary conditions derived here are the zero-value conditions for a normal component of momentum flux and tangential components of vorticity flux, at an interface.
文摘The objective of this paper is to present a new method for designing absorbing or non-reflective boundary conditions (ABC) or (NRBC), illustrated by the case study of the modelling of a solid body in water, specifically the capillary gravity waves generated by its motion at the surface. The study analyses the flow of an inviscid, barotropic, and compressible fluid around the stationary solid body. The dynamic behaviour of the fluid is analysed using a two-dimensional coupled Neumann-Kelvin model extended with capillarity and inertia terms. For computational purposes, it is necessary to truncate the unbounded spatial domain with artificial boundaries and then introduce appropriate absorbing boundary conditions. The propagation of short wavelength waves in a convective fluid medium with significant differences in properties between the interior and the surface of the fluid presents a number of difficulties in the design of these conditions. The results are illustrated numerically and commented upon.
文摘The seismic response analysis of a tailing dam is studied using a fully coupled effective stress approach in conjunction with an advanced multi yield surface plastic constitutive model for tailing material.Strain controlled static and cyclic triaxial tests were carried out to obtain the constitutive model for the tailing material.The tailing materials were collected from the Rampura Agucha tailing dam(Rajasthan State,India).A 2D nonlinear finite element(FE)model was then developed using different boundary conditions from the tailing embankment constructed using the downstream and upstream method of rising using OpenSees software.In first case,the model boundary was fixed in both the X and Y directions,and in the second case,viscous dashpots were introduced for both side and horizontal boundaries.The model was validated with experimental results on tailing material.Analyses were carried out considering five different earthquake motions,which were applied at the base.Comparisons of the different boundary conditions in terms of displacement flow vectors,pore pressure and stress-strain curves during shaking are presented.From the analysis,it was observed that the viscous boundary condition replicates the actual field conditions more accurately than the fixed boundary condition.In addition,it was found that the tailing embankment constructed by the downstream and upstream method of rising is not susceptible to liquefaction and lateral spreading for earthquake motions,even for a magnitude>5.5.
基金funded by National Natural Science Foundation of China (No. 12065015)the Hongliu First-level Discipline Construction Project of Lanzhou University of Technology。
文摘Terahertz(THz) radiation can be generated due to the instability of THz plasma waves in field-effect transistors(FETs). In this work, we discuss the instability of THz plasma waves in the channel of FETs with spin and quantum effects under non-ideal boundary conditions. We obtain a linear dispersion relation by using the hydrodynamic equation, Maxwell equation and spin equation. The influence of source capacitance, drain capacitance, spin effects, quantum effects and channel width on the instability of THz plasma waves under the non-ideal boundary conditions is investigated in great detail. The results of numerical simulation show that the THz plasma wave is unstable when the drain capacitance is smaller than the source capacitance;the oscillation frequency with asymmetric boundary conditions is smaller than that under non-ideal boundary conditions;the instability gain of THz plasma waves becomes lower under non-ideal boundary conditions. This finding provides a new idea for finding efficient THz radiation sources and opens up a new mechanism for the development of THz technology.
基金support of this work by the National Natural Science Foundation of China(No.51405096)the Fundamental Research Funds for the Central Universities(HEUCF210710).
文摘Based on Kirchhoff plate theory and the Rayleigh-Ritz method,the model for free vibration of rectangular plate with rectangular cutouts under arbitrary elastic boundary conditions is established by using the improved Fourier series in combination with the independent coordinate coupling method(ICCM).The effect of the cutout is taken into account by subtracting the energies of the cutouts from the total energies of the whole plate.The vibration displacement function of the hole domain is based on the coordinate system of the hole domain in this method.From the continuity condition of the vibration displacement function at the cutout,the transition matrix between the two coordinate systems is constructed,and the mass and stiffness matrices are completely obtained.As a result,the calculation is simplified and the computational efficiency of the solution is improved.In this paper,numerical examples and modal experiments are presented to validate the effectiveness of the modeling methods,and parameters related to influencing factors of the rectangular plate are analyzed to study the vibration characteristics.
基金MIUR-PRIN project 2017KKJP4X“Innovative numerical methods for evolutionary partial differential equations and applications”.Gabriella Puppo acknowledges also the support of 2019 Ateneo Sapienza research project no.RM11916B51CD40E1.
文摘We address the issue of point value reconstructions from cell averages in the context of third-order finite volume schemes,focusing in particular on the cells close to the boundaries of the domain.In fact,most techniques in the literature rely on the creation of ghost cells outside the boundary and on some form of extrapolation from the inside that,taking into account the boundary conditions,fills the ghost cells with appropriate values,so that a standard reconstruction can be applied also in the boundary cells.In Naumann et al.(Appl.Math.Comput.325:252–270.https://doi.org/10.1016/j.amc.2017.12.041,2018),motivated by the difficulty of choosing appropriate boundary conditions at the internal nodes of a network,a different technique was explored that avoids the use of ghost cells,but instead employs for the boundary cells a different stencil,biased towards the interior of the domain.In this paper,extending that approach,which does not make use of ghost cells,we propose a more accurate reconstruction for the one-dimensional case and a two-dimensional one for Cartesian grids.In several numerical tests,we compare the novel reconstruction with the standard approach using ghost cells.
基金SuppoSed by the NSF of Anhui Provincial Education Depaxtment(KJ2012A265,KJ2012B187)
文摘By using cone theory and the MSnch fixed theorem combined with a monotone iterative technique, we investigate the existence of positive solutions for systems of second- order nonlinear singular differential equations with integral boundary conditions on infinite interval and establish the existence theorem of positive solutions and iterative sequence for approximating the positive solutions. The results in this paper improve some known results.
文摘In this paper, we study a boundary value problem of nonlinear fractional dif- ferential equations of order q （1 〈 q 〈 2） with non-separated integral boundary conditions. Some new existence and uniqueness results are obtained by using some standard fixed point theorems and Leray-Schauder degree theory. Some illustrative examples are also presented. We extend previous results even in the integer case q = 2.
基金supported by National Basic Research Program of China(No.2005CB724105)National Natural Science Foundation of China (No.10477010)National High Technical Research and Development Program of China(No.2007AA04Z141)
文摘Inverse method was used in single crystal superalloy DD6 processing simulation during solidification. Numerical modeling coupled with experiments has been used to estimate the interface heat transfer coefficient （IHTC） between the surface of slab casting and inner mold. Calculated temperature dependent values of IHTC were obtained from a numerical solution. The calculated temperatures agreed well with the measurement of cooling profile.
基金The Hi-tech Research and Development Program (863) of China under contract No.2007AA09Z117the National Key Technology R&D Program under contract No.2011BAC03B02+1 种基金the National Natural Science Fund of China under contract No.40976001the National Marine Renewable Energy Program under contract Nos GHME2010ZC08,No.GHME 2010ZC11 and No.GHME2010ZC01
文摘We developed a Global Ocean Circulation and Tide Model （GOCTM） with coarse grids in the open deep ocean degrading ‘smoothly’ into the highly resolved China Seas （CS） of refined grids to study the tides and circulation there.GOCTM is based on the framework of the Finite Volume approach for better mass conservation through improved transports across the discrete individual control volume.It also takes a full advantage of the geometric flexibility of unstructured mesh using a realistic global topography including the Arctic Ocean.The CS are given a special focus by refining the unstructured grids,but they are embedded into global domain naturally.Furthermore,GOCTM not only successfully avoids the treatment of the open boundaries,but also optimizes the trade-off between computational cost and model accuracy.Meanwhile,GOCTM is driven by the astronomical tide-generating potential and the secondary tide-generating potential directly,together with the wind stress and heat flux.GOCTM succeeds in reproducing the global eight principal tidal harmonic constants.Particularly,the simulated tidal results in the CS are improved compared to some other regional models with the discrepancy of 3.9 cm for M 2 tide.This idea of GOCTM can also be referred for other regional ocean study.
文摘The numerical dynamic model of Chinese mainland lithosphere's stress and strain field was constructed withelasto-viscous creep constitive relation. The most recent data of the stress field of Chinese mainland and thehorizontal movement velocity of the crust blocks were used as constraint conditions. The values of the boundalsforce were computed by trial-and-error method. The effect of the Qinghai-Xizang Plateau's excess potential energyto the movement of Chinese mainland was studied also in this model. The results of the numerical computing showthat, recent rapid raising of the Qinghai-Xizang Plateau and the generation of normal faults in the southem part ofthe plateau resulted from the convergence of Indian Plate to the Eurasian Plate, and also from the rapid convectivethinning of the lower lithosphere. Horizontal extension of the Qinghai-Xizang Plateau is the main dynamic factorto torm the present tectonic framework of the Chinese mainland. The compressive loads on the eastern boundaryof the model were mainly applied by the subduction of the Pacific Plate. The compression from the PhilippinePlate maybe slight.
文摘A time-dependent finite element method (FEM) is developed to analyze the transient hydroelastic responses of very large floating structures (VLFS) subjected to dynamic loads. The hydrodynamic problem is formulated based on the linear theory of fluid and the structural response is analyzed based on the thin plate theory. The FEM truncates the unbounded fluid domain by introducing an artificial boundary surface, thus defining a finite computational domain. At this boundary surface an impedance boundary conditions are applied so that no wave reflections occur. In the proposed scheme, all of the procedures are processed directly in time domain, which is efficient for nonlinear analyses of structure floating on unbounded fluid. Numerical results indicate acceptable accuracy of the proposed method.
文摘A variable coefficient viscoelastic equation with a time-varying delay in the boundary feedback and acoustic boundary conditions and nonlinear source term is considered.Under suitable assumptions, general decay results of the energy are established via suitable Lyapunov functionals and some properties of the convex functions. Our result is obtained without imposing any restrictive growth assumption on the damping term and the elements of the matrix A and the kernel function g.
基金supported by the Deanship of Scientific Research (DSR), King Abdulaziz University, Jeddah, Saudi Arabia
文摘This article concentrates on the steady magnetohydrodynamic （MHD） flow of viscous nanofluid. The flow is caused by a permeable exponentially stretching surface. An incompressible fluid fills the porous space. A comparative study is made for the nanoparticles namely Copper （Cu）, Silver （Ag）, Alumina （A1203） and Titanium Oxide （TiO2）. Water is treated as a base fluid. Convective type boundary conditions are employed in modeling the heat transfer process. The non-linear partial differential equations governing the flow are reduced to an ordinary differential equation by similarity transformations. The obtained equations are then solved for the development of series solutions. Convergence of the obtained series solutions is explicitly discussed. The effects of different parameters on the velocity and temperature profiles are shown and analyzed through graphs.
基金the grants NSFC 11971021National Key R&D Program of China(No.2018YF645B0204404)NSFC 11501399(R.Du)。
文摘Recent years have witnessed growing interests in solving partial differential equations by deep neural networks,especially in the high-dimensional case.Unlike classical numerical methods,such as finite difference method and finite element method,the enforcement of boundary conditions in deep neural networks is highly nontrivial.One general strategy is to use the penalty method.In the work,we conduct a comparison study for elliptic problems with four different boundary conditions,i.e.,Dirichlet,Neumann,Robin,and periodic boundary conditions,using two representative methods:deep Galerkin method and deep Ritz method.In the former,the PDE residual is minimized in the least-squares sense while the corresponding variational problem is minimized in the latter.Therefore,it is reasonably expected that deep Galerkin method works better for smooth solutions while deep Ritz method works better for low-regularity solutions.However,by a number of examples,we observe that deep Ritz method can outperform deep Galerkin method with a clear dependence of dimensionality even for smooth solutions and deep Galerkin method can also outperform deep Ritz method for low-regularity solutions.Besides,in some cases,when the boundary condition can be implemented in an exact manner,we find that such a strategy not only provides a better approximate solution but also facilitates the training process.
文摘Anti-ram bollards used in perimeter protection are tested to meet performance requirements of established standards such as the US Department of State Specification SD-STD-02.01. Under these standards, tests are conducted in prescribed conditions that should be representative of the service installation. In actual project, conditions encountered on site may vary from the test environment and it would be expensive and time consuming to validate each deviation with a physical test. High-fidelity physics-based (HFPB) finite element modeling can provide precise simulations of the behavior of anti-ram bollards. This paper presents the use of HFPB finite element modeling, using LS-DYNA, in an actual project to evaluate the performance of an anti-ram bollard design subjected to various boundary conditions representing the physical conditions encountered on site. The study shows that boundary conditions can have a significant influence on the performance of the anti-ram bollards. This suggests that anti-ram bollards must be designed and engineered according to actual conditions that are found on site. It also shows that HFPB modeling can be an effective tool that supplements physical testing of anti-ram bollards.
