A Riccati solution for Burgers' equation
Navier-Stokes equation approximation for nonhomogeneous case obtained by relating Burgers equation to Riccati equation through similarity transformation
SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Navier-Stokes equation approximation for nonhomogeneous case obtained by relating Burgers equation to Riccati equation through similarity transformation
Forward Euler's explicit, finite-difference formula of extrapolation, is used as a predictor and a convex formula as a corrector to integrate differential equations numerically. An application has been made to Burger's equation.
Consider the initial boundary value problem for Burgers' equation. It is shown that its solutions converge, in time, to a unique steady state. The speed of the convergence depends on the boundary conditions and can be exponentially slow. Methods to speed up the rate of convergence are also discussed.
The properties of the Karhunen-Loeve expansion of a strongly inhomogeneous random process are examined with emphasis on applications to turbulent flow fields. The ability of the KL expansion to represent functions that have both slow and rapid variations in a relatively small number of expansion terms is tested on a one-dimensional model based on the forced Burgers' equation. The rate of the convergence of the expansion is evaluated, and its dependence on the Reynolds number is determined. It is shown that the KL eigenfunctions possess wall boundary layers attached to outer structures that are independent of the Reynolds number (at high Reynolds numbers). It is also shown that the spectrum of eigenvalues is broad at large Reynolds numbers, requiring many terms to represent higher-order derivatives of the function.
In order to obtain the physically relevant discontinuous numerical solution, the steady-state inviscid Burgers' equation is singularly perturbed through the addition of a small amount of viscosity. A 'cell-centered' finite-difference scheme is proposed which employs two points for the inviscid part and four points for the viscous one. While difficulties are experienced in the capture of interior layers centered at node points, computational results for interior layers centered between node points, and for boundary layers, exhibit accurate nonoscillatory solutions whose discontinuities are captured in one cell on both coarse and fine grids.
A stabilization problem for Burgers' equation is considered. Using linearization, various controllers are constructed which minimize certain weighted energy functionals. These controllers produce the desired degree of stability for the closed-loop nonlinear system. A numerical scheme for computing the feedback gain functional is developed and several numerical experiments are performed to show the theoretical results.
Several variations of the TVD scheme, ENO scheme, FCT scheme, and geometrical schemes, such as MUSCL and PPM, are considered. A comparative study of these schemes as applied to the Burgers' equation is presented. The objective is to assess their performance for problems involving formation and propagation of shocks, shock collisions, and expansion of discontinuities.
The Burgers equation with a small viscosity term, initial and periodic boundary conditions is resolved using a spatial approximation constructed from an orthonormal basis of wavelets. The algorithm is directly derived from the notions of multiresolution analysis and tree algorithms. Before the numerical algorithm is described these notions are first recalled. The method uses extensively the localization properties of the wavelets in the physical and Fourier spaces. Moreover, the authors take advantage of the fact that the involved linear operators have constant coefficients. Finally, the algorithm can be considered as a time marching version of the tree algorithm. The most important point is that an adaptive version of the algorithm exists: it allows one to reduce in a significant way the number of degrees of freedom required for a good computation of the solution. Numerical results and description of the different elements of the algorithm are provided in combination with different mathematical comments on the method and some comparison with more classical numerical algorithms.
A systematic evaluation is conducted of all extant numerical schemes for nonlinear scalar transport problems, and several advanced shock-capturing schemes are used to solve the nonlinear Burgers' equation in order to characterize their ability to resolve the sharp discontinuity, expansion zone, and propagation and collision features of shocks. For discontinuous functions, the Warming-Beam scheme generates preshock wiggles, while the Lax-Wendroff scheme generates postshock ones. Such limiters as the MUSCL or the superbee are more compressive than minimod or monotonic limiters. The performance of such TVD schemes as the upwind, the symmetric, and the Roe-Sweby, resemble each other.
The problem of using reduced order dynamic compensators to control a class of nonlinear parabolic distributed parameter systems was considered. Concentration was on a system with unbounded input and output operators governed by Burgers' equation. A linearized model was used to compute low-order-finite-dimensional control laws by minimizing certain energy functionals. Then these laws were applied to the nonlinear model. Standard approaches to this problem employ model/controller reduction techniques in conjunction with linear quadratic Gaussian (LQG) theory. The approach used is based on the finite dimensional Bernstein/Hyland optimal projection theory which yields a fixed-finite-order controller.
A method for active fluid flow control based on control theory is discussed. Dynamic programming and fixed point successive approximations are used to accommodate the nonlinear control problem. The long-term goal of this project is to establish an effective method applicable to complex flows such as turbulence and jets. However, in this report, the method is applied to stochastic Burgers equation as an intermediate step towards this goal. Numerical results are compared with those obtained by gradient search methods.
The well-known analytical solution of Burgers' equation is extended to curvilinear coordinate systems in three dimensions by a method that is much simpler and more suitable to practical applications than that previously used. The results obtained are applied to incompressible flow with cylindrical symmetry, and also to the decay of an initially linearly increasing wind.
Sonic boom propagation through the atmosphere is modeled with an augmented Burgers’ equation which includes nonlinearity and loss mechanisms. This work details an updated discretization of the governing equations which is fully conservative and duality preserving. Adjoint equations, for all the mechanisms involved, are re-derived and implemented using adjoint consistent discretizations. Computation of loudness metrics is performed using digital filters. The updated implementation is demonstrated and compared against the previous formulation for selected cases, and the differences are documented and discussed. The improved discretization results in faster mesh convergence of the loudness metrics and substantially de-creases runtime. In addition, the adjoint solutions provide mesh-converged gradients which are free from spurious oscillations.
Plastic deformation of crystalline materials with isotropic particle attractions proceeds by the creation and migration of dislocations under the influence of external forces. If dislocations are produced and migrated under the action of local forces, then material shape change can occur without the application of surface forces. We investigate how particles with variable diameters can be embedded in colloidal monolayers to produce dislocations on demand. Additionally, we find in simulation that when embedded clusters of variable diameter particles are taken through multiple cycles of swelling and shrinking, large cumulative plastic slip is produced by the creation and biased motion of dislocation pairs in the solid for embedded clusters of particular geometries. In this way, dislocations emitted by these clusters (biased “dislocation emitters”) can be used to reshape colloidal matter. Our results are also applicable to larger-scale swarms of robotic particles that organize into dense ordered two-dimensional (2D) arrangements. We conclude with a discussion of how dislocations fulfill for colloids the role sought by “metamodules” in lattice robotics research and show how successive applications of shear as a unit operation can produce shape change through slicing and swirling.
The study applies mathematical methods of control theory to the problem of control of fluid flow with the long-range objective of developing effective methods for the control of turbulent flows. Model problems are employed through the formalism and language of control theory to present the procedure of how to cast the problem of controlling turbulence into a problem in optimal control theory. Methods of calculus of variations through the adjoint state and gradient algorithms are used to present a suboptimal control and feedback procedure for stationary and time-dependent problems. Two types of controls are investigated: distributed and boundary controls. Several cases of both controls are numerically simulated to investigate the performances of the control algorithm. Most cases considered show significant reductions of the costs to be minimized. The dependence of the control algorithm on the time-descretization method is discussed.
In this article, we make a few remarks on the role that attractors and inertial manifolds play in fluid mechanics problems. We then describe the role of incremental unknowns for approximating attractors and inertial manifolds when finite difference multigrid discretizations are used. The relation with direct numerical simulation and large eddy simulation is also mentioned.
Explore the source record for details and available documents.
The linear stability of the symmetric Burgers vortex to a class of two-dimensional perturbations is demonstrated. The linear stability of the axially symmetric Burgers vortex to two-dimensional low-Reynolds number disturbances is first discussed. The analysis describes the connection of the normal modes with a classical set of special functions. An expansion in Reynolds number gives the variation of the temporal eigenvalues of the normal modes. Then, Burgers-type solutions are sought where the radially inward-straining field is not axially symmetric. A rational perturbation procedure is used to obtain a new class of nonsymmetric Burgers vortices valid when the Reynolds number and the asymmetry in the external strain are small. Finally, the perturbation solutions are numerically extended to larger values of the parameters.