CN101944151B - Wall boundary simulation method in molecular dynamics simulation - Google Patents
Wall boundary simulation method in molecular dynamics simulation Download PDFInfo
- Publication number
- CN101944151B CN101944151B CN2010102996396A CN201010299639A CN101944151B CN 101944151 B CN101944151 B CN 101944151B CN 2010102996396 A CN2010102996396 A CN 2010102996396A CN 201010299639 A CN201010299639 A CN 201010299639A CN 101944151 B CN101944151 B CN 101944151B
- Authority
- CN
- China
- Prior art keywords
- particle
- wall
- coordinate
- speed
- reflection
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Images
Abstract
The invention discloses a wall boundary simulation method in molecular dynamics simulation. The method integrates the advantages of a total reflection method and a random reflection method, can truly reflect the motion state of a particle after colliding with a wall and is also called as a semi-reflection method. The method comprises the following steps of: (1) defining a data structure of the particle in the molecular dynamics simulation; (2) determining the current position A (x0, y0, z0, u0, v0, w0) of the particle in a channel; (3) determining the motion direction of the particle according to the current position of the particle and information of acting force by other particles of the particle, and judging whether the particle can move to other positions beyond the wall boundary in a next step or not according to the current position and the acting force information of the particle; and (4) if the particle moves to a position beyond the wall boundary, reflecting the particle through the wall to a next reflection position B (x, y, z, u, v, w), wherein if the wall exists in the x direction, y is equal to y', z is equal to z', x is equal to 2*x' (which is the coordinate of the wall in the x direction), and the speed direction of the particle is random. In the simulation of the invention, physical conditions of wall constant temperature and wall roughness can be reflected, the position of the reflected particle can be more truly reflected, and the method also can be applied to an irregular nano channel to simplify the calculating method of a velocity vector.
Description
Technical field
The present invention relates to the method in a kind of molecular dynamics simulation, especially the method for the movement locus of simulation particle under wall boundary condition.
Background technology
The dynamics that present increasing scientist is flowed with fluid and nanochannel inner fluid in the molecular dynamics method simulation parallel flat.Along with going deep into of research; Scientists begins the non-straight nanochannel is studied and adopted the molecular dynamics method that this passage is simulated, as: people such as Xi-Jun Fan adopt the non equilibrium molecular dynamics method to simulate the flow characteristics of simple fluid in the periodicity jet pipe; For the research roughness flows and the boundary slip The properties to the microchannel endosome, people such as Cao's handle sun have studied argon gas flowing in zigzag platinum passage; People such as J.Castillo-Tejas simulate at the flow process in 4: 1: 4 contraction-extended channels linear polymerization chain and Newtonian fluid respectively.When using the molecular dynamics method to simulate between two flat boards the flow characteristics with the nanochannel inner fluid, in the former document virtual thermal wall (Virtualthermal wall) situation lower wall surface simulation of boundary condition there are two kinds of methods, a kind of is Allen; M.P.and D.J.Tildesley, the total reflection method that in document Computer Simulation of Liquids., proposes, total reflection method is regarded wall as smooth mirror surface; Direct reflection takes place in particle on wall; Adopt the movement locus of total reflection method simulation particle under wall boundary condition referring to Fig. 1 (a), when particle when the A point moves to B ', owing to the existence of virtual wall; Particle can not occur in the position of B '; Particle can move to the B point after wall generation total reflection, B and B ' point is about the wall symmetry, remains unchanged with the velocity of wall parallel coordinates direction; Remain unchanged with the velocity size of wall vertical coordinate direction, direction is reverse; Though total reflection method can reflect the position of particle after launching truly; But total reflection method can not simulate particle at the movement locus of deciding under the rough condition of wall temperature and wall; And if in the non-straight nanochannel; If adopt total reflection method to come the track under the wall boundary condition of simulation particle, then also have certain degree of difficulty aspect the computing velocity vector; Another kind is Raraport; D.C. in document The art of molecular dynamics simulation., the random reflected method has been proposed; The random reflected method is corresponding with casual emission physically, adopts the movement locus of random reflected method simulation particle under wall boundary condition referring to Fig. 1 (b), and particle reflexes to B ' does vertical some B from vertical line to wall; Velocity magnitude confirms that by wall surface temperature direction at random.Though the random reflected method can simulate the physical condition that wall is rough and decide wall temperature; But the random reflected method since the position after the particle reflection in the unicity of wall vertical direction; Therefore in the process of simulation particle movement locus, there is certain untrue property, can not truly reflects the position after particle reflects.
Therefore, how the motion state of real simulation particle under wall boundary condition becomes the technical matters that the researchist in this area is badly in need of solving.
Summary of the invention
Above-mentioned deficiency to the prior art existence; It is a kind of in molecular dynamics simulation that the present invention provides; When particle is in wall boundary condition following time, comprehensively total reflection method and random reflected method advantage, the wall boundary condition simulation method of the motion state of real embodiment particle under wall boundary condition.
In order to realize the foregoing invention purpose, the present invention takes following technical scheme: the analogy method on wall border in the molecular dynamics simulation is characterized in that said analogy method comprises following steps:
A) data structure of particle in the definition molecular dynamics simulation; Particle is represented the motion state of particle in the molecular dynamics simulation passage by six information in three-dimensional model, said six information are respectively x direction coordinate, y direction coordinate, z direction coordinate, x direction speed, y direction speed and z direction speed; Wherein x direction coordinate is the wall horizontal ordinate, and y direction coordinate is the wall ordinate, and the z direction is the direction vertical with the XY plane;
B) confirm the current location A (x of particle in passage
0, y
0, z
0, u
0, v
0, w
0), x wherein
0The current position coordinates of expression particle on the x direction, y
0The current position coordinates of expression particle on the y direction, z
0The current position coordinates of expression particle on the z direction, u
0The current location speed of expression particle on the x direction, v
0The current location speed of expression particle on the y direction, w
0The current location speed of expression particle on the z direction;
C) receive the acting force information of other particle according to particle current location and particle, confirm the direction of motion of particle, and according to the particle current location and receive force information to judge particle whether next step can move to the wall border with external position; If particle does not move to the wall border with external position, then particle can move along the direction of motion of confirming, arrive next step assigned address B ' (x '; Y ', z ', u '; V ', w '), wherein x ' is the assigned address coordinate of particle on the x direction; The assigned address coordinate of y ' expression particle on the y direction, the assigned address coordinate of z ' expression particle on the z direction, the assigned address speed of u ' expression particle on the x direction; The assigned address speed of v ' expression particle on the y direction, the assigned address speed of w ' expression particle on the z direction;
D) if particle movement is arrived the wall border with external position; Particle can not arrive next step appointed positions B '; Particle is according to the temperature of wall; Reflection arrives next step reflection position B (x, y, z through wall; U; V, w), when there is wall in the x direction; Y=y ' wherein then; Z=z ', the coordinate-x ' of x=2* wall x direction, the particle's velocity direction is at random for
; When there is wall in the y direction; X=x ' wherein then; Z=z '; Coordinate-the y ' of y=2* wall y direction, the particle's velocity direction is at random for
; When there is wall in the z direction; X=x ' wherein then; Y=y '; Coordinate-the z ' of z=2* wall z direction, the particle's velocity direction at random for
.
In the above-mentioned formula; (document The art of moleculardynamics simulation (D.C.Rapaport is adopted in obtaining of this random number to a for the random number between-1～1; Second Edition, Cambridge University Press, 2004) in program calculate; Do not describe in detail at this), k
BBe Boltzmann constant, T is a Current Temperatures, and m is a quality; X representes the reflection position coordinate of particle on the x direction; Y representes the reflection position coordinate of particle on the y direction, and z representes the reflection position coordinate of particle on the z direction, and u is the reflection position speed of particle on the x direction; V is the reflection position speed on the particle y direction, and w is the reflection position speed of particle on the z direction.
The analogy method on wall border in the molecular dynamics simulation of the present invention, its cardinal principle are on the direct reflection basis, to have considered the influence of wall roughness to particle rapidity, and therefore, the inventive method can be called the half reflection method again.
Generally, the macroscopic quantity of statistical system can promptly be passed through particle respectively at x in the inventive method through the motion state of particle; Y, coordinate on the z direction and particle be respectively at x, y, and the speed of z direction reflects the motion state of particle; When the particle movement state changes, x, y; Coordinate on the z direction and particle be respectively at x, y, and the speed of z direction all can change.When adopting the half reflection method particle to be in movement locus under the wall boundary condition to simulate, it had both considered the randomness of particle, promptly particle after running into wall and reflecting, its reflection position speed (u; V w) can be relevant according to the temperature of wall, and its reflection position x, y; (x, y z) then do not meet the x of next step predetermined assigned address of wall, y to z direction coordinate with particle; Z direction coordinate (x ', y ', z ') relevant.Therefore this method is not only considered the temperature of wall when next step position of simulation particle, considered that also particle reflects through wall after, its position must with particle relation is not arranged through the position of reflection.
The beneficial effect that the present invention produces:
(1) the inventive method not only can embody permanent wall temperature condition and the rough physical condition of wall in molecular dynamics simulation, and can in molecular dynamics simulation, truly reflect the position after the particle reflection.It has solved in the random reflected method can not reflect the problem of position afterwards by the actual response particle.
(2) the inventive method can also be simplified the velocity computing method in irregular nanochannel, thereby is convenient to calculate.Solved that the total reflection boundary method can not be simulated the problem of the real motion state of particle under permanent wall temperature condition and the rough condition of wall and in irregular nanochannel, total reflection method can not be simplified the problem of velocity computing method.
(3) the inventive method can be used for the motion conditions of simulation particle after outside various different passages move to virtual wall border, thereby for research particle motion state in different passages, technical support is provided.
Description of drawings
Fig. 1 is the particle movement track synoptic diagram in the molecular dynamics simulation; Wherein (a) is for adopting the movement locus synoptic diagram of total reflection method simulation particle under wall boundary condition; Wherein (b) is for adopting the movement locus synoptic diagram of random reflected method simulation particle under wall boundary condition;
Fig. 2 is that the particle in the molecular dynamics simulation reflects synoptic diagram under two kinds of boundary conditions: (a) half reflection particle trajectory motion synoptic diagram wherein; (b) random reflected particle trajectory motion synoptic diagram;
Fig. 3 is the variable cross section access diagram in the molecular dynamics simulation; Wherein (a) is three-dimensional variable cross section access diagram; (b) be three-dimensional variable cross section passage front elevational schematic;
Fig. 4 is that population density distributes with the comparison diagram to the wall distance in the molecular dynamics simulation;
Fig. 5 is a potential energy change curve in time in the inventive method simulation variable cross section nanochannel.
Embodiment
Below in conjunction with accompanying drawing and embodiment the present invention is done explanation in further detail.
A kind of new wall boundary condition disposal route---half reflection method that combines total reflection method and random reflected method advantage is proposed among this paper.When particle A point moves to B ', like Fig. 1 (a), because the existence of wall; Particle arrives the B point through wall reflection back, and coordinate position is identical with total reflection method particle reflection back, but velocity magnitude is confirmed by wall surface temperature; Direction at random, promptly definite method of velocity is identical with the random reflected method.When simulation non-straight micro-nano passage, such as non-45 ° of ramp ways, if wall adopts total reflection boundary condition treatment method, after the particle reflection, velocity define certain degree of difficulty.The random reflected method has implied two conditions with respect to total reflection method: (1) wall constant temperature; (2) wall is rough, and is more approaching with truth.But the random reflected method has certain untrue property when the simulation nanochannel.As shown in Figure 2, the x direction is a flow direction, when regional A is arrived in particle movement; If adopt the random reflected method; See Fig. 2 (b), particle B appears at B ' after being in the optional position process coordinate adjusting among the A, can not go out the coordinate position of particle movement by actual response.If but adopted half reflection boundary condition analogy method, like Fig. 2 (a), particle B would appear at the optional position among the regional A, can appear at the corresponding B ' point arbitrarily of B after employing half reflection method is regulated through coordinate, thereby more approaching with the real motion position of particle.
Present embodiment is through in the molecular dynamics simulation process; Adopt random reflected method, total reflection method and three kinds of wall boundary condition analogy methods of half reflection to simulate the argon process fluid flow in the parallel flat respectively, thereby further specify the difference of the inventive method and existing random reflected method and total reflection method.Practical implementation is following:
Adopt in the ar atmo simulation NVT system, with speed revised law control system temperature, interaction between atoms adopts the 12-6LJ Potential Model; The equation of motion adopts the Verlet numerical integrating; Time step is 1fs, and simulated time is 100000fs, and simulation is of a size of 6nm * 6nm * 1.5nm; Population is 648, and density is 1g/cm
3. in the present embodiment, true origin is defined as 0 point, x and y direction adopt periodic boundary condition; (the Z direction is meant direction vertical with the xy plane in the three-dimensional channel, and particle only possibly exceed the wall border of z direction, because x and y direction adopt periodic boundary condition when the z direction; So all be unlimited the extension in this embodiment; Do not have wall) be 0 o'clock, the xy plane is a lower wall surface, when z is 1.5nm; The xy plane is a upper wall surface, adopts three kinds of different wall boundary condition analogy methods to obtain the motion state of particle below respectively:
Random reflected method: x=x '; Y=y '; Z=0;
Total reflection method: x=x '; Y=y '; Coordinate-the z ' of z=2* wall z direction; U=u '; V=v '; W=-w '
Half reflection method: x=x '; Y=y '; Coordinate-the z ' of z=2* wall z direction;
In the above-mentioned formula, a is the random number between-1～1, k
BBe Boltzmann constant, T is a Current Temperatures, and m is a quality; X representes the reflection position coordinate of particle on x direction coordinate; Y representes the reflection position coordinate of particle on y direction coordinate, and z representes the reflection position coordinate of particle on z direction coordinate, and u is the reflection position speed of particle on x direction coordinate; V is the reflection position speed of particle on y direction coordinate, and w is the reflection position speed of particle on z direction coordinate." * " among this paper is multiplication sign, and "-" is minus sign.
From the simulation of adopting random reflected method, total reflection method and the half reflection method situation of being rebuffed in flowing of the argon fluid the parallel flat under given conditions, when adopting the simulation of random reflected method, referring to Fig. 1 (b) particle A (x
0, y
0, z
0, u
0, v
0, w
0) when in motion process, exceeding the simulated domain border, particle can not arrive next step assigned address B ' (x ', y ', z ', u ', v ', w ') along given line, is vertical some B of vertical line, x=x ' and can reflex to B ' along particle to wall; Y=y '; The coordinate of z=wall z direction; Velocity magnitude is definite by wall surface temperature,
Direction at random; Total reflection method is analog when adopting, referring to Fig. 1 (a), and particle A (x
0, y
0, z
0, u
0, v
0, w
0) when in motion process, exceeding the simulated domain border, particle can not arrive next step assigned address B ' (x ', y ', z ', u ', v ', w ') along given line; Particle can move to B (x, y, z, u, v, w) point, B (x after wall generation total reflection; Y, z, u, v, w) and B ' (x ', y ', z '; U ', v ', w ') point is about the wall symmetry, remains unchanged with the velocity of wall parallel coordinates direction, remains unchanged with the velocity magnitude of wall vertical coordinate direction, and velocity reversal is reverse; When adopting the simulation of half reflection method, particle A (x
0, y
0, z
0, u
0, v
0, w
0) when in motion process, exceeding the simulated domain border, particle not can along given line arrive next step assigned address B ' (x ', y ', z ', u ', v '; W '), and can on virtual wall border, reflect, arrive reflection position B (x, y, z, u; V, w), its reflection back position B both with particle do not take place reflex time assigned address B ' (x ', y ', z '; U ', v ', w ') relevant, also relevant with the temperature of the wall of particle generation reflex time, i.e. x=x '; Y=y '; Coordinate-the z ' of z=2* wall z direction;
Velocity reversal after the particle reflection at random, thereby obtain the position after the particle reflection.Because the speed after the particle reflection and the speed component u of three directions; V is relevant with w; Therefore; Velocity magnitude after the particle reflection is the synthesizing of velocity reversal of three speed components for the direction of
speed; And speed component has and just has negatively, so synthetic velocity reversal also is that all angles all have, thereby the velocity reversal after the particle reflection is at random.In this instructions; The coordinate of the coordinate of the coordinate of wall x direction, wall y direction and wall z direction is meant the coordinate of wall in three-dimensional cartesian coordinate system; Because the initial point of three-dimensional system of coordinate can be 0 point in actual application; Also can not be 0 point, therefore, this coordinate in three-dimensional system of coordinate of wall be relevant with being provided with of three-dimensional system of coordinate initial point.
2) after simulation finishes, added up the situation that population density distributes under three kinds of wall boundary condition analogy method situation, shown in accompanying drawing 4.
Comparative result explanation: adopt half reflection boundary condition simulation and the result 100% who adopts the simulation of total reflection boundary condition simulation to obtain to coincide, have very strong practicality.Thereby show that the inventive method when the motion state of simulation particle under wall boundary condition, not only can embody permanent wall temperature condition and the rough physical condition of wall, and can in molecular dynamics simulation, truly reflect the position after the particle reflection.
Half reflection boundary method and random reflected boundary method and the total reflection boundary method similarities and differences are seen table 1 and table 2:
Table 1 half reflection boundary method and random reflected boundary method and total reflection boundary method characteristic comparison sheet
Table 2 half reflection boundary method and random reflected boundary method and total reflection boundary method identical point
| The random reflected boundary condition | The total reflection boundary condition | |
| Reflection back particle position | Different | Identical |
| The size of reflection back particle rapidity | Identical | Different |
| The direction of reflection back particle rapidity | Identical | Different |
From table 1 and table 2, can find out; Half reflection boundary condition analogy method had both combined to consider in the simulation of random reflected boundary condition the advantage of wall surface temperature and smooth degree; Also combine the variation of particle physical location in the simulation of total reflection boundary condition; Therefore, it is simulated with the total reflection boundary condition with the simulation of random reflected boundary condition and compares, and more can reflect the variation of the motion state after being rebuffed of particle.
Embodiment 2: adopt half reflection boundary condition analogy method, adopt the flow process of tygon fluid in the molecular dynamics simulation variable cross section passage.Concrete steps are following:
1) adopt the present invention's simulation in the variable cross section channel architecture, to carry out, the flow process of tygon molecule is referring to Fig. 3.With speed revised law control system temperature, the equation of motion adopts the Verlet numerical integrating, and time step is 1fs, and simulated time is 320000fs, and density is 0.2588g/cm
3
2) after simulation finishes, added up the potential energy of different time dot system, shown in accompanying drawing 5.Potential energy is potential energy, is one of sign physical quantity of system's macrostate.When adopting random device commonly used, find to move to when surpassing xy simultaneously at next step with optional position, xz plane passages outside when particle, all can only be after particle is launched in 16 somes appearance fixing, this and actual conditions are not inconsistent; And when adopting total reflection method commonly used, when next step moves to above border, passage inclined-plane when particle, be difficult to confirm the velocity after particle reflects; At this moment; Adopting the inventive method (half reflection method) is a kind of authentic and valid mode, because adopt method of the present invention, after particle is launched; Can appear at the optional position in the passage really, and be easy to calculate velocity.
Analog result explanation: the present invention can be applied in the process of molecular dynamics simulation instance, and the motion state that reflects can the real simulated particle takes place after running into virtual wall, is a kind of authentic and valid numerical simulation calculation method.
The inventive method can be used for the motion state after the multiple passage of simulation particle runs into wall, when the real motion state of simulation particle in passage, both need consider the proper motion state of particle; It is the motion state that particle is not run into virtual wall; (analogy method of this normal condition has been ripe prior art, does not also do detailed elaboration at this patent) need be considered the improper motion state of particle again; It is particle motion state after running into virtual wall; In the simulation to this motion state, prior art has proposed total reflection method and random reflected method, but because these two kinds of defectives that method self is intrinsic; Thereby can not effectively simulate the effective status after particle is rebuffed really; Thereby the applicant proposes " half reflection method ", and this method has overcome the defective of total reflection method and random reflected method, can be authentic and valid simulate the motion state after particle is rebuffed.The proposition of this method be effectively the statistics particle in the position at each quarter, velocity magnitude, and further count the temperature of system, macroscopical physical parameters such as kinetic energy, potential energy, thereby be the motion state technical supports of the different particles of research in different passages.
The inventive method can be simulated the motion state after different particles are rebuffed in different passages; Such as: the argon process fluid flow; Other situation such as tygon process fluid flow in hydrone flow process and the various passage; Need to prove that especially after coordinate transformation became cylindrical coordinates or polar coordinates, the inventive method stood good if analog channel is when being circular channel.
Claims (1)
1. the analogy method on wall border in the molecular dynamics simulation is characterized in that, said analogy method comprises following steps:
A) data structure of particle in the definition molecular dynamics simulation; Particle is represented the motion state of particle in the molecular dynamics simulation passage by six information in three-dimensional model, said six information are respectively x direction coordinate, y direction coordinate, z direction coordinate, x direction speed, y direction speed and z direction speed; Wherein x direction coordinate is the wall horizontal ordinate, and y direction coordinate is the wall ordinate, and the z direction is the direction vertical with the xy plane;
B) confirm the current location A (x of particle in passage
0, y
0, z
0, u
0, v
0, w
0), x wherein
0The current position coordinates of expression particle on the x direction, y
0The current position coordinates of expression particle on the y direction, z
0The current position coordinates of expression particle on the z direction, u
0The current location speed of expression particle on the x direction, v
0The current location speed of expression particle on the y direction, w
0The current location speed of expression particle on the z direction;
C) receive the acting force information of other particle according to particle current location and particle, confirm the direction of motion of particle, and according to the particle current location and receive force information to judge particle whether next step can move to the wall border with external position; If particle does not move to the wall border with external position, then particle can move along the direction of motion of confirming, arrive next step assigned address B ' (x '; Y ', z ', u '; V ', w '), wherein x ' is the assigned address coordinate of particle on the x direction; The assigned address coordinate of y ' expression particle on the y direction, the assigned address coordinate of z ' expression particle on the z direction, the assigned address speed of u ' expression particle on the x direction; The assigned address speed of v ' expression particle on the y direction, the assigned address speed of w ' expression particle on the z direction;
D) if particle movement is arrived the wall border with external position; Particle can not arrive next step appointed positions B '; Particle is according to the temperature of wall; Reflection arrives next step reflection position B (x, y, z through wall; U; V, w), when there is wall in the x direction; Y=y ' wherein then; Z=z ', the coordinate-x ' of x=2* wall x direction, the particle's velocity direction is at random for
; When there is wall in the y direction; X=x ' wherein then; Z=z '; Coordinate-the y ' of y=2* wall y direction, the particle's velocity direction is at random for
; When there is wall in the z direction; X=x ' wherein then; Y=y '; Coordinate-the z ' of z=2* wall z direction, the particle's velocity direction is at random for
;
In the above-mentioned formula, ɑ is the random number between-1～1, k
BBe Boltzmann constant, T is a Current Temperatures, and m is a quality; X representes the reflection position coordinate of particle on the x direction; Y representes the reflection position coordinate of particle on the y direction, and z representes the reflection position coordinate of particle on the z direction, and u is the reflection position speed of particle on the x direction; V is the reflection position speed of particle on the y direction, and w is the reflection position speed of particle on the z direction.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2010102996396A CN101944151B (en) | 2010-09-30 | 2010-09-30 | Wall boundary simulation method in molecular dynamics simulation |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2010102996396A CN101944151B (en) | 2010-09-30 | 2010-09-30 | Wall boundary simulation method in molecular dynamics simulation |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN101944151A CN101944151A (en) | 2011-01-12 |
| CN101944151B true CN101944151B (en) | 2012-06-27 |
Family
ID=43436141
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN2010102996396A Expired - Fee Related CN101944151B (en) | 2010-09-30 | 2010-09-30 | Wall boundary simulation method in molecular dynamics simulation |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN101944151B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US10140683B2 (en) | 2011-12-08 | 2018-11-27 | Five3 Genomics, Llc | Distributed system providing dynamic indexing and visualization of genomic data |
| CN109883649A (en) * | 2019-03-14 | 2019-06-14 | 北京石油化工学院 | The method for studying nano-fluid flow behavior in nanochannel by simulation model |
| CN112820356B (en) * | 2021-01-29 | 2023-03-17 | 广东工业大学 | Method for quickly applying molecular dynamics boundary conditions based on geometric boundary operation |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101329772A (en) * | 2008-07-21 | 2008-12-24 | 北京理工大学 | Emulation modelling method interacted with movable object and water based on SPH |
| CN101813555A (en) * | 2010-04-29 | 2010-08-25 | 浙江工业大学 | Soft abrasive flow field testing method based on level set |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4742387B2 (en) * | 2006-04-20 | 2011-08-10 | 独立行政法人海洋研究開発機構 | Simulation method, simulation program, and simulation apparatus |
-
2010
- 2010-09-30 CN CN2010102996396A patent/CN101944151B/en not_active Expired - Fee Related
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101329772A (en) * | 2008-07-21 | 2008-12-24 | 北京理工大学 | Emulation modelling method interacted with movable object and water based on SPH |
| CN101813555A (en) * | 2010-04-29 | 2010-08-25 | 浙江工业大学 | Soft abrasive flow field testing method based on level set |
Non-Patent Citations (2)
| Title |
|---|
| 康细洋,朱恂,廖强.格子波尔兹曼方法模拟外掠圆管的流动与换热.《化学工程》.2010,第38卷(第4期),21-24. * |
| 罗筱薇,石泳,朱恂,廖强.微纳米尺度上含DNA分子流体传输规律的多尺度算法研究进展.《机械工程学报》.2010,第46卷(第14期),143-153. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN101944151A (en) | 2011-01-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Santarelli et al. | Direct numerical simulations of spherical bubbles in vertical turbulent channel flow | |
| Einaudi et al. | Formation of the slow solar wind in a coronal streamer | |
| CN101944151B (en) | Wall boundary simulation method in molecular dynamics simulation | |
| Huber et al. | Vibration-enhanced droplet motion modes: Simulations of rocking, ratcheting, ratcheting with breakup, and ejection | |
| Carlson et al. | A hybrid CFD-DSMC method of modeling continuum-rarefied flows | |
| Banari et al. | The simulation of turbulent particle‐laden channel flow by the Lattice Boltzmann method | |
| Wang et al. | A velocity corrected unresolved CFD-DEM coupled method to reproduce wake effects at moderate Reynolds number | |
| Liu et al. | Particle velocity distribution function around a single bubble in gas-solid fluidized beds | |
| Padgett et al. | Survey of shock-wave structures of smooth-particle granular flows | |
| Ruan et al. | Thin liquid film resulting from a distributed source on a vertical wall | |
| Sun et al. | A Hybrid Continuum/Particle Approach for Micro‐Scale Gas Flows | |
| Karthik | Computational Fluid Dynamics For Biomimetics: Springtails Are Not a Drag | |
| Margha et al. | Dynamic Transition of Unsteady Supersonic Flow From Mach to Regular Reflection Over a Moving Wedge | |
| Zhen et al. | Lattice Boltzmann simulation of flow past a circular cylinder near a moving wall | |
| Pei et al. | Using numerical simulation to optimize the design of gas turbine flowmeter sensor | |
| Brouwers | Statistical models of large scale turbulent flow | |
| Salewski et al. | Effects of aerodynamic particle interaction in turbulent non-dilute particle-laden flow | |
| Aggarwal et al. | Generation of a liquid sheet by an oblique impingement of interacting jets: A Numerical Investigation | |
| McDonald et al. | Use of the Gaussian moment closure for the modelling of continuum and micron-scale flows with moving boundaries | |
| Abdulkadhim | Solid-particles deposition through a turbulent impinging jet using lattice Boltzmann method | |
| Mele et al. | Optimization by cfd analyses of riblet distribution over a transonic civil aircraft configuration | |
| Bangalee et al. | Numerical investigation of flow over a cube-shaped bluff body: An aerodynamic approach | |
| de Zárate et al. | Shear-induced Casimir forces in liquid layers | |
| Kloss et al. | Interaction of a shock wave with a boundary layer in a micro channel | |
| LALGUDI KANNAN | Aerodynamic study of body-gas interaction of a rough surface |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20120627 Termination date: 20140930 |
|
| EXPY | Termination of patent right or utility model |