
Browse by Category
 All solutions
 ALE / Parameterized Geometry (1)
 Electromagnetics (1)
 Error Messages (82)
 Fluid Dynamics (2)
 General (12)
 Geometry (15)
 Import (16)
 Installation (61)
 Mesh (11)
 Multiphysics (1)
 Physics (11)
 Postprocessing (1)
 Solver (22)
 Structural Mechanics (3)
 User Models (1)
Solution Number:  1118 
Title:  Resolving timedependent waves 
Platform:  All Platforms 
Applies to:  COMSOL Multiphysics, RF Module, MEMS Module, Structural Mechanics Module, Acoustics Module, Wave Optics Module 
Versions:  All versions 
Categories:  General, Mesh, Solver 
Keywords: 
Problem Description
I am solving a timedependent wavetype problem, what mesh and solver settings should I adjust for better accuracy and runtimes?
Solution
The accuracy of a timedependent model solving for electromagnetic, acoustic, structural, or any other kind of wavetype model is limited by how well the finite element mesh resolves the waves in space and by the how well the time steps resolve the temporal variations. This article describes how to control the mesh and the solver settings for such problems.
Before starting your modeling, you should know two things. First, the maximum frequency content of your excitation. Define this as a Global Parameter, for example: f_max. Second, the material properties in all domains. This information governs what size mesh you will need in all domains, since the material properties govern the shortest possible wavelength. If there are any material or boundary condition nonlinearities, consider that these may lead to frequency doubling or other phenomena that may lead to much higher frequency content. If you anticipate any step changes in the loads or boundary conditions, first read: Knowledge Base 1245: Improving runtimes for transient models with step changes in loads
Once you know the maximum frequencies that you expect to model, and the material properties, you can build and mesh your model. By default, second order, (quadratic) elements are used in the electromagnetic, acoustic, and structural transient interfaces. At a bare minimum, the mesh must fulfill the Nyquist criterion, so at least two elements per local wavelength in each domain. However, this will usually be of unacceptably low accuracy. Practically speaking, at least five elements per local wavelength should be used, but very rarely more than ten elements per wavelength. Set up a Usercontrolled mesh with different Size features applied to the different material domains. The Free Triangular (in 2D models) and Free Tetrahedral (for 3D models) meshes are preferred. Also define a Global Parameter (for example: N) to control the number of elements per local wavelength and a set of Variables to define the local wavelengths in different domains, as demonstrated in the screenshot below.
Manually setting the element sizes
By default, the TimeDependent Solver will continuously adjust the time step in order to fulfill your specified tolerances. The Relative tolerance field in the Time Dependent settings (shown below) controls how the solver will adjust the timestep. Smaller numbers will lead to smaller timesteps and higher accuracy. Note that changing the number of output Times in the Time Dependent node controls only the output times, but has little effect on the time steps actually taken by the solver.
The output timesteps and relative tolerance settings
However, since you already know your maximum local mesh size and maximum frequency, it is more computationally efficient to manually set the timestep. The time step should resolve the wave equally well in time as the mesh does in space. Longer time steps will not make optimal use of the mesh, and any shorter time steps will lead to longer solution times with no considerable improvements to the results. The relationship between wavespeed, c, maximum mesh size, h, and time step length, Δt, is known as the CFL number:
CFL = cΔt/h
We have already manually defined the maximum mesh size to be 1/N of the local wavelength (as described above), so this can rewritten in terms of the frequency f, and N the number of elements per local wavelength:
CFL = fNΔt
or, rearranging:
Δt = CFL/Nf
With the default second order, quadratic, mesh elements, the CFL number should be less that 0.2, and a value of 0.1 proves to be near optimal, it is helpful to define this CFL number as a Global Parameter as well.
Implement the manual timestep size in the TimeDependent Solver settings, usually found at Study > Solver Configurations > Solution > TimeDependent Solver, Time Stepping settings. (If these settings are not already in the Study branch, then rightclick on Study and select Show Default Solver.) Set the Time Stepping Method: to Generalized alpha, set the Steps taken by solver: to Manual and set the Time Step to:
CFL/(N*f_max)
The Generalized alpha solver is preferred for the reasons outlined in: Knowledge Base 1062: BDF, Generalized Alpha, and RungeKutta Methods
Manually setting the time step size
For an example of these mesh and solver modifications, see the Gaussian Explosion in the Application Library of the Acoustics Module.
For an example of modeling a nonlinear material within the RF Module, see the Second Harmonic Generation of a Gaussian Beam. The same example is also available within the Wave Optics Module: Second Harmonic Generation of a Gaussian Beam
See also:
Knowledge Base 1240: Manually Setting the Scaling of Variables
Note: If you are solving model using any Time Explicit formulation, such as the Pressure Acoustics, Time Explicit, Convected Wave Equation, Time Explicit or the Electromagnetic Waves, Time Explicit physics interfaces, this knowledgebase does not apply.
Disclaimer
COMSOL makes every reasonable effort to verify the information you view on this page. Resources and documents are provided for your information only, and COMSOL makes no explicit or implied claims to their validity. COMSOL does not assume any legal liability for the accuracy of the data disclosed. Any trademarks referenced in this document are the property of their respective owners. Consult your product manuals for complete trademark details.