Multigrid Simulation Stability via Intermediate Grid Buffer
Find Innovative SolutionsGenerate Solutions
Solution Overview
Problem
Current subgridding methods in FDTD techniques for simulating electromagnetic fields face numerical instability issues at the interface between coarse and fine grids, which complicates the accurate modeling of complex geometries and high-frequency wave phenomena, particularly in applications like mobile communication devices, where accuracy is critical for safety and performance.
Innovation Solution
The method involves defining a computational domain with a first grid, a second grid refinement, and an intermediate grid, where the intermediate grid's perimeter is wider than the second grid's by a specific distance, allowing for temporal and spatial coupling to stabilize the simulation, thereby reducing late-time instability and improving accuracy without the need for extensive computational resources.
Engineering Contradictions & Design Principles
Engineering Contradiction Analysis
1Productivity
If subgridding methods are used to reduce computational resources, then productivity is improved, but reliability deteriorates due to numerical instability at grid interfaces
Solution Approach 1:
The patent introduces an intermediate grid between the coarse and fine grids to act as a buffer zone. This intermediate grid prevents direct interaction between the disparate grid resolutions, thereby eliminating the numerical instability that occurs at the interface of traditional subgridding methods while still allowing computational efficiency gains from using coarser grids in most of the domain.
2Measurement precision
If fine grid refinement is applied throughout the entire domain, then measurement precision is improved, but use of energy increases due to higher computational requirements
Solution Approach 1:
The patent applies fine grid refinement only in specific local regions where high precision is required, rather than uniformly across the entire computational domain. This allows the simulation to achieve high measurement precision in critical areas while significantly reducing the overall computational energy requirements by using coarser grids in regions where high precision is not necessary.
Data Source
Figure 1A~1B
Figure 2A
Figure 2B~3
AI summary
A method for assessing wave propagation arising in a physical system by obtaining a numerical approximation of the physical system to be simulated, the method comprising defining a computational domain comprising a first grid having a plurality of first-grid cells and a first-grid time step, a second grid being a refinement of the first grid and having a plurality of second-grid cells and a second-grid time step, and an intermediate grid having a plurality of intermediate-grid cells of equal size to the first-grid cells and having the second-grid time step, wherein each cell has one or more solution points at which values representing a physical quantity of the physical system to be simulated may be obtained; performing a first update procedure to obtain values representing the physical quantity after a said first-grid time step, the values being obtained for at least one solution point of at least one first-grid cell; performing a coupling procedure using values from the first grid to obtain values at a perimeter of the second grid; and performing a second update procedure, using the values at the perimeter of the second grid, to obtain values representing the physical quantity after a said second-grid time step, the values being obtained for at least one solution point of at least one second-grid cell; wherein the coupling procedure comprises performing temporal coupling by using values from the first grid to obtain values at a perimeter of the intermediate grid; performing an intermediate update procedure to obtain values representing the physical quantity after the said second-grid time step of the second update procedure, the values being obtained for at least one solution point of at least one intermediate-grid cell; and performing spatial coupling by using values from the intermediate grid to obtain the values at the perimeter of the second grid used in the second update procedure; wherein defining the computational domain comprises defining the perimeter of the intermediate grid to be wider than the perimeter of the second grid on at least one side of the second grid by a distance equal to at least two times the size of one of the first-grid cells in that dimension.