1
$\begingroup$

Some equations (such as the non-linear schrödinger equation for pulse propagation) are more easily solved in the spectral form, but still need a representation in the temporal domain to calculate physical parameters such as the energy of the electric field. This results usually in the following steps (with $E$ the electric field in the temporal domain, and $\hat{E}$ the electric field in the spectral domain): $$E_0\underbrace{\rightarrow}_{\text{FFT}}\hat{E}_0\underbrace{\rightarrow}_{\text{Propagation}}\hat{E}_1\underbrace{\rightarrow}_{\text{IFFT}}E_1$$ which requires two different "grids". Is such a calculation possible using FEM, and if yes, how (after requiring two different grids)?

Background to this question is that I have additional equations coupled to that problem which preferably can be solved using FEM and do not require the transformations from the temporal domain to the spectral domain and back (after they only depend on $E_x$, not on $\hat{E}_x$). For this coupling I'd like to avoid having to have two different grids (one for the non-linear schrödinger equation, one for the FEM-covered equation) which I have to interpolate data in between.

$\endgroup$
5
  • $\begingroup$ If you use a DFT matrix $W$, could you form a coupled system of equations (a block matrix expression) by substituting $\hat{E_x}=WE_x$? en.wikipedia.org/wiki/DFT_matrix $\endgroup$ – Charlie S Mar 17 at 14:26
  • $\begingroup$ I think that you could use the same grid, that's the spatial discretization. What changes is the time discretization, isn't it? $\endgroup$ – nicoguaro Mar 18 at 15:17
  • 1
    $\begingroup$ People do this sort of thing all the time in the geodynamo community. $\endgroup$ – Wolfgang Bangerth Mar 18 at 16:29
  • $\begingroup$ @WolfgangBangerth: Do you have examples for that, by chance? I'm completely unfamiliar with this community. $\endgroup$ – arc_lupus Mar 19 at 9:03
  • 1
    $\begingroup$ I'm not familiar with the literature in that arena either, but you could start by looking at the Rayleigh and Calypso codes and the literatures that surround them. Both are available from geodynamics.org. I bet you can also find papers by Gary Glatzmeier on the topic. $\endgroup$ – Wolfgang Bangerth Mar 19 at 13:37

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Browse other questions tagged or ask your own question.