To effectively minimize the electromagnetic field response in the total field solution, we propose a numerical modeling method for the two-dimensional （2D） time- domain transient electromagnetic secondary field of t...To effectively minimize the electromagnetic field response in the total field solution, we propose a numerical modeling method for the two-dimensional （2D） time- domain transient electromagnetic secondary field of the line source based on the DuFort- Frankel finite-difference method. In the proposed method, we included the treatment of the earth-air boundary conductivity, calculated the normalized partial derivative of the induced electromotive force （Emf）, and determined the forward time step. By extending upward the earth-air interface to the air grid nodes and the zero-value boundary conditions, not only we have a method that is more efficient but also simpler than the total field solution. We computed and analyzed the homogeneous half-space model and the fiat layered model with high precision--the maximum relative error is less than 0.01% between our method and the analytical method--and the solution speed is roughly three times faster than the total-field solution. Lastly, we used the model of a thin body embedded in a homogeneous half-space at different delay times to depict the downward and upward spreading characteristics of the induced eddy current, and the physical interaction processes between the electromagnetic field and the underground low-resistivity body.展开更多
This study focuses on the factors that may affect the feasibility of performing elliptical anisotropy analysis on azimuthal PP and PS-wave data in HTI media, with the aim of using the modeling results as guidance in r...This study focuses on the factors that may affect the feasibility of performing elliptical anisotropy analysis on azimuthal PP and PS-wave data in HTI media, with the aim of using the modeling results as guidance in real seismic data application. Our results reveal that there is an offset limitation for both PP- and PS-waves in elliptical anisotropy fitting, and that PS-waves show a wider applicable offset range and larger observable azimuthal anisotropy than PP-waves. The major axis of the elliptical fit to the amplitudes of the R-component is perpendicular to the fracture strike, which is opposite to that in PP-wave analysis. The azimuthal interval travel time of PS-waves shows a nearly elliptical distribution and the major axis of the fit ellipse is perpendicular to the fracture strike, which is same as that in PP-wave analysis. For data within the applicable offset range, the anisotropic magnitude obtained from amplitude and travel time attributes of PP- and PS-waves exhibits a dependence on fracture density, and the major to minor axis ratio of the fit ellipse may be used to infer the relative distribution of fracture densities.展开更多
A three-dimensional numerical model of laser-arc hybrid plasma for aluminum alloy fillet joints is developed in this study.This mod-el accounts for the geometric complexity of fillet joints,the physical properties of ...A three-dimensional numerical model of laser-arc hybrid plasma for aluminum alloy fillet joints is developed in this study.This mod-el accounts for the geometric complexity of fillet joints,the physical properties of shielding gases with varying He-Ar ratios,and the coupling between arc plasma and laser-induced metal plume.The accuracy of the model is validated using a high-speed camera.The effects of varying He contents in the shielding gas on both the temperature and flow velocity of hybrid plasma,as well as the distribu-tion of laser-induced metal vapor mass,were investigated separately.The maximum temperature and size of arc plasma decrease as the He volume ratio increases,the arc distribution becomes more concentrated,and its flow velocity initially decreases and then sharply increases.At high helium content,both the flow velocity of hybrid plasma and metal vapor are high,the metal vapor is con-centrated on the right side of keyhole,and its flow appears chaotic.The flow state of arc plasma is most stable when the shielding gas consists of 50%He+50%Ar.展开更多
To improve the accuracy of the conventional finite-difference method, finitedifference numerical modeling methods of any even-order accuracy are recommended. We introduce any even-order accuracy difference schemes of ...To improve the accuracy of the conventional finite-difference method, finitedifference numerical modeling methods of any even-order accuracy are recommended. We introduce any even-order accuracy difference schemes of any-order derivatives derived from Taylor series expansion. Then, a finite-difference numerical modeling method with any evenorder accuracy is utilized to simulate seismic wave propagation in two-phase anisotropic media. Results indicate that modeling accuracy improves with the increase of difference accuracy order number. It is essential to find the optimal order number, grid size, and time step to balance modeling precision and computational complexity. Four kinds of waves, static mode in the source point, SV wave cusps, reflection and transmission waves are observed in two-phase anisotropic media through modeling.展开更多
Haihua Islands is a large artificial island in Danzhou, Hainan. The construction of Haihua Islands changes the hYdrodynamic environment of Yangpu waters, and further affects its morphological change. Delft3D is used t...Haihua Islands is a large artificial island in Danzhou, Hainan. The construction of Haihua Islands changes the hYdrodynamic environment of Yangpu waters, and further affects its morphological change. Delft3D is used to set up a two dimensional nested hydrodynamic and sediment model for Yangpu waters in this paper, and this paper focuses on simulating the velocity and morphological change due to the construction of Haihua Islands after the verification of the model. The seabed deposition is small because of low suspended sediment concentration and less sand source near Yangpu waters. The bed level erodes in the south area of Xiaochan Reef and the Yangpu channel due to the velocity increase in the area.展开更多
In this paper, a heavy sea fog episode that occurred over the Yellow Sea on 9 March 2005 is investigated. The sea fog patch, with a spatial scale of several hundred kilometers at its mature stage, reduced visibility a...In this paper, a heavy sea fog episode that occurred over the Yellow Sea on 9 March 2005 is investigated. The sea fog patch, with a spatial scale of several hundred kilometers at its mature stage, reduced visibility along the Shandong Peninsula coast to 100 m or much less at some sites. Satellite images, surface observations and soundings at islands and coasts, and analyses from the Japan Meteorology Agency （JMA） axe used to describe and analyze this event. The analysis indicates that this sea fog can be categorized as advection cooling fog. The main features of this sea fog including fog area and its movement axe reasonably reproduced by the Fifth-generation Pennsylvania State University/National Center for Atmospheric Research Mesoscale Model （MM5）. Model results suggest that the formation and evolution of this event can be outlined as： （1） southerly warm/moist advection of low-level air resulted in a strong sea-surface-based inversion with a thickness of about 600 m; （2） when the inversion moved from the warmer East Sea to the colder Yellow Sea, a thermal internal boundary layer （TIBL） gradually formed at the base of the inversion while the sea fog grew in response to cooling and moistening by turbulence mixing; （3） the sea fog developed as the TIBL moved northward and （4） strong northerly cold and dry wind destroyed the TIBL and dissipated the sea fog. The principal findings of this study axe that sea fog forms in response to relatively persistent southerly waxm/moist wind and a cold sea surface, and that turbulence mixing by wind shear is the primary mechanism for the cooling and moistening the marine layer. In addition, the study of sensitivity experiments indicates that deterministic numerical modeling offers a promising approach to the prediction of sea fog over the Yellow Sea but it may be more efficient to consider ensemble numerical modeling because of the extreme sensitivity to model input.展开更多
Severe shield jamming events have been reported during excavation of Uluabat tunnel through adverse geological conditions, which resulted in several stoppages at advancing a single shielded tunnel boring machine(TBM)....Severe shield jamming events have been reported during excavation of Uluabat tunnel through adverse geological conditions, which resulted in several stoppages at advancing a single shielded tunnel boring machine(TBM). To study the jamming mechanism, three-dimensional(3D) simulation of the machine and surrounding ground was implemented using the finite difference code FLAC3D. Numerical analyses were performed for three sections along the tunnel with a higher risk for entrapment due to the combination of overburden and geological conditions. The computational results including longitudinal displacement contours and ground pressure profiles around the shield allow a better understanding of ground behavior within the excavation. Furthermore, they allow realistically assessing the impact of adverse geological conditions on shield jamming. The calculated thrust forces, which are required to move the machine forward, are in good agreement with field observations and measurements. It also proves that the numerical analysis can effectively be used for evaluating the effect of adverse geological environment on TBM entrapments and can be applied to prediction of loads on the shield and preestimating of the required thrust force during excavation through adverse ground conditions.展开更多
A coupled discrete-continuum simulation incorporating a 3D aspect and non-circular particles was performed to analyze soil-pile interactions during pile penetration in sand.A self-developed non-circular particle numer...A coupled discrete-continuum simulation incorporating a 3D aspect and non-circular particles was performed to analyze soil-pile interactions during pile penetration in sand.A self-developed non-circular particle numerical simulation program was used which considered sand near the pile as interacted particles using a discrete element method;the sand away from the pile was simulated as a continuous medium exhibiting linear elastic behaviors.The domain analyzed was divided into two zones.Contact forces at the interface between the two zones were obtained from a discrete zone and applied to the continuum boundaries as nodal forces,while the interface velocities were obtained from the continuum zone and applied to the discrete boundaries.We show that the coupled discrete-continuum simulation can give a microscopic description of the pile penetration process without losing the discrete nature of the zone concerned,and may significantly improve computational efficiency.展开更多
As the depth of excavation increases,rockburst becomes one of the most serious geological hazards damaging equipment and facilities and even causing fatalities in mining and civil engineering.This has forced researche...As the depth of excavation increases,rockburst becomes one of the most serious geological hazards damaging equipment and facilities and even causing fatalities in mining and civil engineering.This has forced researchers worldwide to identify different methods to investigate rockburst-related problems.However,some problems,such as the mechanisms and the prediction of rockbursts,continue to be studied because rockburst is a very complicated phenomenon influenced by the uncertainty and complexity in geological conditions,in situ stresses,induced stresses,etc.Numerical modeling is a widely used method for investigating rockbursts.To date,great achievements have been made owing to the rapid development of information technology(IT)and computer equipment.Hence,it is necessary and meaningful to conduct a review of the current state of the studies for rockburst numerical modeling.In this paper,the categories and the origin of different numerical approaches employed in modeling rockbursts are reviewed and the current usage of various numerical modeling approaches is investigated by a literature research.Later,a state-of-the-art review is implemented to investigate the application of numerical modeling in the mechanism study,and prediction and prevention of rockbursts.The main achievements and problems are highlighted.Finally,this paper discusses the limitations and the future research of numerical modeling for rockbursts.An approach is proposed to provide researchers with a systematic and reasonable numerical modeling framework.展开更多
This article provides the application of the high-order, staggered-grid, finite-difference scheme to model elastic wave propagation in 3-D isotropic media. Here, we use second-order, tempo- ral- and high-order spatial...This article provides the application of the high-order, staggered-grid, finite-difference scheme to model elastic wave propagation in 3-D isotropic media. Here, we use second-order, tempo- ral- and high-order spatial finite-difference formulations with a staggered grid for discretization of the 3-D elastic wave equations of motion. The set of absorbing boundary conditions based on paraxial approximations of 3-D elastic wave equations are applied to the numerical boundaries. The trial re- sults for the salt model show that the numerical dispersion is decreased to a minimum extent, the accuracy high and diffracted waves abundant. It also shows that this method can be used for modeling wave propagation in complex media with the lateral variation of velocity.展开更多
To model the damage process of masonry walls under blast loading, a dynamic continuum damage material model is constructed for brick and mortar separately. The degradation of both the stiffness and strength are govern...To model the damage process of masonry walls under blast loading, a dynamic continuum damage material model is constructed for brick and mortar separately. The degradation of both the stiffness and strength are governed by a damage variable. By using the proposed material model, damage and fragmentation of a typical masonry wall under blast loading at different scaled distances is calculated. The hazard level of the masonry wall to blast loading is evaluated by analyzing the numerical results.展开更多
Estimation of the rock mass modulus of deformation(Em)is one of the most important design parameters in designing many structures in and on rock.This parameter can be obtained by in situ tests,empirical relations betw...Estimation of the rock mass modulus of deformation(Em)is one of the most important design parameters in designing many structures in and on rock.This parameter can be obtained by in situ tests,empirical relations between deformation modulus and rock mass classifcation,and estimating from laboratory tests results.In this paper,a back analysis calculation is performed to present an equation for estimation of the rock mass modulus of deformation using genetic programming(GP)and numerical modeling.A database of 40,960 datasets,including vertical stress(rz),horizontal to vertical stresses ratio(k),Poisson’s ratio(m),radius of circular tunnel(r)and wall displacement of circular tunnel on the horizontal diameter(d)for input parameters and modulus of deformation for output,was established.The selected parameters are easy to determine and rock mass modulus of deformation can be obtained from instrumentation data of any size circular galleries.The resulting RMSE of 0.86 and correlation coeffcient of97%of the proposed equation demonstrated the capability of the computer program(CP)generated by GP.展开更多
Widespread magmatism, metamorphic core complexes(MCCs), and significant lithospheric thinning occurred during the Mesozoic in the North China Craton(NCC). It has been suggested that the coeval exhumation of MCCs with ...Widespread magmatism, metamorphic core complexes(MCCs), and significant lithospheric thinning occurred during the Mesozoic in the North China Craton(NCC). It has been suggested that the coeval exhumation of MCCs with uniform northwest-southeast shear senses and magmatism probably resulted from a decratonization event during the retreat of the paleo-Pacific Plate. Here we used two-dimensional finite element thermomechanical numerical models to investigate critical parameters controlling the formation of MCCs under far-field extensional stress. We observed three end-member deformation modes: the MCC mode, the symmetric-dome mode, and the pure-shear mode. The MCC mode requires a Moho temperature of ≥700 ℃ and an extensional strain rate of ≥5 × 10^(-16)s^(-1), implying that the lithosphere had already thinned when the MCC was formed in the Mesozoic. Considering that the widespread MCCs have the same northwest-southeast extension direction in the NCC, we suggest that the MCCs are surface expressions of both large-scale extension and craton destruction and that rollback of the paleo-Pacific slab might be the common driving force.展开更多
Finite-difference methods with high-order accuracy have been utilized to improve the precision of numerical solution for partial differential equations. However, the computation cost generally increases linearly with ...Finite-difference methods with high-order accuracy have been utilized to improve the precision of numerical solution for partial differential equations. However, the computation cost generally increases linearly with increased order of accuracy. Upon examination of the finite-difference formulas for the first-order and second-order derivatives, and the staggered finite-difference formulas for the first-order derivative, we examine the variation of finite-difference coefficients with accuracy order and note that there exist some very small coefficients. With the order increasing, the number of these small coefficients increases, however, the values decrease sharply. An error analysis demonstrates that omitting these small coefficients not only maintain approximately the same level of accuracy of finite difference but also reduce computational cost significantly. Moreover, it is easier to truncate for the high-order finite-difference formulas than for the pseudospectral for- mulas. Thus this study proposes a truncated high-order finite-difference method, and then demonstrates the efficiency and applicability of the method with some numerical examples.展开更多
The mild-slope equation is familiar to coastal engineers as it can effectively describe wave propagation in nearshore regions. However, its computational method in Cartesian coordinates often renders the model inaccur...The mild-slope equation is familiar to coastal engineers as it can effectively describe wave propagation in nearshore regions. However, its computational method in Cartesian coordinates often renders the model inaccurate in areas with irregular shorelines, such as estuaries and harbors. Based on the hyperbolic mild-slope equation in Cartesian coordinates, the numerical model in orthogonal curvilinear coordinates is developed. The transformed model is discretized by the finite difference method and solved by the ADI method with space-staggered grids. The numerical predictions in curvilinear co- ordinates show good agreemenl with the data obtained in three typical physical expedments, which demonstrates that the present model can be used to simulate wave propagation, for normal incidence and oblique incidence, in domains with complicated topography and boundary conditions.展开更多
This paper describes the interaction between deep-seated landslides and man-made structures such as dams, penstocks, viaducts, and tunnels. Selected case studies are reported first with the intent to gain insights int...This paper describes the interaction between deep-seated landslides and man-made structures such as dams, penstocks, viaducts, and tunnels. Selected case studies are reported first with the intent to gain insights into the complexities associated with the interaction of these structures with deep-seated landslides(generally referred to as deep-seated gravity slope deformations, DSGSDs). The main features, which characterize these landslides, are mentioned together with the interaction problems encountered in each case. Given the main objective of this paper, the numerical modeling methods adopted are outlined as means for increase in the understanding of the interaction problems being investigated. With the above in mind, the attention moves to an important and unique case history dealing with the interaction of a large-size twin-tunnel excavated with an earth pressure balance(EPB)tunnel boring machine(TBM) and a deep-seated landslide, which was reactivated due to the stress changes induced by tunnel excavation in landslide shear zone. The geological and geotechnical conditions are described together with the available monitoring data on the landslide movements, based on the advanced and conventional monitoring tools used. Numerical modeling is illustrated as an aid to back-analyze the monitored surface and subsurface deformations and to assist in finding the appropriate engineering solution for putting the tunnel into service and as a follow-up means for future understanding and control of the interaction problems. The simulation is based on a novel time-dependent model representing the landslide behavior.展开更多
In order to reasonably simulate tidal currents around small structures such as piles in a large-scale model domain, a 2-D hydrodynamic integrated model for Bohai Sea is established with the finite element method. The ...In order to reasonably simulate tidal currents around small structures such as piles in a large-scale model domain, a 2-D hydrodynamic integrated model for Bohai Sea is established with the finite element method. The grid can be discretionarily refined as a non-structure triangle or quadrilateral so that piers can be treated as one or several impermeable elements with an area of 20 to 30 km^2 in a model domain over 85 700 km^2. The computational results of tidal levels and horizontal velocities are in good agreement with the field data. Based on the computed results by the model, the layout of an open 105 DWT liquefied natural gas （LNG）terminal in Caofeidian, Bohal Sea is effectively and reasonably optimized. It can be concluded that the model is suitable and reasonable for direct simulation of tidal currents around small structures in projects.展开更多
Virtual simulation is an economical and efficient method in mechanical system design. Numerical modeling of a spar platform, tethered by a mooring cable with a spherical joint is developed for the dynamic simulation o...Virtual simulation is an economical and efficient method in mechanical system design. Numerical modeling of a spar platform, tethered by a mooring cable with a spherical joint is developed for the dynamic simulation of the floating structure in ocean. The geometry modeling of the spar is created using finite element methods. The submerged part of the spar bears the buoyancy, hydrodynamic drag force, and effect of the added mass and Froude-Krylov force. Strip theory is used to sum up the forces acting on the elements. The geometry modeling of the cable is established based on the lumped-mass-and-spring modeling through which the cable is divided into 10 elements. A new element-fixed local frame is used, which is created by the element orientation vector and relative velocity of the fluid, to express the loads acting on the cable. The bottom of the cable is fixed on the seabed by spring forces, while the top of the cable is connected to the bottom of the spar platform by a modified spherical joint. This system suffers the propagating wave and current in the X-direction and the linear wave theory is applied for setting of the propagating wave. Based on the numerical modeling, the displacement-load relationships are analyzed, and the simulation results of the numerical modeling are compared with those by the commercial simulation code, Proteus DS. The comparison indicates that the numerical modeling of the spar platform tethered by a mooring cable is well developed, which provides an instruction for the optimization of a floating structure tethered by a mooring cable system.展开更多
The Qingchengzi ore field is an important gold-polymetallic center of the North China Cra-ton.It has been recognized that the gold deposits in Qingchengzi were controlled by structures like litho-logical interfaces an...The Qingchengzi ore field is an important gold-polymetallic center of the North China Cra-ton.It has been recognized that the gold deposits in Qingchengzi were controlled by structures like litho-logical interfaces and fractures along mechanically weak bedding and foliation planes,but it still remains poorly understood how the structures affected the localization of the gold deposits.Finite element based numerical modeling was used to reproduce the deformation process of the Baiyun gold deposit during the mineralization period.Paleoproterozoic schist and marble are widely exposed in Qingchengzi,and a large part of the Baiyun gold ores occurs along the interfaces between the schist and the marble.The modeling results suggest that the mechanical contrast between the schist and the marble may be a major reason why the stress was localized along their lithological interfaces under a compressional stress regime.Two parts of their lithological interfaces were identified to be easily stress-localized and first fractured:the interface between the schist and its underlying marble at shallower levels and the one between the schist and its overlying marble at deeper levels.Stress concentration in these two parts is independent on the dipping angle and direction of the interfaces.Therefore,mineralizing fluids may have been concentrated into these two parts.The first one is consistent with the present ore bodies of the Baiyun gold deposit,and the second one could be considered for deep prospecting.These findings also provide implications for the structural controls of lithological interfaces on the mineralization in other gold deposits of this region.展开更多
The finite volume method （FVM） has many advantages in 2-D shallow water numerical simulation. In this study, the finite volume method is used with unstructured triangular grids to simulate the tidal currents. The Ro...The finite volume method （FVM） has many advantages in 2-D shallow water numerical simulation. In this study, the finite volume method is used with unstructured triangular grids to simulate the tidal currents. The Roe scheme is applied in the calculation of the intercell numerical flux, and the MUSCL method is introduced to improve its accuracy. The time integral is a two-step scheme of forecast and revision. For the verification of the present method, the Stoker＇s problem is calculated and the result is compared with the mathematically analytic solutions. The comparison indicates that the method is feasible. A sea area of a port is used as an example to test the method established here. The result shows that the present computational method is satisfactory, and it could be applied to the engineering fields.展开更多
基金supported by the National High Technology Research and Development Program (863 Program)(2009AA06Z108)
文摘To effectively minimize the electromagnetic field response in the total field solution, we propose a numerical modeling method for the two-dimensional （2D） time- domain transient electromagnetic secondary field of the line source based on the DuFort- Frankel finite-difference method. In the proposed method, we included the treatment of the earth-air boundary conductivity, calculated the normalized partial derivative of the induced electromotive force （Emf）, and determined the forward time step. By extending upward the earth-air interface to the air grid nodes and the zero-value boundary conditions, not only we have a method that is more efficient but also simpler than the total field solution. We computed and analyzed the homogeneous half-space model and the fiat layered model with high precision--the maximum relative error is less than 0.01% between our method and the analytical method--and the solution speed is roughly three times faster than the total-field solution. Lastly, we used the model of a thin body embedded in a homogeneous half-space at different delay times to depict the downward and upward spreading characteristics of the induced eddy current, and the physical interaction processes between the electromagnetic field and the underground low-resistivity body.
文摘This study focuses on the factors that may affect the feasibility of performing elliptical anisotropy analysis on azimuthal PP and PS-wave data in HTI media, with the aim of using the modeling results as guidance in real seismic data application. Our results reveal that there is an offset limitation for both PP- and PS-waves in elliptical anisotropy fitting, and that PS-waves show a wider applicable offset range and larger observable azimuthal anisotropy than PP-waves. The major axis of the elliptical fit to the amplitudes of the R-component is perpendicular to the fracture strike, which is opposite to that in PP-wave analysis. The azimuthal interval travel time of PS-waves shows a nearly elliptical distribution and the major axis of the fit ellipse is perpendicular to the fracture strike, which is same as that in PP-wave analysis. For data within the applicable offset range, the anisotropic magnitude obtained from amplitude and travel time attributes of PP- and PS-waves exhibits a dependence on fracture density, and the major to minor axis ratio of the fit ellipse may be used to infer the relative distribution of fracture densities.
基金supported by the National Natural Science Foundation of China(Grant No.52375340,51975263,52405366).
文摘A three-dimensional numerical model of laser-arc hybrid plasma for aluminum alloy fillet joints is developed in this study.This mod-el accounts for the geometric complexity of fillet joints,the physical properties of shielding gases with varying He-Ar ratios,and the coupling between arc plasma and laser-induced metal plume.The accuracy of the model is validated using a high-speed camera.The effects of varying He contents in the shielding gas on both the temperature and flow velocity of hybrid plasma,as well as the distribu-tion of laser-induced metal vapor mass,were investigated separately.The maximum temperature and size of arc plasma decrease as the He volume ratio increases,the arc distribution becomes more concentrated,and its flow velocity initially decreases and then sharply increases.At high helium content,both the flow velocity of hybrid plasma and metal vapor are high,the metal vapor is con-centrated on the right side of keyhole,and its flow appears chaotic.The flow state of arc plasma is most stable when the shielding gas consists of 50%He+50%Ar.
文摘To improve the accuracy of the conventional finite-difference method, finitedifference numerical modeling methods of any even-order accuracy are recommended. We introduce any even-order accuracy difference schemes of any-order derivatives derived from Taylor series expansion. Then, a finite-difference numerical modeling method with any evenorder accuracy is utilized to simulate seismic wave propagation in two-phase anisotropic media. Results indicate that modeling accuracy improves with the increase of difference accuracy order number. It is essential to find the optimal order number, grid size, and time step to balance modeling precision and computational complexity. Four kinds of waves, static mode in the source point, SV wave cusps, reflection and transmission waves are observed in two-phase anisotropic media through modeling.
文摘Haihua Islands is a large artificial island in Danzhou, Hainan. The construction of Haihua Islands changes the hYdrodynamic environment of Yangpu waters, and further affects its morphological change. Delft3D is used to set up a two dimensional nested hydrodynamic and sediment model for Yangpu waters in this paper, and this paper focuses on simulating the velocity and morphological change due to the construction of Haihua Islands after the verification of the model. The seabed deposition is small because of low suspended sediment concentration and less sand source near Yangpu waters. The bed level erodes in the south area of Xiaochan Reef and the Yangpu channel due to the velocity increase in the area.
文摘In this paper, a heavy sea fog episode that occurred over the Yellow Sea on 9 March 2005 is investigated. The sea fog patch, with a spatial scale of several hundred kilometers at its mature stage, reduced visibility along the Shandong Peninsula coast to 100 m or much less at some sites. Satellite images, surface observations and soundings at islands and coasts, and analyses from the Japan Meteorology Agency （JMA） axe used to describe and analyze this event. The analysis indicates that this sea fog can be categorized as advection cooling fog. The main features of this sea fog including fog area and its movement axe reasonably reproduced by the Fifth-generation Pennsylvania State University/National Center for Atmospheric Research Mesoscale Model （MM5）. Model results suggest that the formation and evolution of this event can be outlined as： （1） southerly warm/moist advection of low-level air resulted in a strong sea-surface-based inversion with a thickness of about 600 m; （2） when the inversion moved from the warmer East Sea to the colder Yellow Sea, a thermal internal boundary layer （TIBL） gradually formed at the base of the inversion while the sea fog grew in response to cooling and moistening by turbulence mixing; （3） the sea fog developed as the TIBL moved northward and （4） strong northerly cold and dry wind destroyed the TIBL and dissipated the sea fog. The principal findings of this study axe that sea fog forms in response to relatively persistent southerly waxm/moist wind and a cold sea surface, and that turbulence mixing by wind shear is the primary mechanism for the cooling and moistening the marine layer. In addition, the study of sensitivity experiments indicates that deterministic numerical modeling offers a promising approach to the prediction of sea fog over the Yellow Sea but it may be more efficient to consider ensemble numerical modeling because of the extreme sensitivity to model input.
基金Alexander von Humboldt-Foundation (AvH) for the financial support as a research fellowthe financial support of the Scientific and Technological Research Council of Turkey (TüB_ITAK) under Project No. MAG-114M568
文摘Severe shield jamming events have been reported during excavation of Uluabat tunnel through adverse geological conditions, which resulted in several stoppages at advancing a single shielded tunnel boring machine(TBM). To study the jamming mechanism, three-dimensional(3D) simulation of the machine and surrounding ground was implemented using the finite difference code FLAC3D. Numerical analyses were performed for three sections along the tunnel with a higher risk for entrapment due to the combination of overburden and geological conditions. The computational results including longitudinal displacement contours and ground pressure profiles around the shield allow a better understanding of ground behavior within the excavation. Furthermore, they allow realistically assessing the impact of adverse geological conditions on shield jamming. The calculated thrust forces, which are required to move the machine forward, are in good agreement with field observations and measurements. It also proves that the numerical analysis can effectively be used for evaluating the effect of adverse geological environment on TBM entrapments and can be applied to prediction of loads on the shield and preestimating of the required thrust force during excavation through adverse ground conditions.
基金Project (No.90815008) supported by the National Natural Science Foundation of China
文摘A coupled discrete-continuum simulation incorporating a 3D aspect and non-circular particles was performed to analyze soil-pile interactions during pile penetration in sand.A self-developed non-circular particle numerical simulation program was used which considered sand near the pile as interacted particles using a discrete element method;the sand away from the pile was simulated as a continuous medium exhibiting linear elastic behaviors.The domain analyzed was divided into two zones.Contact forces at the interface between the two zones were obtained from a discrete zone and applied to the continuum boundaries as nodal forces,while the interface velocities were obtained from the continuum zone and applied to the discrete boundaries.We show that the coupled discrete-continuum simulation can give a microscopic description of the pile penetration process without losing the discrete nature of the zone concerned,and may significantly improve computational efficiency.
基金The authors gratefully acknowledge financial support from the China Scholarship Council(Grant No.201808370185).
文摘As the depth of excavation increases,rockburst becomes one of the most serious geological hazards damaging equipment and facilities and even causing fatalities in mining and civil engineering.This has forced researchers worldwide to identify different methods to investigate rockburst-related problems.However,some problems,such as the mechanisms and the prediction of rockbursts,continue to be studied because rockburst is a very complicated phenomenon influenced by the uncertainty and complexity in geological conditions,in situ stresses,induced stresses,etc.Numerical modeling is a widely used method for investigating rockbursts.To date,great achievements have been made owing to the rapid development of information technology(IT)and computer equipment.Hence,it is necessary and meaningful to conduct a review of the current state of the studies for rockburst numerical modeling.In this paper,the categories and the origin of different numerical approaches employed in modeling rockbursts are reviewed and the current usage of various numerical modeling approaches is investigated by a literature research.Later,a state-of-the-art review is implemented to investigate the application of numerical modeling in the mechanism study,and prediction and prevention of rockbursts.The main achievements and problems are highlighted.Finally,this paper discusses the limitations and the future research of numerical modeling for rockbursts.An approach is proposed to provide researchers with a systematic and reasonable numerical modeling framework.
文摘This article provides the application of the high-order, staggered-grid, finite-difference scheme to model elastic wave propagation in 3-D isotropic media. Here, we use second-order, tempo- ral- and high-order spatial finite-difference formulations with a staggered grid for discretization of the 3-D elastic wave equations of motion. The set of absorbing boundary conditions based on paraxial approximations of 3-D elastic wave equations are applied to the numerical boundaries. The trial re- sults for the salt model show that the numerical dispersion is decreased to a minimum extent, the accuracy high and diffracted waves abundant. It also shows that this method can be used for modeling wave propagation in complex media with the lateral variation of velocity.
文摘To model the damage process of masonry walls under blast loading, a dynamic continuum damage material model is constructed for brick and mortar separately. The degradation of both the stiffness and strength are governed by a damage variable. By using the proposed material model, damage and fragmentation of a typical masonry wall under blast loading at different scaled distances is calculated. The hazard level of the masonry wall to blast loading is evaluated by analyzing the numerical results.
文摘Estimation of the rock mass modulus of deformation(Em)is one of the most important design parameters in designing many structures in and on rock.This parameter can be obtained by in situ tests,empirical relations between deformation modulus and rock mass classifcation,and estimating from laboratory tests results.In this paper,a back analysis calculation is performed to present an equation for estimation of the rock mass modulus of deformation using genetic programming(GP)and numerical modeling.A database of 40,960 datasets,including vertical stress(rz),horizontal to vertical stresses ratio(k),Poisson’s ratio(m),radius of circular tunnel(r)and wall displacement of circular tunnel on the horizontal diameter(d)for input parameters and modulus of deformation for output,was established.The selected parameters are easy to determine and rock mass modulus of deformation can be obtained from instrumentation data of any size circular galleries.The resulting RMSE of 0.86 and correlation coeffcient of97%of the proposed equation demonstrated the capability of the computer program(CP)generated by GP.
基金supported by the National Natural Science Foundation of China(Grant No.41774112)。
文摘Widespread magmatism, metamorphic core complexes(MCCs), and significant lithospheric thinning occurred during the Mesozoic in the North China Craton(NCC). It has been suggested that the coeval exhumation of MCCs with uniform northwest-southeast shear senses and magmatism probably resulted from a decratonization event during the retreat of the paleo-Pacific Plate. Here we used two-dimensional finite element thermomechanical numerical models to investigate critical parameters controlling the formation of MCCs under far-field extensional stress. We observed three end-member deformation modes: the MCC mode, the symmetric-dome mode, and the pure-shear mode. The MCC mode requires a Moho temperature of ≥700 ℃ and an extensional strain rate of ≥5 × 10^(-16)s^(-1), implying that the lithosphere had already thinned when the MCC was formed in the Mesozoic. Considering that the widespread MCCs have the same northwest-southeast extension direction in the NCC, we suggest that the MCCs are surface expressions of both large-scale extension and craton destruction and that rollback of the paleo-Pacific slab might be the common driving force.
基金supported by China Scholarship Council and partially by the National "863" Program of China under contract No. 2007AA06Z218.
文摘Finite-difference methods with high-order accuracy have been utilized to improve the precision of numerical solution for partial differential equations. However, the computation cost generally increases linearly with increased order of accuracy. Upon examination of the finite-difference formulas for the first-order and second-order derivatives, and the staggered finite-difference formulas for the first-order derivative, we examine the variation of finite-difference coefficients with accuracy order and note that there exist some very small coefficients. With the order increasing, the number of these small coefficients increases, however, the values decrease sharply. An error analysis demonstrates that omitting these small coefficients not only maintain approximately the same level of accuracy of finite difference but also reduce computational cost significantly. Moreover, it is easier to truncate for the high-order finite-difference formulas than for the pseudospectral for- mulas. Thus this study proposes a truncated high-order finite-difference method, and then demonstrates the efficiency and applicability of the method with some numerical examples.
基金supported by the National Basic Research Program of China ( Grant No.2006CB403302)the National Natural Science Foundation of China (Grant Nos .50839001 and 50709004)the Scientific Research Foundation of the Higher Education Institutions of Liaoning Province (Grant No.2006T018)
文摘The mild-slope equation is familiar to coastal engineers as it can effectively describe wave propagation in nearshore regions. However, its computational method in Cartesian coordinates often renders the model inaccurate in areas with irregular shorelines, such as estuaries and harbors. Based on the hyperbolic mild-slope equation in Cartesian coordinates, the numerical model in orthogonal curvilinear coordinates is developed. The transformed model is discretized by the finite difference method and solved by the ADI method with space-staggered grids. The numerical predictions in curvilinear co- ordinates show good agreemenl with the data obtained in three typical physical expedments, which demonstrates that the present model can be used to simulate wave propagation, for normal incidence and oblique incidence, in domains with complicated topography and boundary conditions.
基金support of Spea Ingegneria Europea SpA and Società Autostrade per l’Italia SpA
文摘This paper describes the interaction between deep-seated landslides and man-made structures such as dams, penstocks, viaducts, and tunnels. Selected case studies are reported first with the intent to gain insights into the complexities associated with the interaction of these structures with deep-seated landslides(generally referred to as deep-seated gravity slope deformations, DSGSDs). The main features, which characterize these landslides, are mentioned together with the interaction problems encountered in each case. Given the main objective of this paper, the numerical modeling methods adopted are outlined as means for increase in the understanding of the interaction problems being investigated. With the above in mind, the attention moves to an important and unique case history dealing with the interaction of a large-size twin-tunnel excavated with an earth pressure balance(EPB)tunnel boring machine(TBM) and a deep-seated landslide, which was reactivated due to the stress changes induced by tunnel excavation in landslide shear zone. The geological and geotechnical conditions are described together with the available monitoring data on the landslide movements, based on the advanced and conventional monitoring tools used. Numerical modeling is illustrated as an aid to back-analyze the monitored surface and subsurface deformations and to assist in finding the appropriate engineering solution for putting the tunnel into service and as a follow-up means for future understanding and control of the interaction problems. The simulation is based on a novel time-dependent model representing the landslide behavior.
文摘In order to reasonably simulate tidal currents around small structures such as piles in a large-scale model domain, a 2-D hydrodynamic integrated model for Bohai Sea is established with the finite element method. The grid can be discretionarily refined as a non-structure triangle or quadrilateral so that piers can be treated as one or several impermeable elements with an area of 20 to 30 km^2 in a model domain over 85 700 km^2. The computational results of tidal levels and horizontal velocities are in good agreement with the field data. Based on the computed results by the model, the layout of an open 105 DWT liquefied natural gas （LNG）terminal in Caofeidian, Bohal Sea is effectively and reasonably optimized. It can be concluded that the model is suitable and reasonable for direct simulation of tidal currents around small structures in projects.
基金Supported by Human Resources Development Program of Korea Institute of Energy Technology Evaluation and Planning(KETEP)Ministry of Trade,Industry and Energy of Korea(Grant No.20134030200290)
文摘Virtual simulation is an economical and efficient method in mechanical system design. Numerical modeling of a spar platform, tethered by a mooring cable with a spherical joint is developed for the dynamic simulation of the floating structure in ocean. The geometry modeling of the spar is created using finite element methods. The submerged part of the spar bears the buoyancy, hydrodynamic drag force, and effect of the added mass and Froude-Krylov force. Strip theory is used to sum up the forces acting on the elements. The geometry modeling of the cable is established based on the lumped-mass-and-spring modeling through which the cable is divided into 10 elements. A new element-fixed local frame is used, which is created by the element orientation vector and relative velocity of the fluid, to express the loads acting on the cable. The bottom of the cable is fixed on the seabed by spring forces, while the top of the cable is connected to the bottom of the spar platform by a modified spherical joint. This system suffers the propagating wave and current in the X-direction and the linear wave theory is applied for setting of the propagating wave. Based on the numerical modeling, the displacement-load relationships are analyzed, and the simulation results of the numerical modeling are compared with those by the commercial simulation code, Proteus DS. The comparison indicates that the numerical modeling of the spar platform tethered by a mooring cable is well developed, which provides an instruction for the optimization of a floating structure tethered by a mooring cable system.
基金The work was financially funded by the National Key R&D Program of China(No.2018YFC0603802)the Basic Research Fund for Central Research Institutes(No.JYYWF20180602)the National Natural Science Foundation of China(No.41822206).
文摘The Qingchengzi ore field is an important gold-polymetallic center of the North China Cra-ton.It has been recognized that the gold deposits in Qingchengzi were controlled by structures like litho-logical interfaces and fractures along mechanically weak bedding and foliation planes,but it still remains poorly understood how the structures affected the localization of the gold deposits.Finite element based numerical modeling was used to reproduce the deformation process of the Baiyun gold deposit during the mineralization period.Paleoproterozoic schist and marble are widely exposed in Qingchengzi,and a large part of the Baiyun gold ores occurs along the interfaces between the schist and the marble.The modeling results suggest that the mechanical contrast between the schist and the marble may be a major reason why the stress was localized along their lithological interfaces under a compressional stress regime.Two parts of their lithological interfaces were identified to be easily stress-localized and first fractured:the interface between the schist and its underlying marble at shallower levels and the one between the schist and its overlying marble at deeper levels.Stress concentration in these two parts is independent on the dipping angle and direction of the interfaces.Therefore,mineralizing fluids may have been concentrated into these two parts.The first one is consistent with the present ore bodies of the Baiyun gold deposit,and the second one could be considered for deep prospecting.These findings also provide implications for the structural controls of lithological interfaces on the mineralization in other gold deposits of this region.
基金This paper was supported bythe Natural Science Foundation of Shandong Province (Grant No.y2004f13)
文摘The finite volume method （FVM） has many advantages in 2-D shallow water numerical simulation. In this study, the finite volume method is used with unstructured triangular grids to simulate the tidal currents. The Roe scheme is applied in the calculation of the intercell numerical flux, and the MUSCL method is introduced to improve its accuracy. The time integral is a two-step scheme of forecast and revision. For the verification of the present method, the Stoker＇s problem is calculated and the result is compared with the mathematically analytic solutions. The comparison indicates that the method is feasible. A sea area of a port is used as an example to test the method established here. The result shows that the present computational method is satisfactory, and it could be applied to the engineering fields.
