On accurate boundary conditions for a shape sensitivity equation method

This paper studies the application of the continuous sensitivity equation method (CSEM) for the Navier–Stokes equations in the particular case of shape parameters. Boundary conditions for shape parameters involve flow derivatives at the boundary. Thus, accurate flow gradients are critical to the success of the CSEM. A new approach is presented to extract accurate flow derivatives at the boundary. High order Taylor series expansions are used on layered patches in conjunction with a constrained least‐squares procedure to evaluate accurate first and second derivatives of the flow variables at the boundary, required for Dirichlet and Neumann sensitivity boundary conditions. The flow and sensitivity fields are solved using an adaptive finite‐element method. The proposed methodology is first verified on a problem with a closed form solution obtained by the Method of Manufactured Solutions. The ability of the proposed method to provide accurate sensitivity fields for realistic problems is then demonstrated. The flow and sensitivity fields for a NACA 0012 airfoil are used for fast evaluation of the nearby flow over an airfoil of different thickness (NACA 0015). Copyright © 2005 John Wiley & Sons, Ltd.     

Evaluation of the drag force by integrating the energy dissipation rate in stokes flow for 2D domains using the FEM

The total drag force on the surface of a body, which is the sum of the form drag and the skin friction drag in a 2D domain, is numerically evaluated by integrating the energy dissipation rate in the whole domain for an incompressible Stokes fluid. The finite element method is used to calculate both the energy dissipation rate in the whole domain as well as the drag on the boundary of the body. The evaluation of the drag and the energy dissipation rate are post-processing operations which are carried out after the velocity field and the pressure field for the flow over a particular profile have been obtained. The results obtained for the flow over three different but constant area profiles—a circle, an ellipse and a cross-section of a prolate spheroid—with uniform inlet velocity are presented and it is shown that the total drag force times the velocity is equal to the total energy dissipation rate in the entire finite flow domain. Hence, by calculating the energy dissipation rate in the domain with unit velocity specified at the far-field boundary enclosing the domain, the drag force on the boundary of the body can be obtained.     

Indirect boundary element method for unsteady linearized flow over prolate and oblate spheroids and hemispheroidal protuberances

The indirect boundary element method was used to study the hydrodynamics of oscillatory viscous flow over prolate and oblate spheroids, and over hemispheroidal bodies hinged to a plate. Analytic techniques, such as spheroidal coordinates, method of images, and series representations, were used to make the numerical methods more efficient. A novel method for computing the hydrodynamic torque was used, since for oscillatory flow the torque cannot be computed directly from the weightings. Instead, a Green's function for torque was derived to compute the torque indirectly from the weightings. For full spheroids, the method was checked by comparing the results to exact solutions at low and high frequencies, and to results computed using the singularity method. For hemispheroids hinged to a plate, the method for low frequencies was checked by comparing the results to previous results, and to exact solutions at high frequencies. Copyright © 2004 John Wiley & Sons, Ltd.     

A projection method for solving incompressible viscous flows on domains with moving boundaries

In this paper, a projection method is presented for solving the flow problems in domains with moving boundaries. In order to track the movement of the domain boundaries, arbitrary‐Lagrangian–Eulerian (ALE) co‐ordinates are used. The unsteady incompressible Navier–Stokes equations on the ALE co‐ordinates are solved by using a projection method developed in this paper. This projection method is based on the Bell's Godunov‐projection method. However, substantial changes are made so that this algorithm is capable of solving the ALE form of incompressible Navier–Stokes equations. Multi‐block structured grids are used to discretize the flow domains. The grid velocity is not explicitly computed; instead the volume change is used to account for the effect of grid movement. A new method is also proposed to compute the freestream capturing metrics so that the geometric conservation law (GCL) can be satisfied exactly in this algorithm. This projection method is also parallelized so that the state of the art high performance computers can be used to match the computation cost associated with the moving grid calculations. Several test cases are solved to verify the performance of this moving‐grid projection method. Copyright © 2004 John Wiley Sons, Ltd.     

Temperature factor effect on the structure of the separated flow within a supersonic gas stream

The effect of the temperature factor, that is, the ratio of the body temperature to the freestream stagnation temperature, on the structure of the separated flow formed in the presence of a concave corner in a supersonic stream is studied. The strong influence of the temperature factor on the separation zone length and the flow-generated aerodynamic characteristics is established. It is shown that for fairly large deflection angles this flow cannot be described by free interaction, or triple deck, theory.     

Boundary control of a hyperbolic system with time lags

This paper considers an optimal boundary control problem fora hyperbolic system in which constant time lags appear in thestate equation and in the boundary condition simultaneously.Making use of Lion's scheme, necessary and sufficient conditionsof optimality for the Neumann problem are derived.     

The Enhancement of Efficiency of High-order Harmonic in Intense Laser Field Based on Asymptotic Boundary Conditions and Symplectic Algorithm

Asymptotic boundary condition (ABC) of laser-atom interaction presented recently is applied to transform the initial value problem of the time-dependent Schrödinger equation (TDSE) in infinite space into the initial and boundary value problem in the finite space, and then the TDSE is discretized into linear canonical equations by substituting the symmetry difference quotient for the 2-order partial derivative. The canonical equation is solved by symplectic algorithm. The ground state and the equal weight coherent superposition of the ground state and the first excited state have been taken as the initial conditions, respectively, while we calculate the population of bound states, the evolution of average distance and the high-order harmonic generation (HHG). The conversion efficiency of HHG can be enhanced by initial coherent superposition state and moderate laser intensities     

A mixed hp FEM for the approximation of fourth‐order singularly perturbed problems on smooth domains

We consider fourth‐order singularly perturbed problems posed on smooth domains and the approximation of their solution by a mixed Finite Element Method on the so‐called Spectral Boundary Layer Mesh. We show that the method converges uniformly, with respect to the singular perturbation parameter, at an exponential rate when the error is measured in the energy norm. Numerical examples illustrate our theoretical findings.