spx_density_summation Module

左志华
2022-06-17
96 statements
Source File

针对实粒子进行密度求和：
1. 连续密度法；
2. 正则密度求和法；
3. 密度求和法。

Uses

Contents

Variables

dsm_obj

Abstract Interfaces

density_summation_fcn

Derived Types

density_summation

Subroutines

init denormalized_density_summation normalized_density_summation continuous_density_summation

Variables

Type	Visibility	Attributes		Name		Initial
type(density_summation),	public		::	dsm_obj			密度求和对象

Abstract Interfaces

abstract interface

private subroutine density_summation_fcn(particle, pairs, n)

Density summation
密度求和

Arguments

Type	Intent	Optional	Attributes		Name
type(particle_type),	intent(inout)			::	particle	particle 区域
type(nnps_pairs_type),	intent(in)			::	pairs	particle pairs 粒子对
integer,	intent(in)			::	n	实弹性粒子数, numbers(1)

Derived Types

type, private :: density_summation

密度求和

Components

Type	Visibility	Attributes		Name		Initial
procedure(density_summation_fcn),	public,	pointer, nopass	::	density			密度求和函数指针

Type-Bound Procedures

procedure , public , :: init Subroutine

Subroutines

private subroutine init(self)

初始化

Arguments

Type	Intent	Optional	Attributes		Name
class(density_summation),	intent(inout)			::	self

private subroutine denormalized_density_summation(particle, pairs, n)

对所有实粒子粒子进行非正则化密度求和

Arguments

Type	Intent	Optional	Attributes		Name
type(particle_type),	intent(inout)			::	particle	区域
type(nnps_pairs_type),	intent(in)			::	pairs	粒子对
integer,	intent(in)			::	n	实弹性粒子数 $N_{real}$

private pure subroutine normalized_density_summation(particle, pairs, n)

Normalized density summation
对所有粒子进行正则化密度求和法

Arguments

Type	Intent	Optional	Attributes		Name
type(particle_type),	intent(inout)			::	particle	particle 区域
type(nnps_pairs_type),	intent(in)			::	pairs	particle pairs 粒子对
integer,	intent(in)			::	n	实弹性粒子数, numbers(1) Read more…

private pure subroutine continuous_density_summation(particle, pairs, n)

Continuous density summation
基于连续性方程的连续密度法，计算密度的增量以在时间积分算法中更新密度值

Arguments

Type	Intent	Optional	Attributes		Name
type(particle_type),	intent(inout)			::	particle	particle 区域
type(nnps_pairs_type),	intent(in)			::	pairs	particle pairs 粒子对
integer,	intent(in)			::	n	实弹性粒子数, numbers(1)