Based on the stochastic AMR model, this paper constructs man-made earthquake catalogues to investigate the property of parameter estimation of the model. Then the stochastic AMR model is applied to the study of severa...Based on the stochastic AMR model, this paper constructs man-made earthquake catalogues to investigate the property of parameter estimation of the model. Then the stochastic AMR model is applied to the study of several strong earthquakes in China and New Zealand. Akaikes AIC criterion is used to discriminate whether an accelerating mode of earthquake activity precedes those events or not. Finally, regional accelerating seismic activity and possible prediction approach for future strong earthquakes are discussed.展开更多
Quasi-likelihood nonlinear models (QLNM) include generalized linear models as a special case. Under some regularity conditions, the rate of the strong consistency of the maximum quasi-likelihood estimation (MQLE) ...Quasi-likelihood nonlinear models (QLNM) include generalized linear models as a special case. Under some regularity conditions, the rate of the strong consistency of the maximum quasi-likelihood estimation (MQLE) is obtained in QLNM. In an important case, this rate is O(n-^1/2(loglogn)^1/2), which is just the rate of LIL of partial sums for i.i.d variables, and thus cannot be improved anymore.展开更多
In this article, we are concerned with the strong solutions for the incompress- ible fluid models of Korteweg type in a bounded domain Ω СR^3. We prove the existence and uniqueness of local strong solutions to the i...In this article, we are concerned with the strong solutions for the incompress- ible fluid models of Korteweg type in a bounded domain Ω СR^3. We prove the existence and uniqueness of local strong solutions to the initial boundary value problem. We point out that in this article we allow the existence of initial vacuum provided initial data satisfy a compatibility condition.展开更多
It is well known that quantitative estimation of slip distributions on fault plane is one of the most important issues for earthquake source inversion related to the fault rupture process. The characteristics of slip ...It is well known that quantitative estimation of slip distributions on fault plane is one of the most important issues for earthquake source inversion related to the fault rupture process. The characteristics of slip distribution on the main fault play a fundamental role to control strong ground motion pattern. A large amount of works have also suggested that variable slip models inverted from longer period ground motion recordings are relevant for the prediction of higher frequency ground motions. Zhang et al. (Chin J Geophys 56:1412-1417, 2013) and Wang et al. (Chin J Geophys 56:1408-1411,2013) published their source inversions for the fault rupturing process soon after the April 20, 2013 Lushan earthquake in Sichuan, China. In this study, first, we synthesize two forward source slip models: the value of maximum slip, fault dimension, size, and dimension of major asperities, and comer wave number obtained from Wang's model is adopted to constrain the gen- eration of k-2 model and crack model. Next, both inverted and synthetic slip models are employed to simulate the ground motions for the Lushan earthquake based on the stochastic finite-fault method. In addition, for a comparison purpose, a stochastic slip model and another k-2 model (k 2 model II) with 2 times value of comer wave number of the original k-2 model (k 2 model I) are also employed for simulation for Lushan event. The simulated results characterized by Modified Mer- calli Intensity (MMI) show that the source slip models based on the inverted and synthetic slip distributions could capture many basic features associated with the ground motion patterns. Moreover, the simulated MMI distributions reflect the rupture directivity effect and the influence of the shallow velocity structure well. On the other hand, the simulated MMI bystochastic slip model and k 2 model II is apparently higher than observed intensity. By contrast, our simulation results show that the higher frequency ground motion is sensitive to the degree of slip roughness; therefore, we suggest that, for realistic ground- motion simulations due to future earthquake, it is imperative to properly estimate the slip roughness distribution.展开更多
The black hole model will be excluded by a very strong radial magnetic field near the Galactic Center which has been detected in 2013. Following it, the explosion mechanism, for both supernova and the hot big bang of ...The black hole model will be excluded by a very strong radial magnetic field near the Galactic Center which has been detected in 2013. Following it, the explosion mechanism, for both supernova and the hot big bang of the Universe, driven by magnetic monopoles is proposed in this paper.展开更多
We study the cosmic constraint to the wCDM (cold dark matter with a constant equation of state w) model via 118 strong gravitational lensing systems which are compiled from SLA CS, BELLS, LSD and SL2S surveys, where...We study the cosmic constraint to the wCDM (cold dark matter with a constant equation of state w) model via 118 strong gravitational lensing systems which are compiled from SLA CS, BELLS, LSD and SL2S surveys, where the ratio between two angular diameter distances Dobs =DA(Zl, Zs ) / D A ( O, Zs ) is taken as a cosmic observable. To obtain this ratio, we adopt two strong tensing models: one is the singular isothermal sphere model (SIS) and the other one is the power-law density profile (PLP) model. Via the Markov chain Monte Carlo method, the posterior distribution of the cosmological model parameters space is obtained. The results show that the cosmological model parameters are not sensitive to the parameterized forms of the power-law index γ. Furthermore, the PLP model gives a relatively tighter constraint to the cosmological parameters than that of the SIS model. The predicted value of Ωm = 0.31+0.44 -0.24 by the SIS model is compatible with that obtained by P1anck2015: Ωm = 0.313 ± 0.013. However, the value of Ωm =0.15+0.13 -0.11 based on the PLP model is smaller and has 1.25σ tension with that obtained by Planck2015.展开更多
Western Sichuan and its vicinity are located in the juncture of three big active blocks, namely, the Chuandian (Sichuan and Yunnan), the Bayan Har and the South China blocks, on the eastern margin of the Qinghai-Xiz...Western Sichuan and its vicinity are located in the juncture of three big active blocks, namely, the Chuandian (Sichuan and Yunnan), the Bayan Har and the South China blocks, on the eastern margin of the Qinghai-Xizang(Tibet) Plateau. Many groups of active faults that are capable of generating earthquakes are developed there. Because there exist lateral secondary active faults, the Chuandian block can be further divided into the central Yumlan and northwestern Sichuan sub-blocks; while the Longmenshan sub-block can be divided on the east end of the Bayan Har block. Joint exploration of deep crustal structure shows that there exist low-velocity and high-conductivity layers in the crust of the Chuandian and Bayan Har blocks, which are one of the important factors that make the upper crust prone to earthquake. The results of geological study and modern GPS observation show that blocks of different orders all have SE- or SSE-trending sliding, clockwise rotation and upwelling movement; but there are some differences in amplitude. This paper has also given the geological or GPS slip rates of main active fault zones and discussed the main scientific problems still existing now.展开更多
The notion of equivalent dispersion of strong seismic waves is further illuminated and a frequency-dependentstochastic model of the Fourier phase of strong earthquake ground motion is presented for use in statistical ...The notion of equivalent dispersion of strong seismic waves is further illuminated and a frequency-dependentstochastic model of the Fourier phase of strong earthquake ground motion is presented for use in statistical estimates of the Fourier phase. The non-stationary features of design ground motions synthesized using the phasemodel in combination with the empirical attenuation model of amplitude spectra of strong ground motion are ableto be consistent with those of observed records in the sense of statistics.展开更多
Using finite differences and entropy inequalities, the global existence of weak solutions to a multidimensional parabolic strongly coupled prey-predator model is obtained. The nonnegativity of the solutions is also sh...Using finite differences and entropy inequalities, the global existence of weak solutions to a multidimensional parabolic strongly coupled prey-predator model is obtained. The nonnegativity of the solutions is also shown.展开更多
An analytical variational method for the ground state of the biased quantum Rabi model in the ultra-strong coupling regime is presented. This analytical variational method can be obtained by a unitary transformation o...An analytical variational method for the ground state of the biased quantum Rabi model in the ultra-strong coupling regime is presented. This analytical variational method can be obtained by a unitary transformation or alternatively by assuming the form of the ground state wave function. The key of the method is to introduce a variational parameter λ,which can be determined by minimizing the energy functional. Using this method, we calculate the physical observables with high accuracy in comparison with the numerical exact ones. Our method evidently improves over the widely used general rotating-wave approximation(GRWA) in both qualitative and quantitative aspects.展开更多
Seismic hazard assessment based on urban active faults can provide scientific bases for city planning and project construction, while numerical simulation of strong ground motion is an important method for seismic haz...Seismic hazard assessment based on urban active faults can provide scientific bases for city planning and project construction, while numerical simulation of strong ground motion is an important method for seismic hazard prediction and assessment. A 3-D physical model in conformity with real strata configuration of (mainly) the Quaternary is a prerequisite to ensure the reliability of the simulation results. In this paper, we give a detailed account of the technical scheme and process for creating a 3-D physical model in Kunming basin. The data used are synthe- sized from seismogeological data, borehole data, topographic data, digital elevation mode (DEM) data, seismic exploration results and wave velocity measurements. Strafigraphic division is based mainly on shear wave velocity, with strata sequence taken into consideration. The model construction is finally accomplished with ArcGIS and many relevant programming techniques via layer-by-layer stacking (in depth direction) of the adjacent medium interfaces (meshes). Meanwhile, a database of 3-D physical models is set up, which provides model data and parameters for strong ground motion simulation. Some processing methods and significant issues are also addressed in the paper in accordance with different types of exploration and experimental data.展开更多
Abstract Accurate simulation of seismic wave propaga- tion in complex geological structures is of particular interest nowadays. However conventional methods may fail to simulate realistic wavefields in environments wi...Abstract Accurate simulation of seismic wave propaga- tion in complex geological structures is of particular interest nowadays. However conventional methods may fail to simulate realistic wavefields in environments with great and rapid structural changes, due for instance to the presence of shadow zones, diffractions and/or edge effects. Different methods, developed to improve seismic model- ing, are typically tested on synthetic configurations against analytical solutions for simple canonical problems or ref- erence methods, or via direct comparison with real data acquired in situ. Such approaches have limitations,especially if the propagation occurs in a complex envi- ronment with strong-contrast reflectors and surface irreg- ularities, as it can be difficult to determine the method which gives the best approximation of the "real" solution, or to interpret the results obtained without an a priori knowledge of the geologic environment. An alternative approach for seismics consists in comparing the synthetic data with high-quality data collected in laboratory experi- ments under controlled conditions for a known configuration. In contrast with numerical experiments, laboratory data possess many of the characteristics of field data, as real waves propagate through models with no numerical approximations. We thus present a comparison of laboratory-scaled measurements of 3D zero-offset wave reflection of broadband pulses from a strong topographic environment immersed in a water tank with numerical data simulated by means of a spectral-element method and a discretized Kirchhoff integral method. The results indicate a good quantitative fit in terms of time arrivals and acceptable fit in amplitudes for all datasets.展开更多
This article is concerned with the estimating problem of semiparametric varyingcoefficient partially linear regression models. By combining the local polynomial and least squares procedures Fan and Huang (2005) prop...This article is concerned with the estimating problem of semiparametric varyingcoefficient partially linear regression models. By combining the local polynomial and least squares procedures Fan and Huang (2005) proposed a profile least squares estimator for the parametric component and established its asymptotic normality. We further show that the profile least squares estimator can achieve the law of iterated logarithm. Moreover, we study the estimators of the functions characterizing the non-linear part as well as the error variance. The strong convergence rate and the law of iterated logarithm are derived for them, respectively.展开更多
Based on the stochastic AMR model, this paper constructs man-made earthquake catalogues to investigate the property of parameter estimation of the model. Then the stochastic AMR model is applied to the study of severa...Based on the stochastic AMR model, this paper constructs man-made earthquake catalogues to investigate the property of parameter estimation of the model. Then the stochastic AMR model is applied to the study of several strong earthquakes in China and New Zealand. Akaikes AIC criterion is used to discriminate whether an accelerating mode of earthquake activity precedes those events or not. Finally, regional accelerating seismic activity and possible prediction approach for future strong earthquakes are discussed.展开更多
In order to accurately simulate strong three-dimensional (3-D) free surface flows and sediment transport, the fully 3- D non-hydrostatic pressure models are developed based on the incompressible Navier-Stokes equati...In order to accurately simulate strong three-dimensional (3-D) free surface flows and sediment transport, the fully 3- D non-hydrostatic pressure models are developed based on the incompressible Navier-Stokes equations and convection-diffusion equation of sediment concentration with the mixing triangle and quadrilateral grids. The governing equations are discretized with the unstructured finite volume method in order to provide conservation properties of mass and momentum, and flexibility with practical application. It is shown that it is first-order accurate on nonuniform plane two-dimensional (2-D) grids and second-order accurate on uniform plane grids. A third-order approximation of the vertical velocity at the top-layer is applied. In such a way, free surface zero stress boundary condition is satisfied maturely, and very few vertical layers are needed to give an accurate solution even for complex discontinuous flow and short wave simulation. The model is applied to four examples to simulate strong 3-D free surface flows and sediment transport where non-hydrostatic pressures have a considerable effect on the velocity field. The newly developed model is verified against analytical solutions with an excellent agreement.展开更多
基金National Natural Science Foundation of China (4007401340134010)Chinese Joint Seismological Science Foundation (042002) and the project during the Tenth Five-year Plan.
文摘Based on the stochastic AMR model, this paper constructs man-made earthquake catalogues to investigate the property of parameter estimation of the model. Then the stochastic AMR model is applied to the study of several strong earthquakes in China and New Zealand. Akaikes AIC criterion is used to discriminate whether an accelerating mode of earthquake activity precedes those events or not. Finally, regional accelerating seismic activity and possible prediction approach for future strong earthquakes are discussed.
基金Supported by the National Natural Sciences Foundation of China (10761011)Mathematical Tianyuan Fund of National Natural Science Fundation of China(10626048)
文摘Quasi-likelihood nonlinear models (QLNM) include generalized linear models as a special case. Under some regularity conditions, the rate of the strong consistency of the maximum quasi-likelihood estimation (MQLE) is obtained in QLNM. In an important case, this rate is O(n-^1/2(loglogn)^1/2), which is just the rate of LIL of partial sums for i.i.d variables, and thus cannot be improved anymore.
基金Supported by NSF (10531020) of Chinathe Programof 985 Innovation Engineering on Information in Xiamen University (2004-2007) and NCETXMU
文摘In this article, we are concerned with the strong solutions for the incompress- ible fluid models of Korteweg type in a bounded domain Ω СR^3. We prove the existence and uniqueness of local strong solutions to the initial boundary value problem. We point out that in this article we allow the existence of initial vacuum provided initial data satisfy a compatibility condition.
基金supported by the CAS/SAFEA International Partnership Program for Creative Research Teams(KZZD-EW-TZ-19)
文摘It is well known that quantitative estimation of slip distributions on fault plane is one of the most important issues for earthquake source inversion related to the fault rupture process. The characteristics of slip distribution on the main fault play a fundamental role to control strong ground motion pattern. A large amount of works have also suggested that variable slip models inverted from longer period ground motion recordings are relevant for the prediction of higher frequency ground motions. Zhang et al. (Chin J Geophys 56:1412-1417, 2013) and Wang et al. (Chin J Geophys 56:1408-1411,2013) published their source inversions for the fault rupturing process soon after the April 20, 2013 Lushan earthquake in Sichuan, China. In this study, first, we synthesize two forward source slip models: the value of maximum slip, fault dimension, size, and dimension of major asperities, and comer wave number obtained from Wang's model is adopted to constrain the gen- eration of k-2 model and crack model. Next, both inverted and synthetic slip models are employed to simulate the ground motions for the Lushan earthquake based on the stochastic finite-fault method. In addition, for a comparison purpose, a stochastic slip model and another k-2 model (k 2 model II) with 2 times value of comer wave number of the original k-2 model (k 2 model I) are also employed for simulation for Lushan event. The simulated results characterized by Modified Mer- calli Intensity (MMI) show that the source slip models based on the inverted and synthetic slip distributions could capture many basic features associated with the ground motion patterns. Moreover, the simulated MMI distributions reflect the rupture directivity effect and the influence of the shallow velocity structure well. On the other hand, the simulated MMI bystochastic slip model and k 2 model II is apparently higher than observed intensity. By contrast, our simulation results show that the higher frequency ground motion is sensitive to the degree of slip roughness; therefore, we suggest that, for realistic ground- motion simulations due to future earthquake, it is imperative to properly estimate the slip roughness distribution.
文摘The black hole model will be excluded by a very strong radial magnetic field near the Galactic Center which has been detected in 2013. Following it, the explosion mechanism, for both supernova and the hot big bang of the Universe, driven by magnetic monopoles is proposed in this paper.
基金Supported by the National Natural Science Foundation of China under Grant No 11275035
文摘We study the cosmic constraint to the wCDM (cold dark matter with a constant equation of state w) model via 118 strong gravitational lensing systems which are compiled from SLA CS, BELLS, LSD and SL2S surveys, where the ratio between two angular diameter distances Dobs =DA(Zl, Zs ) / D A ( O, Zs ) is taken as a cosmic observable. To obtain this ratio, we adopt two strong tensing models: one is the singular isothermal sphere model (SIS) and the other one is the power-law density profile (PLP) model. Via the Markov chain Monte Carlo method, the posterior distribution of the cosmological model parameters space is obtained. The results show that the cosmological model parameters are not sensitive to the parameterized forms of the power-law index γ. Furthermore, the PLP model gives a relatively tighter constraint to the cosmological parameters than that of the SIS model. The predicted value of Ωm = 0.31+0.44 -0.24 by the SIS model is compatible with that obtained by P1anck2015: Ωm = 0.313 ± 0.013. However, the value of Ωm =0.15+0.13 -0.11 based on the PLP model is smaller and has 1.25σ tension with that obtained by Planck2015.
文摘Western Sichuan and its vicinity are located in the juncture of three big active blocks, namely, the Chuandian (Sichuan and Yunnan), the Bayan Har and the South China blocks, on the eastern margin of the Qinghai-Xizang(Tibet) Plateau. Many groups of active faults that are capable of generating earthquakes are developed there. Because there exist lateral secondary active faults, the Chuandian block can be further divided into the central Yumlan and northwestern Sichuan sub-blocks; while the Longmenshan sub-block can be divided on the east end of the Bayan Har block. Joint exploration of deep crustal structure shows that there exist low-velocity and high-conductivity layers in the crust of the Chuandian and Bayan Har blocks, which are one of the important factors that make the upper crust prone to earthquake. The results of geological study and modern GPS observation show that blocks of different orders all have SE- or SSE-trending sliding, clockwise rotation and upwelling movement; but there are some differences in amplitude. This paper has also given the geological or GPS slip rates of main active fault zones and discussed the main scientific problems still existing now.
文摘The notion of equivalent dispersion of strong seismic waves is further illuminated and a frequency-dependentstochastic model of the Fourier phase of strong earthquake ground motion is presented for use in statistical estimates of the Fourier phase. The non-stationary features of design ground motions synthesized using the phasemodel in combination with the empirical attenuation model of amplitude spectra of strong ground motion are ableto be consistent with those of observed records in the sense of statistics.
基金supported by the National Natural Science Foundation of China (Nos. 10701024, 10601011)
文摘Using finite differences and entropy inequalities, the global existence of weak solutions to a multidimensional parabolic strongly coupled prey-predator model is obtained. The nonnegativity of the solutions is also shown.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.11674139,11604009,and 11704025)the Program for Changjiang Scholars and Innovative Research Team in University,China(Grant No.IRT-16R35)+1 种基金the Fundamental Research Funds for the Central Universities,Chinathe financial support of the Future and Emerging Technologies(FET)programme within the Seventh Framework Programme for Research of the European Commission,under FET-Open Grant No.618083(CNTQC)
文摘An analytical variational method for the ground state of the biased quantum Rabi model in the ultra-strong coupling regime is presented. This analytical variational method can be obtained by a unitary transformation or alternatively by assuming the form of the ground state wave function. The key of the method is to introduce a variational parameter λ,which can be determined by minimizing the energy functional. Using this method, we calculate the physical observables with high accuracy in comparison with the numerical exact ones. Our method evidently improves over the widely used general rotating-wave approximation(GRWA) in both qualitative and quantitative aspects.
基金Urban Active Fault Detection Project sponsored by the National Development andReform Commission of China(2004-1138)National Natural Science Foundation of China(40604005)+1 种基金Joint Seismological Science Foundation of China(A07027)Key Project during the 10th Five-year Plan Period from Earthquake Administration of China(J105B-16).
文摘Seismic hazard assessment based on urban active faults can provide scientific bases for city planning and project construction, while numerical simulation of strong ground motion is an important method for seismic hazard prediction and assessment. A 3-D physical model in conformity with real strata configuration of (mainly) the Quaternary is a prerequisite to ensure the reliability of the simulation results. In this paper, we give a detailed account of the technical scheme and process for creating a 3-D physical model in Kunming basin. The data used are synthe- sized from seismogeological data, borehole data, topographic data, digital elevation mode (DEM) data, seismic exploration results and wave velocity measurements. Strafigraphic division is based mainly on shear wave velocity, with strata sequence taken into consideration. The model construction is finally accomplished with ArcGIS and many relevant programming techniques via layer-by-layer stacking (in depth direction) of the adjacent medium interfaces (meshes). Meanwhile, a database of 3-D physical models is set up, which provides model data and parameters for strong ground motion simulation. Some processing methods and significant issues are also addressed in the paper in accordance with different types of exploration and experimental data.
基金the INSIS Institute of the French CNRS,Aix-Marseille Universitythe Carnot Star Institute,the VISTA Projectthe Norwegian Research Council through the ROSE Project for financial support
文摘Abstract Accurate simulation of seismic wave propaga- tion in complex geological structures is of particular interest nowadays. However conventional methods may fail to simulate realistic wavefields in environments with great and rapid structural changes, due for instance to the presence of shadow zones, diffractions and/or edge effects. Different methods, developed to improve seismic model- ing, are typically tested on synthetic configurations against analytical solutions for simple canonical problems or ref- erence methods, or via direct comparison with real data acquired in situ. Such approaches have limitations,especially if the propagation occurs in a complex envi- ronment with strong-contrast reflectors and surface irreg- ularities, as it can be difficult to determine the method which gives the best approximation of the "real" solution, or to interpret the results obtained without an a priori knowledge of the geologic environment. An alternative approach for seismics consists in comparing the synthetic data with high-quality data collected in laboratory experi- ments under controlled conditions for a known configuration. In contrast with numerical experiments, laboratory data possess many of the characteristics of field data, as real waves propagate through models with no numerical approximations. We thus present a comparison of laboratory-scaled measurements of 3D zero-offset wave reflection of broadband pulses from a strong topographic environment immersed in a water tank with numerical data simulated by means of a spectral-element method and a discretized Kirchhoff integral method. The results indicate a good quantitative fit in terms of time arrivals and acceptable fit in amplitudes for all datasets.
基金supported by the National Natural Science Funds for Distinguished Young Scholar (70825004)National Natural Science Foundation of China (NSFC) (10731010 and 10628104)+3 种基金the National Basic Research Program (2007CB814902)Creative Research Groups of China (10721101)Leading Academic Discipline Program, the 10th five year plan of 211 Project for Shanghai University of Finance and Economics211 Project for Shanghai University of Financeand Economics (the 3rd phase)
文摘This article is concerned with the estimating problem of semiparametric varyingcoefficient partially linear regression models. By combining the local polynomial and least squares procedures Fan and Huang (2005) proposed a profile least squares estimator for the parametric component and established its asymptotic normality. We further show that the profile least squares estimator can achieve the law of iterated logarithm. Moreover, we study the estimators of the functions characterizing the non-linear part as well as the error variance. The strong convergence rate and the law of iterated logarithm are derived for them, respectively.
基金National Natural Science Foundation of China (40074013, 40134010), Chinese Joint Seismological Science Foundation (042002) and the project during the Tenth Five-year Plan.
文摘Based on the stochastic AMR model, this paper constructs man-made earthquake catalogues to investigate the property of parameter estimation of the model. Then the stochastic AMR model is applied to the study of several strong earthquakes in China and New Zealand. Akaikes AIC criterion is used to discriminate whether an accelerating mode of earthquake activity precedes those events or not. Finally, regional accelerating seismic activity and possible prediction approach for future strong earthquakes are discussed.
基金financially supported by the Science and Technology Project of the Ministry of Transport(Grant No.2013328352570)
文摘In order to accurately simulate strong three-dimensional (3-D) free surface flows and sediment transport, the fully 3- D non-hydrostatic pressure models are developed based on the incompressible Navier-Stokes equations and convection-diffusion equation of sediment concentration with the mixing triangle and quadrilateral grids. The governing equations are discretized with the unstructured finite volume method in order to provide conservation properties of mass and momentum, and flexibility with practical application. It is shown that it is first-order accurate on nonuniform plane two-dimensional (2-D) grids and second-order accurate on uniform plane grids. A third-order approximation of the vertical velocity at the top-layer is applied. In such a way, free surface zero stress boundary condition is satisfied maturely, and very few vertical layers are needed to give an accurate solution even for complex discontinuous flow and short wave simulation. The model is applied to four examples to simulate strong 3-D free surface flows and sediment transport where non-hydrostatic pressures have a considerable effect on the velocity field. The newly developed model is verified against analytical solutions with an excellent agreement.