ME-FDTD Simulation of Anisotropic Magnetized Plasma

Resolve Bottlenecks,
Find Innovative Solutions
Generate Solutions

Solution Overview

Problem

Traditional finite-difference time-domain (FDTD) methods, such as the FDTD (2,2) method, have limited numerical accuracy and stability when simulating electromagnetic wave propagation in anisotropic magnetized plasma, leading to increased memory usage and reduced computational efficiency due to numerical dispersion and anisotropy errors.

Innovation Solution

The proposed method employs a matrix exponential time-domain finite difference (ME-FDTD) method with fourth-order accuracy in both time and space, using symplectic discretization and matrix exponential techniques to improve numerical calculation accuracy and stability, specifically processing Maxwell and polarization current density equations to obtain iterative equations for electric and magnetic field intensities in anisotropic magnetized plasma.

Engineering Contradictions & Design Principles

VSEngineering Contradiction Analysis

1Ease of manufacture

If traditional FDTD (2,2) method is used, then the method is simple and easy to implement, but the numerical accuracy is limited and memory usage increases

Engineering Contradiction:
Improveease of implementationVSAvoidnumerical accuracy
Core Design Contradiction:
Ease of manufactureVSMeasurement precision

Solution Approach 1:

The patent changes the numerical order parameters of the FDTD method from second-order (2,2) to fourth-order (4,4) in both time and space domains. This parameter change improves numerical accuracy by reducing numerical dispersion and anisotropy errors, while maintaining the explicit solving characteristics of the traditional method. The fourth-order accurate difference equations are derived through systematic expansion of the Maxwell equations.

Inventive Principle:
Principle #35Parameter changes

2Measurement precision

If finer Yee grids are used to reduce numerical dispersion and anisotropy errors, then the numerical accuracy improves, but the memory usage increases and computational efficiency decreases

Engineering Contradiction:
Improvenumerical accuracyVSAvoidcomputational efficiency
Core Design Contradiction:
Measurement precisionVSProductivity

Solution Approach 1:

The patent changes the discretization order parameters from second-order to fourth-order in both time and space. This allows achieving the same numerical accuracy with coarser grids, thereby reducing memory usage and improving computational efficiency. The fourth-order accurate difference equations maintain stability under less strict CFL conditions compared to higher-order methods.

Inventive Principle:
Principle #35Parameter changes

3Measurement precision

If FDTD (2,4) method with fourth-order spatial accuracy is used, then the numerical calculation accuracy improves, but long-term error accumulation occurs and stricter CFL conditions are required

Engineering Contradiction:
Improvenumerical calculation accuracyVSAvoidlong-term stability
Core Design Contradiction:
Measurement precisionVSReliability

Solution Approach 1:

The patent replaces the traditional FDTD time-stepping mechanism with a symplectic integration algorithm based on Hamiltonian mechanics. This substitution preserves the symplectic structure of the discretized difference equation, ensuring long-term stability and energy conservation characteristics. The symplectic algorithm avoids long-term error accumulation while maintaining fourth-order accuracy in both time and space.

Inventive Principle:
Principle #28Mechanics substitution (Replace mechanical system)

4Reliability

If symplectic FDTD (4,4) method is used, then the long-term stability and energy conservation are improved, but the numerical discretization form only holds when split matrix Uα=0 which limits applicability

Engineering Contradiction:
Improvelong-term stabilityVSAvoidapplicability to complex dispersion models
Core Design Contradiction:
ReliabilityVSAdaptability or versatility

Solution Approach 1:

The patent segments the Maxwell equations into electric field equations and polarization current density equations, allowing independent treatment of each component. This segmentation enables the symplectic algorithm to handle complex dispersion models with non-zero split matrices by separately discretizing the field components and their corresponding constitutive relations, thereby extending applicability to general anisotropic magnetized plasma.

Inventive Principle:
Principle #1Segmentation

Data Source

PatentUS20240353458A1Method for processing anisotropic magnetized plasma medium and system thereof
Publication Date: 2024.10.24 ANHUI UNIV
  • US20240353458A1 patent drawing
  • US20240353458A1 patent drawing
  • US20240353458A1 patent drawing

AI summary

A method for processing anisotropic magnetized plasma medium and a system thereof are provided. The method includes obtaining the Maxwell equation and the polarization current density equation based on the electromagnetic characteristics of anisotropic magnetized plasma; processing the Maxwell equation and the polarization current density equation to obtain the electric field intensity, the magnetic field intensity, and the polarization current density after processed; based on the electric field intensity, the magnetic field intensity, and the polarization current density after processed, numerical iterative equations for electric field intensity, magnetic field intensity, and polarization current density in anisotropic magnetized plasma medium are obtained; a numerical modeling simulation electromagnetic model is used to determine the electromagnetic characteristics of the electromagnetic model. The use of the present disclosure to simulate the propagation of electromagnetic waves in anisotropic magnetized plasma medium has higher numerical calculation accuracy.