An essentially conservative adaptive space time conservation element and solution element (CE/SE) method is pro- posed for the effective simulation of shock-induced instability with low computational cost. Its imple...An essentially conservative adaptive space time conservation element and solution element (CE/SE) method is pro- posed for the effective simulation of shock-induced instability with low computational cost. Its implementation is based on redefined conservation elements (CEs) and solution elements (SEs), optimized interpolations and a Courant number insensitive CE/SE scheme. This approach is used in two applications, the Woodward double Mach reflection and a two- component Richtmyer-Meshkov instability experiment. This scheme reveals the essential features of the investigated cases, captures small unstable structures, and yields a solution that is consistent with the results from experiments or other high order methods.展开更多
In order to analyze the mechanism of continental strong shocks from the angle of crustal movement using the data of repeated geodetic survey, this paper has proposed a physical method; it analyzes the mechanism of den...In order to analyze the mechanism of continental strong shocks from the angle of crustal movement using the data of repeated geodetic survey, this paper has proposed a physical method; it analyzes the mechanism of density change due to the occurrence of strong shocks by use of the physical quantity that reflects the time change of crustal density. (1) The general theory of the time change of density in the earth’ s interior and the theory of the time change of single layer density have been introduced, and an algorithm of stepwise iteration has been proposed; (2) The effect of the change of single layer density caused by fault dislocation has been analyzed in brief; (3) The characteristics of the time change of crustal density in the south of the seismogenic region before the 1996 Lijiang earthquake with M_L =7.0 have been studied; (4) The precursor model or causal mechanism of strong shocks possibly existing in the time change of crustal density has been investigated preliminarily.展开更多
An enhanced differential transform method (EDTM), which introduces the Pad@ technique into the standard differential transform method (DTM), is proposed. The enhanced method is applied to the analytic treatment of...An enhanced differential transform method (EDTM), which introduces the Pad@ technique into the standard differential transform method (DTM), is proposed. The enhanced method is applied to the analytic treatment of the shock wave. It accelerates the convergence of the series solution and provides an exact Dower series solution.展开更多
A new method is presented to analyze multi-degree-of-freedom (MDOF) dynamic systems subjected to an external shock excitation. A two-degree-of-freedom theoretical system with linear characteristics is exemplified to i...A new method is presented to analyze multi-degree-of-freedom (MDOF) dynamic systems subjected to an external shock excitation. A two-degree-of-freedom theoretical system with linear characteristics is exemplified to illustrate the procedure of this method. The equations of motion of the dynamic system are established via matrix method. The dynamic responses of the dynamic system under an external shock excitation of a half-sine type are obtained by MATLAB and ANSYS. It is proved that the new method is helpful to analyze MDOF dynamic systems.展开更多
In this paper, the CE/SE method is developed to simulate the two- and three-dimensional flow-field of Pulse Detonation Engine (PDE). The conservation equations with stiff source terms for chemical reaction are solved ...In this paper, the CE/SE method is developed to simulate the two- and three-dimensional flow-field of Pulse Detonation Engine (PDE). The conservation equations with stiff source terms for chemical reaction are solved in two steps. The detailed analysis of computational results of a PDE with a single detonation tube and a PDE with five detonation tubes are given in this paper. Complex wave systems are observed inside and outside a PDE. For a PDE with 5 detonation tubes, there is a big bow shock produced from a number of little shocks near the open ends of tubes. A lot of vortexes interact with shocks and a large expansion wave propagates forward and backward with respect to the PDE in a semi-oval shape.展开更多
Culture shock is one of the common experiences by students and people abroad. They may experience the feeling of anxiety, disorientation, confusion, and even hostility because of misunderstanding to each other."C...Culture shock is one of the common experiences by students and people abroad. They may experience the feeling of anxiety, disorientation, confusion, and even hostility because of misunderstanding to each other."Culture shock"refers to the phenomenon in which people engage in the interaction feel depressed, discomfortable, etc. Culture shock may be caused by the conflict of various value systems when students and people enter a new culture and other value system. Culture shock is also a consequence of the influence of the negative life event. The changes in one's life and lack of support from family members,friends will make them more vulnerable and cause great discomfort in the new environment.展开更多
Many problems at the forefront of theoretical astrophysics require a treatment of dynamical fluid behavior. We present an efficient high-resolution shock-capturing hydrodynamic scheme designed to study such phenomena....Many problems at the forefront of theoretical astrophysics require a treatment of dynamical fluid behavior. We present an efficient high-resolution shock-capturing hydrodynamic scheme designed to study such phenomena. We have implemented a weighted, essentially non-oscillatory (WENO) scheme to fifth order accuracy in space, lILLE approximate Riemann solver is used for the flux computation at cell interface, which does not require spectral decomposition into characteristic waves and so is computationally friendly. For time integration we apply a third order total variation diminishing (TVD) Runge-Kutta scheme. Extensive testing and comparison with schemes that re.quire characteristic decomposition are carried out demonstrating the ability of our scheme to address challenging open questions in astrophysics.展开更多
With the application of electronic detonators, millisecond blasting is regarded as a signifi cant promising approach to improve the rock fragmentation in deep rock blasting. Thus, it is necessary to investigate the fr...With the application of electronic detonators, millisecond blasting is regarded as a signifi cant promising approach to improve the rock fragmentation in deep rock blasting. Thus, it is necessary to investigate the fracturing mechanisms of short-delay blasting. In this work, a rectangle model with two circle boreholes is modeled as a particles assembly based on the discrete element method to simulate the shock wave interactions induced by millisecond blasting. The rectangle model has a size of 12 × 6 m (L × W) and two blast holes have the same diameter of 12 cm. The shock waves are simplifi ed as time-varying forces applied at the particles of walls of the two boreholes. Among a series of numerical tests in this study, the spacing between two adjacent boreholes and delay time of millisecond blasting are considered as two primary variables, and the decoupling charge with a coeffi cient of 1.5 is taken into account in each case. The results show that stress superposition is not a key factor for improving rock fragmentation (tensile stress interactions rather than compressive stress superposition could aff ect the generation of cracks), whereas collision actions from isolated particles or particles with weakened constraints play a crucial role in creating the fracture network. The delay time has an infl uence on causing cracks in rock blasting, however, whether it works heavily depends on the distance between the two holes.展开更多
High order accurate scheme is highly desirable for Slow computation with shocks. After analysis has been made for the reason of the generation of non-physical oscillations around the shock in numerical computations, a...High order accurate scheme is highly desirable for Slow computation with shocks. After analysis has been made for the reason of the generation of non-physical oscillations around the shock in numerical computations, a third-order, upwind biased, shock capturing scheme was proposed. Also, a new shock fitting method, called pseudo shock fitting method, was suggested, which in principle can be with any order of accuracy. Test cases for one dimensional flows show that the new method is very satisfactory.展开更多
The finite element method is presented to attain the numerical simulation of the residual stresses field in the material treated by laser shock processing. The distribution of residual stresses generated by a single l...The finite element method is presented to attain the numerical simulation of the residual stresses field in the material treated by laser shock processing. The distribution of residual stresses generated by a single laser shock with square and round laser spot is predicted and validated by experimental results. With the Finite Element Method (FEM) model, effects of different overlapping rates and impact sequences on the distribution of residual stresses are simulated. The results indicate that: (1) Overlapping laser shock can increase the compressive residual stresses. However, it is not effective on the growth of plastically affected depth; (2) Overlapping rate should be optimized and selected carefully for the large area treatment. Appropriate overlapping rate is beneficial to obtain a homogeneous residual stress field; (3) The impact sequence has a great effect on the residual stress field. It can greatly attenuate the phenomenon of the “residual stress hole” to obtain a homogeneous residual stress field.展开更多
We are concerned with global solutions of multidimensional(M-D)Riemann problems for nonlinear hyperbolic systems of conservation laws,focusing on their global configurations and structures.We present some recent devel...We are concerned with global solutions of multidimensional(M-D)Riemann problems for nonlinear hyperbolic systems of conservation laws,focusing on their global configurations and structures.We present some recent developments in the rigorous analysis of two-dimensional(2-D)Riemann problems involving transonic shock waves through several prototypes of hyperbolic systems of conservation laws and discuss some further M-D Riemann problems and related problems for nonlinear partial differential equations.In particular,we present four different 2-D Riemann problems through these prototypes of hyperbolic systems and show how these Riemann problems can be reformulated/solved as free boundary problems with transonic shock waves as free boundaries for the corresponding nonlinear conservation laws of mixed elliptic-hyperbolic type and related nonlinear partial differential equations.展开更多
The 'trial function method' ( TFM for short) and a routine way in finding traveling,wave solutions to some nonlinear partial differential equations( PDE for short), wer explained. Two types of evolution equati...The 'trial function method' ( TFM for short) and a routine way in finding traveling,wave solutions to some nonlinear partial differential equations( PDE for short), wer explained. Two types of evolution equations are studied, one is a generalized Burgers or KdV equation, the other is the Fisher equation with special nonlinear forms of its reaction rate term. One can see that this method is simple, fast and allowing further extension.展开更多
The deformation of reeds in automotive hydraulic shock absorbers is analyzed with the finite element method. Combination of different thick reeds mounted on different supports is studied. The computational results sho...The deformation of reeds in automotive hydraulic shock absorbers is analyzed with the finite element method. Combination of different thick reeds mounted on different supports is studied. The computational results show that deformation of the overlapped reeds is not always equal to the sum of deflection of single reed under any conditions. Experimental results prove computational results to be correct and computational method effective. The method of analysis and view of point can provide reference to the design and manufacture of hydraulic shock absorbers using reeds.展开更多
In this paper, we construct the equations of generalized thermoelasicity for a non-homogeneous isotropic hollow cylider with a variable modulus of elasticity and thermal conductivity based on the Lord and Shulman theo...In this paper, we construct the equations of generalized thermoelasicity for a non-homogeneous isotropic hollow cylider with a variable modulus of elasticity and thermal conductivity based on the Lord and Shulman theory. The problem has been solved numerically using the finite element method. Numerical results for the displacement, the temperature, the radial stress, and the hoop stress distributions are illustrated graphically. Comparisons are made between the results predicted by the coupled theory and by the theory of generalized thermoelasticity with one relaxation time in the cases of temperature dependent and independent modulus of elasticity.展开更多
We investigate the forming of gamma-ray burst pulses with a simple onedimensional relativistic shock model. The mechanism is that a "central engine" drives forward the nearby plasma inside the fireball to generate a...We investigate the forming of gamma-ray burst pulses with a simple onedimensional relativistic shock model. The mechanism is that a "central engine" drives forward the nearby plasma inside the fireball to generate a series of pressure waves. We give a relativistic geometric recurrence formula that connects the time when the pressure waves are produced and the time when the corresponding shocks occurred. This relation enables us to relate the pulse magnitude with the observation time. Our analysis shows that the evolution of the pressure waves leads to a fast rise and an exponential decay pulses. In determining the width of the pulses, the acceleration time is more important than that of the deceleration.展开更多
Al/Ni-ZrO2 gradient thennal barrier coatings are made on aluminum substrate using plasma spraying method and one direction thermal shock properties of the coatings are studied in this paper. The results show that pore...Al/Ni-ZrO2 gradient thennal barrier coatings are made on aluminum substrate using plasma spraying method and one direction thermal shock properties of the coatings are studied in this paper. The results show that pores in coatings link to form cracks vertical to coating surface. They go through the whole ZrO2 coating once vertical cracks form. When thermal shock cycles increase, horizontal cracks that result in coatings failure forms in the coatings and interface. And vertical cracks delay appearance of horizontal cracks and enhance thermal shock property of coatings. Failure mechanisms of coating thermal shock are discussed using experiments and finite element method.展开更多
The nonlinear propagation of the DIA (dust ion-acoustic) waves in multi-ion dense plasma system containing degener- ate electrons, both positive and negative ions, arbitrary charged dust grains has been investigated b...The nonlinear propagation of the DIA (dust ion-acoustic) waves in multi-ion dense plasma system containing degener- ate electrons, both positive and negative ions, arbitrary charged dust grains has been investigated by employing the reductive perturbation method. The nonlinear waves (solitary and shock waves) have been observed to be formed in case of both positive and negative charged dust grains from the stationary solution of the Korteweg de-Vries (K-dV) equation and Burger’s equation. The fundamental properties of such nonlinear waves have been theoretically analyzed by comparing system potential for both positive and negative dust grains. It has been shown that the basic features of these waves are significantly modified by the positive and negative ions drift speed and polarities of dust grains. The implications of our results in space and laboratory plasmas are briefly discussed.展开更多
In the paper, the problems of dynamie response of multi-rigidbody systems with external impulsive forces are diseussed, and a set of equations with form of Lagrange method is presented. These equations are easy to be ...In the paper, the problems of dynamie response of multi-rigidbody systems with external impulsive forces are diseussed, and a set of equations with form of Lagrange method is presented. These equations are easy to be caleulated by computer.展开更多
This paper investigates the effect of inflow, outflow and shock waves in a single lane highway traffic flow problem. A constant source term has been introduced to demonstrate the inflow and outflow. The classical Ligh...This paper investigates the effect of inflow, outflow and shock waves in a single lane highway traffic flow problem. A constant source term has been introduced to demonstrate the inflow and outflow. The classical Lighthill Whitham and Richards (LWR) model combined with the Greenshields model is used to obtain analytical and numerical solutions. The model is treated as an IBVP and numerical solutions are presented using Lax Friedrichs scheme. Godunov method is also used to present shock wave analysis. The numerical procedures adopted in this investigation yield results which are very much consistent with real life scenario in terms of traffic density and velocity.展开更多
基金Project supported by the National Natural Science Foundation of China (Grant Nos.10732010,10972010,and 11028206)the Opening Project of State Key Laboratory of Explosion Science and Technology,China (Grant No.KFJJ13-5M)
文摘An essentially conservative adaptive space time conservation element and solution element (CE/SE) method is pro- posed for the effective simulation of shock-induced instability with low computational cost. Its implementation is based on redefined conservation elements (CEs) and solution elements (SEs), optimized interpolations and a Courant number insensitive CE/SE scheme. This approach is used in two applications, the Woodward double Mach reflection and a two- component Richtmyer-Meshkov instability experiment. This scheme reveals the essential features of the investigated cases, captures small unstable structures, and yields a solution that is consistent with the results from experiments or other high order methods.
基金the State Key Basic Research Project(G1998040703),China.
文摘In order to analyze the mechanism of continental strong shocks from the angle of crustal movement using the data of repeated geodetic survey, this paper has proposed a physical method; it analyzes the mechanism of density change due to the occurrence of strong shocks by use of the physical quantity that reflects the time change of crustal density. (1) The general theory of the time change of density in the earth’ s interior and the theory of the time change of single layer density have been introduced, and an algorithm of stepwise iteration has been proposed; (2) The effect of the change of single layer density caused by fault dislocation has been analyzed in brief; (3) The characteristics of the time change of crustal density in the south of the seismogenic region before the 1996 Lijiang earthquake with M_L =7.0 have been studied; (4) The precursor model or causal mechanism of strong shocks possibly existing in the time change of crustal density has been investigated preliminarily.
基金Project supported by the National Natural Science Foundation of China(Nos.50909017,51109031, 50921001,11072053,and 51009022)the Doctoral Foundation of Ministry of Education of China(No.20100041120037)+1 种基金the Fundamental Research Funds for the Central Universities (Nos.DUT12LK52 and DUT12LK34)the Major State Basic Research Development Program of China(973 Program)(Nos.2010CB832704 and 2013CB036101)
文摘An enhanced differential transform method (EDTM), which introduces the Pad@ technique into the standard differential transform method (DTM), is proposed. The enhanced method is applied to the analytic treatment of the shock wave. It accelerates the convergence of the series solution and provides an exact Dower series solution.
文摘A new method is presented to analyze multi-degree-of-freedom (MDOF) dynamic systems subjected to an external shock excitation. A two-degree-of-freedom theoretical system with linear characteristics is exemplified to illustrate the procedure of this method. The equations of motion of the dynamic system are established via matrix method. The dynamic responses of the dynamic system under an external shock excitation of a half-sine type are obtained by MATLAB and ANSYS. It is proved that the new method is helpful to analyze MDOF dynamic systems.
基金The project supported by the National Natural Science Foundation of China(59906005)the Teaching Research Award Program for Outstanding Young Teachers in High Education Institutions of MOE,China
文摘In this paper, the CE/SE method is developed to simulate the two- and three-dimensional flow-field of Pulse Detonation Engine (PDE). The conservation equations with stiff source terms for chemical reaction are solved in two steps. The detailed analysis of computational results of a PDE with a single detonation tube and a PDE with five detonation tubes are given in this paper. Complex wave systems are observed inside and outside a PDE. For a PDE with 5 detonation tubes, there is a big bow shock produced from a number of little shocks near the open ends of tubes. A lot of vortexes interact with shocks and a large expansion wave propagates forward and backward with respect to the PDE in a semi-oval shape.
文摘Culture shock is one of the common experiences by students and people abroad. They may experience the feeling of anxiety, disorientation, confusion, and even hostility because of misunderstanding to each other."Culture shock"refers to the phenomenon in which people engage in the interaction feel depressed, discomfortable, etc. Culture shock may be caused by the conflict of various value systems when students and people enter a new culture and other value system. Culture shock is also a consequence of the influence of the negative life event. The changes in one's life and lack of support from family members,friends will make them more vulnerable and cause great discomfort in the new environment.
基金Supported by the National Natural Science Foundation of China.
文摘Many problems at the forefront of theoretical astrophysics require a treatment of dynamical fluid behavior. We present an efficient high-resolution shock-capturing hydrodynamic scheme designed to study such phenomena. We have implemented a weighted, essentially non-oscillatory (WENO) scheme to fifth order accuracy in space, lILLE approximate Riemann solver is used for the flux computation at cell interface, which does not require spectral decomposition into characteristic waves and so is computationally friendly. For time integration we apply a third order total variation diminishing (TVD) Runge-Kutta scheme. Extensive testing and comparison with schemes that re.quire characteristic decomposition are carried out demonstrating the ability of our scheme to address challenging open questions in astrophysics.
基金National Science Foundation for Young Scientists of China under Grant No.51709176National Natural Science Foundation of China under Grant No.51979170+2 种基金Key Project of Hebei Natural Science Foundation under Grant No.F2019210243Hebei Province Science Foundation for Young Scientists under Grant No.E2018210046Open Project of State Key Laboratory of Advanced Electromagnetic Engineering and Technology under Grant No.AEET 2019KF005
文摘With the application of electronic detonators, millisecond blasting is regarded as a signifi cant promising approach to improve the rock fragmentation in deep rock blasting. Thus, it is necessary to investigate the fracturing mechanisms of short-delay blasting. In this work, a rectangle model with two circle boreholes is modeled as a particles assembly based on the discrete element method to simulate the shock wave interactions induced by millisecond blasting. The rectangle model has a size of 12 × 6 m (L × W) and two blast holes have the same diameter of 12 cm. The shock waves are simplifi ed as time-varying forces applied at the particles of walls of the two boreholes. Among a series of numerical tests in this study, the spacing between two adjacent boreholes and delay time of millisecond blasting are considered as two primary variables, and the decoupling charge with a coeffi cient of 1.5 is taken into account in each case. The results show that stress superposition is not a key factor for improving rock fragmentation (tensile stress interactions rather than compressive stress superposition could aff ect the generation of cracks), whereas collision actions from isolated particles or particles with weakened constraints play a crucial role in creating the fracture network. The delay time has an infl uence on causing cracks in rock blasting, however, whether it works heavily depends on the distance between the two holes.
文摘High order accurate scheme is highly desirable for Slow computation with shocks. After analysis has been made for the reason of the generation of non-physical oscillations around the shock in numerical computations, a third-order, upwind biased, shock capturing scheme was proposed. Also, a new shock fitting method, called pseudo shock fitting method, was suggested, which in principle can be with any order of accuracy. Test cases for one dimensional flows show that the new method is very satisfactory.
文摘The finite element method is presented to attain the numerical simulation of the residual stresses field in the material treated by laser shock processing. The distribution of residual stresses generated by a single laser shock with square and round laser spot is predicted and validated by experimental results. With the Finite Element Method (FEM) model, effects of different overlapping rates and impact sequences on the distribution of residual stresses are simulated. The results indicate that: (1) Overlapping laser shock can increase the compressive residual stresses. However, it is not effective on the growth of plastically affected depth; (2) Overlapping rate should be optimized and selected carefully for the large area treatment. Appropriate overlapping rate is beneficial to obtain a homogeneous residual stress field; (3) The impact sequence has a great effect on the residual stress field. It can greatly attenuate the phenomenon of the “residual stress hole” to obtain a homogeneous residual stress field.
基金The research of Gui-Qiang G.Chen was supported in part by the UK Engineering and Physical Sciences Research Council Awards EP/L015811/1,EP/V008854/1,EP/V051121/1the Royal Society-Wolfson Research Merit Award WM090014.
文摘We are concerned with global solutions of multidimensional(M-D)Riemann problems for nonlinear hyperbolic systems of conservation laws,focusing on their global configurations and structures.We present some recent developments in the rigorous analysis of two-dimensional(2-D)Riemann problems involving transonic shock waves through several prototypes of hyperbolic systems of conservation laws and discuss some further M-D Riemann problems and related problems for nonlinear partial differential equations.In particular,we present four different 2-D Riemann problems through these prototypes of hyperbolic systems and show how these Riemann problems can be reformulated/solved as free boundary problems with transonic shock waves as free boundaries for the corresponding nonlinear conservation laws of mixed elliptic-hyperbolic type and related nonlinear partial differential equations.
文摘The 'trial function method' ( TFM for short) and a routine way in finding traveling,wave solutions to some nonlinear partial differential equations( PDE for short), wer explained. Two types of evolution equations are studied, one is a generalized Burgers or KdV equation, the other is the Fisher equation with special nonlinear forms of its reaction rate term. One can see that this method is simple, fast and allowing further extension.
基金This project is supported by Provincial Natural Science Foundation of Liaoning (972028) .Manuscript received on September 16, 19
文摘The deformation of reeds in automotive hydraulic shock absorbers is analyzed with the finite element method. Combination of different thick reeds mounted on different supports is studied. The computational results show that deformation of the overlapped reeds is not always equal to the sum of deflection of single reed under any conditions. Experimental results prove computational results to be correct and computational method effective. The method of analysis and view of point can provide reference to the design and manufacture of hydraulic shock absorbers using reeds.
文摘In this paper, we construct the equations of generalized thermoelasicity for a non-homogeneous isotropic hollow cylider with a variable modulus of elasticity and thermal conductivity based on the Lord and Shulman theory. The problem has been solved numerically using the finite element method. Numerical results for the displacement, the temperature, the radial stress, and the hoop stress distributions are illustrated graphically. Comparisons are made between the results predicted by the coupled theory and by the theory of generalized thermoelasticity with one relaxation time in the cases of temperature dependent and independent modulus of elasticity.
基金Supported by the National Natural Science Foundation of China.
文摘We investigate the forming of gamma-ray burst pulses with a simple onedimensional relativistic shock model. The mechanism is that a "central engine" drives forward the nearby plasma inside the fireball to generate a series of pressure waves. We give a relativistic geometric recurrence formula that connects the time when the pressure waves are produced and the time when the corresponding shocks occurred. This relation enables us to relate the pulse magnitude with the observation time. Our analysis shows that the evolution of the pressure waves leads to a fast rise and an exponential decay pulses. In determining the width of the pulses, the acceleration time is more important than that of the deceleration.
文摘Al/Ni-ZrO2 gradient thennal barrier coatings are made on aluminum substrate using plasma spraying method and one direction thermal shock properties of the coatings are studied in this paper. The results show that pores in coatings link to form cracks vertical to coating surface. They go through the whole ZrO2 coating once vertical cracks form. When thermal shock cycles increase, horizontal cracks that result in coatings failure forms in the coatings and interface. And vertical cracks delay appearance of horizontal cracks and enhance thermal shock property of coatings. Failure mechanisms of coating thermal shock are discussed using experiments and finite element method.
文摘The nonlinear propagation of the DIA (dust ion-acoustic) waves in multi-ion dense plasma system containing degener- ate electrons, both positive and negative ions, arbitrary charged dust grains has been investigated by employing the reductive perturbation method. The nonlinear waves (solitary and shock waves) have been observed to be formed in case of both positive and negative charged dust grains from the stationary solution of the Korteweg de-Vries (K-dV) equation and Burger’s equation. The fundamental properties of such nonlinear waves have been theoretically analyzed by comparing system potential for both positive and negative dust grains. It has been shown that the basic features of these waves are significantly modified by the positive and negative ions drift speed and polarities of dust grains. The implications of our results in space and laboratory plasmas are briefly discussed.
文摘In the paper, the problems of dynamie response of multi-rigidbody systems with external impulsive forces are diseussed, and a set of equations with form of Lagrange method is presented. These equations are easy to be caleulated by computer.
文摘This paper investigates the effect of inflow, outflow and shock waves in a single lane highway traffic flow problem. A constant source term has been introduced to demonstrate the inflow and outflow. The classical Lighthill Whitham and Richards (LWR) model combined with the Greenshields model is used to obtain analytical and numerical solutions. The model is treated as an IBVP and numerical solutions are presented using Lax Friedrichs scheme. Godunov method is also used to present shock wave analysis. The numerical procedures adopted in this investigation yield results which are very much consistent with real life scenario in terms of traffic density and velocity.