Earth medium is not completely elastic, with its viscosity resulting in attenuation and dispersion of seismic waves. Most viscoelastic numerical simulations are based on the finite-difference and finite-element method...Earth medium is not completely elastic, with its viscosity resulting in attenuation and dispersion of seismic waves. Most viscoelastic numerical simulations are based on the finite-difference and finite-element methods. Targeted at viscoelastic numerical modeling for multilayered media, the constant-Q acoustic wave equation is transformed into the corresponding wave integral representation with its Green's function accounting for viscoelastic coefficients. An efficient alternative for full-waveform solution to the integral equation is proposed in this article by extending conventional frequency-domain boundary element methods to viscoelastic media. The viscoelastic boundary element method enjoys a distinct characteristic of the explicit use of boundary continuity conditions of displacement and traction, leading to a semi-analytical solution with sufficient accuracy for simulating the viscoelastic effect across irregular interfaces. Numerical experiments to study the viscoelastic absorption of different Q values demonstrate the accuracy and applicability of the method.展开更多
An analytical model with essential parameters given by a two-phase numerical model is utilized to study the net boundary layer current and sediment transport under skewed asymmetric oscillatory sheet flows. The analyt...An analytical model with essential parameters given by a two-phase numerical model is utilized to study the net boundary layer current and sediment transport under skewed asymmetric oscillatory sheet flows. The analytical model is the first instantaneous type model that can consider phase-lag and asymmetric boundary layer development. The two-phase model supplies the essential phase-lead, instantaneous erosion depth and boundary layer development for the analytical model to enhance the understanding of velocity skewness and acceleration skewness in sediment flux and transport rate. The sediment transport difference between onshore and offshore stages caused by velocity skewness or acceleration skewness is shown to illustrate the determination of net sediment transport by the analytical model. In previous studies about sediment transport in skewed asymmetric sheet flows, the generation of net sediment transport is mainly concluded to the phase-lag effect.However, the phase-lag effect is shown important but not enough for the net sediment transport, while the skewed asymmetric boundary layer development generated net boundary layer current and mobile bed effect are key important in the transport process.展开更多
In this work, an enhanced treatment of the solid boundaries is proposed for smoothed particle hydrodynamics with implicit time integration scheme (Implicit SPH). Three types of virtual particles, i.e., boundary part...In this work, an enhanced treatment of the solid boundaries is proposed for smoothed particle hydrodynamics with implicit time integration scheme (Implicit SPH). Three types of virtual particles, i.e., boundary particles, image particles and mirror particles, are used to impose boundary conditions. Boundary particles are fixed on the solid boundary, and each boundary particle is associated with two fixed image particles inside the fluid domain and two fixed mirror particles outside the fluid domain. The image particles take the flow properties through fluid particles with moving least squares (MLS) interpolation and the properties of mirror particles can be obtained by the corresponding image particles. A repulsive force is also applied for boundary particles to prevent fluid particles from unphysical penetra- tion through solid boundaries. The new boundary treatment method has been validated with five numerical examples. All the numerical results show that Implicit SPH with this new boundary-treatment method can obtain accurate results for non-Newtonian fluids as well as Newtonian fluids, and this method is suitable for complex solid boundaries and can be easily extended to 3D problems.展开更多
In consideration of the electroosmotic flow in a slit microchannel, the con-stitutive relationship of the Eyring fluid model is utilized. Navier's slip condition is used as the boundary condition. The governing equat...In consideration of the electroosmotic flow in a slit microchannel, the con-stitutive relationship of the Eyring fluid model is utilized. Navier's slip condition is used as the boundary condition. The governing equations are solved analytically, yielding the velocity distribution. The approximate expressions of the velocity distribution are also given and discussed. Furthermore, the effects of the dimensionless parameters, the electrokinetic parameter, and the slip length on the flow are studied numerically, and appropriate conclusions are drawn.展开更多
Traditional 3D printing is based on stereolithography or standard tessellation language models,which contain many redundant data and have low precision.This paper proposes a slicing and support structure generation al...Traditional 3D printing is based on stereolithography or standard tessellation language models,which contain many redundant data and have low precision.This paper proposes a slicing and support structure generation algorithm for 3D printing directly on boundary representation(B-rep)models.First,surface slicing is performed by efficiently computing the intersection curves between the faces of the B-rep models and each slicing plane.Then,the normals of the B-rep models are used to detect where the support structures should be located and the support structures are generated.Experimental results show the efficiency and stability of our algorithm.展开更多
This work presents the application of the recently developed “Fifth-Order Comprehensive Adjoint Sensitivity Analysis Methodology for Nonlinear Systems (5<sup>th</sup>-CASAM-N)” to a simplified Bernoulli ...This work presents the application of the recently developed “Fifth-Order Comprehensive Adjoint Sensitivity Analysis Methodology for Nonlinear Systems (5<sup>th</sup>-CASAM-N)” to a simplified Bernoulli model. The 5<sup>th</sup>-CASAM-N builds upon and incorporates all of the lower-order (i.e., the first-, second-, third-, and fourth-order) adjoint sensitivities analysis methodologies. The Bernoulli model comprises a nonlinear model response, uncertain model parameters, uncertain model domain boundaries and uncertain model boundary conditions, admitting closed-form explicit expressions for the response sensitivities of all orders. Illustrating the specific mechanisms and advantages of applying the 5<sup>th</sup>-CASAM-N for the computation of the response sensitivities with respect to the uncertain parameters and boundaries reveals that the 5<sup>th</sup>-CASAM-N provides a fundamental step towards overcoming the curse of dimensionality in sensitivity and uncertainty analysis.展开更多
In this paper,the evaluations of metal ablation processes under high temperature,i.e.,the Al plate ablated by a laser and a heat carrier and the reactor pressure vessel ablated by a core melt,are studied by a novel pe...In this paper,the evaluations of metal ablation processes under high temperature,i.e.,the Al plate ablated by a laser and a heat carrier and the reactor pressure vessel ablated by a core melt,are studied by a novel peridynamic method.Above all,the peridynamic formulation for the heat conduction problem is obtained by Taylor’s expansion technique.Then,a simple and efficient moving boundary model in the peridynamic framework is proposed to handle the variable geometries,in which the ablated states of material points are described by an additional scalar field.Next,due to the automatic non-interpenetration properties of peridynamic method,a contact algorithm is established to determine the contact relationship between the ablated system and the additional heat carrier.In addition,the corresponding computational procedure is listed in detail.Finally,several numerical examples are carried out and the results verify the validity and accuracy of the present method.展开更多
A boundary layer model was developed to predict the capture of inclusions by steel-slag interface in a turbulent fluid flow,which is based on the detailed analysis of inclusion trajectories.The effective boundary laye...A boundary layer model was developed to predict the capture of inclusions by steel-slag interface in a turbulent fluid flow,which is based on the detailed analysis of inclusion trajectories.The effective boundary layer for inclusion removal was proposed by a statistical method.It is noticed that the capture of inclusions by steel-slag interface is not only dependent on the diameter of inclusions but also related to the local turbulent conditions.In high turbulent flow fields,the transport of inclusions is mainly dominated by the turbulent flow,and thus,the effective boundary layer thickness is mainly affected by the level of turbulent kinetic energy and is almost independent of the inclusion diameter.The inertia of inclusions gradually takes over the stochastic effect of turbulent flow,and the effect of inclusion diameter on effective boundary layer thickness becomes more noticeable with the decrease in the level of turbulent kinetic energy.Besides,the effective boundary layer thickness is more susceptible to the inclusion diameter for larger inclusions due to its greater inertia under the same turbulent condition while it principally depends on the level of turbulent kinetic energy for smaller inclusions.As the characteristic velocity increases,the time for inclusions transport and interaction with steel-slag interface decreases,and thus,the effective boundary layer thickness decreases.Moreover,the graphical user interface was developed by using the cubic spline interpolation for ease of coupling the current boundary layer model with the macro-scale model of a turbulent fluid flow in the metallurgical vessel.展开更多
The application of grooved surface structure is an emerging and effective means in turbulence flow control.However,for a realistic configuration,the global flow field described directly by simple application of massiv...The application of grooved surface structure is an emerging and effective means in turbulence flow control.However,for a realistic configuration,the global flow field described directly by simple application of massive grids makes it unfeasible to simulate.In this paper,a boundary surrogate model reproducing the effect of microscopic near-wall region is proposed to improve computational efficiency.The surrogate model trained with Lattice Boltzmann Method(LBM)considering the rarefied effect based on real micro/nanoflow field is new among literature,which accurately shows flow characteristics of the micro/nano structure.With this approach,numerical simulations via Reynolds-averaged Navier Stokes equations with modified wall boundary condition are performed in subsonic and transonic flow.The results show that micro/nano grooved surface structure has the effect of delaying transition from laminar to turbulence,thus reducing the skin friction significantly.Analysis of turbulence intensity and turbulence kinetic energy shows that the near-wall flow field of grooved airfoil is more stable compared with that of the smooth airfoil.The reducing rate of maximum turbulent intensity reaches 13.39%.The paper shows a perspective for further application of micro/nano groove structure to turbulence flow control in aircraft design by providing an accurate and efficient simulation method.展开更多
An urban boundary layer model (UBLM) is improved by incorporating the effect of buildings with a sectional drag coefficient and a height-distributed canopy drag length scale. The improved UBLM is applied to simulate...An urban boundary layer model (UBLM) is improved by incorporating the effect of buildings with a sectional drag coefficient and a height-distributed canopy drag length scale. The improved UBLM is applied to simulate the wind fields over three typical urban blocks over the Beijing area with different height-towidth ratios. For comparisons, the wind fields over the same blocks are simulated by an urban sub-domain scale model resolving the buildings explicitly. The wind fields simulated from the two different methods are in good agreement. Then, two-dimensional building morphological characteristics and urban canopy parameters for Beijing are derived from detailed building height data. Finally, experiements are conducted to investigate the effect of buildings on the wind field in Beijing using the improved UBLM.展开更多
In fractured reservoirs characterized by low matrix permeability,fracture networks control the main fluid flow paths.However,in layered reservoirs,the vertical extension of fractures is often restricted to single laye...In fractured reservoirs characterized by low matrix permeability,fracture networks control the main fluid flow paths.However,in layered reservoirs,the vertical extension of fractures is often restricted to single layers.In this study,we explored the effect of changing marl/shale thickness on fracture extension using comprehensive field data and numerical modeling.The field data were sampled from coastal exposures of Liassic limestone-marl/shale alternations in Wales and Somerset(Bristol Channel Basin,UK).The vertical fracture traces of more than 4000 fractures were mapped in detail.Six sections were selected to represent a variety of layer thicknesses.Besides the field data also thin sections were analyzed.Numerical models of fracture extension in a two-layer limestone-marl system were based on field data and laboratory measurements of Young's moduli.The modeled principal stress magnitude σ3 along the lithological contact was used as an indication for fracture extension through marls.Field data exhibit good correlation(R^2=0.76) between fracture extension and marl thickness,the thicker the marl layer the fewer fractures propagate through.The model results show that almost no tensile stress reaches the top of the marl layer when the marls are thicker than 30 cm.For marls that are less than 20 cm,the propagation of stress is more dependent on the stiffness of the marls.The higher the contrast between limestone and marl stiffness the lower the stress that is transmitted into the marl layer.In both model experiments and field data the critical marl thickness for fracture extension is ca.15-20 cm.This quantification of critical marl thicknesses can be used to improve predictions of fracture networks and permeability in layered rocks.Up-or downsampling methods often ignore spatially continuous impermeable layers with thicknesses that are under the detection limit of seismic data.However,ignoring these layers can lead to overestimates of the overall permeability.Therefore,the understanding of how fractures propagate and terminate through impermeable layers will help to improve the characterization of conventional reservoirs.展开更多
The fluid-structure coupling finite element model and the boundary element model of a complex gearbox's housing are built based on the theory of fluid-structure coupling finite element method and boundary element met...The fluid-structure coupling finite element model and the boundary element model of a complex gearbox's housing are built based on the theory of fluid-structure coupling finite element method and boundary element method. At the same time, the exciting forces of the housing are analyzed and applied to the finite element models. Firstly, vibration of the housing is calculated by the fluid-structure coupling finite element model; secondly, the calculated result is verified by the experiment; finally, sound radiation of the housing is calculated by the boundary element. According to the calculated results, the housing adds some ribs not only to increase the strength, but also to reduce the sound radiation of the housing. At last, the sound radiation of the modified housing is calculated, which shows that the sound radiation of the modified housing with ribs is lower.展开更多
A characteristic like model is presented to calculate the two dimensional flow in open channels with finite or infinite dimensions. The comparisons between the computational results and experimental data show the mo...A characteristic like model is presented to calculate the two dimensional flow in open channels with finite or infinite dimensions. The comparisons between the computational results and experimental data show the model to be robust in simulating two dimensional unsteady flows. The comparisons also show the effectiveness of the open boundary condition in simulation of wave propagation process in an infinite domain with computational domain of finite dimension.展开更多
This work presents the “Second-Order Comprehensive Adjoint Sensitivity Analysis Methodology (2<sup>nd</sup>-CASAM)” for the efficient and exact computation of 1<sup>st</sup>- and 2<sup>...This work presents the “Second-Order Comprehensive Adjoint Sensitivity Analysis Methodology (2<sup>nd</sup>-CASAM)” for the efficient and exact computation of 1<sup>st</sup>- and 2<sup>nd</sup>-order response sensitivities to uncertain parameters and domain boundaries of linear systems. The model’s response (<em>i.e.</em>, model result of interest) is a generic nonlinear function of the model’s forward and adjoint state functions, and also depends on the imprecisely known boundaries and model parameters. In the practically important particular case when the response is a scalar-valued functional of the forward and adjoint state functions characterizing a model comprising N parameters, the 2<sup>nd</sup>-CASAM requires a single large-scale computation using the First-Level Adjoint Sensitivity System (1<sup>st</sup>-LASS) for obtaining all of the first-order response sensitivities, and at most N large-scale computations using the Second-Level Adjoint Sensitivity System (2<sup>nd</sup>-LASS) for obtaining exactly all of the second-order response sensitivities. In contradistinction, forward other methods would require (<em>N</em>2/2 + 3 <em>N</em>/2) large-scale computations for obtaining all of the first- and second-order sensitivities. This work also shows that constructing and solving the 2<sup>nd</sup>-LASS requires very little additional effort beyond the construction of the 1<sup>st</sup>-LASS needed for computing the first-order sensitivities. Solving the equations underlying the 1<sup>st</sup>-LASS and 2<sup>nd</sup>-LASS requires the same computational solvers as needed for solving (<em>i.e.</em>, “inverting”) either the forward or the adjoint linear operators underlying the initial model. Therefore, the same computer software and “solvers” used for solving the original system of equations can also be used for solving the 1<sup>st</sup>-LASS and the 2<sup>nd</sup>-LASS. Since neither the 1<sup>st</sup>-LASS nor the 2<sup>nd</sup>-LASS involves any differentials of the operators underlying the original system, the 1<sup>st</sup>-LASS is designated as a “<u>first-level</u>” (as opposed to a “first-order”) adjoint sensitivity system, while the 2<sup>nd</sup>-LASS is designated as a “<u>second-level</u>” (rather than a “second-order”) adjoint sensitivity system. Mixed second-order response sensitivities involving boundary parameters may arise from all source terms of the 2<sup>nd</sup>-LASS that involve the imprecisely known boundary parameters. Notably, the 2<sup>nd</sup>-LASS encompasses an automatic, inherent, and independent “solution verification” mechanism of the correctness and accuracy of the 2nd-level adjoint functions needed for the efficient and exact computation of the second-order sensitivities.展开更多
Although the fracture behavior of sea urchin spines has been extensively investigated,there is as yet a lack of quantitative estimation on the effect of growth rings on the fracture properties of sea urchin spines.In ...Although the fracture behavior of sea urchin spines has been extensively investigated,there is as yet a lack of quantitative estimation on the effect of growth rings on the fracture properties of sea urchin spines.In sea urchin spines,much denser pores present in growth rings rather than porous layers.The tensile strength and fracture toughness of sea urchin spine samples with different numbers of growth rings are measured by the Boundary Effect Model(BEM).The experimental results of single-edge notched three-point bending tests indicate that the BEM is an appropriate method to estimate the fracture toughness of the present porous sea urchin spines,and the number of growth rings plays an important role in the fracture properties of spines.Specifically,the tensile strength and fracture toughness of sea urchin spines can be significantly improved with the increase in the number of growth rings,and their fracture toughness can even reach a relatively high value compared with some other porous materials with an identical porosity.The present research findings are expected to provide a fundamental insight into the design of high-performance bionic materials with a highly porous structure.展开更多
A one dimensional model is developed for defective gap mode(DGM)with two types of boundary conditions:conducting mesh and conducting sleeve.For a periodically modulated system without defect,the normalized width of...A one dimensional model is developed for defective gap mode(DGM)with two types of boundary conditions:conducting mesh and conducting sleeve.For a periodically modulated system without defect,the normalized width of spectral gaps equals to the modulation factor,which is consistent with previous studies.For a periodic system with local defects introduced by the boundary conditions,it shows that the conducting-mesh-induced DGM is always well confined by spectral gaps while the conducting-sleeve-induced DGM is not.The defect location can be a useful tool to dynamically control the frequency and spatial periodicity of DGM inside spectral gaps.This controllability can be potentially applied to the interaction between gap eigenmodes and energetic particles in fusion plasmas,and optical microcavities and waveguides in photonic crystals.展开更多
A highly efficient absorbing boundary condition suitable for use in the finitedtherence time-domain (FDTD) modelling of acoustic fields is presented in this paper. The new method seeks a least square esthoate of a tra...A highly efficient absorbing boundary condition suitable for use in the finitedtherence time-domain (FDTD) modelling of acoustic fields is presented in this paper. The new method seeks a least square esthoate of a transfer matrix for field components near truncating boundaries by matrir pseud-inversion. The proposed absorbing boundary is considerably more effective than most ekisting ones. The method is also computationally econondcal and robust.The performance of the new method is shown by numerical experiments on a point-source radiation problem, a wedge dimaction problem, and a scattering problem in which a plane wave is scattered by a circular cylinder.展开更多
Using a modified 3D random representative volume(RV)finite element model,the effects of model dimensions(impact region and interval between impact and representative regions),model shapes(rectangular,square,and c...Using a modified 3D random representative volume(RV)finite element model,the effects of model dimensions(impact region and interval between impact and representative regions),model shapes(rectangular,square,and circular),and peening-induced thermal softening on resultant critical quantities(residual stress,Almen intensity,coverage,and arc height)after shot peening are systematically examined.A new quantity,i.e.,the interval between impact and representative regions,is introduced and its optimal value is first determined to eliminate any boundary effect on shot peening results.Then,model dimensions are respectively assessed for all model shapes to reflect the actual shot peening process,based on which shape-independent critical shot peening quantities are obtained.Further,it is found that thermal softening of the target material due to shot peening leads to variances of the surface residual stress and arc height,demonstrating the necessity of considering the thermal effect in a constitutive material model of shot peeing.Our study clarifies some of the finite element modeling aspects and lays the ground for accurate modeling of the SP process.展开更多
The causes and variability of a heavy haze episode in the Beijing region was analyzed. During the episode, the PM2.5 concentration reached a peak value of 450 μg/kg on January 18, 2013 and rapidly decreased to 100μg...The causes and variability of a heavy haze episode in the Beijing region was analyzed. During the episode, the PM2.5 concentration reached a peak value of 450 μg/kg on January 18, 2013 and rapidly decreased to 100μg/kg on January 19, 2013, characterizing a large variability in a very short period. This strong vari- ability provides a good opportunity to study the causes of the haze formation. The in situ measurements (including surface meteorological data and vertical structures of the winds, temperature, humidity, and planetary boundary layer (PBL)) together with a chemical/dynamical regional model (WRF-Chem) were used for the analysis. In order to understand the rapid variability of the PM2.5 concentration in the episode, the correlation between the measured meteorological data (including wind speed, PBL height, relative humidity, etc.) and the measured particle concentration (PM2.5 concentration) was studied. In addition, two sensitive model experiments were performed to study the effect of individual contribution from local emissions and regional surrounding emissions to the heavy haze formation. The results suggest that there were two major meteorological factors in controlling the variability of the PM2.5 concentration, namely, surface wind speed and PBL height. During high wind periods, the horizontal transport of aerosol particles played an important role, and the heavy haze was formed when the wind speeds were very weak (less than 1 m/s). Under weak wind conditions, the horizontal transport of aerosol particles was also weak, and the vertical mixing of aerosol particles played an important role. As a result, the PBL height was a major factor in controlling the variability of the PM2.5 concentration. Under the shallow PBL height, aerosol particles were strongly confined near the surface, producing a high surface PM2.5 concentration. The sensitivity model study suggests that the local emissions (emissions from the Beijing region only) were the major cause for the heavy haze events. With only local emissions, the calculated peak value of the PM2.s concentration was 350μg/kg, which accounted for 78% of the measured peak value (450 μg/kg). In contrast, without the local emissions, the calculated peak value of the PM2.5 concentration was only 100 μg/kg, which accounted for 22% of the measured peak value.展开更多
A coupled numerical method for the direct numerical simulation of particle-fluid systems is formulated and implemented, resolving an order of magnitude smaller than particle size. The particle motion is described by t...A coupled numerical method for the direct numerical simulation of particle-fluid systems is formulated and implemented, resolving an order of magnitude smaller than particle size. The particle motion is described by the time-driven hard-sphere model, while the hydrodynamic equations governing fluid flow are solved by the lattice Boltzmann method (LBM), Particle-fluid coupling is realized by an immersed boundary method (IBM), which considers the effect of boundary on surrounding fluid as a restoring force added to the governing equations of the fluid. The proposed scheme is validated in the classical flow-around-cylinder simulations, and preliminary application of this scheme to fluidization is reported, demonstrating it to be a promising computational strategy for better understanding complex behavior in particle-fluid systems.展开更多
基金supported by the National Natural Science Foundation of China (No. 41130418)the Strategic Leading Science and Technology Programme (Class B) of the Chinese Academy of Sciences (No. XDB10010400)
文摘Earth medium is not completely elastic, with its viscosity resulting in attenuation and dispersion of seismic waves. Most viscoelastic numerical simulations are based on the finite-difference and finite-element methods. Targeted at viscoelastic numerical modeling for multilayered media, the constant-Q acoustic wave equation is transformed into the corresponding wave integral representation with its Green's function accounting for viscoelastic coefficients. An efficient alternative for full-waveform solution to the integral equation is proposed in this article by extending conventional frequency-domain boundary element methods to viscoelastic media. The viscoelastic boundary element method enjoys a distinct characteristic of the explicit use of boundary continuity conditions of displacement and traction, leading to a semi-analytical solution with sufficient accuracy for simulating the viscoelastic effect across irregular interfaces. Numerical experiments to study the viscoelastic absorption of different Q values demonstrate the accuracy and applicability of the method.
基金The National Natural Science Foundation of China under contract Nos 51609244 and 51779258
文摘An analytical model with essential parameters given by a two-phase numerical model is utilized to study the net boundary layer current and sediment transport under skewed asymmetric oscillatory sheet flows. The analytical model is the first instantaneous type model that can consider phase-lag and asymmetric boundary layer development. The two-phase model supplies the essential phase-lead, instantaneous erosion depth and boundary layer development for the analytical model to enhance the understanding of velocity skewness and acceleration skewness in sediment flux and transport rate. The sediment transport difference between onshore and offshore stages caused by velocity skewness or acceleration skewness is shown to illustrate the determination of net sediment transport by the analytical model. In previous studies about sediment transport in skewed asymmetric sheet flows, the generation of net sediment transport is mainly concluded to the phase-lag effect.However, the phase-lag effect is shown important but not enough for the net sediment transport, while the skewed asymmetric boundary layer development generated net boundary layer current and mobile bed effect are key important in the transport process.
基金supported by the National Natural Science Foundation of China(51276192)
文摘In this work, an enhanced treatment of the solid boundaries is proposed for smoothed particle hydrodynamics with implicit time integration scheme (Implicit SPH). Three types of virtual particles, i.e., boundary particles, image particles and mirror particles, are used to impose boundary conditions. Boundary particles are fixed on the solid boundary, and each boundary particle is associated with two fixed image particles inside the fluid domain and two fixed mirror particles outside the fluid domain. The image particles take the flow properties through fluid particles with moving least squares (MLS) interpolation and the properties of mirror particles can be obtained by the corresponding image particles. A repulsive force is also applied for boundary particles to prevent fluid particles from unphysical penetra- tion through solid boundaries. The new boundary treatment method has been validated with five numerical examples. All the numerical results show that Implicit SPH with this new boundary-treatment method can obtain accurate results for non-Newtonian fluids as well as Newtonian fluids, and this method is suitable for complex solid boundaries and can be easily extended to 3D problems.
基金Project supported by the National Natural Science Foundation of China(Nos.11102102 and 91130017)the Independent Innovation Foundation of Shandong University(No.2013ZRYQ002)
文摘In consideration of the electroosmotic flow in a slit microchannel, the con-stitutive relationship of the Eyring fluid model is utilized. Navier's slip condition is used as the boundary condition. The governing equations are solved analytically, yielding the velocity distribution. The approximate expressions of the velocity distribution are also given and discussed. Furthermore, the effects of the dimensionless parameters, the electrokinetic parameter, and the slip length on the flow are studied numerically, and appropriate conclusions are drawn.
基金This work is was funded by National Natural Science Foundation of China under Grant No.61672307.
文摘Traditional 3D printing is based on stereolithography or standard tessellation language models,which contain many redundant data and have low precision.This paper proposes a slicing and support structure generation algorithm for 3D printing directly on boundary representation(B-rep)models.First,surface slicing is performed by efficiently computing the intersection curves between the faces of the B-rep models and each slicing plane.Then,the normals of the B-rep models are used to detect where the support structures should be located and the support structures are generated.Experimental results show the efficiency and stability of our algorithm.
文摘This work presents the application of the recently developed “Fifth-Order Comprehensive Adjoint Sensitivity Analysis Methodology for Nonlinear Systems (5<sup>th</sup>-CASAM-N)” to a simplified Bernoulli model. The 5<sup>th</sup>-CASAM-N builds upon and incorporates all of the lower-order (i.e., the first-, second-, third-, and fourth-order) adjoint sensitivities analysis methodologies. The Bernoulli model comprises a nonlinear model response, uncertain model parameters, uncertain model domain boundaries and uncertain model boundary conditions, admitting closed-form explicit expressions for the response sensitivities of all orders. Illustrating the specific mechanisms and advantages of applying the 5<sup>th</sup>-CASAM-N for the computation of the response sensitivities with respect to the uncertain parameters and boundaries reveals that the 5<sup>th</sup>-CASAM-N provides a fundamental step towards overcoming the curse of dimensionality in sensitivity and uncertainty analysis.
基金supported by the National Natural Science Foundation of China(No.12102416).
文摘In this paper,the evaluations of metal ablation processes under high temperature,i.e.,the Al plate ablated by a laser and a heat carrier and the reactor pressure vessel ablated by a core melt,are studied by a novel peridynamic method.Above all,the peridynamic formulation for the heat conduction problem is obtained by Taylor’s expansion technique.Then,a simple and efficient moving boundary model in the peridynamic framework is proposed to handle the variable geometries,in which the ablated states of material points are described by an additional scalar field.Next,due to the automatic non-interpenetration properties of peridynamic method,a contact algorithm is established to determine the contact relationship between the ablated system and the additional heat carrier.In addition,the corresponding computational procedure is listed in detail.Finally,several numerical examples are carried out and the results verify the validity and accuracy of the present method.
基金the National Natural Science Foundation of China(Grant Nos.51904025 and U22A20171)the Fundamental Research Funds for the Central Universities(Grant No.FRF-IDRY-20-011)+1 种基金National Postdoctoral Program for Innovative Talents(Grant No.BX20190030)the High Steel Center(HSC)at North China University of Technology and University of Science and Technology Beijing,China.
文摘A boundary layer model was developed to predict the capture of inclusions by steel-slag interface in a turbulent fluid flow,which is based on the detailed analysis of inclusion trajectories.The effective boundary layer for inclusion removal was proposed by a statistical method.It is noticed that the capture of inclusions by steel-slag interface is not only dependent on the diameter of inclusions but also related to the local turbulent conditions.In high turbulent flow fields,the transport of inclusions is mainly dominated by the turbulent flow,and thus,the effective boundary layer thickness is mainly affected by the level of turbulent kinetic energy and is almost independent of the inclusion diameter.The inertia of inclusions gradually takes over the stochastic effect of turbulent flow,and the effect of inclusion diameter on effective boundary layer thickness becomes more noticeable with the decrease in the level of turbulent kinetic energy.Besides,the effective boundary layer thickness is more susceptible to the inclusion diameter for larger inclusions due to its greater inertia under the same turbulent condition while it principally depends on the level of turbulent kinetic energy for smaller inclusions.As the characteristic velocity increases,the time for inclusions transport and interaction with steel-slag interface decreases,and thus,the effective boundary layer thickness decreases.Moreover,the graphical user interface was developed by using the cubic spline interpolation for ease of coupling the current boundary layer model with the macro-scale model of a turbulent fluid flow in the metallurgical vessel.
基金funded by AECC Commercial Aircraft Engine in the project of Laminar Flow Design and Turbulent Drag Reduction of Compressor Blade Profile and Experimental Verification。
文摘The application of grooved surface structure is an emerging and effective means in turbulence flow control.However,for a realistic configuration,the global flow field described directly by simple application of massive grids makes it unfeasible to simulate.In this paper,a boundary surrogate model reproducing the effect of microscopic near-wall region is proposed to improve computational efficiency.The surrogate model trained with Lattice Boltzmann Method(LBM)considering the rarefied effect based on real micro/nanoflow field is new among literature,which accurately shows flow characteristics of the micro/nano structure.With this approach,numerical simulations via Reynolds-averaged Navier Stokes equations with modified wall boundary condition are performed in subsonic and transonic flow.The results show that micro/nano grooved surface structure has the effect of delaying transition from laminar to turbulence,thus reducing the skin friction significantly.Analysis of turbulence intensity and turbulence kinetic energy shows that the near-wall flow field of grooved airfoil is more stable compared with that of the smooth airfoil.The reducing rate of maximum turbulent intensity reaches 13.39%.The paper shows a perspective for further application of micro/nano groove structure to turbulence flow control in aircraft design by providing an accurate and efficient simulation method.
基金funded by National Nat-ural Science Foundation of China(Grants Nos.40505002,40652001,and 40775015)Beijing Natural Science Foun-dation(Grant No.8051002)+1 种基金Beijing New Star Project of Science and Technology(Grant No.2005A03)the Ministry of Science and Technology of China(Grant Nos.2008BAC37B04,2006BAJ02A01,and GYHY200906035)
文摘An urban boundary layer model (UBLM) is improved by incorporating the effect of buildings with a sectional drag coefficient and a height-distributed canopy drag length scale. The improved UBLM is applied to simulate the wind fields over three typical urban blocks over the Beijing area with different height-towidth ratios. For comparisons, the wind fields over the same blocks are simulated by an urban sub-domain scale model resolving the buildings explicitly. The wind fields simulated from the two different methods are in good agreement. Then, two-dimensional building morphological characteristics and urban canopy parameters for Beijing are derived from detailed building height data. Finally, experiements are conducted to investigate the effect of buildings on the wind field in Beijing using the improved UBLM.
基金supported by the Deutsche Forschungsgemeinschaft(DFG,grant PH 189/2-1)
文摘In fractured reservoirs characterized by low matrix permeability,fracture networks control the main fluid flow paths.However,in layered reservoirs,the vertical extension of fractures is often restricted to single layers.In this study,we explored the effect of changing marl/shale thickness on fracture extension using comprehensive field data and numerical modeling.The field data were sampled from coastal exposures of Liassic limestone-marl/shale alternations in Wales and Somerset(Bristol Channel Basin,UK).The vertical fracture traces of more than 4000 fractures were mapped in detail.Six sections were selected to represent a variety of layer thicknesses.Besides the field data also thin sections were analyzed.Numerical models of fracture extension in a two-layer limestone-marl system were based on field data and laboratory measurements of Young's moduli.The modeled principal stress magnitude σ3 along the lithological contact was used as an indication for fracture extension through marls.Field data exhibit good correlation(R^2=0.76) between fracture extension and marl thickness,the thicker the marl layer the fewer fractures propagate through.The model results show that almost no tensile stress reaches the top of the marl layer when the marls are thicker than 30 cm.For marls that are less than 20 cm,the propagation of stress is more dependent on the stiffness of the marls.The higher the contrast between limestone and marl stiffness the lower the stress that is transmitted into the marl layer.In both model experiments and field data the critical marl thickness for fracture extension is ca.15-20 cm.This quantification of critical marl thicknesses can be used to improve predictions of fracture networks and permeability in layered rocks.Up-or downsampling methods often ignore spatially continuous impermeable layers with thicknesses that are under the detection limit of seismic data.However,ignoring these layers can lead to overestimates of the overall permeability.Therefore,the understanding of how fractures propagate and terminate through impermeable layers will help to improve the characterization of conventional reservoirs.
文摘The fluid-structure coupling finite element model and the boundary element model of a complex gearbox's housing are built based on the theory of fluid-structure coupling finite element method and boundary element method. At the same time, the exciting forces of the housing are analyzed and applied to the finite element models. Firstly, vibration of the housing is calculated by the fluid-structure coupling finite element model; secondly, the calculated result is verified by the experiment; finally, sound radiation of the housing is calculated by the boundary element. According to the calculated results, the housing adds some ribs not only to increase the strength, but also to reduce the sound radiation of the housing. At last, the sound radiation of the modified housing is calculated, which shows that the sound radiation of the modified housing with ribs is lower.
文摘A characteristic like model is presented to calculate the two dimensional flow in open channels with finite or infinite dimensions. The comparisons between the computational results and experimental data show the model to be robust in simulating two dimensional unsteady flows. The comparisons also show the effectiveness of the open boundary condition in simulation of wave propagation process in an infinite domain with computational domain of finite dimension.
文摘This work presents the “Second-Order Comprehensive Adjoint Sensitivity Analysis Methodology (2<sup>nd</sup>-CASAM)” for the efficient and exact computation of 1<sup>st</sup>- and 2<sup>nd</sup>-order response sensitivities to uncertain parameters and domain boundaries of linear systems. The model’s response (<em>i.e.</em>, model result of interest) is a generic nonlinear function of the model’s forward and adjoint state functions, and also depends on the imprecisely known boundaries and model parameters. In the practically important particular case when the response is a scalar-valued functional of the forward and adjoint state functions characterizing a model comprising N parameters, the 2<sup>nd</sup>-CASAM requires a single large-scale computation using the First-Level Adjoint Sensitivity System (1<sup>st</sup>-LASS) for obtaining all of the first-order response sensitivities, and at most N large-scale computations using the Second-Level Adjoint Sensitivity System (2<sup>nd</sup>-LASS) for obtaining exactly all of the second-order response sensitivities. In contradistinction, forward other methods would require (<em>N</em>2/2 + 3 <em>N</em>/2) large-scale computations for obtaining all of the first- and second-order sensitivities. This work also shows that constructing and solving the 2<sup>nd</sup>-LASS requires very little additional effort beyond the construction of the 1<sup>st</sup>-LASS needed for computing the first-order sensitivities. Solving the equations underlying the 1<sup>st</sup>-LASS and 2<sup>nd</sup>-LASS requires the same computational solvers as needed for solving (<em>i.e.</em>, “inverting”) either the forward or the adjoint linear operators underlying the initial model. Therefore, the same computer software and “solvers” used for solving the original system of equations can also be used for solving the 1<sup>st</sup>-LASS and the 2<sup>nd</sup>-LASS. Since neither the 1<sup>st</sup>-LASS nor the 2<sup>nd</sup>-LASS involves any differentials of the operators underlying the original system, the 1<sup>st</sup>-LASS is designated as a “<u>first-level</u>” (as opposed to a “first-order”) adjoint sensitivity system, while the 2<sup>nd</sup>-LASS is designated as a “<u>second-level</u>” (rather than a “second-order”) adjoint sensitivity system. Mixed second-order response sensitivities involving boundary parameters may arise from all source terms of the 2<sup>nd</sup>-LASS that involve the imprecisely known boundary parameters. Notably, the 2<sup>nd</sup>-LASS encompasses an automatic, inherent, and independent “solution verification” mechanism of the correctness and accuracy of the 2nd-level adjoint functions needed for the efficient and exact computation of the second-order sensitivities.
基金This work was supported by the National Natural Science Foundation of China[Grant No.51902043]the Fundamental Research Funds for the Central Universities[Grant No.N2102007,and N2102002]。
文摘Although the fracture behavior of sea urchin spines has been extensively investigated,there is as yet a lack of quantitative estimation on the effect of growth rings on the fracture properties of sea urchin spines.In sea urchin spines,much denser pores present in growth rings rather than porous layers.The tensile strength and fracture toughness of sea urchin spine samples with different numbers of growth rings are measured by the Boundary Effect Model(BEM).The experimental results of single-edge notched three-point bending tests indicate that the BEM is an appropriate method to estimate the fracture toughness of the present porous sea urchin spines,and the number of growth rings plays an important role in the fracture properties of spines.Specifically,the tensile strength and fracture toughness of sea urchin spines can be significantly improved with the increase in the number of growth rings,and their fracture toughness can even reach a relatively high value compared with some other porous materials with an identical porosity.The present research findings are expected to provide a fundamental insight into the design of high-performance bionic materials with a highly porous structure.
基金supported by National Natural Science Foundation of China(No.11405271)
文摘A one dimensional model is developed for defective gap mode(DGM)with two types of boundary conditions:conducting mesh and conducting sleeve.For a periodically modulated system without defect,the normalized width of spectral gaps equals to the modulation factor,which is consistent with previous studies.For a periodic system with local defects introduced by the boundary conditions,it shows that the conducting-mesh-induced DGM is always well confined by spectral gaps while the conducting-sleeve-induced DGM is not.The defect location can be a useful tool to dynamically control the frequency and spatial periodicity of DGM inside spectral gaps.This controllability can be potentially applied to the interaction between gap eigenmodes and energetic particles in fusion plasmas,and optical microcavities and waveguides in photonic crystals.
文摘A highly efficient absorbing boundary condition suitable for use in the finitedtherence time-domain (FDTD) modelling of acoustic fields is presented in this paper. The new method seeks a least square esthoate of a transfer matrix for field components near truncating boundaries by matrir pseud-inversion. The proposed absorbing boundary is considerably more effective than most ekisting ones. The method is also computationally econondcal and robust.The performance of the new method is shown by numerical experiments on a point-source radiation problem, a wedge dimaction problem, and a scattering problem in which a plane wave is scattered by a circular cylinder.
基金the financial support from China Scholarship Council (CSC) (No. 201406025083)National Natural Science Foundation of China (NSFC) (Nos. 51305012 and 51675024)+3 种基金Aviation Science Fund of China (No. 2014ZB51)financial support from NSFC (No. 51375031)financial support from NSFC (No. 51628101)National Sciences and Engineering Research Council (NSERC) Discovery grant (No. RGPIN 418469-2012)
文摘Using a modified 3D random representative volume(RV)finite element model,the effects of model dimensions(impact region and interval between impact and representative regions),model shapes(rectangular,square,and circular),and peening-induced thermal softening on resultant critical quantities(residual stress,Almen intensity,coverage,and arc height)after shot peening are systematically examined.A new quantity,i.e.,the interval between impact and representative regions,is introduced and its optimal value is first determined to eliminate any boundary effect on shot peening results.Then,model dimensions are respectively assessed for all model shapes to reflect the actual shot peening process,based on which shape-independent critical shot peening quantities are obtained.Further,it is found that thermal softening of the target material due to shot peening leads to variances of the surface residual stress and arc height,demonstrating the necessity of considering the thermal effect in a constitutive material model of shot peeing.Our study clarifies some of the finite element modeling aspects and lays the ground for accurate modeling of the SP process.
基金supported by National Science Foundation of China(NSFC) under Grant No.41205100,41275186, 41175007National Department Public Benefit Research Foundation of China under Grant No.GYHY201306065The National Center for Atmospheric Research,U.S.is sponsored by the National Science Foundation,U.S
文摘The causes and variability of a heavy haze episode in the Beijing region was analyzed. During the episode, the PM2.5 concentration reached a peak value of 450 μg/kg on January 18, 2013 and rapidly decreased to 100μg/kg on January 19, 2013, characterizing a large variability in a very short period. This strong vari- ability provides a good opportunity to study the causes of the haze formation. The in situ measurements (including surface meteorological data and vertical structures of the winds, temperature, humidity, and planetary boundary layer (PBL)) together with a chemical/dynamical regional model (WRF-Chem) were used for the analysis. In order to understand the rapid variability of the PM2.5 concentration in the episode, the correlation between the measured meteorological data (including wind speed, PBL height, relative humidity, etc.) and the measured particle concentration (PM2.5 concentration) was studied. In addition, two sensitive model experiments were performed to study the effect of individual contribution from local emissions and regional surrounding emissions to the heavy haze formation. The results suggest that there were two major meteorological factors in controlling the variability of the PM2.5 concentration, namely, surface wind speed and PBL height. During high wind periods, the horizontal transport of aerosol particles played an important role, and the heavy haze was formed when the wind speeds were very weak (less than 1 m/s). Under weak wind conditions, the horizontal transport of aerosol particles was also weak, and the vertical mixing of aerosol particles played an important role. As a result, the PBL height was a major factor in controlling the variability of the PM2.5 concentration. Under the shallow PBL height, aerosol particles were strongly confined near the surface, producing a high surface PM2.5 concentration. The sensitivity model study suggests that the local emissions (emissions from the Beijing region only) were the major cause for the heavy haze events. With only local emissions, the calculated peak value of the PM2.s concentration was 350μg/kg, which accounted for 78% of the measured peak value (450 μg/kg). In contrast, without the local emissions, the calculated peak value of the PM2.5 concentration was only 100 μg/kg, which accounted for 22% of the measured peak value.
基金sponsored by Ministry of Finance under the grant ZDYZ2008-2National Key Science and Technology Project under the grant 2008ZX05014-003-006HZthe Chinese Academy of Sciences under the grant KGCX2-YW-124
文摘A coupled numerical method for the direct numerical simulation of particle-fluid systems is formulated and implemented, resolving an order of magnitude smaller than particle size. The particle motion is described by the time-driven hard-sphere model, while the hydrodynamic equations governing fluid flow are solved by the lattice Boltzmann method (LBM), Particle-fluid coupling is realized by an immersed boundary method (IBM), which considers the effect of boundary on surrounding fluid as a restoring force added to the governing equations of the fluid. The proposed scheme is validated in the classical flow-around-cylinder simulations, and preliminary application of this scheme to fluidization is reported, demonstrating it to be a promising computational strategy for better understanding complex behavior in particle-fluid systems.