A general theorem for the Stokes flow over a plane boundary with mixed stick-slip boundary conditions is established. This is done by using a representation for the velocity and pressure fields in the three-dimensiona...A general theorem for the Stokes flow over a plane boundary with mixed stick-slip boundary conditions is established. This is done by using a representation for the velocity and pressure fields in the three-dimensional Stokes flow in terms of a biharmonic function and a harmonic function. The earlier theorem for the Stokes flow due to fundamental singularities before a no-slip plane boundary is shown to be a special case of the present theorem. Furthermore, in terms of the Stokes stream function, a corollary of the theorem is also derived, providing a solution to the problem of the axisymmetric Stokes flow along a rigid plane with stick-slip boundary conditions. The formulae for the drag and torque exerted by the fluid on the boundary are established. An illustrative example is given.展开更多
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.展开更多
We consider a strongly non-linear degenerate parabolic-hyperbolic problem with p(x)-Laplacian diffusion flux function. We propose an entropy formulation and prove the existence of an entropy solution.
Let?denote a smooth,bounded domain in R^(N)(N≥2).Suppose that g is a nondecreasing C^(1)positive function and assume that b(x)is continuous and nonnegative inΩ,and that it may be singular on■Ω.In this paper,we pro...Let?denote a smooth,bounded domain in R^(N)(N≥2).Suppose that g is a nondecreasing C^(1)positive function and assume that b(x)is continuous and nonnegative inΩ,and that it may be singular on■Ω.In this paper,we provide sufficient and necessary conditions on the existence of boundary blow-up solutions to the p-Laplacian problem△_(p)u=b(x)g(u)for x∈Ω,u(x)→+∞as dist(x,■Ω)→0.The estimates of such solutions are also investigated.Moreover,when b has strong singularity,the nonexistence of boundary blow-up(radial)solutions and infinitely many radial solutions are also considered.展开更多
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.展开更多
We study the initial-boundary value problem of the Navier-Stokes equations for incompressible fluids in a general domain in R^n with compact and smooth boundary, subject to the kinematic and vorticity boundary conditi...We study the initial-boundary value problem of the Navier-Stokes equations for incompressible fluids in a general domain in R^n with compact and smooth boundary, subject to the kinematic and vorticity boundary conditions on the non-flat boundary. We observe that, under the nonhomogeneous boundary conditions, the pressure p can be still recovered by solving the Neumann problem for the Poisson equation. Then we establish the well-posedness of the unsteady Stokes equations and employ the solution to reduce our initial-boundary value problem into an initial-boundary value problem with absolute boundary conditions. Based on this, we first establish the well-posedness for an appropriate local linearized problem with the absolute boundary conditions and the initial condition (without the incompressibility condition), which establishes a velocity mapping. Then we develop apriori estimates for the velocity mapping, especially involving the Sobolev norm for the time-derivative of the mapping to deal with the complicated boundary conditions, which leads to the existence of the fixed point of the mapping and the existence of solutions to our initial-boundary value problem. Finally, we establish that, when the viscosity coefficient tends zero, the strong solutions of the initial-boundary value problem in R^n(n ≥ 3) with nonhomogeneous vorticity boundary condition converge in L^2 to the corresponding Euler equations satisfying the kinematic condition.展开更多
In this paper, first we calculate finite-difference coefficients of implicit finite- difference methods (IFDM) for the first and second-order derivatives on normal grids and first- order derivatives on staggered gri...In this paper, first we calculate finite-difference coefficients of implicit finite- difference methods (IFDM) for the first and second-order derivatives on normal grids and first- order derivatives on staggered grids and find that small coefficients of high-order IFDMs exist. Dispersion analysis demonstrates that omitting these small coefficients can retain approximately the same order accuracy but greatly reduce computational costs. Then, we introduce a mirrorimage symmetric boundary condition to improve IFDMs accuracy and stability and adopt the hybrid absorbing boundary condition (ABC) to reduce unwanted reflections from the model boundary. Last, we give elastic wave modeling examples for homogeneous and heterogeneous models to demonstrate the advantages of the proposed scheme.展开更多
Material identification is critical for understanding the relationship between mechanical properties and the associated mechanical functions.However,material identification is a challenging task,especially when the ch...Material identification is critical for understanding the relationship between mechanical properties and the associated mechanical functions.However,material identification is a challenging task,especially when the characteristic of the material is highly nonlinear in nature,as is common in biological tissue.In this work,we identify unknown material properties in continuum solid mechanics via physics-informed neural networks(PINNs).To improve the accuracy and efficiency of PINNs,we develop efficient strategies to nonuniformly sample observational data.We also investigate different approaches to enforce Dirichlet-type boundary conditions(BCs)as soft or hard constraints.Finally,we apply the proposed methods to a diverse set of time-dependent and time-independent solid mechanic examples that span linear elastic and hyperelastic material space.The estimated material parameters achieve relative errors of less than 1%.As such,this work is relevant to diverse applications,including optimizing structural integrity and developing novel materials.展开更多
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.展开更多
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.展开更多
Boundary conditions for the classical solution of the Terzaghi one-dimensional consolidation equation conflict with the equation's initial condition. As such, the classical initial-boundary value problem for the Terz...Boundary conditions for the classical solution of the Terzaghi one-dimensional consolidation equation conflict with the equation's initial condition. As such, the classical initial-boundary value problem for the Terzaghi one-dimensional consolidation equation is not well-posed. Moreover, the classical boundary conditions of the equation can only be applied to problems with either perfectly pervious or perfectly impervious boundaries. General boundary conditions are proposed to overcome these shortcomings and thus transfer the solution of the Terzaghi one-dimensional consolidation equation to a well-posed initial boundary value problem. The solution for proposed general boundary conditions is validated by comparing it to the classical solution. The actual field drainage conditions can be simulated by adjusting the values of parameters b and c given in the proposed general botmdary conditions. For relatively high coefficient of consolidation, just one term in series expansions is enough to obtain results with acceptable accuracy.展开更多
文摘A general theorem for the Stokes flow over a plane boundary with mixed stick-slip boundary conditions is established. This is done by using a representation for the velocity and pressure fields in the three-dimensional Stokes flow in terms of a biharmonic function and a harmonic function. The earlier theorem for the Stokes flow due to fundamental singularities before a no-slip plane boundary is shown to be a special case of the present theorem. Furthermore, in terms of the Stokes stream function, a corollary of the theorem is also derived, providing a solution to the problem of the axisymmetric Stokes flow along a rigid plane with stick-slip boundary conditions. The formulae for the drag and torque exerted by the fluid on the boundary are established. An illustrative example is given.
基金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.
文摘We consider a strongly non-linear degenerate parabolic-hyperbolic problem with p(x)-Laplacian diffusion flux function. We propose an entropy formulation and prove the existence of an entropy solution.
基金supported by the Beijing Natural Science Foundation(1212003)。
文摘Let?denote a smooth,bounded domain in R^(N)(N≥2).Suppose that g is a nondecreasing C^(1)positive function and assume that b(x)is continuous and nonnegative inΩ,and that it may be singular on■Ω.In this paper,we provide sufficient and necessary conditions on the existence of boundary blow-up solutions to the p-Laplacian problem△_(p)u=b(x)g(u)for x∈Ω,u(x)→+∞as dist(x,■Ω)→0.The estimates of such solutions are also investigated.Moreover,when b has strong singularity,the nonexistence of boundary blow-up(radial)solutions and infinitely many radial solutions are also considered.
文摘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.
基金supported in part by the National Science Foundation under Grants DMS-0807551, DMS-0720925, and DMS-0505473the Natural Science Foundationof China (10728101)supported in part by EPSRC grant EP/F029578/1
文摘We study the initial-boundary value problem of the Navier-Stokes equations for incompressible fluids in a general domain in R^n with compact and smooth boundary, subject to the kinematic and vorticity boundary conditions on the non-flat boundary. We observe that, under the nonhomogeneous boundary conditions, the pressure p can be still recovered by solving the Neumann problem for the Poisson equation. Then we establish the well-posedness of the unsteady Stokes equations and employ the solution to reduce our initial-boundary value problem into an initial-boundary value problem with absolute boundary conditions. Based on this, we first establish the well-posedness for an appropriate local linearized problem with the absolute boundary conditions and the initial condition (without the incompressibility condition), which establishes a velocity mapping. Then we develop apriori estimates for the velocity mapping, especially involving the Sobolev norm for the time-derivative of the mapping to deal with the complicated boundary conditions, which leads to the existence of the fixed point of the mapping and the existence of solutions to our initial-boundary value problem. Finally, we establish that, when the viscosity coefficient tends zero, the strong solutions of the initial-boundary value problem in R^n(n ≥ 3) with nonhomogeneous vorticity boundary condition converge in L^2 to the corresponding Euler equations satisfying the kinematic condition.
基金supported by the National Natural Science Foundation of China(NSFC)(Grant No. 41074100)the Program for New Century Excellent Talents in University of Ministry of Education of China(Grant No. NCET-10-0812)
文摘In this paper, first we calculate finite-difference coefficients of implicit finite- difference methods (IFDM) for the first and second-order derivatives on normal grids and first- order derivatives on staggered grids and find that small coefficients of high-order IFDMs exist. Dispersion analysis demonstrates that omitting these small coefficients can retain approximately the same order accuracy but greatly reduce computational costs. Then, we introduce a mirrorimage symmetric boundary condition to improve IFDMs accuracy and stability and adopt the hybrid absorbing boundary condition (ABC) to reduce unwanted reflections from the model boundary. Last, we give elastic wave modeling examples for homogeneous and heterogeneous models to demonstrate the advantages of the proposed scheme.
基金funded by the Cora Topolewski Cardiac Research Fund at the Children’s Hospital of Philadelphia(CHOP)the Pediatric Valve Center Frontier Program at CHOP+4 种基金the Additional Ventures Single Ventricle Research Fund Expansion Awardthe National Institutes of Health(USA)supported by the program(Nos.NHLBI T32 HL007915 and NIH R01 HL153166)supported by the program(No.NIH R01 HL153166)supported by the U.S.Department of Energy(No.DE-SC0022953)。
文摘Material identification is critical for understanding the relationship between mechanical properties and the associated mechanical functions.However,material identification is a challenging task,especially when the characteristic of the material is highly nonlinear in nature,as is common in biological tissue.In this work,we identify unknown material properties in continuum solid mechanics via physics-informed neural networks(PINNs).To improve the accuracy and efficiency of PINNs,we develop efficient strategies to nonuniformly sample observational data.We also investigate different approaches to enforce Dirichlet-type boundary conditions(BCs)as soft or hard constraints.Finally,we apply the proposed methods to a diverse set of time-dependent and time-independent solid mechanic examples that span linear elastic and hyperelastic material space.The estimated material parameters achieve relative errors of less than 1%.As such,this work is relevant to diverse applications,including optimizing structural integrity and developing novel materials.
基金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.
基金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.
基金Foundation item: Project(50608038) supported by the National Natural Science Foundation of China
文摘Boundary conditions for the classical solution of the Terzaghi one-dimensional consolidation equation conflict with the equation's initial condition. As such, the classical initial-boundary value problem for the Terzaghi one-dimensional consolidation equation is not well-posed. Moreover, the classical boundary conditions of the equation can only be applied to problems with either perfectly pervious or perfectly impervious boundaries. General boundary conditions are proposed to overcome these shortcomings and thus transfer the solution of the Terzaghi one-dimensional consolidation equation to a well-posed initial boundary value problem. The solution for proposed general boundary conditions is validated by comparing it to the classical solution. The actual field drainage conditions can be simulated by adjusting the values of parameters b and c given in the proposed general botmdary conditions. For relatively high coefficient of consolidation, just one term in series expansions is enough to obtain results with acceptable accuracy.