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 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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
In this paper the characteristics of tidal flow and seasonal variation of seidment content in the Hangzhou Bay and their affecting factors are studied. Field investigations and data analysis indicate that the sediment...In this paper the characteristics of tidal flow and seasonal variation of seidment content in the Hangzhou Bay and their affecting factors are studied. Field investigations and data analysis indicate that the sediment movement is mainly influenced by the Yangtze estuary and the sediment of the Yangtze estuary is induced by wind wave and tidal flow. Owing to the variation of dynamic conditions, the instantaneous sediment content is controlled by tidal flow, wind wave, depth of water and tidal range synthetically. A sediment content relationship formula is established with related factors. A non-equilibrium 2-dimensional numerical model of suspended sediment transportation is set up, and the finite element method is applied. The computation results of the model is in accordance with field data.展开更多
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.展开更多
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.展开更多
The development of functional relationships between the observed deposition rate and the experimental conditions is an important step toward understanding and optimizing low-pressure chemical vapor deposition(LPCVD)or...The development of functional relationships between the observed deposition rate and the experimental conditions is an important step toward understanding and optimizing low-pressure chemical vapor deposition(LPCVD)or low-pressure chemical vapor infiltration(LPCVI).In the field of ceramic matrix composites(CMCs),methyltrichlorosilane(CH3 SiCl3,MTS)is the most widely used source gas system for SiC,because stoichiometric SiC deposit can be facilitated at 900°C–1300°C.However,the reliability and accuracy of existing numerical models for these processing conditions are rarely reported.In this study,a comprehensive transport model was coupled with gas-phase and surface kinetics.The resulting gas-phase kinetics was confirmed via the measured concentration of gaseous species.The relationship between deposition rate and 24 gaseous species has been effectively evaluated by combining the special superiority of the novel extreme machine learning method and the conventional sticking coefficient method.Surface kinetics were then proposed and shown to reproduce the experimental results.The proposed simulation strategy can be used for different material systems.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
A process-based 3D numerical model for surfzone hydrodynamics and beach evolution was established. Comparisons between the experimental data and model results proved that the model could effectively describe the hydro...A process-based 3D numerical model for surfzone hydrodynamics and beach evolution was established. Comparisons between the experimental data and model results proved that the model could effectively describe the hydrodynamics, sediment transport feature and sandbar migration process in the surfzone with satisfactory precision. A series of numerical simulations on the wave breaking and shoaling up to a barred beach were carried out based on the model system. Analyzed from the model results, the wave-induced current system in the surfzone consists of two major processes, which are the phase-averaged undertow caused by wave breaking and the net drift caused by both of the nonlinear wave motion and surface roller effect. When storm waves come to the barred beach, the strong offshore undertow along the beach suppresses the onshore net drift, making the initial sandbar migrate to the seaside. Under the condition of calm wave environment, both the undertow and net drift flow to the shoreline at the offshore side of the sandbar, and then push the initial sandbar to the shoreline. The consideration of surface roller has significant impact on the modeling results of the sandbar migration. As the roller transfer rate increases, the sandbar moves onshore especially under the storm wave condition.展开更多
基金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.
基金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.
文摘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.
文摘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.
文摘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.
基金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 work was financially supported by the funds for Youth Scholar of the Ministry of Communications,P.R.China
文摘In this paper the characteristics of tidal flow and seasonal variation of seidment content in the Hangzhou Bay and their affecting factors are studied. Field investigations and data analysis indicate that the sediment movement is mainly influenced by the Yangtze estuary and the sediment of the Yangtze estuary is induced by wind wave and tidal flow. Owing to the variation of dynamic conditions, the instantaneous sediment content is controlled by tidal flow, wind wave, depth of water and tidal range synthetically. A sediment content relationship formula is established with related factors. A non-equilibrium 2-dimensional numerical model of suspended sediment transportation is set up, and the finite element method is applied. The computation results of the model is in accordance with field data.
基金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.
文摘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.
基金the National Key R&D Program of China(Grants No.2017YFB0703200)National Natural Science Foundation of China(Grants Nos.51702100,51972268)China Postdoctoral Science Foundation(Grants No.2018M643075)for financial support。
文摘The development of functional relationships between the observed deposition rate and the experimental conditions is an important step toward understanding and optimizing low-pressure chemical vapor deposition(LPCVD)or low-pressure chemical vapor infiltration(LPCVI).In the field of ceramic matrix composites(CMCs),methyltrichlorosilane(CH3 SiCl3,MTS)is the most widely used source gas system for SiC,because stoichiometric SiC deposit can be facilitated at 900°C–1300°C.However,the reliability and accuracy of existing numerical models for these processing conditions are rarely reported.In this study,a comprehensive transport model was coupled with gas-phase and surface kinetics.The resulting gas-phase kinetics was confirmed via the measured concentration of gaseous species.The relationship between deposition rate and 24 gaseous species has been effectively evaluated by combining the special superiority of the novel extreme machine learning method and the conventional sticking coefficient method.Surface kinetics were then proposed and shown to reproduce the experimental results.The proposed simulation strategy can be used for different material systems.
文摘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.
文摘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 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.
基金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.
基金financially supported by the National Key Research and Development Program of China(Grant No.2016YFC0402603)the National Natural Science Foundation of China(Grant Nos.51779112,51509119,and 51609029)+2 种基金the Project of Tianjin Natural Science Foundation(Grant No.16JCQNJC06900)the Open Project of State Key Laboratory of Hydrology-Water Resources and Hydraulic Engineering(Grant No.2014492211)the Fundamental Research Funds for the Central Public Welfare Research Institutes(Grant Nos.TKS170101and TKS170202)
文摘A process-based 3D numerical model for surfzone hydrodynamics and beach evolution was established. Comparisons between the experimental data and model results proved that the model could effectively describe the hydrodynamics, sediment transport feature and sandbar migration process in the surfzone with satisfactory precision. A series of numerical simulations on the wave breaking and shoaling up to a barred beach were carried out based on the model system. Analyzed from the model results, the wave-induced current system in the surfzone consists of two major processes, which are the phase-averaged undertow caused by wave breaking and the net drift caused by both of the nonlinear wave motion and surface roller effect. When storm waves come to the barred beach, the strong offshore undertow along the beach suppresses the onshore net drift, making the initial sandbar migrate to the seaside. Under the condition of calm wave environment, both the undertow and net drift flow to the shoreline at the offshore side of the sandbar, and then push the initial sandbar to the shoreline. The consideration of surface roller has significant impact on the modeling results of the sandbar migration. As the roller transfer rate increases, the sandbar moves onshore especially under the storm wave condition.
