Time Step Criteria for Solving Unsteady Engineering Field Problems

Author(s):  
Rabi Mohtar ◽  
Larry Segerlina

Abstract Parabolic equations govern a variety of time-dependent problems in science and engineering. Applying numerical methods such as the finite element or the finite difference method in the space domain changes the partial differential equation to a system ordinary differential equations (ODE’s). Another numerical method is needed to solve the ODE’s in the time domain. The following paper presents time step estimates to solve the system of ODE’s that satisfy both stability and accuracy criteria. These time steps are related to the Froude and Courant Numbers. Suggestions as to which time scheme is the best to use is also presented.

2013 ◽  
Vol 347-350 ◽  
pp. 1758-1762
Author(s):  
Lei Zhang ◽  
Tong Bin Yu ◽  
De Xin Qu ◽  
Xiao Gang Xie

The microstrip circuit is mostly analyzed in transform domain, because its equivalent circuit equation is often nonlinear differential equation, which is easily analyzed in transform domain relatively, but hardly did in time domain, so the analysis of microstrip circuit is a hard work in time domain. In this paper, the FDTD method is used to analyze the microstrip circuit in time domain, by transforming the nonlinear differential equation into time domain iterative equation, selecting suitable time step, and having an iterative computing, the time domain numerical solution can be solved. The FDTD method analyzing the microstrip circuit provides a new way of thought for analyzing microstrip circuit in time domain.


2017 ◽  
Vol 2017 ◽  
pp. 1-12
Author(s):  
Jui-Hsiang Kao

This research develops an Advance-Tracing Boundary Element Method in the time domain to calculate the waves that radiate from an immersed obstacle moving with random acceleration. The moving velocity of the immersed obstacle is multifrequency and is projected along the normal direction of every element on the obstacle. The projected normal velocity of every element is presented by the Fourier series and includes the advance-tracing time, which is equal to a quarter period of the moving velocity. The moving velocity is treated as a known boundary condition. The computing scheme is based on the boundary integral equation in the time domain, and the approach process is carried forward in a loop from the first time step to the last. At each time step, the radiated pressure on each element is updated until obtaining a convergent result. The Advance-Tracing Boundary Element Method is suitable for calculating the radiating problem from an arbitrary obstacle moving with random acceleration in the time domain and can be widely applied to the shape design of an immersed obstacle in order to attain security and confidentiality.


Geophysics ◽  
1974 ◽  
Vol 39 (3) ◽  
pp. 355-355
Author(s):  
Shri Krishna Singh

In this paper Verma obtains a time‐domain solution by inverting the frequency‐domain solution given by Wait (1952). However, it has been recently pointed out by Singh (1973a) (see also Wait, 1973) that there is an error in the quasi‐static solution of Wait. Wait neglected the axially symmetric inducted electric current in the cylinder giving rise to a secondary transverse magnetic field outside (the n=0 term in the scattered wavefield). Singh (1973a) has shown that this term dominates. [It should be noted that Wait in his other works on the cylinder retains this term (e.g., Wait, 1959).] It is clear that this term would be dominant in the time‐domain also. This has been shown by Singh (1972, 1973b). Since the theoretical solution given by Verma in the paper under discussion is incomplete, his interpretation schemes are meaningless.


Author(s):  
Omar Ramadan

Purpose The purpose of this paper is to derive a unified formulation for incorporating different dispersive models into the explicit and implicit finite difference time domain (FDTD) simulations. Design/methodology/approach In this paper, dispersive integro-differential equation (IDE) FDTD formulation is presented. The resultant IDE is written in the discrete time domain by applying the trapezoidal recursive convolution and central finite differences schemes. In addition, unconditionally stable implicit split-step (SS) FDTD implementation is also discussed. Findings It is found that the time step stability limit of the explicit IDE-FDTD formulation maintains the conventional Courant–Friedrichs–Lewy (CFL) constraint but with additional stability limits related to the dispersive model parameters. In addition, the CFL stability limit can be removed by incorporating the implicit SS scheme into the IDE-FDTD formulation, but this is traded for degradation in the accuracy of the formulation. Research limitations/implications The stability of the explicit FDTD scheme is bounded not only by the CFL limit but also by additional condition related to the dispersive material parameters. In addition, it is observed that implicit JE-IDE FDTD implementation decreases as the time step exceeds the CFL limit. Practical implications Based on the presented formulation, a single dispersive FDTD code can be written for implementing different dispersive models such as Debye, Drude, Lorentz, critical point and the quadratic complex rational function. Originality/value The proposed formulation not only unifies the FDTD implementation of the frequently used dispersive models with the minimal storage requirements but also can be incorporated with the implicit SS scheme to remove the CFL time step stability constraint.


Author(s):  
Wei Qiu ◽  
Hongxuan Peng

Based on the panel-free method, large-amplitude motions of floating offshore structures have been computed by solving the body-exact problem in the time domain using the exact geometry. The body boundary condition is imposed on the instantaneous wetted surface exactly at each time step. The free surface boundary is assumed linear so that the time-domain Green function can be applied. The instantaneous wetted surface is obtained by trimming the entire NURBS surfaces of a floating structure. At each time step, Gaussian points are automatically distributed on the instantaneous wetted surface. The velocity potentials and velocities are computed accurately on the body surface by solving the desingularized integral equations. Nonlinear Froude-Krylov forces are computed on the instantaneous wetted surface under the incident wave profile. Validation studies have been carried out for a Floating Production Storage and Offloading (FPSO) vessel. Computed results were compared with experimental results and solutions by the panel method.


2014 ◽  
Vol 14 (02) ◽  
pp. 1350065 ◽  
Author(s):  
K. LIU ◽  
N. ZHANG ◽  
H. XIA ◽  
G. DE ROECK

The interaction between a bridge and a train moving on the bridge is a coupled dynamic problem. The equations of motion of the bridge and the vehicle are coupled by the time dependent contact forces. At each time step, the motion of the bridge influences the forces transferred to the vehicle and this, in turn, changes the forces acting on the bridge. In this paper, a comparison of three different time domain solution algorithms for the coupled equation of motion of the train–bridge system is presented. Guidelines are given for a good choice of the time step.


2012 ◽  
Vol 1 (2) ◽  
pp. 26
Author(s):  
Y. Liu ◽  
B. Ravelo ◽  
J. Ben Hadj Slama

This paper is devoted on the application of the computational method for calculating the transient electromagnetic (EM) near-field (NF) radiated by electronic structures from the frequency-dependent data for the arbitrary wave form perturbations i(t). The method proposed is based on the fast Fourier transform (FFT). The different steps illustrating the principle of the method is described. It is composed of three successive steps: the synchronization of the input excitation spectrum I(f) and the given frequency data H0(f), the convolution of the two inputs data and then, the determination of the time-domain emissions H(t). The feasibility of the method is verified with standard EM 3D simulations. In addition to this method, an extraction technique of the time-dependent z-transversal EM NF component Xz(t) from the frequency-dependent x- and y- longitudinal components Hx(f) and Hy(f) is also presented. This technique is based on the conjugation of the plane wave spectrum (PWS) transform and FFT. The feasibility of the method is verified with a set of dipole radiations. The method introduced in this paper is particularly useful for the investigation of time-domain emissions for EMC applications by considering transient EM interferences (EMIs).


1992 ◽  
Vol 114 (3) ◽  
pp. 175-184 ◽  
Author(s):  
Y. Li ◽  
A. Kareem

The wave forces computed at the displaced position of offshore structures may introduce additional drift forces. This contribution is particularly significant for compliant offshore structures that are configured by design to experience large excursions under the environmental load effects, e.g., tension leg platform. In a random sea environment, this feature can be included in the time domain analysis by synthesizing drag and diffraction forces through a summation of a large number of harmonics with an appropriate phase relationship that reflects the platform displaced position. This approach is not only limited to the time domain analysis, but the superposition of a large number of trigonometric terms in such an analysis requires a considerable computational effort. This paper presents a computationally efficient procedure in both the time and frequency domains that permits inclusion of the time-dependent drift forces, introduced by the platform displacement, in terms of linear and nonlinear feedback contributions. These time-dependent feedback forces are expressed in terms of the applied wave loads by linear and quadratic transformations. It is demonstrated that the results obtained by this approach exhibit good agreement with the procedure based on the summation of trigonometric functions.


Sign in / Sign up

Export Citation Format

Share Document