This paper presents a method to characterize asphalt pavement macrotexture using the gray-tone difference matrix (GTDM)and discusses the potentials of the GTDM indicators for skid resistance evaluation.There are 37 ...This paper presents a method to characterize asphalt pavement macrotexture using the gray-tone difference matrix (GTDM)and discusses the potentials of the GTDM indicators for skid resistance evaluation.There are 37 field sites included in the data collection,which cover 6 types of asphalt pavement surfaces. The mean profile depth derived from 3-D macrotexture measurements (MPD3 ) has a significant relationship with the mean texture depth (MTD ),which can be described by a logarithm model with R2 of 0.962.There is no significant linear relationship between the friction coefficient at a speed of 60 km/h (DFT60 )and macrotexture indicators.A nonlinear model with British pendulum number (BPN ) incorporated can relate DFT60 to MTD or indicator fcon .A comparison with MTD shows that GTDM-based fcon has a potential to be a macrotexture indicator for skid resistance evaluation,which describes the general height difference and the average local height difference of pavement macrotexture. A relatively high fcon is helpful for improving asphalt pavement skid resistance.展开更多
A family of high_order accuracy explicit difference schemes for solving 2_dimension parabolic P.D.E. are constructed. Th e stability condition is r=Δt/Δx 2=Δt/Δy 2【1/2 and the truncation err or is O(Δt 3+Δx...A family of high_order accuracy explicit difference schemes for solving 2_dimension parabolic P.D.E. are constructed. Th e stability condition is r=Δt/Δx 2=Δt/Δy 2【1/2 and the truncation err or is O(Δt 3+Δx 4).展开更多
The error equation of a rotating inertial navigation system was introduced. The effect of the system's main error source (constant drift of gyro and zero bias of accelerometer) under rotating conditions for the sy...The error equation of a rotating inertial navigation system was introduced. The effect of the system's main error source (constant drift of gyro and zero bias of accelerometer) under rotating conditions for the system was analyzed. Validity of theoretical analysis was shown via simulation, and that provides a theoretical foundation for a rotating strap-down inertial navigation system during actual experimentation and application.展开更多
A modified one-dimensional transient hygrothermal model for multilayer wall was proposed using air humidity ratio and temperature as the driving potentials.The solution for the governing equations was obtained numeric...A modified one-dimensional transient hygrothermal model for multilayer wall was proposed using air humidity ratio and temperature as the driving potentials.The solution for the governing equations was obtained numerically by implementing the finite-difference scheme.To evaluate the accuracy of the model,a test system was built up to measure relative humidity and temperature within a porous wall and compare with the prediction of the model.The prediction results have good agreement with the experimental results.For the interface close to indoor side,the maximum deviation of temperature between calculated and test data is 1.87 K,and the average deviation is 0.95 K;the maximum deviation of relative humidity is 11.4%,and the average deviation is 5.7%.For the interface close to outdoor side,the maximum deviation of temperature between prediction and measurement is 1.78 K,and the average deviation is 1.1 K;the maximum deviation of relative humidity is 9.9%,and the average deviation is 4.2%.展开更多
We present a highly efficient lattice Boltzmann model for simulating compressible flows. This model is based on the combination of an appropriate finite difference scheme, a 16-discrete-velocity model [Kataoka and Tsu...We present a highly efficient lattice Boltzmann model for simulating compressible flows. This model is based on the combination of an appropriate finite difference scheme, a 16-discrete-velocity model [Kataoka and Tsutahara, Phys. Rev. E 69 (2004) 035701(R)] and reasonable dispersion and dissipation terms. The dispersion term effectively reduces the oscillation at the discontinuity and enhances numerical precision. The dissipation term makes the new model more easily meet with the yon Neumann stability condition. This model works for both high-speed and low-speed flows with arbitrary specific-heat-ratio. With the new model simulation results for the well-known benchmark problems get a high accuracy compared with the analytic or experimental ones. The used benchmark tests include (i) Shock tubes such as the Sod, Lax, Sjogreen, Colella explosion wave, and collision of two strong shocks, (ii) Regular and Mach shock reflections, and (iii) Shock wave reaction on cylindrical bubble problems. With a more realistic equation of state or free-energy functional, the new model has the potential tostudy the complex procedure of shock wave reaction on porous materials.展开更多
A three dimensional numerical model of nonlinear wave action on a quasi-ellipse caisson in a time domain was developed in this paper. Navier-Stokes equations were solved by the finite difference method, and the volume...A three dimensional numerical model of nonlinear wave action on a quasi-ellipse caisson in a time domain was developed in this paper. Navier-Stokes equations were solved by the finite difference method, and the volume of fluid (VOF) method was employed to trace the free surface. The partial cell method was used to deal with the irregular boundary typical of this type of problem during first-time wave interaction with the structure, and a satisfactory result was obtained. The numerical model was verified and used to investigate the effects of the relative wave height H/d, relative caisson width kD, and relative length-width ratio B/D on the wave forces of the quasi-ellipse caisson. It was shown that the relative wave height H/d has a significant effect on the wave forces of the caisson. Compared with the non-dimensional inline wave force, the relative length-width ratio BID was shown to have significant influence on the non-dimensional transverse wave force.展开更多
In this paper,using natural earthquake P-wave arrival time data recorded by the seismic network in the surrounding area of Madoi,the three-dimensional fine P-wave crustal velocity structure at depths above 60 km in th...In this paper,using natural earthquake P-wave arrival time data recorded by the seismic network in the surrounding area of Madoi,the three-dimensional fine P-wave crustal velocity structure at depths above 60 km in the epicenter of the Madoi Ms7.4 earthquake was inverted using the double-difference seismic tomography method.On the basis of the relocation of the source of the aftershock sequence,we summarized the strip-shaped distribution characteristics along the strike of the Jiangcuo fault,revealing the significant heterogeneity of the crustal velocity structure in the source area.Research has found that most of the Madoi Ms7.4 aftershocks were located in the weak area of the high-speed anomaly in the upper crust.The focal depth changed with the velocity structure,showing obvious fluctuation and segmentation characteristics.There was a good correspondence between the spatial distribution and the velocity structure.The high-velocity bodies of the upper crust in the hypocenter area provided a medium environment for earthquake rupture,the low-velocity bodies of the middle crust formed the deep material,and the migration channel and the undulating shape of the high-speed body in the lower crust corroborated the strong pushing action in the region.The results confirmed that under the continuous promotion of tectonic stress in the Madoi area,the high-speed body of the Jiangcuo fault blocked the migration of weak materials in the middle crust.When the stress accumulation exceeded the limit,the Madoi Ms7.4 earthquake occurred.Meanwhile,the nonuniform velocity structure near the fault plane determined the location of the main shock and the spatiotemporal distribution of the aftershock sequence.展开更多
A class of two-level high-order accuracy explicit difference scheme for solving 3-D parabolic P.D.E is constructed. Its truncation error is (Δt2+Δx4) and the stability condition is r=Δt/Δx2=Δt/Δy2=Δt/Δz2≤1/6.
Based on the VHF lightning locating system,a three-dimensional-space cell-gridded approach is used to extract the lighting channel and calculate the length of the channel.Through clustering of the located radiation so...Based on the VHF lightning locating system,a three-dimensional-space cell-gridded approach is used to extract the lighting channel and calculate the length of the channel.Through clustering of the located radiation sources and then extracting the lightning channel,it can accurately obtain the length of the channel.To validate the feasibility of the approach,a simulation experiment is designed,and it shows the length error is no more than 10%.The relationship between the NO production of per unit arc length and atmospheric pressure obtained in laboratory is applied to the NOX production of per unit flash length at different altitudes in this paper.The channel length and the NOX production of 11 negative cloud-to-ground flashes and 59 intracloud flashes in an isolated thunderstorm in the northeastern Qinghai-Tibet Plateau are calculated.The results show that the average channel lengths of per cloud-to-ground and intracloud flash are 28.9 and 22.3 km respectively;the average NOX productions of per cloud-to-ground and intracloud flash are 1.89×1025 and 0.42×1025 molecules,respectively.展开更多
We study the effect of "ghost forces" for a quasicontinuum method in three dimension with a planar interface. "Ghost forces" are the inconsistency of the quasicontinuum method across the interface between the atom...We study the effect of "ghost forces" for a quasicontinuum method in three dimension with a planar interface. "Ghost forces" are the inconsistency of the quasicontinuum method across the interface between the atomistic region and the continuum region. Numerical results suggest that "ghost forces" may lead to a negilible error on the solution, while lead to a finite size error on the gradient of the solution. The error has a layer-like profile, and the interfacial layer width is of O(ε). The error in certain component of the displacement gradient decays algebraically from O(1) to O(ε) away from the interface. A surrogate model is proposed and analyzed, which suggests the same scenario for the effect of "ghost forces". Our analysis is based on the explicit solution of the surrogate model.展开更多
基金The National Natural Science Foundation of China(No.50908004,51178013)
文摘This paper presents a method to characterize asphalt pavement macrotexture using the gray-tone difference matrix (GTDM)and discusses the potentials of the GTDM indicators for skid resistance evaluation.There are 37 field sites included in the data collection,which cover 6 types of asphalt pavement surfaces. The mean profile depth derived from 3-D macrotexture measurements (MPD3 ) has a significant relationship with the mean texture depth (MTD ),which can be described by a logarithm model with R2 of 0.962.There is no significant linear relationship between the friction coefficient at a speed of 60 km/h (DFT60 )and macrotexture indicators.A nonlinear model with British pendulum number (BPN ) incorporated can relate DFT60 to MTD or indicator fcon .A comparison with MTD shows that GTDM-based fcon has a potential to be a macrotexture indicator for skid resistance evaluation,which describes the general height difference and the average local height difference of pavement macrotexture. A relatively high fcon is helpful for improving asphalt pavement skid resistance.
文摘A family of high_order accuracy explicit difference schemes for solving 2_dimension parabolic P.D.E. are constructed. Th e stability condition is r=Δt/Δx 2=Δt/Δy 2【1/2 and the truncation err or is O(Δt 3+Δx 4).
基金the Nature Science Foundation of China under Grant No.60604019 and No.6075001
文摘The error equation of a rotating inertial navigation system was introduced. The effect of the system's main error source (constant drift of gyro and zero bias of accelerometer) under rotating conditions for the system was analyzed. Validity of theoretical analysis was shown via simulation, and that provides a theoretical foundation for a rotating strap-down inertial navigation system during actual experimentation and application.
基金Project(51078127) supported by the National Natural Science Foundation of ChinaProject(JJ201109091631) supported by the Foundation for Young Scientists of Jiangxi Education Department, China
文摘A modified one-dimensional transient hygrothermal model for multilayer wall was proposed using air humidity ratio and temperature as the driving potentials.The solution for the governing equations was obtained numerically by implementing the finite-difference scheme.To evaluate the accuracy of the model,a test system was built up to measure relative humidity and temperature within a porous wall and compare with the prediction of the model.The prediction results have good agreement with the experimental results.For the interface close to indoor side,the maximum deviation of temperature between calculated and test data is 1.87 K,and the average deviation is 0.95 K;the maximum deviation of relative humidity is 11.4%,and the average deviation is 5.7%.For the interface close to outdoor side,the maximum deviation of temperature between prediction and measurement is 1.78 K,and the average deviation is 1.1 K;the maximum deviation of relative humidity is 9.9%,and the average deviation is 4.2%.
基金Supported by the Science Foundations of LCP and CAEP under Grant Nos.2009A0102005 and 2009B0101012the National Basic Research Program (973 Program) under Grant No.2007CB815105the National Natural Science Foundation under Grant Nos.10775018,10702010,and 10775088
文摘We present a highly efficient lattice Boltzmann model for simulating compressible flows. This model is based on the combination of an appropriate finite difference scheme, a 16-discrete-velocity model [Kataoka and Tsutahara, Phys. Rev. E 69 (2004) 035701(R)] and reasonable dispersion and dissipation terms. The dispersion term effectively reduces the oscillation at the discontinuity and enhances numerical precision. The dissipation term makes the new model more easily meet with the yon Neumann stability condition. This model works for both high-speed and low-speed flows with arbitrary specific-heat-ratio. With the new model simulation results for the well-known benchmark problems get a high accuracy compared with the analytic or experimental ones. The used benchmark tests include (i) Shock tubes such as the Sod, Lax, Sjogreen, Colella explosion wave, and collision of two strong shocks, (ii) Regular and Mach shock reflections, and (iii) Shock wave reaction on cylindrical bubble problems. With a more realistic equation of state or free-energy functional, the new model has the potential tostudy the complex procedure of shock wave reaction on porous materials.
基金Supported by the Science and Technology Program on Transportation Construction in Western China Ministry of Communications under Grant No.2004-328-832-51
文摘A three dimensional numerical model of nonlinear wave action on a quasi-ellipse caisson in a time domain was developed in this paper. Navier-Stokes equations were solved by the finite difference method, and the volume of fluid (VOF) method was employed to trace the free surface. The partial cell method was used to deal with the irregular boundary typical of this type of problem during first-time wave interaction with the structure, and a satisfactory result was obtained. The numerical model was verified and used to investigate the effects of the relative wave height H/d, relative caisson width kD, and relative length-width ratio B/D on the wave forces of the quasi-ellipse caisson. It was shown that the relative wave height H/d has a significant effect on the wave forces of the caisson. Compared with the non-dimensional inline wave force, the relative length-width ratio BID was shown to have significant influence on the non-dimensional transverse wave force.
基金support:Seismic Regime Tracking Project of CEA (2023010123)Combination Project with Monitoring,Prediction and Scientific Research of Earthquake Technology,CEA (3JH-202302019).
文摘In this paper,using natural earthquake P-wave arrival time data recorded by the seismic network in the surrounding area of Madoi,the three-dimensional fine P-wave crustal velocity structure at depths above 60 km in the epicenter of the Madoi Ms7.4 earthquake was inverted using the double-difference seismic tomography method.On the basis of the relocation of the source of the aftershock sequence,we summarized the strip-shaped distribution characteristics along the strike of the Jiangcuo fault,revealing the significant heterogeneity of the crustal velocity structure in the source area.Research has found that most of the Madoi Ms7.4 aftershocks were located in the weak area of the high-speed anomaly in the upper crust.The focal depth changed with the velocity structure,showing obvious fluctuation and segmentation characteristics.There was a good correspondence between the spatial distribution and the velocity structure.The high-velocity bodies of the upper crust in the hypocenter area provided a medium environment for earthquake rupture,the low-velocity bodies of the middle crust formed the deep material,and the migration channel and the undulating shape of the high-speed body in the lower crust corroborated the strong pushing action in the region.The results confirmed that under the continuous promotion of tectonic stress in the Madoi area,the high-speed body of the Jiangcuo fault blocked the migration of weak materials in the middle crust.When the stress accumulation exceeded the limit,the Madoi Ms7.4 earthquake occurred.Meanwhile,the nonuniform velocity structure near the fault plane determined the location of the main shock and the spatiotemporal distribution of the aftershock sequence.
文摘A class of two-level high-order accuracy explicit difference scheme for solving 3-D parabolic P.D.E is constructed. Its truncation error is (Δt2+Δx4) and the stability condition is r=Δt/Δx2=Δt/Δy2=Δt/Δz2≤1/6.
基金supported by the National Natural Science Foundation of Chi-na(Grant Nos.41075002,40775004,41030960)R&D Special Fund for Public Welfare Industry(Meteorology)(Grant No.GYHY201006005-03)
文摘Based on the VHF lightning locating system,a three-dimensional-space cell-gridded approach is used to extract the lighting channel and calculate the length of the channel.Through clustering of the located radiation sources and then extracting the lightning channel,it can accurately obtain the length of the channel.To validate the feasibility of the approach,a simulation experiment is designed,and it shows the length error is no more than 10%.The relationship between the NO production of per unit arc length and atmospheric pressure obtained in laboratory is applied to the NOX production of per unit flash length at different altitudes in this paper.The channel length and the NOX production of 11 negative cloud-to-ground flashes and 59 intracloud flashes in an isolated thunderstorm in the northeastern Qinghai-Tibet Plateau are calculated.The results show that the average channel lengths of per cloud-to-ground and intracloud flash are 28.9 and 22.3 km respectively;the average NOX productions of per cloud-to-ground and intracloud flash are 1.89×1025 and 0.42×1025 molecules,respectively.
基金supported by National Natural Science Foundation of China(Grant Nos.1093201191230203 and 11021101)
文摘We study the effect of "ghost forces" for a quasicontinuum method in three dimension with a planar interface. "Ghost forces" are the inconsistency of the quasicontinuum method across the interface between the atomistic region and the continuum region. Numerical results suggest that "ghost forces" may lead to a negilible error on the solution, while lead to a finite size error on the gradient of the solution. The error has a layer-like profile, and the interfacial layer width is of O(ε). The error in certain component of the displacement gradient decays algebraically from O(1) to O(ε) away from the interface. A surrogate model is proposed and analyzed, which suggests the same scenario for the effect of "ghost forces". Our analysis is based on the explicit solution of the surrogate model.