- NSF-PAR ID:
- 10379976
- Date Published:
- Journal Name:
- Inverse Problems
- Volume:
- 38
- Issue:
- 11
- ISSN:
- 0266-5611
- Page Range / eLocation ID:
- 114001
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
We propose a new iterative scheme to compute the numerical solution to an over-determined boundary value problem for a general quasilinear elliptic PDE. The main idea is to repeatedly solve its linearization by using the quasi-reversibility method with a suitable Carleman weight function. The presence of the Carleman weight function allows us to employ a Carleman estimate to prove the convergence of the sequence generated by the iterative scheme above to the desired solution. The convergence of the iteration is fast at an exponential rate without the need of an initial good guess. We apply this method to compute solutions to some general quasilinear elliptic equations and a large class of first-order Hamilton-Jacobi equations. Numerical results are presented.more » « less
-
Abstract We consider the three-dimensional Euler equations in a domain with a free boundary with no surface tension. We assume that
is such that in an arbitrarily small neighborhood of the free boundary, and we use the Lagrangian approach to derive an a priori estimate that can be used to prove local-in-time existence and uniqueness of solutions under the Rayleigh–Taylor stability condition. -
Abstract High fidelity near-wall velocity measurements in wall bounded fluid flows continue to pose a challenge and the resulting limitations on available experimental data cloud our understanding of the near-wall velocity behavior in turbulent boundary layers. One of the challenges is the spatial averaging and limited spatial resolution inherent to cross-correlation-based particle image velocimetry (PIV) methods. To circumvent this difficulty, we implement an explicit no-slip boundary condition in a wavelet-based optical flow velocimetry (wOFV) method. It is found that the no-slip boundary condition on the velocity field can be implemented in wOFV by transforming the constraint to the wavelet domain through a series of algebraic linear transformations, which are formulated in terms of the known wavelet filter matrices, and then satisfying the resulting constraint on the wavelet coefficients using constrained optimization for the optical flow functional minimization. The developed method is then used to study the classical problem of a turbulent channel flow using synthetic data from a direct numerical simulation (DNS) and experimental particle image data from a zero pressure gradient, high Reynolds number turbulent boundary layer. The results obtained by successfully implementing the no-slip boundary condition are compared to velocity measurements from wOFV without the no-slip condition and to a commercial PIV code, using the velocity from the DNS as ground truth. It is found that wOFV with the no-slip condition successfully resolves the near-wall profile with enhanced accuracy compared to the other velocimetry methods, as well as other derived quantities such as wall shear and turbulent intensity, without sacrificing accuracy away from the wall, leading to state of the art measurements in the
region of the turbulent boundary layer when applied to experimental particle images. -
To walk over constrained environments, bipedal robots must meet concise control objectives of speed and foot placement. The decisions made at the current step need to factor in their effects over a time horizon. Such step-to-step control is formulated as a two-point boundary value problem (2-BVP). As the dimensionality of the biped increases, it becomes increasingly difficult to solve this 2-BVP in real-time. The common method to use a simple linearized model for real-time planning followed by mapping on the high dimensional model cannot capture the nonlinearities and leads to potentially poor performance for fast walking speeds. In this paper, we present a framework for real-time control based on using partial feedback linearization (PFL) for model reduction, followed by a data-driven approach to find a quadratic polynomial model for the 2-BVP. This simple step-tostep model along with constraints is then used to formulate and solve a quadratically constrained quadratic program to generate real-time control commands. We demonstrate the efficacy of the approach in simulation on a 5-link biped following a reference velocity profile and on a terrain with ditches.more » « less
-
We consider the linear third order (in time) PDE known as the SMGTJ-equation, defined on a bounded domain, under the action of either Dirichlet or Neumann boundary control
. Optimal interior and boundary regularity results were given in [\begin{document}$ g $\end{document} 1 ], after [41 ], when , which, moreover, in the canonical case\begin{document}$ g \in L^2(0, T;L^2(\Gamma)) \equiv L^2(\Sigma) $\end{document} , were expressed by the well-known explicit representation formulae of the wave equation in terms of cosine/sine operators [\begin{document}$ \gamma = 0 $\end{document} 19 ], [17 ], [24 ,Vol Ⅱ]. The interior or boundary regularity theory is however the same, whether or\begin{document}$ \gamma = 0 $\end{document} , since\begin{document}$ 0 \neq \gamma \in L^{\infty}(\Omega) $\end{document} is responsible only for lower order terms. Here we exploit such cosine operator based-explicit representation formulae to provide optimal interior and boundary regularity results with\begin{document}$ \gamma \neq 0 $\end{document} "smoother" than\begin{document}$ g $\end{document} , qualitatively by one unit, two units, etc. in the Dirichlet boundary case. To this end, we invoke the corresponding results for wave equations, as in [\begin{document}$ L^2(\Sigma) $\end{document} 17 ]. Similarly for the Neumann boundary case, by invoking the corresponding results for the wave equation as in [22 ], [23 ], [37 ] for control smoother than , and [\begin{document}$ L^2(0, T;L^2(\Gamma)) $\end{document} 44 ] for control less regular in space than . In addition, we provide optimal interior and boundary regularity results when the SMGTJ equation is subject to interior point control, by invoking the corresponding wave equations results [\begin{document}$ L^2(\Gamma) $\end{document} 42 ], [24 ,Section 9.8.2].